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

We present ultraviolet (UV) spectroscopy and photometry of four Type Ia supernovae (SNe 2004dt, 2004ef, 2005M, and 2005cf) obtained with the UV prism of the Advanced Camera for Surveys on the Hubble Space Telescope. This dataset provides unique spectral time series down to 2000 Angstrom. Significant diversity is seen in the near maximum-light spectra (~ 2000--3500 Angstrom) for this small sample. The corresponding photometric data, together with archival data from Swift Ultraviolet/Optical Telescope observations, provide further evidence of increased dispersion in the UV emission with respect to the optical. The peak luminosities measured in uvw1/F250W are found to correlate with the B-band light-curve shape parameter dm15(B), but with much larger scatter relative to the correlation in the broad-band B band (e.g., ~0.4 mag versus ~0.2 mag for those with 0.8 < dm15 < 1.7 mag). SN 2004dt is found as an outlier of this correlation (at > 3 sigma), being brighter than normal SNe Ia such as SN 2005cf by ~0.9 mag and ~2.0 mag in the uvw1/F250W and uvm2/F220W filters, respectively. We show that different progenitor metallicity or line-expansion velocities alone cannot explain such a large discrepancy. Viewing-angle effects, such as due to an asymmetric explosion, may have a significant influence on the flux emitted in the UV region. Detailed modeling is needed to disentangle and quantify the above effects.

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.

We report the discovery of a supernova (SN) with the highest apparent energy output to date and conclude that it represents an extreme example of the Type IIn subclass. The SN, which was discovered behind the Large Magellanic Cloud at z = 0.289 by the SuperMACHO microlensing survey, peaked at M_R = -21.5 mag and only declined by 2.9 mag over 4.7 years after the peak. Over this period, SN 2003ma had an integrated bolometric luminosity of 4 x 10^51 ergs, more than any other SN to date. The radiated energy is close to the limit allowed by conventional core-collapse explosions. Optical spectra reveal that SN 2003ma has persistent single-peaked intermediate-width hydrogen lines, a signature of interaction between the SN and a dense circumstellar medium. The light curves show further evidence for circumstellar interaction, including a long plateau with a shape very similar to the classic SN IIn 1988Z -- however, SN 2003ma is ten times more luminous at all epochs. The fast velocity measured for the intermediate-width H_alpha component (~6000 km/s) points towards an extremely energetic explosion (> 10^52 ergs), which imparts a faster blast-wave speed to the post-shock material and a higher luminosity from the interaction than is observed in typical SNe IIn. Mid-infrared observations of SN 2003ma suggest an infrared light echo is produced by normal interstellar dust at a distance ~0.5 pc from the SN.

We present 2323 High-Amplitude \delta-Scuti (HADS) candidates discovered in the Large Magellanic Cloud (LMC) by the SuperMACHO survey (Rest et al. 2005). Frequency analyses of these candidates reveal that several are multimode pulsators, including 119 whose largest amplitude of pulsation is in the fundamental (F) mode and 19 whose largest amplitude of pulsation is in the first overtone (FO) mode. Using Fourier decomposition of the HADS light curves, we find that the period-luminosity (PL) relation defined by the FO pulsators does not show a clear separation from the PL-relation defined by the F pulsators. This differs from other instability strip pulsators such as type c RR Lyrae. We also present evidence for a larger amplitude, subluminous population of HADS similar to that observed in Fornax (Poretti et al. 2008).

We present ugriz light curves for 146 spectroscopically confirmed or spectroscopically probable Type Ia supernovae from the 2005 season of the SDSS-II Supernova survey. The light curves have been constructed using a photometric technique that we call scene modelling, which is described in detail here; the major feature is that supernova brightnesses are extracted from a stack of images without spatial resampling or convolution of the image data. This procedure produces accurate photometry along with accurate estimates of the statistical uncertainty, and can be used to derive photometry taken with multiple telescopes. We discuss various tests of this technique that demonstrate its capabilities. We also describe the methodology used for the calibration of the photometry, and present calibrated magnitudes and fluxes for all of the spectroscopic SNe Ia from the 2005 season.

We present measurements of the Hubble diagram for 103 Type Ia supernovae (SNe) with redshifts 0.04 < z < 0.42, discovered during the first season (Fall 2005) of the Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey. These data fill in the redshift "desert" between low- and high-redshift SN Ia surveys. We combine the SDSS-II measurements with new distance estimates for published SN data from the ESSENCE survey, the Supernova Legacy Survey, the Hubble Space Telescope, and a compilation of nearby SN Ia measurements. Combining the SN Hubble diagram with measurements of Baryon Acoustic Oscillations from the SDSS Luminous Red Galaxy sample and with CMB temperature anisotropy measurements from WMAP, we estimate the cosmological parameters w and Omega_M, assuming a spatially flat cosmological model (FwCDM) with constant dark energy equation of state parameter, w. For the FwCDM model and the combined sample of 288 SNe Ia, we find w = -0.76 +- 0.07(stat) +- 0.11(syst), Omega_M = 0.306 +- 0.019(stat) +- 0.023(syst) using MLCS2k2 and w = -0.96 +- 0.06(stat) +- 0.12(syst), Omega_M = 0.265 +- 0.016(stat) +- 0.025(syst) using the SALT-II fitter. We trace the discrepancy between these results to a difference in the rest-frame UV model combined with a different luminosity correction from color variations; these differences mostly affect the distance estimates for the SNLS and HST supernovae. We present detailed discussions of systematic errors for both light-curve methods and find that they both show data-model discrepancies in rest-frame $U$-band. For the SALT-II approach, we also see strong evidence for redshift-dependence of the color-luminosity parameter (beta). Restricting the analysis to the 136 SNe Ia in the Nearby+SDSS-II samples, we find much better agreement between the two analysis methods but with larger uncertainties.

We describe a general analysis package for supernova (SN) light curves, called SNANA, that contains a simulation, light curve fitter, and cosmology fitter. The software is designed with the primary goal of using SNe Ia as distance indicators for the determination of cosmological parameters, but it can also be used to study efficiencies for analyses of SN rates, estimate contamination from non-Ia SNe, and optimize future surveys. Several SN models are available within the same software architecture, allowing technical features such as K-corrections to be consistently used among multiple models, and thus making it easier to make detailed comparisons between models. New and improved light-curve models can be easily added. The software works with arbitrary surveys and telescopes and has already been used by several collaborations, leading to more robust and easy-to-use code. This software is not intended as a final product release, but rather it is designed to undergo continual improvements from the community as more is learned about SNe. Below we give an overview of the SNANA capabilities, as well as some of its limitations. Interested users can find software downloads and more detailed information from the manuals at http://www.sdss.org/supernova/SNANA.html .

We present a measurement of the rate of type Ia supernovae (SNe Ia) from the first of three seasons of data from the SDSS-II Supernova Survey. For this measurement, we include 17 SNe Ia at redshift $z\le0.12$. Assuming a flat cosmology with $\Omega_m = 0.3=1-\Omega_\Lambda$, we find a volumetric SN Ia rate of $[2.93^{+0.17}_{-0.04}({\rm systematic})^{+0.90}_{-0.71}({\rm statistical})] \times 10^{-5} {\rm SNe} {\rm Mpc}^{-3} h_{70}^3 {\rm year}^{-1}$, at a volume-weighted mean redshift of 0.09. This result is consistent with previous measurements of the SN Ia rate in a similar redshift range. The systematic errors are well controlled, resulting in the most precise measurement of the SN Ia rate in this redshift range. We use a maximum likelihood method to fit SN rate models to the SDSS-II Supernova Survey data in combination with other rate measurements, thereby constraining models for the redshift-evolution of the SN Ia rate. Fitting the combined data to a simple power-law evolution of the volumetric SN Ia rate, $r_V \propto (1+z)^{\beta}$, we obtain a value of $\beta = 1.5 \pm 0.6$, i.e. the SN Ia rate is determined to be an increasing function of redshift at the $\sim 2.5 \sigma$ level. Fitting the results to a model in which the volumetric SN rate, $r_V=A\rho(t)+B\dot \rho(t)$, where $\rho(t)$ is the stellar mass density and $\dot \rho(t)$ is the star formation rate, we find $A = (2.8 \pm 1.2) \times 10^{-14} \mathrm{SNe} \mathrm{M}_{\sun}^{-1} \mathrm{year}^{-1}$, $B = (9.3^{+3.4}_{-3.1})\times 10^{-4} \mathrm{SNe} \mathrm{M}_{\sun}^{-1}$.

We report the discovery and orbit determination of 14 trans-Neptunian objects (TNOs) from the ESSENCE Supernova Survey difference imaging dataset. Two additional objects discovered in a similar search of the SDSS-II Supernova Survey database were recovered in this effort. ESSENCE repeatedly observed fields far from the Solar System ecliptic (-21 deg < beta < -5 deg), reaching limiting magnitudes per observation of I~23.1 and R~23.7. We examine several of the newly detected objects in detail, including 2003 UC_414 which orbits entirely between Uranus and Neptune and lies very close to a dynamical region that would make it stable for the lifetime of the Solar System. 2003 SS_422 and 2007 TA_418 have high eccentricities and large perihelia, making them candidate members of an outer class of trans-Neptunian objects. We also report a new member of the ''extended'' or ''detached'' scattered disk, 2004 VN_112, and verify the stability of its orbit using numerical simulations. This object would have been visible to ESSENCE for only ~2% of its orbit, suggesting a vast number of similar objects across the sky. We emphasize that off-ecliptic surveys are optimal for uncovering the diversity of such objects, which in turn will constrain the history of gravitational influences that shaped our early Solar System.

Using effective temperature and metallicity derived from SDSS spectra for ~60,000 F and G type main sequence stars (0.2<g-r<0.6), we develop polynomial models for estimating these parameters from the SDSS u-g and g-r colors. We apply this method to SDSS photometric data for about 2 million F/G stars and measure the unbiased metallicity distribution for a complete volume-limited sample of stars at distances between 500 pc and 8 kpc. The metallicity distribution can be exquisitely modeled using two components with a spatially varying number ratio, that correspond to disk and halo. The two components also possess the kinematics expected for disk and halo stars. The metallicity of the halo component is spatially invariant, while the median disk metallicity smoothly decreases with distance from the Galactic plane from -0.6 at 500 pc to -0.8 beyond several kpc. The absence of a correlation between metallicity and kinematics for disk stars is in a conflict with the traditional decomposition in terms of thin and thick disks. We detect coherent substructures in the kinematics--metallicity space, such as the Monoceros stream, which rotates faster than the LSR, and has a median metallicity of [Fe/H]=-0.96, with an rms scatter of only ~0.15 dex. We extrapolate our results to the performance expected from the Large Synoptic Survey Telescope (LSST) and estimate that the LSST will obtain metallicity measurements accurate to 0.2 dex or better, with proper motion measurements accurate to ~0.2 mas/yr, for about 200 million F/G dwarf stars within a distance limit of ~100 kpc (g<23.5). [abridged]

We present multiepoch spectra of 13 high-redshift Type Ia supernovae (SNe Ia) drawn from the literature, the ESSENCE and SNLS projects, and our own separate dedicated program on the ESO Very Large Telescope. We use the Supernova Identification (SNID) code of Blondin & Tonry to determine the spectral ages in the supernova rest frame. Comparison with the observed elapsed time yields an apparent aging rate consistent with the 1/(1+z) factor (where z is the redshift) expected in a homogeneous, isotropic, expanding universe. These measurements thus confirm the expansion hypothesis, while unambiguously excluding models that predict no time dilation, such as Zwicky's "tired light" hypothesis. We also test for power-law dependencies of the aging rate on redshift. The best-fit exponent for these models is consistent with the expected 1/(1+z) factor.

This paper presents spectroscopy of supernovae discovered in the first season of the Sloan Digital Sky Survey-II Supernova Survey. This program searches for and measures multi-band light curves of supernovae in the redshift range z = 0.05 - 0.4, complementing existing surveys at lower and higher redshifts. Our goal is to better characterize the supernova population, with a particular focus on SNe Ia, improving their utility as cosmological distance indicators and as probes of dark energy. Our supernova spectroscopy program features rapid-response observations using telescopes of a range of apertures, and provides confirmation of the supernova and host-galaxy types as well as precise redshifts. We describe here the target identification and prioritization, data reduction, redshift measurement, and classification of 129 SNe Ia, 16 spectroscopically probable SNe Ia, 7 SNe Ib/c, and 11 SNe II from the first season. We also describe our efforts to measure and remove the substantial host galaxy contamination existing in the majority of our SN spectra.

We present a new public archive of light-motion curves in Sloan Digital Sky Survey (SDSS) Stripe 82, covering 99 deg in right ascension from RA = 20.7 h to 3.3 h and spanning 2.52 deg in declination from Dec = -1.26 to 1.26 deg, for a total sky area of ~249 sq deg. Stripe 82 has been repeatedly monitored in the u, g, r, i and z bands over a seven-year baseline. Objects are cross-matched between runs, taking into account the effects of any proper motion. The resulting catalogue contains almost 4 million light-motion curves of stellar objects and galaxies. The photometry are recalibrated to correct for varying photometric zeropoints, achieving ~20 mmag and ~30 mmag root-mean-square (RMS) accuracy down to 18 mag in the g, r, i and z bands for point sources and extended sources, respectively. The astrometry are recalibrated to correct for inherent systematic errors in the SDSS astrometric solutions, achieving ~32 mas and ~35 mas RMS accuracy down to 18 mag for point sources and extended sources, respectively. For each light-motion curve, 229 photometric and astrometric quantities are derived and stored in a higher-level catalogue. On the photometric side, these include mean exponential and PSF magnitudes along with uncertainties, RMS scatter, chi^2 per degree of freedom, various magnitude distribution percentiles, object type (stellar or galaxy), and eclipse, Stetson and Vidrih variability indices. On the astrometric side, these quantities include mean positions, proper motions as well as their uncertainties and chi^2 per degree of freedom. The here presented light-motion curve catalogue is complete down to r~21.5 and is at present the deepest large-area photometric and astrometric variability catalogue available.

A 2.5 x 100 degree region along the celestial equator (Stripe 82) has been imaged repeatedly from 1998 to 2005 by the Sloan Digital Sky Survey. A new catalogue of ~4 million light-motion curves, together with over 200 derived statistical quantities, for objects in Stripe 82 brighter than r~21.5 has been constructed by combining these data by Bramich et al. (2007). This catalogue is at present the deepest catalogue of its kind. Extracting the ~130000 objects with highest signal-to-noise ratio proper motions, we build a reduced proper motion diagram to illustrate the scientific promise of the catalogue. In this diagram disk and halo subdwarfs are well-separated from the cool white dwarf sequence. Our sample of 1049 cool white dwarf candidates includes at least 8 and possibly 21 new ultracool H-rich white dwarfs (T_eff < 4000K) and one new ultracool He-rich white dwarf candidate identified from their SDSS optical and UKIDSS infrared photometry. At least 10 new halo white dwarfs are also identified from their kinematics.

The Sloan Digital Sky Survey-II Supernova Survey has identified a large number of new transient sources in a 300 sq. deg. region along the celestial equator during its first two seasons of a three-season campaign. Multi-band (ugriz) light curves were measured for most of the sources, which include solar system objects, Galactic variable stars, active galactic nuclei, supernovae (SNe), and other astronomical transients. The imaging survey is augmented by an extensive spectroscopic follow-up program to identify SNe, measure their redshifts, and study the physical conditions of the explosions and their environment through spectroscopic diagnostics. During the survey, light curves are rapidly evaluated to provide an initial photometric type of the SNe, and a selected sample of sources are targeted for spectroscopic observations. In the first two seasons, 476 sources were selected for spectroscopic observations, of which 403 were identified as SNe. For the Type Ia SNe, the main driver for the Survey, our photometric typing and targeting efficiency is 90%. Only 6% of the photometric SN Ia candidates were spectroscopically classified as non-SN Ia instead, and the remaining 4% resulted in low signal-to-noise, unclassified spectra. This paper describes the search algorithm and the software, and the real-time processing of the SDSS imaging data. We also present the details of the supernova candidate selection procedures and strategies for follow-up spectroscopic and imaging observations of the discovered sources.

We present the first large-scale effort of creating composite spectra of high-redshift type Ia supernovae (SNe Ia) and comparing them to low-redshift counterparts. Through the ESSENCE project, we have obtained 107 spectra of 88 high-redshift SNe Ia with excellent light-curve information. In addition, we have obtained 397 spectra of low-redshift SNe through a multiple-decade effort at Lick and Keck Observatories, and we have used 45 UV spectra obtained by HST/IUE. The low-redshift spectra act as a control sample when comparing to the ESSENCE spectra. In all instances, the ESSENCE and Lick composite spectra appear very similar. The addition of galaxy light to the Lick composite spectra allows a nearly perfect match of the overall spectral-energy distribution with the ESSENCE composite spectra, indicating that the high-redshift SNe are more contaminated with host-galaxy light than their low-redshift counterparts. This is caused by observing objects at all redshifts with the same slit width, which corresponds to different projected distances. After correcting for the galaxy-light contamination, subtle differences in the spectra remain. We have estimated the systematic errors when using current spectral templates for K-corrections to be ~0.02 mag. The variance in the composite spectra give an estimate of the intrinsic variance in low-redshift maximum-light SN spectra of ~3% in the optical and growing toward the UV. The difference between the maximum light low and high-redshift spectra constrain SN evolution between our samples to be < 10% in the rest-frame optical.

The Sloan Digital Sky Survey-II (SDSS-II) has embarked on a multi-year project to identify and measure light curves for intermediate-redshift (0.05 < z < 0.35) Type Ia supernovae (SNe Ia) using repeated five-band (ugriz) imaging over an area of 300 sq. deg. The survey region is a stripe 2.5 degrees wide centered on the celestial equator in the Southern Galactic Cap that has been imaged numerous times in earlier years, enabling construction of a deep reference image for discovery of new objects. Supernova imaging observations are being acquired between 1 September and 30 November of 2005-7. During the first two seasons, each region was imaged on average every five nights. Spectroscopic follow-up observations to determine supernova type and redshift are carried out on a large number of telescopes. In its first two three-month seasons, the survey has discovered and measured light curves for 327 spectroscopically confirmed SNe Ia, 30 probable SNe Ia, 14 confirmed SNe Ib/c, 32 confirmed SNe II, plus a large number of photometrically identified SNe Ia, 94 of which have host-galaxy spectra taken so far. This paper provides an overview of the project and briefly describes the observations completed during the first two seasons of operation.

The first cosmological results from the ESSENCE supernova survey (Wood-Vasey et al. 2007) are extended to a wider range of cosmological models including dynamical dark energy and non-standard cosmological models. We fold in a greater number of external data sets such as the recent Higher-z release of high-redshift supernovae (Riess et al. 2007) as well as several complementary cosmological probes. Model comparison statistics such as the Bayesian and Akaike information criteria are applied to gauge the worth of models. These statistics favor models that give a good fit with fewer parameters. Based on this analysis, the preferred cosmological model is the flat cosmological constant model, where the expansion history of the universe can be adequately described with only one free parameter describing the energy content of the universe. Among the more exotic models that provide good fits to the data, we note a preference for models whose best-fit parameters reduce them to the cosmological constant model.

We present extensive ugrizYHJK photometry and optical spectroscopy of SN 2005gj obtained by the SDSS-II and CSP Supernova Projects, which give excellent coverage during the first 150 days after the time of explosion. These data show that SN 2005gj is the second clear case, after SN 2002ic, of a thermonuclear explosion in a dense circumstellar environment. Both the presence of singly and doubly ionized iron-peak elements (FeIII and weak SII, SiII) near maximum light as well as the spectral evolution show that SN 2002ic-like events are Type Ia explosions. Independent evidence comes from the exponential decay in luminosity of SN 2005gj, pointing to an exponential density distribution of the ejecta. The interaction of the supernova ejecta with the dense circumstellar medium is stronger than in SN 2002ic: (1) the supernova lines are weaker; (2) the Balmer emission lines are more luminous; and (3) the bolometric luminosity is higher close to maximum light. The velocity evolution of the Halpha components suggest that the CSM around SN 2005gj is clumpy and it has a flatter density distribution compared with the steady wind solution, in agreement with SN 2002ic. An early X-ray observation with Chandra gives an upper-limit on the mass loss rate from the companion of < 2x10^{-4} Msun/yr.

We quantify the variability of faint unresolved optical sources using a catalog based on multiple SDSS imaging observations. The catalog covers SDSS Stripe 82, and contains 58 million photometric observations in the SDSS ugriz system for 1.4 million unresolved sources. In each photometric bandpass we compute various low-order lightcurve statistics and use them to select and study variable sources. We find that 2% of unresolved optical sources brighter than g=20.5 appear variable at the 0.05 mag level (rms) simultaneously in the g and r bands. The majority (2/3) of these variable sources are low-redshift (<2) quasars, although they represent only 2% of all sources in the adopted flux-limited sample. We find that at least 90% of quasars are variable at the 0.03 mag level (rms) and confirm that variability is as good a method for finding low-redshift quasars as is the UV excess color selection (at high Galactic latitudes). We analyze the distribution of lightcurve skewness for quasars and find that is centered on zero. We find that about 1/4 of the variable stars are RR Lyrae stars, and that only 0.5% of stars from the main stellar locus are variable at the 0.05 mag level. The distribution of lightcurve skewness in the g-r vs. u-g color-color diagram on the main stellar locus is found to be bimodal (with one mode consistent with Algol-like behavior). Using over six hundred RR Lyrae stars, we demonstrate rich halo substructure out to distances of 100 kpc. We extrapolate these results to expected performance by the Large Synoptic Survey Telescope and estimate that it will obtain well-sampled 2% accurate, multi-color lightcurves for ~2 million low-redshift quasars, and will discover at least 50 million variable stars.

We present extensive u'g'r'i'BVRIYJHKs photometry and optical spectroscopy of SN 2005hk. These data reveal that SN 2005hk was nearly identical in its observed properties to SN 2002cx, which has been called ``the most peculiar known type Ia supernova.'' Both supernovae exhibited high ionization SN 1991T-like pre-maximum spectra, yet low peak luminosities like SN 1991bg. The spectra reveal that SN 2005hk, like SN 2002cx, exhibited expansion velocities that were roughly half those of typical type Ia supernovae. The R and I light curves of both supernovae were also peculiar in not displaying the secondary maximum observed for normal type Ia supernovae. Our YJH photometry of SN 2005hk reveals the same peculiarity in the near-infrared. By combining our optical and near-infrared photometry of SN 2005hk with published ultraviolet light curves obtained with the Swift satellite, we are able to construct a bolometric light curve from ~10 days before to ~60 days after B maximum. The shape and unusually low peak luminosity of this light curve, plus the low expansion velocities and absence of a secondary maximum at red and near-infrared wavelengths, are all in reasonable agreement with model calculations of a 3D deflagration which produces ~0.25 M_sun of 56Ni.

We describe a standard star catalog constructed using multiple SDSS photometric observations (at least four per band, with a median of ten) in the $ugriz$ system. The catalog includes 1.01 million non-variable unresolved objects from the equatorial stripe 82 ($|\delta_{J2000}|<$ 1.266$^\circ$) in the RA range 20h 34m to 4h 00m, and with the corresponding $r$ band (approximately Johnson V band) magnitudes in the range 14--22. The distributions of measurements for individual sources demonstrate that the photometric pipeline correctly estimates random photometric errors, which are below 0.01 mag for stars brighter than (19.5, 20.5, 20.5, 20, 18.5) in $ugriz$, respectively (about twice as good as for individual SDSS runs). Several independent tests of the internal consistency suggest that the spatial variation of photometric zeropoints is not larger than $\sim$0.01 mag (rms). In addition to being the largest available dataset with optical photometry internally consistent at the $\sim$1% level, this catalog provides practical definition of the SDSS photometric system. Using this catalog, we show that photometric zeropoints for SDSS observing runs can be calibrated within nominal uncertainty of 2% even for data obtained through 1 mag thick clouds, and demonstrate the existence of He and H white dwarf sequences using photometric data alone. Based on the properties of this catalog, we conclude that upcoming large-scale optical surveys such as the Large Synoptic Survey Telescope will be capable of delivering robust 1% photometry for billions of sources.

We compare Stetson's photometric standards with measurements listed in a standard star catalog constructed using repeated SDSS imaging observations. The SDSS catalog includes over 700,000 candidate standard stars from the equatorial stripe 82 (|Dec|<1.266 deg) in the RA range 20h 34' to 4h 00', and with the $r$ band magnitudes in the range 14--21. The distributions of measurements for individual sources demonstrate that the SDSS photometric pipeline correctly estimates random photometric errors, which are below 0.01 mag for stars brighter than (19.5, 20.5, 20.5, 20, 18.5) in ugriz, respectively (about twice as good as for individual SDSS runs). We derive mean photometric transformations between the SDSS gri and the BVRI system using 1165 Stetson stars found in the equatorial stripe 82, and then study the spatial variation of the difference in zeropoints between the two catalogs. Using third order polynomials to describe the color terms, we find that photometric measurements for main-sequence stars can be transformed between the two systems with systematic errors smaller than a few millimagnitudes. The spatial variation of photometric zeropoints in the two catalogs typically does not exceed 0.01 magnitude. Consequently, the SDSS Standard Star Catalog for Stripe 82 can be used to calibrate new data in both the SDSS ugriz and the BVRI systems with a similar accuracy.

We describe the implementation and optimization of the ESSENCE supernova survey, which we have undertaken to measure the equation of state parameter of the dark energy. We present a method for optimizing the survey exposure times and cadence to maximize our sensitivity to the dark energy equation of state parameter w=P/rho c^2 for a given fixed amount of telescope time. For our survey on the CTIO 4m telescope, measuring the luminosity distances and redshifts for supernovae at modest redshifts (z~0.5 +- 0.2) is optimal for determining w. We describe the data analysis pipeline based on using reliable and robust image subtraction to find supernovae automatically and in near real-time. Since making cosmological inferences with supernovae relies crucially on accurate measurement of their brightnesses, we describe our efforts to establish a thorough calibration of the CTIO 4m natural photometric system. In its first four years, ESSENCE has discovered and spectroscopically confirmed 102 type Ia SNe, at redshifts from 0.10 to 0.78, identified through an impartial, effective methodology for spectroscopic classification and redshift determination. We present the resulting light curves for the all type Ia supernovae found by ESSENCE and used in our measurement of w, presented in Wood-Vasey et al, 2007.

We present constraints on the dark energy equation-of-state parameter, w=P/(rho c^2), using 60 Type Ia supernovae (SNe Ia) from the ESSENCE supernova survey. We derive a set of constraints on the nature of the dark energy assuming a flat Universe. By including constraints on (Omega_M, w) from baryon acoustic oscillations, we obtain a value for a static equation-of-state parameter w=-1.05^{+0.13}_{-0.12} (stat; 1 sigma) +- 0.11 (sys) and Omega_M=0.274^{+0.033}_{-0.020} (stat; 1 sigma) with a best-fit chi^2/DoF of 0.96. These results are consistent with those reported by the SuperNova Legacy Survey in a similar program measuring supernova distances and redshifts. We evaluate sources of systematic error that afflict supernova observations and present Monte Carlo simulations that explore these effects. Currently, the largest systematic currently with the potential to affect our measurements is the treatment of extinction due to dust in the supernova host galaxies. Combining our set of ESSENCE SNe Ia with the SuperNova Legacy Survey SNe Ia, we obtain a joint constraint of w=-1.07^{+0.09}_{-0.09} (stat; 1 sigma) +- 0.12 (sys), Omega_M=0.267^{+0.028}_{-0.018} (stat; 1 sigma) with a best-fit chi^2/DoF of 0.91. The current SN Ia data are fully consistent with a cosmological constant.

We present a set of 11 type Ia supernova (SN Ia) lightcurves with dense, pre-maximum sampling. These supernovae (SNe), in galaxies behind the Large Magellanic Cloud (LMC), were discovered by the SuperMACHO survey. The SNe span a redshift range of z = 0.11 - 0.35. Our lightcurves contain some of the earliest pre-maximum observations of SNe Ia to date. We also give a functional model that describes the SN Ia lightcurve shape (in our VR-band). Our function uses the "expanding fireball" model of Goldhaber et al. (1998) to describe the rising lightcurve immediately after explosion but constrains it to smoothly join the remainder of the lightcurve. We fit this model to a composite observed VR-band lightcurve of three SNe between redshifts of 0.135 to 0.165. These SNe have not been K-corrected or adjusted to account for reddening. In this redshift range, the observed VR-band most closely matches the rest frame V-band. Using the best fit to our functional description of the lightcurve, we find the time between explosion and observed VR-band maximum to be 17.6+-1.3(stat)+-0.07(sys) rest-frame days for a SN Ia with a VR-band Delta m_{-10} of 0.52mag. For the redshifts sampled, the observed VR-band time-of-maximum brightness should be the same as the rest-frame V-band maximum to within 1.1 rest-frame days.

In principle, the light from historical supernovae could still be visible as scattered-light echoes even centuries later. However, while echoes have been discovered around some nearby extragalactic supernovae well after the explosion, targeted searches have not recovered any echoes in the regions of historical Galactic supernovae. The discovery of echoes can allow us to pinpoint the supernova event both in position and age and, most importantly, allow us to acquire spectra of the echo light to type the supernova centuries after the direct light from the explosion first reached the Earth. Here we report on the discovery of three faint new variable surface brightness complexes with high apparent proper motion pointing back to well-defined positions in the Large Magellanic Cloud (LMC). These positions correspond to three of the six smallest (and likely youngest) previously catalogued supernova remnants, and are believed to be due to thermonuclear (Type Ia) supernovae. Using the distance and proper motions of these echo arcs, we estimate ages of 610 and 410 yr for the echoes #2 and #3.

Using archival data of low-redshift (z < 0.01) Type Ia supernovae (SN Ia) and recent observations of high-redshift (0.16 < z <0.64; Matheson et al. 2005) SN Ia, we study the "uniformity'' of the spectroscopic properties of nearby and distant SN Ia. We find no difference in the measures we describe here. In this paper, we base our analysis solely on line-profile morphology, focusing on measurements of the velocity location of maximum absorption (vabs) and peak emission (vpeak). We find that the evolution of vabs and vpeak for our sample lines (Ca II 3945, Si II 6355, and S II 5454, 5640) is similar for both the low- and high-redshift samples. We find that vabs for the weak S II 5454, 5640 lines, and vpeak for S II 5454, can be used to identify fast-declining [dm15 > 1.7] SN Ia, which are also subluminous. In addition, we give the first direct evidence in two high-z SN Ia spectra of a double-absorption feature in Ca II 3945, an event also observed, though infrequently, in low-redshift SN Ia spectra (6/22 SN Ia in our local sample). We report for the first time the unambiguous and systematic intrinsic blueshift of peak emission of optical P-Cygni line profiles in Type Ia spectra, by as much as 8000 km/s. All the high-z SN Ia analyzed in this paper were discovered and followed up by the ESSENCE collaboration, and are now publicly available.

The proper usage of Type Ia supernovae (SNe Ia) as distance indicators has revolutionized cosmology, and added a new dominant component to the energy density of the Universe, dark energy. Following the discovery and confirmation era, the currently ongoing SNe Ia surveys aim to determine the properties of the dark energy. ESSENCE is a five year ground-based supernova survey aimed at finding and characterizing 200 SNe Ia in the redshift domain z=[0.2-0.8]. The goal of the project is to put constraints on the equation of state parameter, w, of the dark energy with an accuracy of <10%. This paper presents these ongoing efforts in the context of the current developments in observational cosmology.

Characterizing the nature and spatial distribution of the lensing objects that produce the previously measured microlensing optical depth toward the Large Magellanic Cloud (LMC) remains an open problem. We present an appraisal of the ability of the SuperMACHO Project, a next-generation microlensing survey directed toward the LMC, to discriminate between various proposed lensing populations. We consider two scenarios: lensing by a uniform foreground screen of objects and self-lensing by LMC stars. We have carried out extensive simulations, based upon data obtained during the first year of the project, to assess the SuperMACHO survey's ability to discriminate between these two scenarios. We find that the event rate itself shows significant sensitivity to the choice of the LMC luminosity function, limiting the conclusions which can be drawn from the absolute rate. If instead we determine the differential event rate across the LMC, we will decrease the impact of these systematic biases and render our conclusions more robust. With this approach the SuperMACHO Project should be able to distinguish between the two categories of lens populations. This will provide important constraints on the nature of the lensing objects and their contributions to the Galactic dark matter halo.

We present broad-band light curves of nine supernovae ranging in redshift from 0.5 to 0.8. The supernovae were discovered as part of the ESSENCE project, and the light curves are a combination of Cerro Tololo 4-m and Hubble Space Telescope (HST) photometry. On the basis of spectra and/or light-curve fitting, eight of these objects are definitely Type Ia supernovae, while the classification of one is problematic. The ESSENCE project is a five-year endeavor to discover about 200 high-redshift Type Ia supernovae, with the goal of tightly constraining the time average of the equation-of-state parameter [w = p/(rho c^2)] of the "dark energy." To help minimize our systematic errors, all of our ground-based photometry is obtained with the same telescope and instrument. In 2003 the highest-redshift subset of ESSENCE supernovae was selected for detailed study with HST. Here we present the first photometric results of the survey. We find that all but one of the ESSENCE SNe have slowly declining light curves, and the sample is not representative of the low-redshift set of ESSENCE Type Ia supernovae. This is unlikely to be a sign of evolution in the population. We attribute the decline-rate distribution of HST events to a selection bias at the high-redshift edge of our sample and find that such a bias will infect other magnitude-limited SN Ia searches unless appropriate precautions are taken.

In preparation for the Supernova Survey of the Sloan Digital Sky Survey (SDSS) II, a proposed 3-year extension to the SDSS, we have conducted an early engineering and science run during the fall of 2004, which consisted of approximately 20 scheduled nights of repeated imaging of half of the southern equatorial stripe. Transient supernova-like events were detected in near real-time and photometric measurements were made in the five SDSS filter bandpasses with a cadence of ~2 days. Candidate type Ia supernovae (SNe) were pre-selected based on their colors, light curve shape, and the properties of the host galaxy. Follow-up spectroscopic observations were performed with the Astrophysical Research Consortium 3.5m telescope and the 9.2m Hobby-Eberly Telescope to confirm their types and measure the redshifts. The 2004 campaign resulted in 22 spectroscopically confirmed SNe, which includes 16 type Ia, 5 type II, and 1 type Ib/c. These SN Ia will help fill in the sparsely sampled redshift interval of z = 0.05 - 0.35, the so-called 'redshift desert', in the Hubble diagram. Detailed investigation of the spectral properties of these moderate-redshift SNe Ia will also provide a bridge between local SNe and high-redshift objects, and will help us understand the systematics for future cosmological applications that require high photometric precision. Finally, the large survey volume also provides the opportunity to select unusual supernovae for spectroscopic study that are poorly sampled in other surveys. We report on some of the early results from this program and discuss potential future applications.

We present the results of spectroscopic observations of targets discovered during the first two years of the ESSENCE project. The goal of ESSENCE is to use a sample of ~200 Type Ia supernovae (SNe Ia) at moderate redshifts (0.2 < z < 0.8) to place constraints on the equation of state of the Universe. Spectroscopy not only provides the redshifts of the objects, but also confirms that some of the discoveries are indeed SNe Ia. This confirmation is critical to the project, as techniques developed to determine luminosity distances to SNe Ia depend upon the knowledge that the objects at high redshift are the same as the ones at low redshift. We describe the methods of target selection and prioritization, the telescopes and detectors, and the software used to identify objects. The redshifts deduced from spectral matching of high-redshift SNe Ia with low-redshift SNe Ia are consistent with those determined from host-galaxy spectra. We show that the high-redshift SNe Ia match well with low-redshift templates. We include all spectra obtained by the ESSENCE project, including 52 SNe Ia, 5 core-collapse SNe, 12 active galactic nuclei, 19 galaxies, 4 possibly variable stars, and 16 objects with uncertain identifications.

We present the first results from our next-generation microlensing survey, the SuperMACHO project. We are using the CTIO 4m Blanco telescope and the MOSAIC imager to carry out a search for microlensing toward the Large Magellanic Cloud (LMC). We plan to ascertain the nature of the population responsible for the excess microlensing rate seen by the MACHO project. Our observing strategy is optimized to measure the differential microlensing rate across the face of the LMC. We find this derivative to be relatively insensitive to the details of the LMC's internal structure but a strong discriminant between Galactic halo and LMC self lensing. In December 2003 we completed our third year of survey operations. 2003 also marked the first year of real-time microlensing alerts and photometric and spectroscopic followup. We have extracted several dozen microlensing candidates, and we present some preliminary light curves and related information. Similar to the MACHO project, we find SNe behind the LMC to be a significant contaminant - this background has not been completely removed from our current single-color candidate sample. Our follow-up strategy is optimized to discriminate between SNe and true microlensing.

We present a detailed observational strategy for finding Near Earth Objects (NEOs) with the Sloan Digital Sky Survey (SDSS) telescope. We investigate strategies in normal, unbinned mode as well as binning the CCDs 2x2 or 3x3, which affects the sky coverage rate and the limiting apparent magnitude. We present results from 1 month, 3 year and 10 year simulations of such surveys. For each cadence and binning mode, we evaluate the possibility of achieving the Spaceguard goal of detecting 90% of 1 km NEOs (absolute magnitude H <= 18 for an albedo of 0.1). We find that an unbinned survey is most effective at detecting H <= 20 NEOs in our sample. However, a 3x3 binned survey reaches the Spaceguard Goal after only seven years of operation. As the proposed large survey telescopes (PanStarss; LSST) are at least 5-10 years from operation, an SDSS NEO survey could make a significant contribution to the detection and photometric characterization of the NEO population.

We present photometric and spectroscopic observations of 23 high redshift supernovae spanning a range of z=0.34-1.03, 9 of which are unambiguously classified as Type Ia. These supernovae were discovered during the IfA Deep Survey, which began in September 2001 and observed a total of 2.5 square degrees to a depth of approximately m=25-26 in RIZ over 9-17 visits, typically every 1-3 weeks for nearly 5 months, with additional observations continuing until April 2002. We give a brief description of the survey motivations, observational strategy, and reduction process. This sample of 23 high-redshift supernovae includes 15 at z>0.7, doubling the published number of objects at these redshifts, and indicates that the evidence for acceleration of the universe is not due to a systematic effect proportional to redshift. In combination with the recent compilation of Tonry et al. (2003), we calculate cosmological parameter density contours which are consistent with the flat universe indicated by the CMB (Spergel et al. 2003). Adopting the constraint that Omega_total = 1.0, we obtain best-fit values of (Omega_m, Omega_Lambda)=(0.33, 0.67) using 22 SNe from this survey augmented by the literature compilation. We show that using the empty-beam model for gravitational lensing does not eliminate the need for Omega_Lambda > 0. Experience from this survey indicates great potential for similar large-scale surveys while also revealing the limitations of performing surveys for z>1 SNe from the ground.

The potential of the Sloan Digital Sky Survey for wide-field variability studies is illustrated using multi-epoch observations for 3,000,000 point sources observed in 700 deg2 of sky, with time spans ranging from 3 hours to 3 years. These repeated observations of the same sources demonstrate that SDSS delivers ~0.02 mag photometry with well behaved and understood errors. We show that quasars dominate optically faint (r > 18) point sources that are variable on time scales longer than a few months, while for shorter time scales, and at bright magnitudes, most variable sources are stars.

The commissioning year of the Sloan Digital Sky Survey has demonstrated that many cataclysmic variables have been missed in previous surveys with brighter limits. We report the identification of 22 cataclysmic variables, of which 19 are new discoveries and 3 are known systems (SW UMa, BH Lyn and Vir4). A compendium of positions, colors and characteristics of these systems obtained from the SDSS photometry and spectroscopy is presented along with data obtained during follow-up studies with the Apache Point Observatory (APO) and Manastash Ridge Observatory (MRO) telescopes. We have determined orbital periods for 3 of the new systems: two show dwarf nova outbursts, and the third is a likely magnetic system with eclipses of its region of line emission. Based on these results, we expect the completed survey to locate at least 400 new CVs. Most of these will be faint systems with low accretion rates that will provide new constraints on binary evolution models.

We present BVRI photometry of the Type Ia supernovae 1999da, 1999dk, 1999gp, 2000bk, and 2000ce, plus infrared photometry of three of these. These objects exhibit the full range of decline rates of Type Ia supernovae. Combined optical and infrared data show that families of V - infrared color curves can be used to derive the host extinction (A_V) of these objects. Existing data do not yet allow us to construct these loci for all color indices and supernova decline rates, but the V-K color evolution is sufficiently uniform that it allows the determination of host extinction over a wide range of supernova decline rates to an accuracy of roughly +/- 0.1 mag. We introduce a new empirical parameter, the mean I-band flux 20 to 40 days after maximum light, and show how it is directly related to the decline rate.