We present a long temporal baseline variability survey in the Frontier Field MACSJ1149. In this study, we identify active galactic nuclei (AGNs) and other transient sources via their variability using over a decade of Hubble Space Telescope (HST) images for thousands of galaxies in the cluster region and detect significant variability in galaxies extending down to an apparent nuclear magnitude of m$_{i}$ $<$ 26.5. Our analysis utilizes HST images obtained in six different wavelengths from 435 nm to 1.6 microns and covers time scales 12 hours to 12 years apart. We find that $\sim$2% of galaxies in these images are variable with 49 AGN candidates and 4 new supernovae candidates detected. Half of the variables are in the cluster and these are primarily elliptical galaxies displaying variability only in the near-infrared bands. About 20% of the AGN candidates have morphologies and colors consistent with quasars, though most of the variables appear to be dominated by the host galaxy light. The structure function for these sources show a greater amplitude of variability at shorter wavelengths with slopes shallower than typical quasars. We also report a previously unknown Einstein cross identified in this field.

This document presentes general guidelines to investigators proposing ultra-rapid target of opportunity (ToO) observations with the Hubble Space Telescope (HST). Establishing clear plans in advance and communicating with STScI staff, particularly the Program Coordinator, are crucial to minimising the time between triggering a ToO and executing the observations.

We present analytical reconstructions of type Ia supernova (SN Ia) delay time distributions (DTDs) by way of two independent methods: by a Markov chain Monte Carlo best-fit technique comparing the volumetric SN Ia rate history to today's compendium cosmic star-formation history, and secondly through a maximum likelihood analysis of the star formation rate histories of individual galaxies in the GOODS/CANDELS field, in comparison to their resultant SN Ia yields. We adopt a flexible skew-normal DTD model, which could match a wide range of physically motivated DTD forms. We find a family of solutions that are essentially exponential DTDs, similar in shape to the $\beta\approx-1$ power-law DTDs, but with more delayed events (>1 Gyr in age) than prompt events (<1 Gyr). Comparing these solutions to delay time measures separately derived from field galaxies and galaxy clusters, we find the skew-normal solutions can accommodate both without requiring a different DTD form in different environments. These model fits are generally inconsistent with results from single-degenerate binary population synthesis models, and are seemingly supportive of double-degenerate progenitors for most SN Ia events.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

The Wide Field InfraRed Survey Telescope (WFIRST) was the highest ranked large space-based mission of the 2010 New Worlds, New Horizons decadal survey. It is now a NASA mission in formulation with a planned launch in the mid-2020s. A primary mission objective is to precisely constrain the nature of dark energy through multiple probes, including Type Ia supernovae (SNe Ia). Here, we present the first realistic simulations of the WFIRST SN survey based on current hardware specifications and using open-source tools. We simulate SN light curves and spectra as viewed by the WFIRST wide-field channel (WFC) imager and integral-field channel (IFC) spectrometer, respectively. We examine 11 survey strategies with different time allocations between the WFC and IFC, two of which are based upon the strategy described by the WFIRST Science Definition Team, which measures SN distances exclusively from IFC data. We propagate statistical and, crucially, systematic uncertainties to predict the Dark Energy Task Force figure of merit (FoM) for each strategy. Of the strategies investigated, we find the most successful to be WFC-focused. However, further work in constraining systematics is required to fully optimize the use of the IFC. Even without improvements to other cosmological probes, the WFIRST SN survey has the potential to increase the FoM by more than an order of magnitude from the current values. Although the survey strategies presented here have not been fully optimized, these initial investigations are an important step in the development of the final hardware design and implementation of the WFIRST mission.

Time domain and multi-messenger astrophysics are growing and important modes of observational astronomy that will help define astrophysics in the 2020s. Significant effort is being put into developing the components of a follow-up system for dynamically turning survey alerts into data. This system consists of: 1) brokers that will aggregate, classify, and filter alerts; 2) Target Observation Managers (TOMs) for prioritizing targets and managing observations and data; and 3) observatory interfaces, schedulers, and facilities along with data reduction software and science archives. These efforts need continued community support and funding in order to complete and maintain them. Many of the efforts can be community open-source software projects but they will benefit from the leadership of professional software developers. The coordination should be done by institutions that are involved in the follow-up system such as the national observatories (e.g. LSST/Gemini/NOAO Mid-scale/Community Science and Data Center) or a new MMA institute. These tools will help the community to produce the most science from new facilities and will provide new capabilities for all users of the facilities that adopt them.

Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role.

Astrophysical transients have been observed for millennia and have shaped our most basic assumptions about the Universe. In the last century, systematic searches have grown from detecting handfuls of transients per year to over 7000 in 2018 alone. As these searches have matured, we have discovered both large samples of "normal" classes and new, rare classes. Recently, a transient was the first object observed in both gravitational waves and light. Ground-based observatories, including LSST, will discover thousands of transients in the optical, but these facilities will not provide the high-fidelity near-infrared (NIR) photometry and high-resolution imaging of a space-based observatory. WFIRST can fill this gap. With its survey designed to measure the expansion history of the Universe with Type Ia supernovae, WFIRST will also discover and monitor thousands of other transients in the NIR, revealing the physics for these high-energy events. Small-scale GO programs, either as a supplement to the planned survey or as specific target-of-opportunity observations, would significantly expand the scope of transient science that can be studied with WFIRST.

Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed primarily to deliver the best and brightest high-redshift candidates from the first billion years after the Big Bang. RELICS observed 41 massive galaxy clusters with Hubble and Spitzer at 0.4-1.7um and 3.0-5.0um, respectively. We selected 21 clusters based on Planck PSZ2 mass estimates and the other 20 based on observed or inferred lensing strength. Our 188-orbit Hubble Treasury Program obtained the first high-resolution near-infrared images of these clusters to efficiently search for lensed high-redshift galaxies. We observed 46 WFC3/IR pointings (~200 arcmin^2) with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered 322 z ~ 6 - 10 candidates, including the brightest known at z ~ 6, and the most distant spatially-resolved lensed arc known at z ~ 10. Spitzer IRAC imaging (945 hours awarded, plus 100 archival) has crucially enabled us to distinguish z ~ 10 candidates from z ~ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. We delivered reduced HST images and catalogs of all clusters to the public via MAST and reduced Spitzer images via IRSA. We have also begun delivering lens models of all clusters, to be completed before the JWST GO call for proposals.

The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

Empirical models of supernova (SN) spectral energy distributions (SEDs) are widely used for SN survey simulations and photometric classifications. The existing library of SED models has excellent optical templates but limited, poorly constrained coverage of ultraviolet (UV) and infrared (IR) wavelengths. However, both regimes are critical for the design and operation of future SN surveys, particularly at IR wavelengths that will be accessible with the James Webb Space Telescope (JWST) and the Wide-Field Infrared Survey Telescope (WFIRST). We create a public repository of improved empirical SED templates using a sampling of Type Ia and core-collapse (CC) photometric light curves to extend the Type Ia parameterized SALT2 model and a set of SN Ib, SN Ic, and SN II SED templates into the UV and near-IR. We apply this new repository of extrapolated SN SED models to examine how future surveys can discriminate between CC and Type Ia SNe at UV and IR wavelengths, and present an open-source software package written in Python, SNSEDextend, that enables a user to generate their own extrapolated SEDs.

Massive foreground galaxy clusters magnify and distort the light of objects behind them, permitting a view into both the extremely distant and intrinsically faint galaxy populations. We present here the z ~ 6 - 8 candidate high-redshift galaxies from the Reionization Lensing Cluster Survey (RELICS), a Hubble and Spitzer Space Telescope survey of 41 massive galaxy clusters spanning an area of ~200 arcmin^2. These clusters were selected to be excellent lenses and we find similar high-redshift sample sizes and magnitude distributions as CLASH. We discover 321 candidate galaxies with photometric redshifts between z ~ 6 to z ~ 8, including extremely bright objects with H-band magnitudes of m_AB ~ 23 mag. As a sample, the observed (lensed) magnitudes of these galaxies are among the brightest known at z> 6, comparable to much wider, blank-field surveys. RELICS demonstrates the efficiency of using strong gravitational lenses to produce high-redshift samples in the epoch of reionization. These brightly observed galaxies are excellent targets for follow-up study with current and future observatories, including the James Webb Space Telescope.

Strong gravitational lensing by galaxy clusters magnifies background galaxies, enhancing our ability to discover statistically significant samples of galaxies at z>6, in order to constrain the high-redshift galaxy luminosity functions. Here, we present the first five lens models out of the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury Program, based on new HST WFC3/IR and ACS imaging of the clusters RXC J0142.9+4438, Abell 2537, Abell 2163, RXC J2211.7-0349, and ACT-CLJ0102-49151. The derived lensing magnification is essential for estimating the intrinsic properties of high-redshift galaxy candidates, and properly accounting for the survey volume. We report on new spectroscopic redshifts of multiply imaged lensed galaxies behind these clusters, which are used as constraints, and detail our strategy to reduce systematic uncertainties due to lack of spectroscopic information. In addition, we quantify the uncertainty on the lensing magnification due to statistical and systematic errors related to the lens modeling process, and find that in all but one cluster, the magnification is constrained to better than 20% in at least 80% of the field of view, including statistical and systematic uncertainties. The five clusters presented in this paper span the range of masses and redshifts of the clusters in the RELICS program. We find that they exhibit similar strong lensing efficiencies to the clusters targeted by the Hubble Frontier Fields within the WFC3/IR field of view. Outputs of the lens models are made available to the community through the Mikulski Archive for Space Telescopes

The Hubble Source Catalog (HSC) combines lists of sources detected on images obtained with the WFPC2, ACS and WFC3 instruments aboard the Hubble Space Telescope (HST) available in the Hubble Legacy Archive. The catalog contains time-domain information with about two million of its sources detected with the same instrument and filter in at least five HST visits. The Hubble Catalog of Variables (HCV) project aims to identify HSC sources showing significant brightness variations. A magnitude-dependent threshold in the median absolute deviation of photometric measurements (an outlier-resistant measure of lightcurve scatter) is adopted as the variability-detection statistic. It is supplemented with a cut in $\chi_{\rm red}^2$ that removes sources with large photometric errors. A pre-processing procedure involving bad image identification, outlier rejection and computation of local magnitude zero-point corrections is applied to HSC lightcurves before computing the variability detection statistic. About 52000 HSC sources are identified as candidate variables, among which 7800 show variability in more than one filter. Visual inspection suggests that $\sim 70\%$ of the candidates detected in multiple filters are true variables while the remaining $\sim 30\%$ are sources with aperture photometry corrupted by blending, imaging artifacts or image processing anomalies. The candidate variables have AB magnitudes in the range 15-27$^{m}$ with the median 22$^{m}$. Among them are the stars in our own and nearby galaxies as well as active galactic nuclei.

The Hubble Catalog of Variables (HCV) project aims to identify the variable sources in the Hubble Source Catalog (HSC), which includes about 92 million objects with over 300 million measurements detected by the WFPC2, ACS and WFC3 cameras on board of the Hubble Space Telescope (HST), by using an automated pipeline containing a set of detection and validation algorithms. All the HSC sources with more than a predefined number of measurements in a single filter/instrument combination are pre-processed to correct systematic effect and to remove the bad measurements. The corrected data are used to compute a number of variability indexes to determine the variability status of each source. The final variable source catalog will contain variables stars, active galactic nuclei (AGNs), supernovae (SNs) or even new types of variables, reaching an unprecedented depth (V$\leq$27 mag). At the end of the project, the first release of the HCV will be available at the Mikulski Archive for Space Telescopes (MAST) and the ESA Hubble Science Archives. The HCV pipeline will be deployed at the Space Telescope Science Institute (STScI) so that an updated HCV may be generated following future releases of HSC.

We present an analysis of 15 Type Ia supernovae (SNe Ia) at redshift z > 1 (9 at 1.5 < z < 2.3) recently discovered in the CANDELS and CLASH Multi-Cycle Treasury programs using WFC3 on the Hubble Space Telescope. We combine these SNe Ia with a new compilation of 1050 SNe Ia, jointly calibrated and corrected for simulated survey biases to produce accurate distance measurements. We present unbiased constraints on the expansion rate at six redshifts in the range 0.07 < z < 1.5 based only on this combined SN Ia sample. The added leverage of our new sample at z > 1.5 leads to a factor of ~3 improvement in the determination of the expansion rate at z = 1.5, reducing its uncertainty to ~20%, a measurement of H(z=1.5)/H0=2.67 (+0.83,-0.52). We then demonstrate that these six measurements alone provide a nearly identical characterization of dark energy as the full SN sample, making them an efficient compression of the SN Ia data. The new sample of SNe Ia at z > 1 usefully distinguishes between alternative cosmological models and unmodeled evolution of the SN Ia distance indicators, placing empirical limits on the latter. Finally, employing a realistic simulation of a potential WFIRST SN survey observing strategy, we forecast optimistic future constraints on the expansion rate from SNe Ia.

A massive galaxy cluster can serve as a magnifying glass for distant stellar populations, with strong gravitational lensing exposing details in the lensed background galaxies that would otherwise be undetectable. The MACS J0416.1-2403 cluster (hereafter MACS0416) is one of the most efficient lenses in the sky, and in 2014 it was observed with high-cadence imaging from the Hubble Space Telescope (HST). Here we describe two unusual transient events that appeared behind MACS0416 in a strongly lensed galaxy at redshift $z = 1.0054 \pm 0.0002$. These transients---designated HFF14Spo-NW and HFF14Spo-SE and collectively nicknamed "Spock"---were faster and fainter than any supernova (SN), but significantly more luminous than a classical nova. They reached peak luminosities of $\sim10^{41}$ erg s$^{-1}$ ($M_{\rm AB} < -14$ mag) in 5 rest-frame days, then faded below detectability in roughly the same time span. Models of the cluster lens suggest that these events may be spatially coincident at the source plane, but are most likely not temporally coincident. We find that HFF14Spo can be explained as a luminous blue variable (LBV), a recurrent nova (RN), or a pair of stellar microlensing events. To distinguish between these hypotheses will require a clarification of the positions of nearby critical curves, along with high-cadence monitoring of the field that could detect new transient episodes in the host galaxy.

We aim to construct an exceptionally deep (V ~< 27) catalog of variable objects in selected Galactic and extragalactic fields visited multiple times by the Hubble Space Telescope (HST). While HST observations of some of these fields were searched for specific types of variables before (most notably, the extragalactic Cepheids), we attempt a systematic study of the population of variable objects of all types at the magnitude range not easily accessible with ground-based telescopes. The variability timescales that can be probed range from hours to years depending on how often a particular field has been visited. For source extraction and cross-matching of sources between visits we rely on the Hubble Source Catalog which includes 10^7 objects detected with WFPC2, ACS, and WFC3 HST instruments. The lightcurves extracted from the HSC are corrected for systematic effects by applying local zero-point corrections and are screened for bad measurements. For each lightcurve we compute variability indices sensitive to a broad range of variability types. The indices characterize the overall lightcurve scatter and smoothness. Candidate variables are selected as having variability index values significantly higher than expected for objects of similar brightness in the given set of observations. The Hubble Catalog of Variables will be released in 2018.

The Proposal Auto-Categorizer and Manager (PACMan) tool was written to respond to concerns on subjective flaws and potential biases in some aspects of the proposal review process for time allocation for the {\it Hubble Space Telescope} (HST), and to partially alleviate some of the anticipated additional workload from the {\it James Webb Space Telescope} (JWST) proposal review. PACMan is essentially a mixed-method Naive Bayesian spam filtering routine, with multiple pools representing scientific categories, that utilizes the Robinson method for combining token (or word) probabilities. PACMan was trained to make similar programmatic decisions in science category sorting, panelist selection, and proposal-to-panelists assignments to those made by individuals and committees in the Science Policies Group (SPG) at Space Telescope Science Institute. Based on training from the previous cycle's proposals, PACMan made the same science category assignments for proposals in Cycle 24 as did the SPG, an average of 87\% of the time. Tests for similar science categorizations, based on training using proposals from additional cycles, show that this accuracy can be further improved, to the $>95\%$ level. This tool will be used to augment or replace key functions in the TAC review processes in future HST and JWST cycles.

In Hubble Space Telescope (HST) imaging taken on 10 November 2014, four images of supernova (SN) "Refsdal" (redshift z=1.49) appeared in an Einstein-cross--like configuration (images S1-S4) around an early-type galaxy in the cluster MACS J1149.5+2223 (z=0.54). Almost all lens models of the cluster have predicted that the SN should reappear within a year in a second host-galaxy image created by the cluster's potential. In HST observations taken on 11 December 2015, we find a new source at the predicted position of the new image of SN Refsdal approximately 8" from the previous images S1-S4. This marks the first time the appearance of a SN at a particular time and location in the sky was successfully predicted in advance! We use these data and the light curve from the first four observed images of SN Refsdal to place constraints on the relative time delay and magnification of the new image (SX), compared to images S1-S4. This enables us, for the first time, to test "blind" lens model predictions of both magnifications and time delays for a lensed SN. We find that the timing and brightness of the new image are consistent with the blind predictions of a fraction of the models. The reappearance illustrates the discriminatory power of this blind test and its utility to uncover sources of systematic uncertainty. From planned HST photometry, we expect to reach a precision of 1-2% on the time delay between S1-S4 and SX.

The Hubble Source Catalog is designed to help optimize science from the Hubble Space Telescope by combining the tens of thousands of visit-based source lists in the Hubble Legacy Archive into a single master catalog. Version 1 of the Hubble Source Catalog includes WFPC2, ACS/WFC, WFC3/UVIS, and WFC3/IR photometric data generated using SExtractor software to produce the individual source lists. The catalog includes roughly 80 million detections of 30 million objects involving 112 different detector/filter combinations, and about 160 thousand HST exposures. Source lists from Data Release 8 of the Hubble Legacy Archive are matched using an algorithm developed by Budavari & Lubow (2012). The mean photometric accuracy for the catalog as a whole is better than 0.10 mag, with relative accuracy as good as 0.02 mag in certain circumstances (e.g., bright isolated stars). The relative astrometric residuals are typically within 10 mas, with a value for the mode (i.e., most common value) of 2.3 mas. The absolute astrometric accuracy is better than $\sim$0.1 arcsec for most sources, but can be much larger for a fraction of fields that could not be matched to the PanSTARRS, SDSS, or 2MASS reference systems. In this paper we describe the database design with emphasis on those aspects that enable the users to fully exploit the catalog while avoiding common misunderstandings and potential pitfalls. We provide usage examples to illustrate some of the science capabilities and data quality characteristics, and briefly discuss plans for future improvements to the Hubble Source Catalog.

We present the first year of Hubble Space Telescope imaging of the unique supernova (SN) 'Refsdal', a gravitationally lensed SN at z=1.488$\pm$0.001 with multiple images behind the galaxy cluster MACS J1149.6+2223. The first four observed images of SN Refsdal (images S1-S4) exhibited a slow rise (over ~150 days) to reach a broad peak brightness around 20 April, 2015. Using a set of light curve templates constructed from SN 1987A-like peculiar Type II SNe, we measure time delays for the four images relative to S1 of 4$\pm$4 (for S2), 2$\pm$5 (S3), and 24$\pm$7 days (S4). The measured magnification ratios relative to S1 are 1.15$\pm$0.05 (S2), 1.01$\pm$0.04 (S3), and 0.34$\pm$0.02 (S4). None of the template light curves fully captures the photometric behavior of SN Refsdal, so we also derive complementary measurements for these parameters using polynomials to represent the intrinsic light curve shape. These more flexible fits deliver fully consistent time delays of 7$\pm$2 (S2), 0.6$\pm$3 (S3), and 27$\pm$8 days (S4). The lensing magnification ratios are similarly consistent, measured as 1.17$\pm$0.02 (S2), 1.00$\pm$0.01 (S3), and 0.38$\pm$0.02 (S4). We compare these measurements against published predictions from lens models, and find that the majority of model predictions are in very good agreement with our measurements. Finally, we discuss avenues for future improvement of time delay measurements -- both for SN Refsdal and for other strongly lensed SNe yet to come.

We have acquired Hubble Space Telescope (HST) and Very Large Telescope near-infrared spectra and images of supernova (SN) Refsdal after its discovery as an Einstein cross in Fall 2014. The HST light curve of SN Refsdal matches the distinctive, slowly rising light curves of SN 1987A-like supernovae (SNe), and we find strong evidence for a broad H-alpha P-Cygni profile in the HST grism spectrum at the redshift (z = 1.49) of the spiral host galaxy. SNe IIn, powered by circumstellar interaction, could provide a good match to the light curve of SN Refsdal, but the spectrum of a SN IIn would not show broad and strong H-alpha absorption. From the grism spectrum, we measure an H-alpha expansion velocity consistent with those of SN 1987A-like SNe at a similar phase. The luminosity, evolution, and Gaussian profile of the H-alpha emission of the WFC3 and X-shooter spectra, separated by ~2.5 months in the rest frame, provide additional evidence that supports the SN 1987A-like classification. In comparison with other examples of SN 1987A-like SNe, SN Refsdal has a blue B-V color and a high luminosity for the assumed range of potential magnifications. If SN Refsdal can be modeled as a scaled version of SN 1987A, we estimate it would have an ejecta mass of 20+-5 solar masses. The evolution of the light curve at late times will provide additional evidence about the potential existence of any substantial circumstellar material (CSM). Using MOSFIRE and X-shooter spectra, we estimate a subsolar host-galaxy metallicity (8.3+-0.1 dex and <8.4 dex, respectively) near the explosion site.

We present two supernovae (SNe) discovered with the Hubble Space Telescope (HST) in the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS), an HST multi-cycle treasury program. We classify both objects as Type Ia SNe and find redshifts of z = 1.80+-0.02 and 2.26 +0.02 -0.10, the latter of which is the highest redshift Type Ia SN yet seen. Using light curve fitting we determine luminosity distances and find that both objects are consistent with a standard Lambda-CDM cosmological model. These SNe were observed using the HST Wide Field Camera 3 infrared detector (WFC3-IR), with imaging in both wide- and medium-band filters. We demonstrate that the classification and redshift estimates are significantly improved by the inclusion of single-epoch medium-band observations. This medium-band imaging approximates a very low resolution spectrum (lambda/delta lambda ~ 100) which can isolate broad spectral absorption features that differentiate Type Ia SNe from their most common core collapse cousins. This medium-band method is also insensitive to dust extinction and (unlike grism spectroscopy) it is not affected by contamination from the SN host galaxy or other nearby sources. As such, it can provide a more efficient - though less precise - alternative to IR spectroscopy for high-z SNe.

Supernova "Refsdal," multiply imaged by cluster MACSJ1149.5+2223, represents a rare opportunity to make a true blind test of model predictions in extragalactic astronomy, on a time scale that is short compared to a human lifetime. In order to take advantage of this event, we produced seven gravitational lens models with five independent methods, based on Hubble Space Telescope (HST) Hubble Frontier Field images, along with extensive spectroscopic follow-up observations by HST, the Very Large and the Keck Telescopes. We compare the model predictions and show that they agree reasonably well with the measured time delays and magnification ratios between the known images, even though these quantities were not used as input. This agreement is encouraging, considering that the models only provide statistical uncertainties, and do not include additional sources of uncertainties such as structure along the line of sight, cosmology, and the mass sheet degeneracy. We then present the model predictions for the other appearances of SN "Refsdal." A future image will reach its peak in the first half of 2016, while another image appeared between 1994 and 2004. The past image would have been too faint to be detected in existing archival images. The future image should be approximately one third as bright as the brightest known image (i.e., H_AB~25.7 mag at peak and H_AB~26.7 mag six months before peak), and thus detectable in single-orbit HST images. We will find out soon whether our predictions are correct.

We present a compilation of UBV RIz light curves of 51 type II supernovae discovered during the course of four different surveys during 1986 to 2003: the Cerro Tololo Supernova Survey, the Calan/Tololo Supernova Program (C&T), the Supernova Optical and Infrared Survey (SOIRS), and the Carnegie Type II Supernova Survey (CATS). The photometry is based on template-subtracted images to eliminate any potential host galaxy light contamination, and calibrated from foreground stars. This work presents these photometric data, studies the color evolution using different bands, and explores the relation between the magnitude at maximum brightness and the brightness decline parameter (s) from maximum light through the end of the recombination phase. This parameter is found to be shallower for redder bands and appears to have the best correlation in the B band. In addition, it also correlates with the plateau duration, being thus shorter (longer) for larger (smaller) s values.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and Cluster Lensing And Supernova survey with Hubble (CLASH) multi-cycle treasury programs with the Hubble Space Telescope (HST) have provided new opportunities to probe the rate of core-collapse supernovae (CCSNe) at high redshift, now extending to $z\approx2.5$. Here we use a sample of approximately 44 CCSNe to determine volumetric rates, $R_{CC}$, in six redshift bins in the range $0.1<z<2.5$. Together with rates from our previous HST program, and rates from the literature, we trace a more complete history of $R_{CC}(z)$, with $R_{CC}=0.72\pm0.06$ yr$^{-1}$ Mpc$^{-3}$ 10$^{-4}$ $h_{70}^{3}$ at $z<0.08$, and increasing to $3.7^{+3.1}_{-1.6}$ yr$^{-1}$ Mpc$^{-3}$ 10$^{-4}$ $h_{70}^{3}$ to $z\approx2.0$. The statistical precision in each bin is several factors better than than the systematic error, with significant contributions from host extinction, and average peak absolute magnitudes of the assumed luminosity functions for CCSN types. Assuming negligible time delays from stellar formation to explosion, we find these composite CCSN rates to be in excellent agreement with cosmic star formation rate density (SFRs) derived largely from dust-corrected rest-frame UV emission, with a scaling factor of $k=0.0091\pm0.0017\,M^{-1}_{\odot}$, and inconsistent (to $>95\%$ confidence) with SFRs from IR luminous galaxies, or with SFR models that include simple evolution in the initial mass function over time. This scaling factor is expected if the fraction of the IMF contributing to CCSN progenitors is in the 8 to 50 $M_{\odot}$ range. It is not supportive, however, of an upper mass limit for progenitors at $<20\,M_{\odot}$.

SN HFF14Tom is a Type Ia Supernova (SN) discovered at z = 1.3457 +- 0.0001 behind the galaxy cluster Abell 2744 (z = 0.308). In a cosmology-independent analysis, we find that HFF14Tom is 0.77 +- 0.15 magnitudes brighter than unlensed Type Ia SNe at similar redshift, implying a lensing magnification of mu_obs = 2.03 +- 0.29. This observed magnification provides a rare opportunity for a direct empirical test of galaxy cluster lens models. Here we test 17 lens models, 13 of which were generated before the SN magnification was known, qualifying as pure "blind tests". The models are collectively fairly accurate: 8 of the models deliver median magnifications that are consistent with the measured mu to within 1-sigma. However, there is a subtle systematic bias: the significant disagreements all involve models overpredicting the magnification. We evaluate possible causes for this mild bias, and find no single physical or methodological explanation to account for it. We do find that model accuracy can be improved to some extent with stringent quality cuts on multiply-imaged systems, such as requiring that a large fraction have spectroscopic redshifts. In addition to testing model accuracies as we have done here, Type Ia SN magnifications could also be used as inputs for future lens models of Abell 2744 and other clusters, providing valuable constraints in regions where traditional strong- and weak-lensing information is unavailable.

In 1964, Refsdal hypothesized that a supernova whose light traversed multiple paths around a strong gravitational lens could be used to measure the rate of cosmic expansion. We report the discovery of such a system. In Hubble Space Telescope imaging, we have found four images of a single supernova forming an Einstein cross configuration around a redshift z=0.54 elliptical galaxy in the MACS J1149.6+2223 cluster. The cluster's gravitational potential also creates multiple images of the z=1.49 spiral supernova host galaxy, and a future appearance of the supernova elsewhere in the cluster field is expected. The magnifications and staggered arrivals of the supernova images probe the cosmic expansion rate, as well as the distribution of matter in the galaxy and cluster lenses.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) was a multi-cycle treasury program on the Hubble Space Telescope (HST) that surveyed a total area of ~0.25 deg^2 with ~900 HST orbits spread across 5 fields over 3 years. Within these survey images we discovered 65 supernovae (SN) of all types, out to z~2.5. We classify ~24 of these as Type Ia SN (SN Ia) based on host-galaxy redshifts and SN photometry (supplemented by grism spectroscopy of 6 SN). Here we present a measurement of the volumetric SN Ia rate as a function of redshift, reaching for the first time beyond z=2 and putting new constraints on SN Ia progenitor models. Our highest redshift bin includes detections of SN that exploded when the universe was only ~3 Gyr old and near the peak of the cosmic star-formation history. This gives the CANDELS high-redshift sample unique leverage for evaluating the fraction of SN Ia that explode promptly after formation (<500 Myr). Combining the CANDELS rates with all available SN Ia rate measurements in the literature we find that this prompt SN Ia fraction is fP=0.53 +0.09 -0.10 (stat) +0.10 -0.26 (sys), consistent with a delay time distribution that follows a simple t^{-1} power law for all times t>40 Myr. However, a mild tension is apparent between ground-based low-z surveys and space-based high-z surveys. In both CANDELS and the sister HST program CLASH, we find a low rate of SN Ia at z>1. This could be a hint that prompt progenitors are in fact relatively rare, accounting for only ~20% of all SN Ia explosions -- though further analysis and larger samples will be needed to examine that suggestion.

We report observations of three gravitationally lensed supernovae (SNe) in the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program. These objects, SN CLO12Car (z = 1.28), SN CLN12Did (z = 0.85), and SN CLA11Tib (z = 1.14), are located behind three different clusters, MACSJ1720.2+3536 (z = 0.391), RXJ1532.9+3021 (z = 0.345), and Abell 383 (z = 0.187), respectively. Each SN was detected in Hubble Space Telescope (HST) optical and infrared images. Based on photometric classification, we find that SNe CLO12Car and CLN12Did are likely to be Type Ia supernovae (SNe Ia), while the classification of SN CLA11Tib is inconclusive. Using multi-color light-curve fits to determine a standardized SN Ia luminosity distance, we infer that SN CLO12Car was approximately 1.0 +/- 0.2 mag brighter than field SNe Ia at a similar redshift and ascribe this to gravitational lens magnification. Similarly, SN CLN12Did is approximately 0.2 +/- 0.2 mag brighter than field SNe Ia. We derive independent estimates of the predicted magnification from CLASH strong+weak lensing maps of the clusters: 0.83 +/- 0.16 mag for SN CLO12Car, 0.28 +/- 0.08 mag for SN CLN12Did, and 0.43 +/- 0.11 mag for SN CLA11Tib. The two SNe Ia provide a new test of the cluster lens model predictions: we find that the magnifications based on the SN Ia brightness and those predicted by the lens maps are consistent. Our results herald the promise of future observations of samples of cluster-lensed SNe Ia (from the ground or space) to help illuminate the dark-matter distribution in clusters of galaxies, through the direct determination of absolute magnifications.

We present the supernova (SN) sample and Type-Ia SN (SN Ia) rates from the Cluster Lensing And Supernova survey with Hubble (CLASH). Using the Advanced Camera for Surveys and the Wide Field Camera 3 on the Hubble Space Telescope (HST), we have imaged 25 galaxy-cluster fields and parallel fields of non-cluster galaxies. We report a sample of 27 SNe discovered in the parallel fields. Of these SNe, ~13 are classified as SN Ia candidates, including four SN Ia candidates at redshifts z > 1.2. We measure volumetric SN Ia rates to redshift 1.8 and add the first upper limit on the SN Ia rate in the range 1.8 < z < 2.4. The results are consistent with the rates measured by the HST/GOODS and Subaru Deep Field SN surveys. We model these results together with previous measurements at z < 1 from the literature. The best-fitting SN Ia delay-time distribution (DTD; the distribution of times that elapse between a short burst of star formation and subsequent SN Ia explosions) is a power law with an index of -1.00 +0.06(0.09) -0.06(0.10) (statistical) +0.12 -0.08 (systematic), where the statistical uncertainty is a result of the 68% and 95% (in parentheses) statistical uncertainties reported for the various SN Ia rates (from this work and from the literature), and the systematic uncertainty reflects the range of possible cosmic star-formation histories. We also test DTD models produced by an assortment of published binary population synthesis (BPS) simulations. The shapes of all BPS double-degenerate DTDs are consistent with the volumetric SN Ia measurements, when the DTD models are scaled up by factors of 3-9. In contrast, all BPS single-degenerate DTDs are ruled out by the measurements at a >99% significance level.

The RCT 1.3-meter telescope, formerly known as the Kitt Peak National Observatory (KPNO) 50-inch telescope, has been refurbished as a fully robotic telescope, using an autonomous scheduler to take full advantage of the observing site without the requirement of a human presence. Here we detail the current configuration of the RCT, and present as a demonstration of its high-priority science goals, the broadband {\it UBVRI} photometric calibration of the optical facility. In summary, we find the linear color transformation and extinction corrections to be consistent with similar optical KPNO facilities, to within a photometric precision of 10% (at $1\sigma$). While there were identified instrumental errors likely adding to the overall uncertainty, associated with since-resolved issues in engineering and maintenance of the robotic facility, a preliminary verification of this calibration gave good indication that the solution is robust, perhaps to a higher precision than this initial calibration implies. The RCT has been executing regular science operations since 2009, and is largely meeting the science requirements set in its acquisition and re-design.

We present the discovery of a Type Ia supernova (SN) at redshift $z = 1.914$ from the CANDELS multi-cycle treasury program on the \textit{Hubble Space Telescope (HST)}. This SN was discovered in the infrared using the Wide-Field Camera 3, and it is the highest-redshift Type Ia SN yet observed. We classify this object as a SN\,Ia by comparing its light curve and spectrum with those of a large sample of Type Ia and core-collapse supernovae (SNe). Its apparent magnitude is consistent with that expected from the $\Lambda$CDM concordance cosmology. We discuss the use of spectral evidence for classification of $z > 1.5$ SNe\,Ia using {\it HST} grism simulations, finding that spectral data alone can frequently rule out SNe\,II, but distinguishing between SNe\,Ia and SNe\,Ib/c can require prohibitively long exposures. In such cases, a quantitative analysis of the light curve may be necessary for classification. Our photometric and spectroscopic classification methods can aid the determination of SN rates and cosmological parameters from the full high-redshift CANDELS SN sample.

We use a sample of 45 core collapse supernovae detected with the Advanced Camera for Surveys on-board the Hubble Space Telescope to derive the core collapse supernova rate in the redshift range 0.1<z<1.3. In redshift bins centered on <z>=0.39, <z>=0.73, and <z>=1.11, we find rates 3.00 {+1.28}{-0.94}{+1.04}{-0.57}, 7.39 {+1.86}{-1.52}{+3.20}{-1.60}, and 9.57 {+3.76}{-2.80}{+4.96}{-2.80}, respectively, given in units yr^{-1} Mpc^{-3} 10^{-4} h_{70}^3. The rates have been corrected for host galaxy extinction, including supernovae missed in highly dust enshrouded environments in infrared bright galaxies. The first errors represent statistical while the second are the estimated systematic errors. We perform a detailed discussion of possible sources of systematic errors and note that these start to dominate over statistical errors at z>0.5, emphasizing the need to better control the systematic effects. For example, a better understanding of the amount of dust extinction in the host galaxies and knowledge of the supernova luminosity function, in particular the fraction of faint M > -15 supernovae, is needed to better constrain the rates. When comparing our results with the core collapse supernova rate based on the star formation rate, we find a good agreement, consistent with the supernova rate following the star formation rate, as expected.

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 report the discovery of a Type Ia supernova (SNIa) at redshift z=1.55 with the infrared detector of the Wide Field Camera 3 (WFC3-IR) on the Hubble Space Telescope (HST). This object was discovered in CANDELS imaging data of the Hubble Ultra Deep Field, and followed as part of the CANDELS+CLASH Supernova project, comprising the SN search components from those two HST multi-cycle treasury programs. This is the highest redshift SNIa with direct spectroscopic evidence for classification. It is also the first SN Ia at z>1 found and followed in the infrared, providing a full light curve in rest-frame optical bands. The classification and redshift are securely defined from a combination of multi-band and multi-epoch photometry of the SN, ground-based spectroscopy of the host galaxy, and WFC3-IR grism spectroscopy of both the SN and host. This object is the first of a projected sample at z>1.5 that will be discovered by the CANDELS and CLASH programs. The full CANDELS+CLASH SN Ia sample will enable unique tests for evolutionary effects that could arise due to differences in SN Ia progenitor systems as a function of redshift. This high-z sample will also allow measurement of the SN Ia rate out to z~2, providing a complementary constraint on SN Ia progenitor models.

In this paper we investigate the properties of SN 2009E, which exploded in a relatively nearby spiral galaxy (NGC 4141) and that is probably the faintest 1987A-like supernova discovered so far. Spectroscopic observations which started about 2 months after the supernova explosion, highlight significant differences between SN 2009E and the prototypical SN 1987A. Modelling the data of SN 2009E allows us to constrain the explosion parameters and the properties of the progenitor star, and compare the inferred estimates with those available for the similar SNe 1987A and 1998A. The light curve of SN 2009E is less luminous than that of SN 1987A and the other members of this class, and the maximum light curve peak is reached at a slightly later epoch than in SN 1987A. Late-time photometric observations suggest that SN 2009E ejected about 0.04 solar masses of 56Ni, which is the smallest 56Ni mass in our sample of 1987A-like events. Modelling the observations with a radiation hydrodynamics code, we infer for SN 2009E a kinetic plus thermal energy of about 0.6 foe, an initial radius of ~7 x 10^12 cm and an ejected mass of ~19 solar masses. The photospheric spectra show a number of narrow (v~1800 km/s) metal lines, with unusually strong Ba II lines. The nebular spectrum displays narrow emission lines of H, Na I, [Ca II] and [O I], with the [O I] feature being relatively strong compared to the [Ca II] doublet. The overall spectroscopic evolution is reminiscent of that of the faint 56Ni-poor type II-plateau supernovae. This suggests that SN 2009E belongs to the low-luminosity, low 56Ni mass, low-energy tail in the distribution of the 1987A-like objects in the same manner as SN 1997D and similar events represent the faint tail in the distribution of physical properties for normal type II-plateau supernovae.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5\sigma\ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5\sigma\ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

Using the Hubble Space Telescope ACS imaging of the GOODS North and South fields during Cycles 11, 12, and 13, we derive empirical constraints on the delay-time distribution function for type Ia supernovae. We extend our previous analysis to the three-year sample of 56 SNe Ia over the range 0.2<z<1.8, using a Markov chain Monte Carlo to determine the best-fit unimodal delay-time distribution function. The test, which ultimately compares the star formation rate density history to the unbinned volumetric SN Ia rate history from the GOODS/HST-SN survey, reveals a SN Ia delay-time distribution that is tightly confined to 3-4 Gyrs (to >95% confidence). This result is difficult to resolve with any intrinsic delay-time distribution function (bimodal or otherwise), in which a substantial fraction (e.g., >10%) of events are ``prompt'', requiring less than approximately 1 Gyr to develop from formation to explosion. The result is, however, strongly motivated by the decline in the number of SNe Ia at z>1.2. Sub-samples of the HST-SN data confined to lower redshifts (z<1) show plausible delay-time distributions that are dominated by prompt events, which is more consistent with results from low-redshift supernova samples and supernova host galaxy properties. Scenarios in which a substantial fraction of z>1.2 supernovae are extraordinarily obscured by dust may partly explain the differences in low-z and high-z results. Other possible resolutions may include environmental dependencies (such as gas-phase metallicity) that affect the progenitor mechanism efficiency, especially in the early universe.

We present a comparative study of the galactic and small scale environments of gamma-ray bursts (GRB) and core collapse supernovae (CCSN). We use a sample of 34 GRB hosts at z<1.2, and a comparison sample of 58 supernova hosts located within the Great Observatories Origins Deep Survey footprint. We fit template spectra to the available photometric data, which span the range 0.45-24 micron, and extract absolute magnitudes, stellar masses and star formation rates from the resulting fits. Our results broadly corroborate previous findings, but offer significant enhancements in spectral coverage and a factor 2-3 increase in sample size. Specifically, we find that CCSN occur frequently in massive spirals (spiral fraction ~50%). In contrast GRBs occur in small, relatively low mass galaxies with high specific and surface star formation rates, and have a spiral fraction of only ~10%. A comparison of the rest frame absolute magnitudes of the GRB and CCSN sample is less conclusive than found in previous work, suggesting that while GRB hosts are typically both smaller and bluer than those of CCSN their total blue light luminosities are only slightly lower. We suggest this is likely due to rapid periods of intensified star formation activity, as indicated by the high specific star formation rates, which both create the GRB progenitors and briefly significantly enhance the host galaxy blue luminosity. Finally, our analysis of local environments of GRBs and CCSN shows that GRBs are highly concentrated on their host light, and further occur in regions of higher absolute surface luminosity than CCSN.

In the next decade Type Ia supernovae (SNe Ia) will be used to test theories predicting changes in the Dark Energy equation of state with time. Ultimately this requires a dedicated space mission like JDEM. SNe Ia are mature cosmological probes --- their limitations are well characterized, and a path to improvement is clear. Dominant systematic errors include photometric calibration, selection effects, reddening, and population-dependent differences. Building on past lessons, well-controlled new surveys are poised to make strides in these areas: the Palomar Transient Factory, Skymapper, La Silla QUEST, Pan-STARRS, the Dark Energy Survey, LSST, and JDEM. They will obviate historical calibrations and selection biases, and allow comparisons via large subsamples. Some systematics follow from our ignorance of SN Ia progenitors, which there is hope of determining with SN Ia rate studies from 0<z<4. Aside from cosmology, SNe Ia regulate galactic and cluster chemical evolution, inform stellar evolution, and are laboratories for extreme physics. Essential probes of SNe Ia in these contexts include spectroscopy from the UV to the IR, X-ray cluster and SN remnant observations, spectropolarimetry, and advanced theoretical studies. While there are an abundance of discovery facilities planned, there is a deficit of follow-up resources. Living in the systematics era demands deep understanding rather than larger statistics. NOAO ReSTAR initiative to build 2-4m telescopes would provide necessary follow-up capability. Finally, to fully exploit LSST, well-matched wide-field spectroscopic capabilities are desirable.

We present a new compilation of Type Ia supernovae (SNe Ia), a new dataset of low-redshift nearby-Hubble-flow SNe and new analysis procedures to work with these heterogeneous compilations. This ``Union'' compilation of 414 SN Ia, which reduces to 307 SNe after selection cuts, includes the recent large samples of SNe Ia from the Supernova Legacy Survey and ESSENCE Survey, the older datasets, as well as the recently extended dataset of distant supernovae observed with HST. A single, consistent and blind analysis procedure is used for all the various SN Ia subsamples, and a new procedure is implemented that consistently weights the heterogeneous data sets and rejects outliers. We present the latest results from this Union compilation and discuss the cosmological constraints from this new compilation and its combination with other cosmological measurements (CMB and BAO). The constraint we obtain from supernovae on the dark energy density is $\Omega_\Lambda= 0.713^{+0.027}_{-0.029} (stat)}^{+0.036}_{-0.039} (sys)}$, for a flat, LCDM Universe. Assuming a constant equation of state parameter, $w$, the combined constraints from SNe, BAO and CMB give $w=-0.969^{+0.059}_{-0.063}(stat)^{+0.063}_{-0.066} (sys)$. While our results are consistent with a cosmological constant, we obtain only relatively weak constraints on a $w$ that varies with redshift. In particular, the current SN data do not yet significantly constrain $w$ at $z>1$. With the addition of our new nearby Hubble-flow SNe Ia, these resulting cosmological constraints are currently the tightest available.

We use the HST ACS imaging of the two GOODS fields during Cycles 11, 12, and 13 to derive the Type Ia supernova rate in four redshift intervals in the range 0.2<z<1.8. Compared to our previous results based on Cycle 11 observations only, we have increased the coverage of the search by a factor 2.7, making the total area searched equivalent to one square degree. The sample now consists of 56 Type Ia supernovae. The rates we derive are consistent with our results based on the Cycle 11 observations. In particular, the few detected supernovae at z>1.4 supports our previous result that there is a drop in the Type Ia supernova rate at high redshift, suggesting a long time delay between the formation of the progenitor star and the explosion of the supernova. If described by a simple one-parameter model, we find a characteristic delay time of 2-3 Gyr. However, a number of recent results based on e.g., low redshift supernova samples and supernova host galaxy properties suggest that the supernova delay time distribution is bimodal. In this model, a major fraction of the Type Ia supernova rate is 'prompt' and follows the star formation rate, while a smaller fraction of the rate has a long delay time, making this channel proportional to mass. While our results are fully consistent with the bimodal model at low redshifts, the low rate we find at z>1.4 appears to contradict these results. Models that corrects for star formation hidden by dust may explain at least part of the differences. Here we discuss this possibility together with other ways to reconcile data and models.

We have discovered 21 new Type Ia supernovae (SNe Ia) with the Hubble Space Telescope (HST) and have used them to trace the history of cosmic expansion over the last 10 billion years. These objects, which include 13 spectroscopically confirmed SNe Ia at z > 1, were discovered during 14 epochs of reimaging of the GOODS fields North and South over two years with the Advanced Camera for Surveys on HST. Together with a recalibration of our previous HST-discovered SNe Ia, the full sample of 23 SNe Ia at z > 1 provides the highest-redshift sample known. Combined with previous SN Ia datasets, we measured H(z) at discrete, uncorrelated epochs, reducing the uncertainty of H(z>1) from 50% to under 20%, strengthening the evidence for a cosmic jerk--the transition from deceleration in the past to acceleration in the present. The unique leverage of the HST high-redshift SNe Ia provides the first meaningful constraint on the dark energy equation-of-state parameter at z >1. The result remains consistent with a cosmological constant (w(z)=-1), and rules out rapidly evolving dark energy (dw/dz >>1). The defining property of dark energy, its negative pressure, appears to be present at z>1, in the epoch preceding acceleration, with ~98% confidence in our primary fit. Moreover, the z>1 sample-averaged spectral energy distribution is consistent with that of the typical SN Ia over the last 10 Gyr, indicating that any spectral evolution of the properties of SNe Ia with redshift is still below our detection threshold.

We present and discuss correlations for optical and near-infrared (5500-10030 angstroms) line intensity measurements at many positions in the Crab Nebula. These correlations suggest the existence of gas produced by a range of nuclear processing, from material in which synthesis ended with the CNO-cycle, to some helium-burning and nitrogen depletion, to regions containing enriched products of oxygen-burning. The latter exhibit a gradual, linear rise of [Ni II] emission with increasing argon enrichment, whereas gas with less nuclear processing shows markedly different [Ni II] emission characteristics, including the highest derived abundances. This suggests two origins for stable, neutron-rich nickel in the nebula: a type of "alpha-rich freezeout" in the more highly processed material, and possibly removal of ions from the neutron star in other regions. In addition, the data indicate that anomalously strong observed [C I] emission comes from broad, low-ionization zones. Although the strongest He I emission could also be enhanced in similar low-ionization gas, correlations between relevant line ratios argue against that interpretation, strengthening the case for an exceptionally high helium mass fraction in some locations.

When massive stars exhaust their fuel they collapse and often produce the extraordinarily bright explosions known as core-collapse supernovae. On occasion, this stellar collapse also powers an even more brilliant relativistic explosion known as a long-duration gamma-ray burst. One would then expect that long gamma-ray bursts and core-collapse supernovae should be found in similar galactic environments. Here we show that this expectation is wrong. We find that the long gamma-ray bursts are far more concentrated on the very brightest regions of their host galaxies than are the core-collapse supernovae. Furthermore, the host galaxies of the long gamma-ray bursts are significantly fainter and more irregular than the hosts of the core-collapse supernovae. Together these results suggest that long-duration gamma-ray bursts are associated with the most massive stars and may be restricted to galaxies of limited chemical evolution. Our results directly imply that long gamma-ray bursts are relatively rare in galaxies such as our own Milky Way.

We present ground-based and Hubble Space Telescope optical observations of the X-ray flash (XRF) 020903, covering 300 days. The afterglow showed a very rapid rise in the first day, followed by a relatively slow decay in the next few days. There was a clear bump in the light curve after ~25 days, accompanied by a drastic change in the spectral energy distribution. The light curve and the spectral energy distribution are naturally interpreted as the emergence -- and subsequent decay -- of a supernova (SN), similar to SN 1998bw. At peak luminosity, the SN is estimated to be 0.8 +/- 0.1 mag fainter than SN1998bw. This argues in favor of the existence of a supernova associated with this X-ray flash. A spectrum obtained 35 days after the burst shows emission lines from the host galaxy. We use this spectrum to put an upper limit on the oxygen abundance of the host at [O/H] < -0.6 dex. We also discuss a possible trend between the softness of several bursts and the early behavior of the optical afterglow, in the sense that XRFs and X-ray rich GRBs seem to have a plateau phase or even a rising light curve. This can be naturally explained in models where XRFs are similar to GRBs but seen off the jet axis.

We present 39 nights of optical photometry, 34 nights of infrared photometry, and 4 nights of optical spectroscopy of the Type Ia SN 1999ac. This supernova was discovered two weeks before maximum light, and observations were begun shortly thereafter. At early times its spectra resembled the unusual SN 1999aa and were characterized by very high velocities in the Ca II H and K lines, but very low velocities in the Si II 6355 A line. The optical photometry showed a slow rise to peak brightness but, quite peculiarly, was followed by a more rapid decline from maximum. Thus, the B- and V-band light curves cannot be characterized by a single stretch factor. We argue that the best measure of the nature of this object is not the decline rate parameter Delta m_15 (B). The B-V colors were unusual from 30 to 90 days after maximum light in that they evolved to bluer values at a much slower rate than normal Type Ia supernovae. The spectra and bolometric light curve indicate that this event was similar to the spectroscopically peculiar slow decliner SN 1999aa.

We present a photometric search for objects with point-source components that are optically variable on timescales of weeks--months in the Hubble Ultra Deep Field (HUDF) to i'(AB)=28.0 mag. The data are split into four sub-stacks of approximately equal exposure times. Objects exhibiting the signature of optical variability are selected by studying the photometric error distribution between the four different epochs, and selecting 622 candidates as 3.0 sigma outliers from the original catalog of 4644 objects. Of these, 45 are visually confirmed as free of contamination from close neighbors or various types of image defects. Four lie within the positional error boxes of Chandra X-ray sources, and two of these are spectroscopically confirmed AGN. The photometric redshift distribution of the selected variable sample is compared to that of field galaxies, and we find that a constant fraction of ~1% of all field objects show variability over the range of 0.1<z<4.5. Combined with other recent HUDF results, as well as those of recent state-of-the-art numerical simulations, we discuss a potential link between the hierarchical merging of galaxies and the growth of AGN.

In this paper, we address whether the growth of supermassive black-holes has kept pace with the process of galaxy assembly. For this purpose, we first searched the Hubble Ultra Deep Field (HUDF) for "tadpole galaxies", which have a knot at one end and an extended tail. They appear dynamically unrelaxed -- presumably early-stage mergers -- and make up ~6% of the field galaxy population. Their redshift distribution follows that of field galaxies, indicating that -- if tadpole galaxies are indeed dynamically young -- the process of galaxy assembly generally kept up with the reservoir of field galaxies as a function of epoch. Next, we present a search for HUDF objects with point-source components that are optically variable (at the >~3.0 sigma level) on timescales of weeks--months. Among 4644 objects to i_AB=28.0 mag (10 sigma), 45 have variable point-like components, which are likely weak AGN. About 1% of all field objects show variability for 0.1 < z < 4.5, and their redshift distribution is similar to that of field galaxies. Hence supermassive black-hole growth in weak AGN likely also kept up with the process of galaxy assembly. However, the faint AGN sample has almost no overlap with the tadpole sample, which was predicted by recent hydrodynamical numerical simulations. This suggests that tadpole galaxies are early-stage mergers, which likely preceded the ``turn-on'' of the AGN component and the onset of visible point-source variability by >~1 Gyr.

The Hubble Ultra Deep Field Survey has not only provided the deepest optical and near infrared views of universe, but has enabled a search for the most distant supernovae to z~2.2. We have found four supernovae by searching spans of integrations of the Ultra Deep Field and the Ultra Deep Field Parallels taken with the Hubble Space Telescope paired with the Advanced Camera for Surveys and the Near Infrared Multi Object Spectrometer. Interestingly, none of these supernovae were at z>1.4, despite the substantially increased sensitivity per unit area to such objects over the Great Observatories Origins Deep Survey. We present the optical photometric data for the four supernovae. We also show that the low frequency of Type Ia supernovae observed at z>1.4 is statistically consistent with current estimates of the global star formation history combined with the non-trivial assembly time of SN Ia progenitors.

We have discovered an extended Lyman alpha plume associated with a compact source at redshift 5.4 in slitless spectroscopic data from the Grism ACS Program for Extragalactic Science (GRAPES) project. The spatial extent of the emission is about 6 x 1.5 kpc (1 x 0.25 arcsec). Combining our grism data and the broadband images from the Hubble UltraDeep Field (UDF) images, we find a Lyman alpha line flux of 2e-17 erg/cm2/s and surface brightness 7e-17 erg/cm2/s/arcsec2. The UDF images show diffuse continuum emission associated with the Lyman alpha plume (hereafter UDF 5225), with three embedded knots. The morphology of UDF 5225 is highly suggestive of a galaxy in assembly. It is moreover possible that the prominent Lyman alpha emission from this object is due to an active nucleus, and that we are seeing the simultaneous growth through accretion of a galaxy and its central black hole. Followup observations at higher spectral resolution could test this hypothesis.

We use a sample of 42 supernovae detected with the Advanced Camera for Surveys on-board the Hubble Space Telescope as part of the Great Observatories Origins Deep Survey to measure the rate of core collapse supernovae to z~0.7 and type Ia supernovae to z~1.6. This significantly increases the redshift range where supernova rates have been estimated from observations. The rate of core collapse supernovae can be used as an independent probe of the cosmic star formation rate. Based on the observations of 17 core collapse supernovae, we measure an increase in the core collapse supernova rate by a factor of 1.6 in the range 0.3<z<0.7, and an overall increase by a factor of 7 to z~0.7 in comparison to the local core collapse supernova rate. The increase in the rate in this redshift range in consistent with recent measurements of the star formation rate derived from UV-luminosity densities and IR datasets. Based on 25 type Ia supernovae, we find a SN Ia rate that is a factor 3-5 higher at z~1 compared to earlier estimates at lower redshifts (z<0.5), implying that the type Ia supernova rate traces a higher star formation rate at redshifts z>1 compared to low redshift. At higher redshift (z>1), we find a suggested decrease in the type Ia rate with redshift. This evolution of the Ia rate with redshift is consistent with a type Ia progenitor model where there is a substantial delay between the formation of the progenitor star and the explosion of the supernova. Assuming that the type Ia progenitor stars have initial main sequence masses 3-8 M_Sun, we find that 5-7% of the available progenitors explode as type Ia supernovae.

We present results from the Hubble Higher-z Supernova Search, the first space-based open field survey for supernovae (SNe). In cooperation with the Great Observatories Origins Deep Survey, we have used the Hubble Space Telescope with the Advanced Camera for Surveys to cover 300 square arcmin in the area of the Chandra Deep Field South and the Hubble Deep Field North on five separate search epochs (separated by 45 day intervals) to a limiting magnitude of z'=26. These deep observations have allowed us to discover 42 SNe in the redshift range 0.2 < z < 1.6. As these data span a large range in redshift, they are ideal for testing the validity of Type Ia supernova progenitor models with the distribution of expected ``delay times,'' from progenitor star formation to SN Ia explosion, and the SN rates these models predict. Through a Bayesian maximum likelihood test, we determine which delay-time models best reproduce the redshift distribution of SNe Ia discovered in this survey. We find that models that require a large fraction of ``prompt'' (less than 2 Gyr) SNe Ia poorly reproduce the observed redshift distribution and are rejected at 95% confidence. We find that Gaussian models best fit the observed data for mean delay times in the range of 3 to 4 Gyr.

We have discovered 16 Type Ia supernovae (SNe Ia) with the Hubble Space Telescope (HST) and have used them to provide the first conclusive evidence for cosmic deceleration that preceded the current epoch of cosmic acceleration. These objects, discovered during the course of the GOODS ACS Treasury program, include 6 of the 7 highest-redshift SNe Ia known, all at z>1.25, and populate the Hubble diagram in unexplored territory. The luminosity distances to these and 170 previous SNe Ia are provided. A purely kinematic interpretation of the SN Ia sample provides evidence at the > 99% confidence level for a transition from deceleration to acceleration or similarly, strong evidence for a cosmic jerk. Using a simple model of the expansion history, the transition between the two epochs is constrained to be at z=0.46 +/- 0.13. The data are consistent with the cosmic concordance model of Omega_M ~ 0.3, Omega_Lambda~0.7 (chi^2_dof=1.06), and are inconsistent with a simple model of evolution or dust as an alternative to dark energy. For a flat Universe with a cosmological constant. When combined with external flat-Universe constraints we find w=-1.02 + 0.13 - 0.19 (and $<-0.76 at the 95% confidence level) for an assumed static equation of state of dark energy, P = w\rho c^2. Joint constraints on both the recent equation of state of dark energy, $w_0$, and its time evolution, dw/dz, are a factor of ~8 more precise than its first estimate and twice as precise as those without the SNe Ia discovered with HST. Our constraints are consistent with the static nature of and value of w expected for a cosmological constant (i.e., w_0 = -1.0, dw/dz = 0), and are inconsistent with very rapid evolution of dark energy. We address consequences of evolving dark energy for the fate of the Universe.

We present near infrared photometry of the Type Ia supernova 1999ee; also, optical and infrared photometry of the Type Ia SNe 2000bh, 2000ca, and 2001ba. For SNe 1999ee and 2000bh we present the first-ever SN photometry at 1.035 microns (the Y-band). We present K-corrections which transform the infrared photometry in the observer's frame to the supernova rest frame. Using our infrared K-corrections and stretch factors derived from optical photometry, we construct JHK templates which can be used to determine the apparent magnitudes at maximum if one has some data in the window -12 to +10 d with respect to T(B_max). Following up previous work on the uniformity of V minus IR loci of Type Ia supernovae of mid-range decline rates, we present unreddened loci for slow decliners. We also discuss evidence for a continuous change of color at a given epoch as a function of decline rate.

We present discoveries of SNe Ia at z > 1 and the photometric diagnostic used to discriminate them from other types of SNe detected during the GOODS Hubble Space Telescope Treasury program with the Advanced Camera for Surveys (ACS). Photometric redshift measurements of the hosts combined with deep f606w, f775w, and f850lp imaging discriminates hydrogen-rich SNe II from SNe I at z > 1 by exploiting the ultraviolet (UV) deficit in the energy distributions of SNe Ia. This sorting allows rapid follow-up of space-based discoveries. Subsequent spectroscopy of 11 GOODS SNe Ia obtained from the ground and with the grism on ACS confirmed the reliability of our photometric screening. We present the highest-redshift spectrum of any supernova published to date, SN Ia 2002fw at z = 1.3 observed near maximum brightness with the ACS grism. The lack of UV flux for true SNe Ia provides an effective tool for our ongoing efforts to build a sample of SNe Ia at 1.2 < z < 1.8 which will extend the useful range of the magnitude-redshift relation of SNe Ia.

We present observations of the first two supernovae discovered with the recently installed Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. The supernovae were found in Wide Field Camera images of the Hubble Deep Field North taken with the F775W, F850LP, and G800L optical elements as part of the ACS guaranteed time observation program. Spectra extracted from the ACS G800L grism exposures confirm that the objects are Type Ia supernovae (SNe Ia) at redshifts z=0.47 and z=0.95. Follow-up HST observations have been conducted with ACS in F775W and F850LP and with NICMOS in the near-infrared F110W bandpass, yielding a total of 9 flux measurements in the 3 bandpasses over a period of 50 days in the observed frame. We discuss many of the important issues in doing accurate photometry with the ACS. We analyze the multi-band light curves using two different fitting methods to calibrate the supernovae luminosities and place them on the SNe Ia Hubble diagram. The resulting distances are consistent with the redshift-distance relation of the accelerating universe model, although evolving intergalactic grey dust remains as a less likely possibility. The relative ease with which these SNe Ia were found, confirmed, and monitored demonstrates the potential ACS holds for revolutionizing the field of high-redshift SNe Ia, and therefore of testing the accelerating universe cosmology and constraining the "epoch of deceleration".

SN 1999aw was discovered during the first campaign of the Nearby Galaxies Supernova Search (NGSS) project. This luminous, slow-declining (Delta m_{15} (B) = 0.81 \pm 0.03) Type Ia supernova was noteworthy in at least two respects. First, it occurred in an extremely low luminosity host galaxy that was not visible in the template images, nor in initial subsequent deep imaging. Secondly, the photometric and spectral properties of this supernova indicate that it very likely was similar to the subclass of Type Ia supernovae whose prototype is SN 1999aa. This paper presents the BVRI and J_{s} HK_{s} lightcurves of SN 1999aw (through 100 days past maximum light), as well as several epochs of optical spectra. From these data we calculate the bolometric light curve, and give estimates of the luminosity at maximum light and the initial ^{56}Ni mass. In addition, we present deep BVI images obtained recently with the Baade 6.5-meter telescope at Las Campanas Observatory which reveal the remarkably low-luminosity host galaxy.

We present UBVRIz lightcurves of the Type Ia SN 1999ee and the Type Ib/c SN 1999ex, both located in the galaxy IC 5179. SN 1999ee has an extremely well sampled lightcurve spanning from 10 days before Bmax through 53 days after peak. Near maximum we find systematic differences ~0.05 mag in photometry measured with two different telescopes, even though the photometry is reduced to the same local standards around the supernova using the specific color terms for each instrumental system. We use models for our bandpasses and spectrophotometry of SN 1999ee to derive magnitude corrections (S-corrections) and remedy this problem. This exercise demonstrates the need of accurately characterizing the instrumental system before great photometric accuracies of Type Ia supernovae can be claimed. It also shows that this effect can have important astrophysical consequences since a small systematic shift of 0.02 mag in the B-V color can introduce a 0.08 mag error in the extinction corrected peak B magnitudes of a supernova and thus lead to biased cosmological parameters. The data for the Type Ib/c SN 1999ex present us with the first ever observed shock breakout of a supernova of this class. These observations show that shock breakout occurred 18 days before Bmax and support the idea that Type Ib/c supernovae are due to core collapse of massive stars rather than thermonuclear disruption of white dwarfs.

We present optical and IR spectroscopy of the first two months of evolution of the Type II SN 1999em. We combine these data with high-quality optical/IR photometry beginning only three days after shock breakout, in order to study the performance of the ``Expanding Photosphere Method'' (EPM) in the determination of distances. With this purpose we develop a technique to measure accurate photospheric velocities by cross-correlating observed and model spectra. The application of this technique to SN 1999em shows that we can reach an average uncertainty of 11% in velocity from an individual spectrum. Our analysis shows that EPM is quite robust to the effects of dust. In particular, the distances derived from the VI filters change by only 7% when the adopted visual extinction in the host galaxy is varied by 0.45 mag. The superb time sampling of the BVIZJHK light-curves of SN 1999em permits us to study the internal consistency of EPM and test the dilution factors computed from atmosphere models for Type II plateau supernovae. We find that, in the first week since explosion, the EPM distances are up to 50% lower than the average, possibly due the presence of circumstellar material. Over the following 65 days, on the other hand, our tests lend strong credence to the atmosphere models, and confirm previous claims that EPM can produce consistent distances without having to craft specific models to each supernova. This is particularly true for the VI filters which yield distances with an internal consistency of 4%. From the whole set of BVIZJHK photometry, we obtain an average distance of 7.5+/-0.5 Mpc, where the quoted uncertainty (7%) is a conservative estimate of the internal precision of the method obtained from the analysis of the first 70 days of the supernova evolution.