We measure the weak-lensing masses and galaxy distributions of four massive galaxy clusters observed during the Science Verification phase of the Dark Energy Survey. This pathfinder study is meant to 1) validate the DECam imager for the task of measuring weak-lensing shapes, and 2) utilize DECam's large field of view to map out the clusters and their environments over 90 arcmin. We conduct a series of rigorous tests on astrometry, photometry, image quality, PSF modeling, and shear measurement accuracy to single out flaws in the data and also to identify the optimal data processing steps and parameters. We find Science Verification data from DECam to be suitable for the lensing analysis described in this paper. The PSF is generally well-behaved, but the modeling is rendered difficult by a flux-dependent PSF width and ellipticity. We employ photometric redshifts to distinguish between foreground and background galaxies, and a red-sequence cluster finder to provide cluster richness estimates and cluster-galaxy distributions. By fitting NFW profiles to the clusters in this study, we determine weak-lensing masses that are in agreement with previous work. For Abell 3261, we provide the first estimates of redshift, weak-lensing mass, and richness. In addition, the cluster-galaxy distributions indicate the presence of filamentary structures attached to 1E 0657-56 and RXC J2248.7-4431, stretching out as far as 1 degree (approximately 20 Mpc), showcasing the potential of DECam and DES for detailed studies of degree-scale features on the sky.

The Dark Energy Camera has captured a large set of images as part of Science Verification (SV) for the Dark Energy Survey. The SV footprint covers a lar ge portion of the outer Large Magellanic Cloud (LMC), providing photometry 1.5 magnitudes fainter than the main sequence turn-off of the oldest LMC stel lar population. We derive geometrical and structural parameters for various stellar populations in the LMC disk. For the distribution of all LMC stars, we find an inclination of $i=-38.14^{\circ}\pm0.08^{\circ}$ (near side in the North) and a position angle for the line of nodes of $\theta_0=129.51^{\circ}\pm0.17^{\circ}$. We find that stars younger than $\sim 4$ Gyr are more centrally concentrated than older stars. Fitting a projected exponential disk shows that the scale radius of the old populations is $R_{>4 Gyr}=1.41\pm0.01$ kpc, while the younger population has $R_{<4 Gyr}=0.72\pm0.01$ kpc. Howe ver, the spatial distribution of the younger population deviates significantly from the projected exponential disk model. The distribution of old stars suggests a large truncation radius of $R_{t}=13.5\pm0.8$ kpc. If this truncation is dominated by the tidal field of the Galaxy, we find that the LMC is $\simeq 24^{+9}_{-6}$ times less massive than the encircled Galactic mass. By measuring the Red Clump peak magnitude and comparing with the best-fit LM C disk model, we find that the LMC disk is warped and thicker in the outer regions north of the LMC centre. Our findings may either be interpreted as a warped and flared disk in the LMC outskirts, or as evidence of a spheroidal halo component

We present DES13S2cmm, the first spectroscopically-confirmed superluminous supernova (SLSN) from the Dark Energy Survey (DES). We briefly discuss the data and search algorithm used to find this event in the first year of DES operations, and outline the spectroscopic data obtained from the European Southern Observatory (ESO) Very Large Telescope to confirm its redshift (z = 0.663 +/- 0.001 based on the host-galaxy emission lines) and likely spectral type (type I). Using this redshift, we find M_U_peak = -21.05 +0.10 -0.09 for the peak, rest-frame U-band absolute magnitude, and find DES13S2cmm to be located in a faint, low metallicity (sub-solar), low stellar-mass host galaxy (log(M/M_sun) = 9.3 +/- 0.3); consistent with what is seen for other SLSNe-I. We compare the bolometric light curve of DES13S2cmm to fourteen similarly well-observed SLSNe-I in the literature and find it possesses one of the slowest declining tails (beyond +30 days rest frame past peak), and is the faintest at peak. Moreover, we find the bolometric light curves of all SLSNe-I studied herein possess a dispersion of only 0.2-0.3 magnitudes between +25 and +30 days after peak (rest frame) depending on redshift range studied; this could be important for 'standardising' such supernovae, as is done with the more common type Ia. We fit the bolometric light curve of DES13S2cmm with two competing models for SLSNe-I - the radioactive decay of 56Ni, and a magnetar - and find that while the magnetar is formally a better fit, neither model provides a compelling match to the data. Although we are unable to conclusively differentiate between these two physical models for this particular SLSN-I, further DES observations of more SLSNe-I should break this degeneracy, especially if the light curves of SLSNe-I can be observed beyond 100 days in the rest frame of the supernova.

We present the combination of optical data from the Science Verification phase of the Dark Energy Survey (DES) with near infrared data from the ESO VISTA Hemisphere Survey (VHS). The deep optical detections from DES are used to extract fluxes and associated errors from the shallower VHS data. Joint 7-band ($grizYJK$) photometric catalogues are produced in a single 3 sq-deg DECam field centred at 02h26m$-$04d36m where the availability of ancillary multi-wavelength photometry and spectroscopy allows us to test the data quality. Dual photometry increases the number of DES galaxies with measured VHS fluxes by a factor of $\sim$4.5 relative to a simple catalogue level matching and results in a $\sim$1.5 mag increase in the 80\% completeness limit of the NIR data. Almost 70\% of DES sources have useful NIR flux measurements in this initial catalogue. Photometric redshifts are estimated for a subset of galaxies with spectroscopic redshifts and initial results, although currently limited by small number statistics, indicate that the VHS data can help reduce the photometric redshift scatter at both $z<0.5$ and $z>1$. We present example DES+VHS colour selection criteria for high redshift Luminous Red Galaxies (LRGs) at $z\sim0.7$ as well as luminous quasars. Using spectroscopic observations in this field we show that the additional VHS fluxes enable a cleaner selection of both populations with $<$10\% contamination from galactic stars in the case of spectroscopically confirmed quasars and $<0.5\%$ contamination from galactic stars in the case of spectroscopically confirmed LRGs. The combined DES+VHS dataset, which will eventually cover almost 5000 sq-deg, will therefore enable a range of new science and be ideally suited for target selection for future wide-field spectroscopic surveys.

We present results from a study of the photometric redshift performance of the Dark Energy Survey (DES), using the early data from a Science Verification (SV) period of observations in late 2012 and early 2013 that provided science-quality images for almost 200 sq.~deg.~at the nominal depth of the survey. We assess the photometric redshift performance using about 15000 galaxies with spectroscopic redshifts available from other surveys. These galaxies are used, in different configurations, as a calibration sample, and photo-$z$'s are obtained and studied using most of the existing photo-$z$ codes. A weighting method in a multi-dimensional color-magnitude space is applied to the spectroscopic sample in order to evaluate the photo-$z$ performance with sets that mimic the full DES photometric sample, which is on average significantly deeper than the calibration sample due to the limited depth of spectroscopic surveys. Empirical photo-$z$ methods using, for instance, Artificial Neural Networks or Random Forests, yield the best performance in the tests, achieving core photo-$z$ resolutions $\sigma_{68} \sim 0.08$. Moreover, the results from most of the codes, including template fitting methods, comfortably meet the DES requirements on photo-$z$ performance, therefore, providing an excellent precedent for future DES data sets.

PreCam, a precursor observational campaign supporting the Dark Energy Survey (DES), is designed to produce a photometric and astrometric catalog of nearly a hundred thousand standard stars within the DES footprint, while the PreCam instrument also serves as a prototype testbed for the Dark Energy Camera (DECam)'s hardware and software. This catalog represents a potential 100-fold increase in Southern Hemisphere photometric standard stars, and therefore will be an important component in the calibration of the Dark Energy Survey. We provide details on the PreCam instrument's design, construction and testing, as well as results from a subset of the 51 nights of PreCam survey observations on the University of Michigan Department of Astronomy's Curtis-Schmidt telescope at Cerro Tololo Inter-American Observatory. We briefly describe the preliminary data processing pipeline that has been developed for PreCam data and the preliminary results of the instrument performance, as well as astrometry and photometry of a sample of stars previously included in other southern sky surveys.

We present an analysis of supernova light curves simulated for the upcoming Dark Energy Survey (DES) supernova search. The simulations employ a code suite that generates and fits realistic light curves in order to obtain distance modulus/redshift pairs that are passed to a cosmology fitter. We investigated several different survey strategies including field selection, supernova selection biases, and photometric redshift measurements. Using the results of this study, we chose a 30 square degree search area in the griz filter set. We forecast 1) that this survey will provide a homogeneous sample of up to 4000 Type Ia supernovae in the redshift range 0.05<z<1.2, and 2) that the increased red efficiency of the DES camera will significantly improve high-redshift color measurements. The redshift of each supernova with an identified host galaxy will be obtained from spectroscopic observations of the host. A supernova spectrum will be obtained for a subset of the sample, which will be utilized for control studies. In addition, we have investigated the use of combined photometric redshifts taking into account data from both the host and supernova. We have investigated and estimated the likely contamination from core-collapse supernovae based on photometric identification, and have found that a Type Ia supernova sample purity of up to 98% is obtainable given specific assumptions. Furthermore, we present systematic uncertainties due to sample purity, photometric calibration, dust extinction priors, filter-centroid shifts, and inter-calibration. We conclude by estimating the uncertainty on the cosmological parameters that will be measured from the DES supernova data.

We present the results of a study of selection criteria to identify Type Ia supernovae photometrically in a simulated mixed sample of Type Ia supernovae and core collapse supernovae. The simulated sample is a mockup of the expected results of the Dark Energy Survey. Fits to the MLCS2k2 and SALT2 Type Ia supernova models are compared and used to help separate the Type Ia supernovae from the core collapse sample. The Dark Energy Task Force Figure of Merit (modified to include core collapse supernovae systematics) is used to discriminate among the various selection criteria. This study of varying selection cuts for Type Ia supernova candidates is the first to evaluate core collapse contamination using the Figure of Merit. Different factors that contribute to the Figure of Merit are detailed. With our analysis methods, both SALT2 and MLCS2k2 Figures of Merit improve with tighter selection cuts and higher purities, peaking at 98% purity.

MSH 11-62 (G291.1-0.9) is a composite supernova remnant for which radio and X-ray observations have identified the remnant shell as well as its central pulsar wind nebula. The observations suggest a relatively young system expanding into a low density region. Here we present a study of MSH 11-62 using observations with the Chandra, XMM-Newton, and Fermi observatories, along with radio observations from the Australia Telescope Compact Array (ATCA). We identify a compact X-ray source that appears to be the putative pulsar that powers the nebula, and show that the X-ray spectrum of the nebula bears the signature of synchrotron losses as particles diffuse into the outer nebula. Using data from the Fermi LAT, we identify gamma-ray emission originating from MSH 11-62. With density constraints from the new X-ray measurements of the remnant, we model the evolution of the composite system in order to constrain the properties of the underlying pulsar and the origin of the gamma-ray emission.

Supernova cosmology without spectroscopic confirmation is an exciting new frontier which we address here with the Bayesian Estimation Applied to Multiple Species (BEAMS) algorithm and the full three years of data from the Sloan Digital Sky Survey II Supernova Survey (SDSS-II SN). BEAMS is a Bayesian framework for using data from multiple species in statistical inference when one has the probability that each data point belongs to a given species, corresponding in this context to different types of supernovae with their probabilities derived from their multi-band lightcurves. We run the BEAMS algorithm on both Gaussian and more realistic SNANA simulations with of order 10^4 supernovae, testing the algorithm against various pitfalls one might expect in the new and somewhat uncharted territory of photometric supernova cosmology. We compare the performance of BEAMS to that of both mock spectroscopic surveys and photometric samples which have been cut using typical selection criteria. The latter typically are either biased due to contamination or have significantly larger contours in the cosmological parameters due to small data-sets. We then apply BEAMS to the 792 SDSS-II photometric supernovae with host spectroscopic redshifts. In this case, BEAMS reduces the area of the (\Omega_m,\Omega_\Lambda) contours by a factor of three relative to the case where only spectroscopically confirmed data are used (297 supernovae). In the case of flatness, the constraints obtained on the matter density applying BEAMS to the photometric SDSS-II data are \Omega_m(BEAMS)=0.194\pm0.07. This illustrates the potential power of BEAMS for future large photometric supernova surveys such as LSST.

We describe a new formalism to fit the parameters $\alpha$ and $\beta$ that are used in the SALT2 model to determine the standard magnitudes of Type Ia supernovae. The new formalism describes the intrinsic scatter in Type Ia supernovae by a covariance matrix in place of the single parameter normally used. We have applied this formalism to the Sloan Digital Sky Survey Supernova Survey (SDSS-II) data and conclude that the data are best described by $\alpha=0.135_{-.017}^{+.033}$ and $\beta=3.19_{-0.24}^{+0.14}$, where the error is dominated by the uncertainty in the form of the intrinsic scatter matrix. Our result depends on the introduction of a more general form for the intrinsic scatter of the distance moduli of Type Ia supernovae than is conventional, resulting in a larger value of $\beta$ and a larger uncertainty than the conventional approach. Although this analysis results in a larger value of $\beta$ and a larger error, the SDSS data differ (at a 98% confidence level) with $\beta=4.1$, the value expected for extinction by the type of dust found in the Milky Way. We have modeled the distribution of supernovae Ia in terms of their color and conclude that there is strong evidence that variation in color is a significant contributor to the scatter of supernovae Ia around their standard candle magnitude.

Large planned photometric surveys will discover hundreds of thousands of supernovae (SNe), outstripping the resources available for spectroscopic follow-up and necessitating the development of purely photometric methods to exploit these events for cosmological study. We present a light-curve fitting technique for SN Ia photometric redshift (photo-z) estimation in which the redshift is determined simultaneously with the other fit parameters. We implement this "LCFIT+Z" technique within the frameworks of the MLCS2k2 and SALT-II light-curve fit methods and determine the precision on the redshift and distance modulus. This method is applied to a spectroscopically confirmed sample of 296 SNe Ia from the SDSS-II Supernova Survey and 37 publicly available SNe Ia from the Supernova Legacy Survey (SNLS). We have also applied the method to a large suite of realistic simulated light curves for existing and planned surveys, including SDSS, SNLS, and LSST. When intrinsic SN color fluctuations are included, the photo-z precision for the simulation is consistent with that in the data. Finally, we compare the LCFIT+Z photo-z precision with previous results using color-based SN photo-z estimates.

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.

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 have undertaken the simulation of hydrodynamic flows with bulk Lorentz factors in the range 10^2--10^6. We discuss the application of an existing relativistic, hydrodynamic primitive-variable recovery algorithm to a study of pulsar winds, and, in particular, the refinement made to admit such ultra-relativistic flows. We show that an iterative quartic root finder breaks down for Lorentz factors above 10^2 and employ an analytic root finder as a solution. We find that the former, which is known to be robust for Lorentz factors up to at least 50, offers a 24% speed advantage. We demonstrate the existence of a simple diagnostic allowing for a hybrid primitives recovery algorithm that includes an automatic, real-time toggle between the iterative and analytical methods. We further determine the accuracy of the iterative and hybrid algorithms for a comprehensive selection of input parameters and demonstrate the latter's capability to elucidate the internal structure of ultra-relativistic plasmas. In particular, we discuss simulations showing that the interaction of a light, ultra-relativistic pulsar wind with a slow, dense ambient medium can give rise to asymmetry reminiscent of the Guitar nebula leading to the formation of a relativistic backflow harboring a series of internal shockwaves. The shockwaves provide thermalized energy that is available for the continued inflation of the PWN bubble. In turn, the bubble enhances the asymmetry, thereby providing positive feedback to the backflow.

We present simulations for the Dark Energy Survey (DES) using a new code suite (SNANA) that generates realistic supernova light curves accounting for atmospheric seeing conditions and intrinsic supernova luminosity variations using MLCS2k2 or SALT2 models. Errors include stat-noise from photo-statistics and sky noise. We applied SNANA to simulate DES supernova observations and employed an MLCS-based fitter to obtain the distance modulus for each simulated light curve. We harnessed the light curves in order to study selection biases for high-redshift supernovae and to constrain the optimal DES observing strategy using the Dark Energy Task Force figure of merit.

The brightest, or first-ranked, galaxies (BCGs) in rich clusters show a very small dispersion in luminosity, making them excellent standard candles. This small dispersion has raised questions about the nature of BCGs. Are they simply the extremes of normal galaxies formed via a stochastic process, or do they belong to a special class of atypical objects? Arguments have been proposed on both sides of the debate. Bhavsar (1989) suggested that the distribution in magnitudes can only be explained by a two-population model. Thus, a new controversy has arisen. Do first-ranked galaxies consist of one or two populations of objects? We examine an older and newer data set and present our results. Two-population models do better than do one-population models. A simple model where a random boost in the magnitude of a fraction of bright normal galaxies forms a class of atypical galaxies best describes the observed distribution of BCG magnitudes. Moreover, the parameters that describe the model and the parameters of the boost have a strong physical basis.