Search

In late 2014, four images of Supernova (SN) "Refsdal," the first known example of a strongly lensed SN with multiple resolved images, were detected in the MACS J1149 galaxy-cluster field. Following the images' discovery, the SN was predicted to reappear within hundreds of days at a new position ~8 arcseconds away in the field. The observed reappearance in late 2015 makes it possible to carry out Refsdal's (1964) original proposal to use a multiply imaged SN to measure the Hubble constant H0, since the time delay between appearances should vary inversely with H0. Moreover, the position, brightness, and timing of the reappearance enable a novel test of the blind predictions of galaxy-cluster models, which are typically constrained only by the positions of multiply imaged galaxies. We have developed a new photometry pipeline that uses DOLPHOT to measure the fluxes of the five images of SN Refsdal from difference images. We apply four separate techniques to perform a blind measurement of the relative time delays and magnification ratios (mu_i/mu_1) between the last image SX and the earlier images S1-S4. We measure the relative time delay of SX-S1 to be 376.0+5.6-5.5 days and the relative magnification to be 0.30+0.05-0.03. This corresponds to a 1.5% precision on the time delay and 17% precision for the magnification ratios, and includes uncertainties due to millilensing and microlensing. In an accompanying paper, we place initial and blind constraints on the value of the Hubble constant., Comment: Published in ApJ. Companion paper presenting H0 constraints published in Science (DOI: 10.1126/science.abh1322)

The gravitationally lensed Supernova Refsdal appeared in multiple images, produced through gravitational lensing by a massive foreground galaxy cluster. After the supernova appeared in 2014, lens models of the galaxy cluster predicted an additional image of the supernova would appear in 2015, which was subsequently observed. We use the time delays between the images to perform a blinded measurement of the expansion rate of the Universe, quantified by the Hubble constant (H0). Using eight cluster lens models, we infer H0 = 64.8 +4.4-4.3 km / s / Mpc, where Mpc is the megaparsec. Using the two models most consistent with the observations, we find H0 = 66.6 +4.1-3.3 km / s / Mpc. The observations are best reproduced by models that assign dark-matter halos to individual galaxies and the overall cluster., Comment: Published in Science on May 11, 2023; this version updated to reflect minor edits to galley proofs. Companion paper presenting time-delay and relative magnification measurements published in ApJ (DOI: 10.3847/1538-4357/ac4ccb)

NASA will launch the Nancy Grace Roman Space Telescope (Roman) in the second half of this decade, which will allow for a generation-defining measurement of dark energy through multiple probes, including Type Ia supernovae (SNe Ia). To improve decisions on survey strategy, we have created the first simulations of realistic Roman images that include artificial SNe Ia injected as point sources in the images. Our analysis combines work done on Roman simulations for weak gravitational lensing studies as well as catalog-level simulations of SN samples. We have created a time series of images over two years containing $\sim$ 1,050 SNe Ia, covering a 1 square degree subarea of a planned 5 square degree deep survey. We have released these images publicly for community use along with input catalogs of all injected sources. We create secondary products from these images by generating coadded images and demonstrating recovery of transient sources using image subtraction. We perform first-use analyses on these images in order to measure galaxy-detection efficiency, point source-detection efficiency, and host-galaxy association biases. The simulated images can be found here: https://roman.ipac.caltech.edu/sims/SN_Survey_Image_sim.html., Comment: 12 pages, 12 figures. Submitted to MNRAS. For simulated images see https://roman.ipac.caltech.edu/sims/SN_Survey_Image_sim.html

Gravitationally lensed supernovae (glSNe) are of interest for time delay cosmology and SN physics. However, glSNe detections are rare, owing to the intrinsic rarity of SN explosions, the necessity of alignment with a foreground lens, and the relatively short window of detectability. We present the Las Cumbres Observatory Lensed Supernova Search, LCOLSS, a targeted survey designed for detecting glSNe in known strong lensing systems. Using cadenced r'-band imaging, LCOLSS targeted 114 galaxy-galaxy lensing systems with high expected SN rates, based on estimated star formation rates. No plausible glSN was detected by LCOLSS during our two year observing program. We carry out an analysis here to measure a detection efficiency for these observations. We then perform Monte Carlo simulations using the predicted supernova rates to determine the expected number of glSN detections. The results of the simulation suggest an expected number of detections and $68\%$ Poisson confidence intervals, $N_{SN} = 0.20, [0,2.1] $, $N_{Ia} = 0.08, [0,2.0]$, $N_{CC} = 0.12, [0,2.0]$, for all SNe, Type Ia SNe, and core-collapse (CC) SNe respectively. These results are broadly consistent with the absence of a detection in our survey. Analysis of the survey strategy can provide insights for future efforts to develop targeted glSN discovery programs, especially considering the large anticipated yields of upcoming surveys., Comment: 19 pages, 10 figures, 5 tables. Accepted by MNRAS

When the light from a distant object passes very near to a foreground galaxy or cluster, gravitational lensing can cause it to appear as multiple images on the sky. If the source is variable, it can be used to constrain the cosmic expansion rate and dark energy models. Achieving these cosmological goals requires many lensed transients with precise time delay measurements. Lensed supernovae (SN) are attractive for this purpose because they have relatively simple photometric behavior, with well-understood light curve shapes and colours $-$ in contrast to the stochastic variation of quasars. Here we report the discovery of a multiply-imaged supernova, AT2016jka ("SN Requiem"). It appeared in an evolved galaxy at $z=1.95$, gravitationally lensed by a foreground galaxy cluster. It is likely a Type Ia supernova $-$ the explosion of a low-mass stellar remnant, whose light curve can be used to measure cosmic distances. In archival Hubble Space Telescope imaging, three lensed images of the supernova are detected with relative time delays of $<$200 days. We predict a fourth image will appear close to the cluster core in the year 2037$\pm$2. Observation of the fourth image could provide a time delay precision of $\approx$7 days, $<1\%$ of the extraordinary 20 year baseline. The SN classification and the predicted reappearance time could be improved with further lens modelling and a comprehensive analysis of systematic uncertainties., Comment: Accepted for publication in a peer-reviewed journal. Main text = 6 pages, 3 figures, 1 table; Full document = 28 pages, 12 figures with Methods, Supplemental Info and references. v2: reformatted; minor corrections in SI

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., Comment: Accepted ApJS; MAST archive will be live concurrent with publication

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., Comment: 20 pages, including 14 figures

While Type Ia Supernovae (SNe Ia) are one of the most mature cosmological probes, the next era promises to be extremely exciting in the number of different ways SNe Ia are used to measure various cosmological parameters. Here we review the experiments in the 2020s that will yield orders of magnitudes more SNe Ia, and the new understandings and capabilities to constrain systematic uncertainties at a level to match these statistics. We then discuss five different cosmological probes with SNe Ia: the conventional Hubble diagram for measuring dark energy properties, the distance ladder for measuring the Hubble constant, peculiar velocities and weak lensing for measuring sigma8 and strong-lens measurements of H0 and other cosmological parameters. For each of these probes, we discuss the experiments that will provide the best measurements and also the SN Ia-related systematics that affect each one., Comment: Science White paper submitted to Astro2020 Decadal Survey. Additional co-authors added

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., Comment: submitted to Astro2020

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., Comment: 29 pages, 6 figures, submitted to ApJ. For reduced images, catalogs, lens models, and more, see relics.stsci.edu

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., Comment: 14 pages, 3 figures, 2 tables

Recently, there have been two landmark discoveries of gravitationally lensed supernovae: the first multiply-imaged SN, "Refsdal", and the first Type Ia SN resolved into multiple images, SN iPTF16geu. Fitting the multiple light curves of such objects can deliver measurements of the lensing time delays, which are the difference in arrival times for the separate images. These measurements provide precise tests of lens models or constraints on the Hubble constant and other cosmological parameters that are independent of the local distance ladder. Over the next decade, accurate time delay measurements will be needed for the tens to hundreds of lensed SNe to be found by wide-field time-domain surveys such as LSST and WFIRST. We have developed an open source software package for simulations and time delay measurements of multiply-imaged SNe, including an improved characterization of the uncertainty caused by microlensing. We describe simulations using the package that suggest a before-peak detection of the leading image enables a more accurate and precise time delay measurement (by ~1 and ~2 days, respectively), when compared to an after-peak detection. We also conclude that fitting the effects of microlensing without an accurate prior often leads to biases in the time delay measurement and over-fitting to the data, but that employing a Gaussian Process Regression (GPR) technique is sufficient for determining the uncertainty due to microlensing., Comment: Accepted 4/3 to ApJ

The most distant galaxies known are at z~10-11, observed 400-500 Myr after the Big Bang. The few z~10-11 candidates discovered to date have been exceptionally small- barely resolved, if at all, by the Hubble Space Telescope. Here we present the discovery of SPT0615-JD, a fortuitous z~10 (z_phot=9.9+/-0.6) galaxy candidate stretched into an arc over ~2.5" by the effects of strong gravitational lensing. Discovered in the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury program and companion S-RELICS Spitzer program, this candidate has a lensed H-band magnitude of 25.7+/-0.1 AB mag. With a magnification of \mu~4-7 estimated from our lens models, the de-lensed intrinsic magnitude is 27.6+/-0.3 AB mag, and the half-light radius is r_e<0.8 kpc, both consistent with other z>9 candidates. The inferred stellar mass (log [M* /M_Sun]=9.7^{+0.7}_{-0.5}) and star formation rate (\log [SFR/M_Sun yr^{-1}]=1.3^{+0.2}_{-0.3}) indicate that this candidate is a typical star-forming galaxy on the z>6 SFR-M* relation. We note that three independent lens models predict two counterimages, at least one of which should be of a similar magnitude to the arc, but these counterimages are not yet detected. Counterimages would not be expected if the arc were at lower redshift. However, the only spectral energy distributions capable of fitting the Hubble and Spitzer photometry well at lower redshifts require unphysical combinations of z~2 galaxy properties. The unprecedented lensed size of this z~10 candidate offers the potential for the James Webb Space Telescope to study the geometric and kinematic properties of a galaxy observed 500 Myr after the Big Bang., Comment: 7 pages, 4 figures. Submitted to ApJ Letters

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, Comment: Accepted to ApJ

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., Comment: 20 pages, 15 figures, submitted to ApJ. v2 edit to Table 1

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., Comment: 14 pages, 5 figures, 7 tables; submitted to ApJ

Galaxy-cluster gravitational lenses can magnify background galaxies by a total factor of up to ~50. Here we report an image of an individual star at redshift z=1.49 (dubbed "MACS J1149 Lensed Star 1 (LS1)") magnified by >2000. A separate image, detected briefly 0.26 arcseconds from LS1, is likely a counterimage of the first star demagnified for multiple years by a >~3 solar-mass object in the cluster. For reasonable assumptions about the lensing system, microlensing fluctuations in the stars' light curves can yield evidence about the mass function of intracluster stars and compact objects, including binary fractions and specific stellar evolution and supernova models. Dark-matter subhalos or massive compact objects may help to account for the two images' long-term brightness ratio., Comment: Accepted by Nature Astronomy on February 23, 2018

Using the power of gravitational lensing magnification by massive galaxy clusters, the Hubble Frontier Fields provide deep views of six patches of the high redshift Universe. The combination of deep Hubble imaging and exceptional lensing strength has revealed the greatest numbers of multiply-imaged galaxies available to constrain models of cluster mass distributions. However, even with $\mathcal{O}(100)$ images per cluster, the uncertainties associated with the reconstructions are not negligible. The goal of this paper is to show the diversity of model magnification predictions. We examine 7 and 9 mass models of Abell 2744 and MACS J0416, respectively, submitted to the Mikulski Archive for Space Telescopes for public distribution in September 2015. The dispersion between model predictions increases from 30% at common low magnifications ($\mu\sim2$) to 70% at rare high magnifications ($\mu\sim40$). MACS J0416 exhibits smaller dispersions than Abell 2744 for $2<\mu<10$. We show that magnification maps based on different lens inversion techniques typically differ from each other by more than their quoted statistical errors. This suggests that some models underestimate the true uncertainties, which are primarily due to various lensing degeneracies. Though the exact mass sheet degeneracy is broken, its generalized counterpart is not broken at least in Abell 2744. Other, local degeneracies are also present in both clusters. Our comparison of models is complementary to the comparison of reconstructions of known synthetic mass distributions. By focusing on observed clusters, we can identify those that are best constrained, and therefore provide the clearest view of the distant Universe., Comment: 17 pages, 18 figures, MNRAS, in press

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.195\%$ 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}$., Comment: 13 pages, 6 figures

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., Comment: This version corrects errors in values reported in the text for x1, c, and mB for both Type Ia SNe (erratum accepted to AJ). This error was only in the text, and did not affect any figures or conclusions

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., Comment: 21 pages, 8 figures, accepted to ApJ

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., Comment: Published in the 6 March 2015 issue of Science; 17 pages, 7 figures, and 3 tables including Supplementary Materials

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., Comment: 20 pages + 10-page appendix, 12 figures + 5 in appendix; Accepted to AJ ; latest version includes updates to discussion and figures, responding to referee and others

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., Comment: ApJ in press

Aims. We investigate the degree of improvement in dark energy constraints that can be achieved by extending Type Ia Supernova (SN Ia) samples to redshifts z > 1.5 with the Hubble Space Telescope (HST), particularly in the ongoing CANDELS and CLASH multi-cycle treasury programs. Methods. Using the popular CPL parametrization of the dark energy, w = w0 +wa(1-a), we generate mock SN Ia samples that can be projected out to higher redshifts. The synthetic datasets thus generated are fitted to the CPL model, and we evaluate the improvement that a high-z sample can add in terms of ameliorating the statistical and systematic uncertainties on cosmological parameters. Results. In an optimistic but still very achievable scenario, we find that extending the HST sample beyond CANDELS+CLASH to reach a total of 28 SN Ia at z > 1.0 could improve the uncertainty in the wa parameter by up to 21%. The corresponding improvement in the figure of merit (FoM) would be as high as 28%. Finally, we consider the use of high-redshift SN Ia samples to detect non-cosmological evolution in SN Ia luminosities with redshift, finding that such tests could be undertaken by future spacebased infrared surveys using the James Webb Space Telescope (JWST)., Comment: 9 pages, 6 figures, accepted for publication on Astronomy & Astrophysics

Past analyses of Type Ia Supernovae (SNe Ia) have identified an irreducible scatter of 5-10% in distance widely attributed to an intrinsic dispersion in luminosity. Another, equally valid, source of this scatter is intrinsic dispersion in color. Misidentification of the true source of this scatter can bias both the retrieved color-luminosity relation and cosmological parameter measurements. The size of this bias depends on the magnitude of the intrinsic color dispersion relative to the distribution of colors that correlate with distance. We produce a realistic simulation of a misattribution of intrinsic scatter, and find a negative bias in the recovered color-luminosity relation, beta, of dbeta -1.0 (~33%) and a positive bias in the equation of state parameter, w, of dw +0.04 (~4%). We re-analyze current published data sets with the assumptions that the distance scatter is predominantly the result of color. Unlike previous analyses, we find that the data are consistent with a Milky Way reddening law R_V=3.1, and that a Milky Way dust model better predicts the asymmetric color-luminosity trends than the conventional luminosity scatter hypothesis. We also determine that accounting for color variation reduces the correlation between various Host galaxy properties and Hubble residuals by ~20%., Comment: 9 pages, 5 figures. Submitted to The Astrophysical Journal

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., Comment: 10 pages, 5 figures; accepted for publication in ApJ; minor changes

We present a candidate for the most distant galaxy known to date with a photometric redshift z = 10.7 +0.6 / -0.4 (95% confidence limits; with z < 9.5 galaxies of known types ruled out at 7.2-sigma). This J-dropout Lyman Break Galaxy, named MACS0647-JD, was discovered as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). We observe three magnified images of this galaxy due to strong gravitational lensing by the galaxy cluster MACSJ0647.7+7015 at z = 0.591. The images are magnified by factors of ~8, 7, and 2, with the brighter two observed at ~26th magnitude AB (~0.15 uJy) in the WFC3/IR F160W filter (~1.4 - 1.7 um) where they are detected at >~ 12-sigma. All three images are also confidently detected at >~ 6-sigma in F140W (~1.2 - 1.6 um), dropping out of detection from 15 lower wavelength HST filters (~0.2 - 1.4 um), and lacking bright detections in Spitzer/IRAC 3.6um and 4.5um imaging (~3.2 - 5.0 um). We rule out a broad range of possible lower redshift interlopers, including some previously published as high redshift candidates. Our high redshift conclusion is more conservative than if we had neglected a Bayesian photometric redshift prior. Given CLASH observations of 17 high mass clusters to date, our discoveries of MACS0647-JD at z ~ 10.8 and MACS1149-JD1 at z ~ 9.6 are consistent with a lensed luminosity function extrapolated from lower redshifts. This would suggest that low luminosity galaxies could have reionized the universe. However given the significant uncertainties based on only two galaxies, we cannot yet rule out the sharp drop off in number counts at z >~ 10 suggested by field searches., Comment: Accepted for publication in the Astrophysical Journal. 23 pages, 18 figures

We present VLT/X-shooter observations of a high redshift, type Ia supernova host galaxy, discovered with HST/WFC3 as part of the CANDELS Supernova project. The galaxy exhibits strong emission lines of Ly{\alpha}, [O II], H{\beta}, [O III], and H{\alpha} at z = 1.54992(+0.00008-0.00004). From the emission-line fluxes and SED fitting of broad-band photometry we rule out AGN activity and characterize the host galaxy as a young, low mass, metal poor, starburst galaxy with low intrinsic extinction and high Ly{\alpha} escape fraction. The host galaxy stands out in terms of the star formation, stellar mass, and metallicity compared to its lower redshift counterparts, mainly because of its high specific star-formation rate. If valid for a larger sample of high-redshift SN Ia host galaxies, such changes in the host galaxy properties with redshift are of interest because of the potential impact on the use of SN Ia as standard candles in cosmology., Comment: 25 pages, 8 figures. Accepted for publication in ApJ

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., Comment: 10 pages, 6 figures, accepted for publication in ApJ

We present Hubble Space Telescope Wide Field Camera 3 slitless grism spectroscopy of 28 emission-line galaxies at z~2, in the GOODS-S region of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). The high sensitivity of these grism observations, with 1-sigma detections of emission lines to f > 2.5x10^{-18} erg/s/cm^2, means that the galaxies in the sample are typically ~7 times less massive (median M_* = 10^{9.5} M_sun) than previously studied z~2 emission-line galaxies. Despite their lower mass, the galaxies have OIII/Hb ratios which are very similar to previously studied z~2 galaxies and much higher than the typical emission-line ratios of local galaxies. The WFC3 grism allows for unique studies of spatial gradients in emission lines, and we stack the two-dimensional spectra of the galaxies for this purpose. In the stacked data the OIII emission line is more spatially concentrated than the Hb emission line with 98.1 confidence. We additionally stack the X-ray data (all sources are individually undetected), and find that the average L(OIII)/L(0.5-10 keV) ratio is intermediate between typical z~0 obscured active galaxies and star-forming galaxies. Together the compactness of the stacked OIII spatial profile and the stacked X-ray data suggest that at least some of these low-mass, low-metallicity galaxies harbor weak active galactic nuclei., Comment: ApJ accepted. 8 pages, 6 figures

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., Comment: Submitted to Astrophysical Journal Supplement Series; Revised version, subsequent to referee report

The young stellar cluster illuminating the W40 H II region, one of the nearest massive star forming regions, has been observed with the ACIS detector on board the Chandra X-ray Observatory. Due to its high obscuration, this is a poorly-studied stellar cluster with only a handful of bright stars visible in the optical band, including three OB stars identified as primary excitation sources. We detect 225 X-ray sources, of which 85% are confidently identified as young stellar members of the region. Two potential distances of the cluster, 260 pc and 600 pc, are used in the paper. Supposing the X-ray luminosity function to be universal, it supports a 600 pc distance as a lower limit for W40 and a total population of at least 600 stars down to 0.1 Mo under the assumption of a coeval population with a uniform obscuration. In fact, there is strong spatial variation in Ks-band-excess disk fraction and non-uniform obscuration due to a dust lane that is identified in absorption in optical, infrared and X-ray. The dust lane is likely part of a ring of material which includes the molecular core within W40. In contrast to the likely ongoing star formation in the dust lane, the molecular core is inactive. The star cluster has a spherical morphology, an isothermal sphere density profile, and mass segregation down to 1.5 Mo. However, other cluster properties, including a \leq{1} Myr age estimate and ongoing star formation, indicate that the cluster is not dynamically relaxed. X-ray diffuse emission and a powerful flare from a young stellar object are also reported., Comment: Accepted for publication in The Astrophysical Journal. 60 pages, 16 figures

The IfA Deep survey uncovered ~130 thermonuclear supernovae (TNSNe, i.e. Type Ia) candidates at redshifts from z=0.1 out to beyond z=1. The TNSN explosion rates derived from these data have been controversial, conflicting with evidence emerging from other surveys. This work revisits the IfA Deep survey to re-evaluate the photometric evidence. Applying the SOFT program to the light curves of all SN candidates, we derive new classification grades and redshift estimates. We find a volumetric rate for z~0.5 that is substantially smaller than the originally published values, bringing the revised IfA Deep rate into good agreement with other surveys. With our improved photometric analysis techniques, we are able to confidently extend the rate measurements to higher redshifts, and we find a steadily increasing TNSN rate, with no indication of a peak out to z=1.05., Comment: Accepted for publication in ApJ. 12 pages, 3 figures. Table 2 to be available in full in ApJ online, or on request to authors

We report results from the Supernova Photometric Classification Challenge (SNPCC), a publicly released mix of simulated supernovae (SNe), with types (Ia, Ibc, and II) selected in proportion to their expected rate. The simulation was realized in the griz filters of the Dark Energy Survey (DES) with realistic observing conditions (sky noise, point-spread function and atmospheric transparency) based on years of recorded conditions at the DES site. Simulations of non-Ia type SNe are based on spectroscopically confirmed light curves that include unpublished non-Ia samples donated from the Carnegie Supernova Project (CSP), the Supernova Legacy Survey (SNLS), and the Sloan Digital Sky Survey-II (SDSS-II). A spectroscopically confirmed subset was provided for training. We challenged scientists to run their classification algorithms and report a type and photo-z for each SN. Participants from 10 groups contributed 13 entries for the sample that included a host-galaxy photo-z for each SN, and 9 entries for the sample that had no redshift information. Several different classification strategies resulted in similar performance, and for all entries the performance was significantly better for the training subset than for the unconfirmed sample. For the spectroscopically unconfirmed subset, the entry with the highest average figure of merit for classifying SNe~Ia has an efficiency of 0.96 and an SN~Ia purity of 0.79. As a public resource for the future development of photometric SN classification and photo-z estimators, we have released updated simulations with improvements based on our experience from the SNPCC, added samples corresponding to the Large Synoptic Survey Telescope (LSST) and the SDSS, and provided the answer keys so that developers can evaluate their own analysis., Comment: accepted by PASP

We present deep 3.6 cm radio continuum observations of the HII region W40 obtained using the Very Large Array in its A and B configurations. We detect a total of 20 compact radio sources in a region of 4 x 4 arcmin, with 11 of them concentrated in a band with 30 arcsec of extent. We also present JHK photometry of the W40 cluster taken with the QUIRC instrument on the University of Hawaii 2.2 meter telescope. These data reveal that 15 of the 20 VLA sources have infrared counterparts, and 10 show radio variability with periods less than 20 days. Based on these combined radio and IR data, we propose that 8 of the radio sources are candidate ultracompact HII regions, 7 are likely to be young stellar objects, and 2 may be shocked interstellar gas., Comment: Accepted for publication in the Astronomical Journal. 10 pages, 4 figures, plus online table in electronic version

Wide field surveys will soon be discovering Type Ia supernovae (SNe) at rates of several thousand per year. Spectroscopic follow-up can only scratch the surface for such enormous samples, so these extensive data sets will only be useful to the extent that they can be characterized by the survey photometry alone. In a companion paper (Rodney and Tonry, 2009) we introduced the SOFT method for analyzing SNe using direct comparison to template light curves, and demonstrated its application for photometric SN classification. In this work we extend the SOFT method to derive estimates of redshift and luminosity distance for Type Ia SNe, using light curves from the SDSS and SNLS surveys as a validation set. Redshifts determined by SOFT using light curves alone are consistent with spectroscopic redshifts, showing a root-mean-square scatter in the residuals of RMS_z=0.051. SOFT can also derive simultaneous redshift and distance estimates, yielding results that are consistent with the currently favored Lambda-CDM cosmological model. When SOFT is given spectroscopic information for SN classification and redshift priors, the RMS scatter in Hubble diagram residuals is 0.18 mags for the SDSS data and 0.28 mags for the SNLS objects. Without access to any spectroscopic information, and even without any redshift priors from host galaxy photometry, SOFT can still measure reliable redshifts and distances, with an increase in the Hubble residuals to 0.37 mags for the combined SDSS and SNLS data set. Using Monte Carlo simulations we predict that SOFT will be able to improve constraints on time-variable dark energy models by a factor of 2-3 with each new generation of large-scale SN surveys., Comment: 20 pages, 7 figures, accepted to ApJ; paper 1 is arXiv:0910.3702

Modern supernova (SN) surveys are now uncovering stellar explosions at rates that far surpass what the world's spectroscopic resources can handle. In order to make full use of these SN datasets, it is necessary to use analysis methods that depend only on the survey photometry. This paper presents two methods for utilizing a set of SN light curve templates to classify SN objects. In the first case we present an updated version of the Bayesian Adaptive Template Matching program (BATM). To address some shortcomings of that strictly Bayesian approach, we introduce a method for Supernova Ontology with Fuzzy Templates (SOFT), which utilizes Fuzzy Set Theory for the definition and combination of SN light curve models. For well-sampled light curves with a modest signal to noise ratio (S/N>10), the SOFT method can correctly separate thermonuclear (Type Ia) SNe from core collapse SNe with 98% accuracy. In addition, the SOFT method has the potential to classify supernovae into sub-types, providing photometric identification of very rare or peculiar explosions. The accuracy and precision of the SOFT method is verified using Monte Carlo simulations as well as real SN light curves from the Sloan Digital Sky Survey and the SuperNova Legacy Survey. In a subsequent paper the SOFT method is extended to address the problem of parameter estimation, providing estimates of redshift, distance, and host galaxy extinction without any spectroscopy., Comment: 26 pages, 12 figures. Accepted to ApJ

The W40 complex is a nearby site of recent massive star formation composed of a dense molecular cloud adjacent to an HII region that contains an embedded OB star cluster. The HII region is beginning to blister out and break free from its envelope of molecular gas, but our line of sight to the central stars is largely obscured by intervening dust. Several bright OB stars in W40 - visible at optical, infrared, or cm wavelengths - are providing the ionizing flux that heats the HII region. The known stellar component of W40 is dominated by a small number of partly or fully embedded OB stars which have been studied at various wavelengths, but the lower mass stellar population remains largely unexamined. Despite its modest optical appearance, at 600pc W40 is one of the nearest massive star forming regions, and with a UV flux of about 1/10th of the Orion Nebula Cluster, this neglected region deserves detailed investigation., Comment: appearing in the Handbook of Star Forming Regions, Vol. II, Bo Reipurth, ed

We observed the young open cluster NGC 2301 for 14 nights in Feb. 2004 using the orthogonal transfer CCD camera (OPTIC). We used PSF shaping techniques ("square stars") during the observations allowing a larger dynamic range (4.5 magnitudes) of high photometric precision results ($\le$2 mmag) to be obtained. These results are better than similar observing campaigns using standard CCD imagers. This paper discusses our observational techniques and presents initial results for the variability statistics found in NGC 2301. Details of the variability statistics as functions of color, variability type, stellar type, and cluster location will appear in paper II.