Search

We introduce a new computational technique for searching for faint moving sources in astronomical images. Starting from a maximum likelihood estimate for the probability of the detection of a source within a series of images, we develop a massively parallel algorithm for searching through candidate asteroid trajectories that utilizes Graphics Processing Units (GPU). This technique can search over 10^10 possible asteroid trajectories in stacks of the order 10-15 4K x 4K images in under a minute using a single consumer grade GPU. We apply this algorithm to data from the 2015 campaign of the High Cadence Transient Survey (HiTS) obtained with the Dark Energy Camera (DECam). We find 39 previously unknown Kuiper Belt Objects in the 150 square degrees of the survey. Comparing these asteroids to an existing model for the inclination distribution of the Kuiper Belt we demonstrate that we recover a KBO population above our detection limit consistent with previous studies. Software used in this analysis is made available as an open source package., Comment: Accepted to AJ

In this paper, we describe the optical imaging data processing pipeline developed for the Subaru Telescope's Hyper Suprime-Cam (HSC) instrument. The HSC Pipeline builds on the prototype pipeline being developed by the Large Synoptic Survey Telescope's Data Management system, adding customizations for HSC, large-scale processing capabilities, and novel algorithms that have since been reincorporated into the LSST codebase. While designed primarily to reduce HSC Subaru Strategic Program (SSP) data, it is also the recommended pipeline for reducing general-observer HSC data. The HSC pipeline includes high level processing steps that generate coadded images and science-ready catalogs as well as low-level detrending and image characterizations., Comment: 39 pages, 21 figures, 2 tables. Submitted to Publications of the Astronomical Society of Japan

The ESSENCE survey discovered 213 Type Ia supernovae at redshifts 0.1 < z < 0.81 between 2002 and 2008. We present their R and I-band photometry, measured from images obtained using the MOSAIC II camera at the CTIO 4 m Blanco telescope, along with rapid-response spectroscopy for each object. We use our spectroscopic follow-up observations to determine an accurate, quantitative classification and a precise redshift. Through an extensive calibration program we have improved the precision of the CTIO Blanco natural photometric system. We use several empirical metrics to measure our internal photometric consistency and our absolute calibration of the survey. We assess the effect of various potential sources of systematic bias on our measured fluxes, and we estimate that the dominant term in the systematic error budget from the photometric calibration on our absolute fluxes is ~1%., Comment: (40 pages, 21 figures, accepted for publication in ApJS)

We present a study of the statistical flare rates of M dwarfs (dMs) with close white dwarf (WD) companions (WD+dM; typical separations < 1 au). Our previous analysis demonstrated that dMs with close WD companions are more magnetically active than their field counterparts. One likely implication of having a close binary companion is increased stellar rotation through disk-disruption, tidal effects, and/or angular momentum exchange; increased stellar rotation has long been associated with an increase in stellar activity. Previous studies show a strong correlation between dMs that are magnetically active (showing H{\alpha} in emission) and the frequency of stellar flare rates. We examine the difference between the flare rates observed in close WD+dM binary systems and field dMs. Our sample consists of a subset of 181 close WD+dM pairs from Morgan et al. (2012) observed in the Sloan Digital Sky Survey Stripe 82, where we obtain multi-epoch observations in the Sloan ugriz-bands. We find an increase in the overall flaring fraction in the close WD+dM pairs (0.09$\pm$0.03%) compared to the field dMs (0.0108$\pm$0.0007%; Kowalski et al. 2009) and a lower flaring fraction for active WD+dMs (0.05$\pm$0.03%) compared to active dMs (0.28$\pm$0.05%; Kowalski et al. 2009). We discuss how our results constrain both the single and binary dM flare rates. Our results also constrain dM multiplicity, our knowledge of the Galactic transient background, and may be important for the habitability of attending planets around dMs with close companions., Comment: 14 pages, 7 figures

The Large Synoptic Survey Telescope (LSST) is a large-aperture, wide-field, ground-based survey system that will image the sky in six optical bands from 320 to 1050 nm, uniformly covering approximately $18,000$deg$^2$ of the sky over 800 times. The LSST is currently under construction on Cerro Pach\'on in Chile, and expected to enter operations in 2022. Once operational, the LSST will explore a wide range of astrophysical questions, from discovering "killer" asteroids to examining the nature of Dark Energy. The LSST will generate on average 15 TB of data per night, and will require a comprehensive Data Management system to reduce the raw data to scientifically useful catalogs and images with minimum human intervention. These reductions will result in a real-time alert stream, and eleven data releases over the 10-year duration of LSST operations. To enable this processing, the LSST project is developing a new, general-purpose, high-performance, scalable, well documented, open source data processing software stack for O/IR surveys. Prototypes of this stack are already capable of processing data from existing cameras (e.g., SDSS, DECam, MegaCam), and form the basis of the Hyper-Suprime Cam (HSC) Survey data reduction pipeline., Comment: 10 pages, 3 figures, to appear in the Proceedings of ADASS XXV

Active galactic nuclei (AGNs) are well-known to exhibit flux variability across a wide range of wavelength regimes, but the precise origin of the variability at different wavelengths remains unclear. To investigate the relatively unexplored near-IR variability of the most luminous AGNs, we conduct a search for variability using well sampled JHKs-band light curves from the 2MASS survey calibration fields. Our sample includes 27 known quasars with an average of 924 epochs of observation over three years, as well as one spectroscopically confirmed blazar (SDSSJ14584479+3720215) with 1972 epochs of data. This is the best-sampled NIR photometric blazar light curve to date, and it exhibits correlated, stochastic variability that we characterize with continuous auto-regressive moving average (CARMA) models. None of the other 26 known quasars had detectable variability in the 2MASS bands above the photometric uncertainty. A blind search of the 2MASS calibration field light curves for AGN candidates based on fitting CARMA(1,0) models (damped-random walk) uncovered only 7 candidates. All 7 were young stellar objects within the {\rho} Ophiuchus star forming region, five with previous X-ray detections. A significant {\gamma}-ray detection (5{\sigma}) for the known blazar using 4.5 years of Fermi photon data is also found. We suggest that strong NIR variability of blazars, such as seen for SDSSJ14584479+3720215, can be used as an efficient method of identifying previously-unidentified {\gamma}-ray blazars, with low contamination from other AGN., Comment: 6 pages, 3 figures, ApJ Accepted

We present the fiducial main sequence stellar locus traced by 10 photometric colors observed by SDSS, 2MASS, and WISE. Median colors are determined using 1,052,793 stars with r-band extinction less than 0.125. We use this locus to measure the dust extinction curve relative to the r-band, which is consistent with previous measurements in the SDSS and 2MASS bands. The WISE band extinction coefficients are larger than predicted by standard extinction models. Using 13 lines of sight, we find variations in the extinction curve in H, Ks, and WISE bandpasses. Relative extinction decreases towards Galactic anti-center, in agreement with prior studies. Relative extinction increases with Galactic latitude, in contrast to previous observations. This indicates a universal mid-IR extinction law does not exist due to variations in dust grain size and chemistry with Galactocentric position. A preliminary search for outliers due to warm circumstellar dust is also presented, using stars with high signal-to-noise in the W3-band. We find 199 such outliers, identified by excess emission in Ks-W3. Inspection of SDSS images for these outliers reveals a large number of contaminants due to nearby galaxies. Six sources appear to be genuine dust candidates, yielding a fraction of systems with infrared excess of 0.12$\pm$0.05%., Comment: 11 pages, 10 figures, MNRAS Accepted. Tables 1 and 2 available online: https://github.com/jradavenport/wise_locus

We present the use of continuous-time autoregressive moving average (CARMA) models as a method for estimating the variability features of a light curve, and in particular its power spectral density (PSD). CARMA models fully account for irregular sampling and measurement errors, making them valuable for quantifying variability, forecasting and interpolating light curves, and for variability-based classification. We show that the PSD of a CARMA model can be expressed as a sum of Lorentzian functions, which makes them extremely flexible and able to model a broad range of PSDs. We present the likelihood function for light curves sampled from CARMA processes, placing them on a statistically rigorous foundation, and we present a Bayesian method to infer the probability distribution of the PSD given the measured lightcurve. Because calculation of the likelihood function scales linearly with the number of data points, CARMA modeling scales to current and future massive time-domain data sets. We conclude by applying our CARMA modeling approach to light curves for an X-ray binary, two AGN, a long-period variable star, and an RR-Lyrae star, in order to illustrate their use, applicability, and interpretation., Comment: 18 pages, 16 figures, in press at ApJ. Revised to match accepted version. Code for utilizing the CARMA models is at https://github.com/bckelly80/carma_pack

This paper describes the data release of the Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey conducted between 2005 and 2007. Light curves, spectra, classifications, and ancillary data are presented for 10,258 variable and transient sources discovered through repeat ugriz imaging of SDSS Stripe 82, a 300 deg2 area along the celestial equator. This data release is comprised of all transient sources brighter than r~22.5 mag with no history of variability prior to 2004. Dedicated spectroscopic observations were performed on a subset of 889 transients, as well as spectra for thousands of transient host galaxies using the SDSS-III BOSS spectrographs. Photometric classifications are provided for the candidates with good multi-color light curves that were not observed spectroscopically. From these observations, 4607 transients are either spectroscopically confirmed, or likely to be, supernovae, making this the largest sample of supernova candidates ever compiled. We present a new method for SN host-galaxy identification and derive host-galaxy properties including stellar masses, star-formation rates, and the average stellar population ages from our SDSS multi-band photometry. We derive SALT2 distance moduli for a total of 1443 SN Ia with spectroscopic redshifts as well as photometric redshifts for a further 677 purely-photometric SN Ia candidates. Using the spectroscopically confirmed subset of the three-year SDSS-II SN Ia sample and assuming a flat Lambda-CDM cosmology, we determine Omega_M = 0.315 +/- 0.093 (statistical error only) and detect a non-zero cosmological constant at 5.7 sigmas., Comment: Submitted to ApJS. Full catalogs and datafiles are available here: http://sdssdp62.fnal.gov/sdsssn/DataRelease/index.html

With upcoming all sky surveys such as LSST poised to generate a deep digital movie of the optical sky, variability-based AGN selection will enable the construction of highly-complete catalogs with minimum contamination. In this study, we generate $g$-band difference images and construct light curves for QSO/AGN candidates listed in SDSS Stripe 82 public catalogs compiled from different methods, including spectroscopy, optical colors, variability, and X-ray detection. Image differencing excels at identifying variable sources embedded in complex or blended emission regions such as Type II AGNs and other low-luminosity AGNs that may be omitted from traditional photometric or spectroscopic catalogs. To separate QSOs/AGNs from other sources using our difference image light curves, we explore several light curve statistics and parameterize optical variability by the characteristic damping timescale ($\tau$) and variability amplitude. By virtue of distinguishable variability parameters of AGNs, we are able to select them with high completeness of 93.4% and efficiency (i.e., purity) of 71.3%. Based on optical variability, we also select highly variable blazar candidates, whose infrared colors are consistent with known blazars. One third of them are also radio detected. With the X-ray selected AGN candidates, we probe the optical variability of X-ray detected optically-extended sources using their difference image light curves for the first time. A combination of optical variability and X-ray detection enables us to select various types of host-dominated AGNs. Contrary to the AGN unification model prediction, two Type II AGN candidates (out of 6) show detectable variability on long-term timescales like typical Type I AGNs. This study will provide a baseline for future optical variability studies of extended sources., Comment: 21 pages, 22 figures, Accepted for publication in The Astrophysical Journal

We describe the construction of a highly reliable sample of approximately 7,000 optically faint periodic variable stars with light curves obtained by the asteroid survey LINEAR across 10,000 sq.deg of northern sky. Majority of these variables have not been cataloged yet. The sample flux limit is several magnitudes fainter than for most other wide-angle surveys; the photometric errors range from ~0.03 mag at $r=15$ to ~0.20 mag at r=18. Light curves include on average 250 data points, collected over about a decade. Using SDSS-based photometric recalibration of the LINEAR data for about 25 million objects, we selected ~200,000 most probable candidate variables and visually confirmed and classified approximately 7,000 periodic variables using phased light curves. The reliability and uniformity of visual classification across eight human classifiers was calibrated and tested using a SDSS Stripe 82 region variable star catalog, and verified using an unsupervised machine learning approach. The resulting sample of periodic LINEAR variables is dominated by 3,900 RR Lyrae stars and 2,700 eclipsing binary stars of all subtypes, and includes small fractions of relatively rare populations such as asymptotic giant branch stars and SX Phoenicis stars. We discuss the distribution of these mostly uncataloged variables in various diagrams constructed with optical-to-infrared SDSS, 2MASS and WISE photometry, and with LINEAR light curve features. An interesting side result is a robust and precise quantitative description of a strong correlation between the light-curve period and color/spectral type for close and contact eclipsing binary stars. These large samples of robustly classified variable stars will enable detailed statistical studies of the Galactic structure and physics of binary and other stars, and we make them publicly available., Comment: 71 pages, 21 figures. For additional machine readable tables and light curves, see http://www.astro.washington.edu/users/ivezic/linear/PaperIII/PLV.html . Accepted by AJ

We present a sample of ~5,000 RR Lyrae stars selected from the recalibrated LINEAR dataset and detected at heliocentric distances between 5 kpc and 30 kpc over ~8,000 deg^2 of sky. The coordinates and light curve properties, such as period and Oosterhoff type, are made publicly available. We find evidence for the Oosterhoff dichotomy among field RR Lyrae stars, with the ratio of the type II and I subsamples of about 1:4. The number density distribution of halo RRab stars as a function of galactocentric distance can be described as an oblate ellipsoid with the axis ratio q=0.63 and with either a single or a double power law with a power-law index in the range -2 to -3. Using a group-finding algorithm EnLink, we detected seven candidate halo groups, only one of which is statistically spurious. Three of these groups are near globular clusters (M53/NGC 5053, M3, M13), and one is near a known halo substructure (Virgo Stellar Stream); the remaining three groups do not seem to be near any known halo substructures or globular clusters, and seem to have a higher ratio of Oosterhoff type II to Oosterhoff type I RRab stars than what is found in the halo. The extended morphology and the position (outside the tidal radius) of some of the groups near globular clusters is suggestive of tidal streams possibly originating from globular clusters. Spectroscopic followup of detected halo groups is encouraged., Comment: 49 pages, 18 figures, one machine-readable table (downloadable from "Other formats/Source"), AJ in press; updated references and Table 1 to match the published version. The animations referenced in the paper may be downloaded from http://www.astro.caltech.edu/~bsesar/media.html

The Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) observed 10 transits of XO-2b over a period of three years. We present measurements which confirm previous estimates of system parameters like the normalized semi-major axis (a/R_{*}), stellar density (\rho_{*}), impact parameter (b) and orbital inclination (i_{orb}). Our errors on system parameters like a/R_{*} and \rho_{*} have improved by ~40% compared to previous best ground-based measurements. Our study of the transit times show no evidence for transit timing variations and we are able to rule out co-planar companions with masses \ge 0.20 \mearth\ in low order mean motion resonance with XO-2b. We also explored the stability of the XO-2 system given various orbital configurations of a hypothetical planet near the 2:1 mean motion resonance. We find that a wide range of orbits (including Earth-mass perturbers) are both dynamically stable and produce observable TTVs. We find that up to 51% of our stable simulations show TTVs that are smaller than the typical transit timing errors (~20 sec) measured for XO-2b, and hence remain undetectable., Comment: 33 pages, 6 figures, 9 tables (submitted to ApJ, Apr 21 2013). arXiv admin note: text overlap with arXiv:1211.6828

We present transit observations of the WASP-2 exoplanet system by the Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) program. Model fitting to these data allows us to improve measurements of the hot-Jupiter exoplanet WASP-2b and its orbital parameters by a factor of ~2 over prior studies; we do not find evidence for transit depth variations. We do find reduced chi^2 values greater than 1.0 in the observed minus computed transit times. A sinusoidal fit to the residuals yields a timing semi-amplitude of 32 seconds and a period of 389 days. However, random rearrangements of the data provide similar quality fits, and we cannot with certainty ascribe the timing variations to mutual exoplanet interactions. This inconclusive result is consistent with the lack of incontrovertible transit timing variations (TTVs) observed in other hot-Jupiter systems. This outcome emphasizes that unique recognition of TTVs requires dense sampling of the libration cycle (e.g. continuous observations from space-based platforms). However, even in systems observed with the Kepler spacecraft, there is a noted lack of transiting companions and TTVs in hot-Jupiter systems. This result is more meaningful, and indicates that hot-Jupiter systems, while they are easily observable from the ground, do not appear to be currently configured in a manner favorable to the detection of TTVs. The future of ground-based TTV studies may reside in resolving secular trends, and/or implementation at extreme quality observing sites to minimize atmospheric red noise., Comment: ApJ Letters, accepted

The Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) observed eleven transits of TrES-3b over two years in order to constrain system parameters and look for transit timing and depth variations. We describe an updated analysis protocol for APOSTLE data, including the reduction pipeline, transit model and Markov Chain Monte Carlo analyzer. Our estimates of the system parameters for TrESb are consistent with previous estimates to within the 2\sigma\ confidence level. We improved the errors (by 10--30%) on system parameters like the orbital inclination ($i_{\text{orb}}$), impact parameter (b) and stellar density (\rho$_{\star}$) compared to previous measurements. The near-grazing nature of the system, and incomplete sampling of some transits, limited our ability to place reliable uncertainties on individual transit depths and hence we do not report strong evidence for variability. Our analysis of the transit timing data show no evidence for transit timing variations and our timing measurements are able to rule out Super-Earth and Gas Giant companions in low order mean motion resonance with TrES-3b., Comment: 30 pages, 8 tables, 5 figures (submitted to ApJ, Sep 2012)

We investigate the use of optical photometric variability to select and identify blazars in large-scale time-domain surveys, in part to aid in the identification of blazar counterparts to the ~30% of gamma-ray sources in the Fermi 2FGL catalog still lacking reliable associations. Using data from the optical LINEAR asteroid survey, we characterize the optical variability of blazars by fitting a damped random walk model to individual light curves with two main model parameters, the characteristic timescales of variability (tau), and driving amplitudes on short timescales (sigma). Imposing cuts on minimum tau and sigma allows for blazar selection with high efficiency E and completeness C. To test the efficacy of this approach, we apply this method to optically variable LINEAR objects that fall within the several-arcminute error ellipses of gamma-ray sources in the Fermi 2FGL catalog. Despite the extreme stellar contamination at the shallow depth of the LINEAR survey, we are able to recover previously-associated optical counterparts to Fermi AGN with E > 88% and C = 88% in Fermi 95% confidence error ellipses having semimajor axis r < 8'. We find that the suggested radio counterpart to Fermi source 2FGL J1649.6+5238 has optical variability consistent with other gamma-ray blazars, and is likely to be the gamma-ray source. Our results suggest that the variability of the non-thermal jet emission in blazars is stochastic in nature, with unique variability properties due to the effects of relativistic beaming. After correcting for beaming, we estimate that the characteristic timescale of blazar variability is ~3 years in the rest-frame of the jet, in contrast with the ~320 day disk flux timescale observed in quasars. The variability-based selection method presented will be useful for blazar identification in time-domain optical surveys, and is also a probe of jet physics., Comment: version resubmitted to ApJ with incorporated referee's comments. 12 pages, 10 figures

We present follow-up observations and analysis of the recently discovered short period low-mass eclipsing binary, SDSS J001641-000925. With an orbital period of 0.19856 days, this system has one of the shortest known periods for an M dwarf binary system. Medium-resolution spectroscopy and multi-band photometry for the system are presented. Markov chain Monte Carlo modeling of the light curves and radial velocities yields estimated masses for the stars of M1 = 0.54 +/- 0.07 Msun and M2 = 0.34 +/- 0.04 Msun, and radii of R1 = 0.68 +/- 0.03 Rsun and R2 = 0.58 +/- 0.03 Rsun respectively. This solution places both components above the critical Roche overfill limit, providing strong evidence that SDSS J001641-000925 is the first verified M-dwarf contact binary system. Within the follow-up spectroscopy we find signatures of non-solid body rotation velocities, which we interpret as evidence for mass transfer or loss within the system. In addition, our photometry samples the system over 9 years, and we find strong evidence for period decay at the rate of dP/dt ~8 s/yr. Both of these signatures raise the intriguing possibility that the system is in over-contact, and actively losing angular momentum, likely through mass loss. This places SDSS J001641-000925 as not just the first M-dwarf over-contact binary, but one of the few systems of any spectral type known to be actively undergoing coalescence. Further study SDSS J001641-000925 is on-going to verify the nature of the system, which may prove to be a unique astrophysical laboratory., Comment: 11 figures, ApJ Accepted

We present the ensemble properties of 31 comets (27 resolved and 4 unresolved) observed by the Sloan Digital Sky Survey (SDSS). This sample of comets represents about 1 comet per 10 million SDSS photometric objects. Five-band (u,g,r,i,z) photometry is used to determine the comets' colors, sizes, surface brightness profiles, and rates of dust production in terms of the Af{\rho} formalism. We find that the cumulative luminosity function for the Jupiter Family Comets in our sample is well fit by a power law of the form N(< H) \propto 10(0.49\pm0.05)H for H < 18, with evidence of a much shallower fit N(< H) \propto 10(0.19\pm0.03)H for the faint (14.5 < H < 18) comets. The resolved comets show an extremely narrow distribution of colors (0.57 \pm 0.05 in g - r for example), which are statistically indistinguishable from that of the Jupiter Trojans. Further, there is no evidence of correlation between color and physical, dynamical, or observational parameters for the observed comets., Comment: 19 pages, 8 tables, 11 figures, to appear in Icarus

We present the first rates of flares from M dwarf stars in both red optical and near infrared (NIR) filters. We have studied ~50,000 M dwarfs from the SDSS Stripe 82 area, and 1,321 M dwarfs from the 2MASS Calibration Scan Point Source Working Database that overlap SDSS imaging fields. We assign photometric spectral types from M0 to M6 using (r-i) and (i-z) colors for every star in our sample. Stripe 82 stars each have 50-100 epochs of data, while 2MASS Calibration stars have ~1900 epochs. From these data we estimate the observed rates and theoretical detection thresholds for flares in eight photometric bands as a function of spectral type. Optical flare rates are found to be in agreement with previous studies, while the frequency per hour of NIR flare detections is found to be more than two orders of magnitude lower. An excess of small amplitude flux increases in all bands exhibits a power-law distribution, which we interpret as the result of flares below our detection thresholds. In order to investigate the recovery efficiency for flares in each filter, we extend a two-component flare model into the NIR. Quiescent M0-M6 spectral templates were used with the model to predict the photometric response of flares from u to Ks. We determine that red optical filters are sensitive to flares with u-band amplitudes >2 mag, and NIR filters to flares with delta u>4.5 mag. Our model predicts that M0 stars have the best color-contrast for J-band detections, but M4-M6 stars should show the highest rate of NIR flares with amplitudes of delta J ~0.01 mag. Characterizing flare rates and photometric variations at longer wavelengths is important for predicting the signatures of M dwarf variability in next-generation surveys, and we discuss their impact on surveys such as LSST., Comment: ApJ in Press. 12 Pages, 17 Figures

We provide a quantitative description and statistical interpretation of the optical continuum variability of quasars. The Sloan Digital Sky Survey (SDSS) has obtained repeated imaging in five UV-to-IR photometric bands for 33,881 spectroscopically confirmed quasars. About 10,000 quasars have an average of 60 observations in each band obtained over a decade along Stripe 82 (S82), whereas the remaining ~25,000 have 2-3 observations due to scan overlaps. The observed time lags span the range from a day to almost 10 years, and constrain quasar variability at rest-frame time lags of up to 4 years, and at rest-frame wavelengths from 1000A to 6000A. We publicly release a user-friendly catalog of quasars from the SDSS Data Release 7 that have been observed at least twice in SDSS or once in both SDSS and the Palomar Observatory Sky Survey, and we use it to analyze the ensemble properties of quasar variability. Based on a damped random walk (DRW) model defined by a characteristic time scale and an asymptotic variability amplitude that scale with the luminosity, black hole mass, and rest wavelength for individual quasars calibrated in S82, we can fully explain the ensemble variability statistics of the non-S82 quasars such as the exponential distribution of large magnitude changes. All available data are consistent with the DRW model as a viable description of the optical continuum variability of quasars on time scales of ~5-2000 days in the rest frame. We use these models to predict the incidence of quasar contamination in transient surveys such as those from PTF and LSST., Comment: 33 pages, 19 figures, replaced with accepted version. Catalog is available at http://www.astro.washington.edu/users/ivezic/macleod/qso_dr7/

We present details of the construction and characterization of the coaddition of the Sloan Digital Sky Survey Stripe 82 \ugriz\ imaging data. This survey consists of 275 deg$^2$ of repeated scanning by the SDSS camera of $2.5\arcdeg$ of $\delta$ over $-50\arcdeg \le \alpha \le 60\arcdeg$ centered on the Celestial Equator. Each piece of sky has $\sim 20$ runs contributing and thus reaches $\sim2$ magnitudes fainter than the SDSS single pass data, i.e. to $r\sim 23.5$ for galaxies. We discuss the image processing of the coaddition, the modeling of the PSF, the calibration, and the production of standard SDSS catalogs. The data have $r$-band median seeing of 1.1\arcsec, and are calibrated to $\le 1%$. Star color-color, number counts, and psf size vs modelled size plots show the modelling of the PSF is good enough for precision 5-band photometry. Structure in the psf-model vs magnitude plot show minor psf mis-modelling that leads to a region where stars are being mis-classified as galaxies, and this is verified using VVDS spectroscopy. As this is a wide area deep survey there are a variety of uses for the data, including galactic structure, photometric redshift computation, cluster finding and cross wavelength measurements, weak lensing cluster mass calibrations, and cosmic shear measurements., Comment: 18 page, 10 figures, submitted to ApJ. Small changes in text to be consistent with revised photo-z catalog of Reis et al. arXiv:1111.6620v2

We describe photometric recalibration of data obtained by the asteroid survey LINEAR. Although LINEAR was designed for astrometric discovery of moving objects, the dataset described here contains over 5 billion photometric measurements for about 25 million objects, mostly stars. We use SDSS data from the overlapping ~10,000 deg^2 of sky to recalibrate LINEAR photometry, and achieve errors of 0.03 mag for sources not limited by photon statistics, with errors of 0.2 mag at r~18. With its 200 observations per object on average, LINEAR data provide time domain information for the brightest 4 magnitudes of SDSS survey. At the same time, LINEAR extends the deepest similar wide-area variability survey, the Northern Sky Variability Survey, by 3 mag. We briefly discuss the properties of about 7,000 visually confirmed periodic variables, dominated by roughly equal fractions of RR Lyrae stars and eclipsing binary stars, and analyze their distribution in optical and infra-red color-color diagrams. The LINEAR dataset is publicly available from the SkyDOT website (http://skydot.lanl.gov)., Comment: Corrected some typos to match the published version, AJ in press

Time-resolved databases with large spatial coverage are quickly becoming a standard tool for all types of astronomical studies. We report preliminary results from our search for stellar flares in the 2MASS calibration fields. A sample of 4343 M dwarfs, spatially matched between the SDSS and the 2MASS calibration fields, each with hundreds to thousands of epochs in near infrared bandpasses, is analyzed using a modified Welch-Stetson index to characterize the variability. A Monte Carlo model was used to assess the noise of the variability index. We find significnat residuals above the noise with power-law slopes of -3.37 and -4.05 for our JH and HKs distributions respectively. This is evidence for flares being observed from M dwarfs in infrared photometry., Comment: 7 Pages, 3 Figures, Submitted to the 16th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy., Comment: 596 pages. Also available at full resolution at http://www.lsst.org/lsst/scibook

We present an improved analysis of halo substructure traced by RR Lyrae stars in the SDSS stripe 82 region. With the addition of SDSS-II data, a revised selection method based on new ugriz light curve templates results in a sample of 483 RR Lyrae stars that is essentially free of contamination. The main result from our first study persists: the spatial distribution of halo stars at galactocentric distances 5--100 kpc is highly inhomogeneous. At least 20% of halo stars within 30 kpc from the Galactic center can be statistically associated with substructure. We present strong direct evidence, based on both RR Lyrae stars and main sequence stars, that the halo stellar number density profile significantly steepens beyond a Galactocentric distance of ~30 kpc, and a larger fraction of the stars are associated with substructure. By using a novel method that simultaneously combines data for RR Lyrae and main sequence stars, and using photometric metallicity estimates for main sequence stars derived from deep co-added u-band data, we measure the metallicity of the Sagittarius dSph tidal stream (trailing arm) towards R.A.2h-3h and Dec~0 deg to be 0.3 dex higher ([Fe/H]=-1.2) than that of surrounding halo field stars. Together with a similar result for another major halo substructure, the Monoceros stream, these results support theoretical predictions that an early forming, smooth inner halo, is metal poor compared to high surface brightness material that have been accreted onto a later-forming outer halo. The mean metallicity of stars in the outer halo that are not associated with detectable clumps may still be more metal-poor than the bulk of inner-halo stars, as has been argued from other data sets., Comment: Submitted to ApJ, 68 pages, 26 figures, supplemental material (light curves, templates, animation) can be downloaded from http://www.astro.washington.edu/bsesar/S82_RRLyr.html

The next decade of survey astronomy has the potential to transform our knowledge of variable stars. Stellar variability underpins our knowledge of the cosmological distance ladder, and provides direct tests of stellar formation and evolution theory. Variable stars can also be used to probe the fundamental physics of gravity and degenerate material in ways that are otherwise impossible in the laboratory. The computational and engineering advances of the past decade have made large-scale, time-domain surveys an immediate reality. Some surveys proposed for the next decade promise to gather more data than in the prior cumulative history of astronomy. The actual implementation of these surveys will have broad implications for the types of science that will be enabled. We examine the design considerations for an optimal time-domain photometric survey dedicated to variable star science, including: observing cadence, wavelength coverage, photometric and astrometric accuracy, single-epoch and cumulative depth, overall sky coverage, and data access by the broader astronomical community. The best surveys must combine aspects from each of these considerations to fully realize the potential for the next decade of time-domain science., Comment: pages; white paper submitted to the Astro2010 Decadal Survey Science Frontier Panel on Stars and Stellar Evolution

We report the discovery of a minor planet (2006 SQ372) on an orbit with a perihelion of 24 AU and a semimajor axis of 796 AU. Dynamical simulations show that this is a transient orbit and is unstable on a timescale of 200 Myrs. Falling near the upper semimajor axis range of the scattered disk and the lower semimajor axis range of the Oort Cloud, previous membership in either class is possible. By modeling the production of similar orbits from the Oort Cloud as well as from the scattered disk, we find that the Oort Cloud produces 16 times as many objects on SQ372-like orbits as the scattered disk. Given this result, we believe this to be the most distant long-period comet ever discovered. Furthermore, our simulation results also indicate that 2000 OO67 has had a similar dynamical history. Unaffected by the "Jupiter-Saturn Barrier," these two objects are most likely long-period comets from the inner Oort Cloud.

Gliese 436 is an M dwarf with a mass of 0.45 Msun and hosts the extrasolar planet GL 436b [3, 6, 7, 2], which is currently the least massive transiting planet with a mass of ~23.17 Mearth [10], and the only planet known to transit an M dwarf. GL 436b represents the first transiting detection of the class of extrasolar planets known as "Hot Neptunes" that have masses within a few times that of Neptune's mass (~17 Mearth) and orbital semimajor axis <0.1 AU about the host star. Unlike most other known transiting extrasolar planets, GL 436b has a high eccentricity (e~0.16). This brings to light a new parameter space for habitability zones of extrasolar planets with host star masses much smaller than typical stars of roughly a solar mass. This unique system is an ideal candidate for orbital perturbation and transit-time variation (TTV) studies to detect smaller, possibly Earth-mass planets in the system. In April 2008 we began a long-term intensive campaign to obtain complete high-precision light curves using the Apache Point Observatory's 3.5-meter telescope, NMSU's 1-meter telescope (located at APO), and Sommers Bausch Observatory's 24" telescope. These light curves are being analyzed together, along with amateur and other professional astronomer observations. Results of our analysis are discussed. Continued measurements over the next few years are needed to determine if additional planets reside in the system, and to study the impact of other manifestations on the light curves, such as star spots and active regions., Comment: 4 pages, 3 figures. To appear in "Proceedings of the 15th Cambridge Workshop on Cool Stars, Stellar Systems and the Sun", 2009, AIP Conference Proceedings vol. 1094, ed. Eric Stempels

We present ground-based observations of the transiting Neptune-mass planet Gl 436b obtained with the 3.5-meter telescope at Apache Point Observatory and other supporting telescopes. Included in this is an observed transit in early 2005, over two years before the earliest reported transit detection. We have compiled all available transit data to date and perform a uniform modeling using the JKTEBOP code. We do not detect any transit timing variations of amplitude greater than ~1 minute over the ~3.3 year baseline. We do however find possible evidence for a self-consistent trend of increasing orbital inclination, transit width, and transit depth, which supports the supposition that Gl 436b is being perturbed by another planet of < 12 Earth Masses in a non-resonant orbit., Comment: 6 pages, 3 figures, 2 tables, Accepted to ApJ Letters

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world., Comment: 57 pages, 32 color figures, version with high-resolution figures available from https://www.lsst.org/overview

I provide an incomplete inventory of the astronomical variability that will be found by next-generation time-domain astronomical surveys. These phenomena span the distance range from near-Earth satellites to the farthest Gamma Ray Bursts. The surveys that detect these transients will issue alerts to the greater astronomical community; this decision process must be extremely robust to avoid a slew of ``false'' alerts, and to maintain the community's trust in the surveys. I review the functionality required of both the surveys and the telescope networks that will be following them up, and the role of VOEvents in this process. Finally, I offer some ideas about object and event classification, which will be explored more thoroughly by other articles in these proceedings., Comment: 4 pages. In Proceedings of Hot-wiring the Transient Universe (HTU) 2007. To be published in Astronomische Nachrichten, March 2008

We have developed an end-to-end photometric data processing pipeline to compare current photometric algorithms commonly used on ground-based imaging data. This testbed is exceedingly adaptable, and enables us to perform many research and development tasks, including image subtraction and co-addition, object detection and measurements, the production of photometric catalogs, and the creation and stocking of database tables with time-series information. This testing has been undertaken to evaluate existing photometry algorithms for consideration by a next-generation image processing pipeline for the Large Synoptic Survey Telescope (LSST). We outline the results of our tests for four packages: The Sloan Digital Sky Survey's (SDSS) Photo package, Daophot and Allframe, DoPhot, and two versions of Source Extractor (SExtractor). The ability of these algorithms to perform point-source photometry, astrometry, shape measurements, star-galaxy separation, and to measure objects at low signal-to-noise is quantified. We also perform a detailed crowded field comparison of Daophot and Allframe, and profile the speed and memory requirements in detail for SExtractor. We find that both Daophot and Photo are able to perform aperture photometry to high enough precision to meet LSST's science requirements, and less adequately at PSF-fitting photometry. Photo performs the best at simultaneous point and extended-source shape and brightness measurements. SExtractor is the fastest algorithm, and recent upgrades in the software yield high-quality centroid and shape measurements with little bias towards faint magnitudes. Allframe yields the best photometric results in crowded fields., Comment: Accepted for publication in PASP. Appendix is included only here in this eprint. 11 figures, 18 tables

AM CVn systems are a rare (about a dozen previously known) class of cataclysmic variables, arguably encompassing the shortest orbital periods (down to about 10 minutes) of any known binaries. Both binary components are thought to be degenerate (or partially so), likely with mass-transfer from a helium-rich donor onto a white dwarf, driven by gravitational radiation. Although rare, AM CVn systems are of high interest as possible SN Ia progenitors, and because they are predicted to be common sources of gravity waves in upcoming experiments such as LISA. We have identified four new AM CVn candidates from the Sloan Digital Sky Survey (SDSS) spectral database. All four show hallmark spectroscopic characteristics of the AM CVn class: each is devoid of hydrogen features, and instead shows a spectrum dominated by helium. All four show double-peaked emission, indicative of helium-dominated accretion disks. Limited time-series CCD photometric follow-on data have been obtained for three of the new candidates from the ARC 3.5m; most notably, a 28.3 minute binary period with sharp, deep eclipses is discovered in one case, SDSS J0926+3624. This is the first confirmed eclipsing AM CVn, and our data allow initial estimates of binary parameters for this ultracompact system. The four new SDSS objects also provide a substantial expansion of the currently critically-small sample of AM CVn systems., Comment: 7 pages, 3 figures, submitted to the Astronomical Journal

We present previously unpublished photometry of three Large Magellanic Cloud (LMC) microlensing events and show that the new photometry confirms the microlensing interpretation of these events. These events were discovered by the MACHO Project alert system and were also recovered by the analysis of the 5.7 year MACHO data set. This new photometry provides a substantial increase in the signal-to-noise ratio over the previously published photometry and in all three cases, the gravitational microlensing interpretation of these events is strengthened. The new data consist of MACHO-Global Microlensing Alert Network (GMAN) follow-up images from the CTIO 0.9 telescope plus difference imaging photometry of the original MACHO data from the 1.3m "Great Melbourne" telescope at Mt. Stromlo. We also combine microlensing light curve fitting with photometry from high resolution HST images of the source stars to provide further confirmation of these events and to show that the microlensing interpretation of event MACHO-LMC-23 is questionable. Finally, we compare our results with the analysis of Belokurov, Evans & Le Du who have attempted to classify candidate microlensing events with a neural network method, and we find that their results are contradicted by the new data and more powerful light curve fitting analysis for each of the four events considered in this paper. The failure of the Belokurov, Evans & Le Du method is likely to be due to their use of a set of insensitive statistics to feed their neural networks., Comment: 29 pages with 8 included postscript figures, accepted by the Astrophysical Journal