Recent studies of Tidal Disruption Events (TDEs) have revealed unexpected correlations between the TDE rate and the large-scale properties of the host galaxies. In this review, we present the host galaxy properties of all TDE candidates known to date and quantify their distributions. We consider throughout the differences between observationally-identified types of TDEs and differences from spectroscopic control samples of galaxies. We focus here on the black hole and stellar masses of TDE host galaxies, their star formation histories and stellar populations, the concentration and morphology of the optical light, the presence of AGN activity, and the extra-galactic environment of the TDE hosts. We summarize the state of several possible explanations for the links between the TDE rate and host galaxy type. We present estimates of the TDE rate for different host galaxy types and quantify the degree to which rate enhancement in some types results in rate suppression in others. We discuss the possibilities for using TDE host galaxies to assist in identifying TDEs in upcoming large transient surveys and possibilities for TDE observations to be used to study their host galaxies.

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 last decade has seen a proliferation of mentoring programs that provide high-school students authentic research experiences. Such programs expose students to front-line research, equip them with basic research skills (including coding skills), and introduce them to scientist role models. Mentors in such programs range from undergraduate students to faculty members. Here, I describe the founding and first two years of operation of the Harvard Science Research Mentoring Program (SRMP). This program specifically recruits advanced graduate students and postdoctoral scholars to serve as mentors. By mentoring high-school students over a long timescale (September to May), early-career scientists gain hands-on experience in the skills required to advise students|skills that are often required of them in future academic positions yet seldom taught by academic institutions. Finally, I invite directors of existing and prospective SRMPs to join the Global SPHERE Network, through which directors of SRMPs around the world can share their experiences, best practices, and questions.

The light curves of Type Ia supernovae are routinely used to constrain cosmology models. Driven by radioactive decay of 56Ni, the light curves steadily decline over time, but >150 days past explosion, the near-infrared portion is poorly characterized. We report a year-long plateau in the near-infrared light curve at 150-500 days, followed by a second decline phase accompanied by a possible appearance of [Fe I] emission lines. This near-infrared plateau contrasts sharply with Type IIP plateaus and requires a new physical mechanism. We suggest a such as masking of the "near-infrared catastrophe," a predicted yet unobserved sharp light-curve decline, by scattering of ultraviolet photons to longer wavelengths. The transition off the plateau could be due to a change in the dominant ionization state of the supernova ejecta. Our results shed new light on the complex radiative transfer processes that take place in Type Ia supernovae and enhance their use as "standard candles."

We present and analyse a new tidal disruption event (TDE), AT2017eqx at redshift z=0.1089, discovered by Pan-STARRS and ATLAS. The position of the transient is consistent with the nucleus of its host galaxy; it peaks at a luminosity of $L \approx 10^{44}$ erg s$^{-1}$; and the spectrum shows a persistent blackbody temperature $T \gtrsim 20,000$ K with broad H I and He II emission. The lines are initially centered at zero velocity, but by 100 days the H I lines disappear while the He II develops a blueshift of $\gtrsim 5,000$ km s$^{-1}$. Both the early- and late-time morphologies have been seen in other TDEs, but the complete transition between them is unprecedented. The evolution can be explained by combining an extended atmosphere, undergoing slow contraction, with a wind in the polar direction becoming visible at late times. Our observations confirm that a lack of hydrogen a TDE spectrum does not indicate a stripped star, while the proposed model implies that much of the diversity in TDEs may be due to the observer viewing angle. Modelling the light curve suggests AT2017eqx resulted from the complete disruption of a solar-mass star by a black hole of $\sim 10^{6.3} M_\odot$. The host is another quiescent, Balmer-strong galaxy, though fainter and less centrally concentrated than most TDE hosts. Radio limits rule out a relativistic jet, while X-ray limits at 500 days are among the deepest for a TDE at this phase.

We present the status of the Dark Energy Spectroscopic Instrument (DESI) and its plans and opportunities for the coming decade. DESI construction and its initial five years of operations are an approved experiment of the US Department of Energy and is summarized here as context for the Astro2020 panel. Beyond 2025, DESI will require new funding to continue operations. We expect that DESI will remain one of the world's best facilities for wide-field spectroscopy throughout the decade. More about the DESI instrument and survey can be found at https://www.desi.lbl.gov.

Astrophysics lies at the crossroads of big datasets (such as the Large Synoptic Survey Telescope and Gaia), open source software to visualize and interpret high dimensional datasets (such as Glue, WorldWide Telescope, and OpenSpace), and uniquely skilled software engineers who bridge data science and research fields. At the same time, more than 4,000 planetariums across the globe immerse millions of visitors in scientific data. We have identified the potential for critical synergy across data, software, hardware, locations, and content that -- if prioritized over the next decade -- will drive discovery in astronomical research. Planetariums can and should be used for the advancement of scientific research. Current facilities such as the Hayden Planetarium in New York City, Adler Planetarium in Chicago, Morrison Planetarium in San Francisco, the Iziko Planetarium and Digital Dome Research Consortium in Cape Town, and Visualization Center C in Norrkoping are already developing software which ingests catalogs of astronomical and multi-disciplinary data critical for exploration research primarily for the purpose of creating scientific storylines for the general public. We propose a transformative model whereby scientists become the audience and explorers in planetariums, utilizing software for their own investigative purposes. In this manner, research benefits from the authentic and unique experience of data immersion contained in an environment bathed in context and equipped for collaboration. Consequently, in this white paper we argue that over the next decade the research astronomy community should partner with planetariums to create visualization-based research opportunities for the field. Realizing this vision will require new investments in software and human capital.

In the upcoming decade cadenced wide-field imaging surveys will increase the number of identified $z<0.3$ Type~Ia supernovae (SNe~Ia) from the hundreds to the hundreds of thousands. The increase in the number density and solid-angle coverage of SNe~Ia, in parallel with improvements in the standardization of their absolute magnitudes, now make them competitive probes of the growth of structure and hence of gravity. The peculiar velocity power spectrum is sensitive to the growth index $\gamma$, which captures the effect of gravity on the linear growth of structure through the relation $f=\Omega_M^\gamma$. We present the first projections for the precision in $\gamma$ for a range of realistic SN peculiar-velocity survey scenarios. In the next decade the peculiar velocities of SNe~Ia in the local $z<0.3$ Universe will provide a measure of $\gamma$ to $\pm 0.01$ precision that can definitively distinguish between General Relativity and leading models of alternative gravity.

During the next decade, gravitational waves will be observed from hundreds of binary inspiral events. When the redshifts of the host galaxies are known, these events can be used as `standard sirens', sensitive to the expansion rate of the Universe. Measurements of the Hubble constant $H_0$ from standard sirens can be done independently of other cosmological probes, and events occurring at $z<0.1$ will allow us to infer $H_0$ independently of cosmological models. The next generation of spectroscopic galaxy surveys will play a crucial role in reducing systematic uncertainties in $H_0$ from standard sirens, particularly for the numerous `dark sirens' which do not have an electromagnetic counterpart. In combination with large spectroscopic data sets, standard sirens with an EM counterpart are expected to constrain $H_0$ to $\sim 1-2\%$ precision within the next decade. This is competitive with the best estimates of $H_0$ obtained to date and will help illuminate the current tension between existing measurements. Information on the galaxies that host the gravitational wave events will also shed light on the origin and evolution of compact object binaries.

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

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.

Type Ia supernovae are understood to arise from the thermonuclear explosion of a carbon-oxygen white dwarf, yet the evolutionary mechanisms leading to such events remain unknown. Many proposed channels, including the classical single-degenerate scenario, invoke a hot, luminous evolutionary phase for the progenitor, in which it is a prodigious source of photoionizing emission. Here, we examine the environment of SN 2014J for evidence of a photoionized nebula in pre- and post-explosion [O III] $\lambda5007$ \AA\ and H$\beta$ images taken with the Hubble Space Telescope. From the absence of any extended emission, we exclude a stable nuclear-burning white dwarf at the location of SN 2014J in the last ~100,000 years, assuming a typical warm interstellar medium (ISM) particle density of 1 cm$^{-3}$. These limits greatly exceed existing X-ray constraints at temperatures typical of known supersoft sources. Significant extreme-UV/soft X-ray emission prior to explosion remains plausible for lower ISM densities (e.g., $n_{\rm ISM}\sim 0.1~\rm{cm}^{-3}$). In this case, however, any putative nebula would be even more extended, allowing deeper follow-up observations to resolve this ambiguity in the near future.

Recent works have studied the late-time light curves of Type Ia supernovae (SNe Ia) when these were older than 500 days past B-band maximum light. Of these, SN 2014J, which exploded in the nearby galaxy M82, was studied with the Advanced Camera for Surveys onboard the Hubble Space Telescope (HST) by Yang et al. Here, I report complementary photometry of SN 2014J taken with the HST Wide Field Camera 3 when it was ~360-1300 days old. My F555W measurements are consistent with the F606W measurements of Yang et al., but the F438W measurements are ~1 mag fainter than their F475W measurements. I corroborate their finding that even though SN 2014J has spatially resolved light echoes, its photometry is not contaminated by an unresolved echo. Finally, I compare the F438W and F555W light curves of SN 2014J to those of the other late-time SNe Ia observed to date and show that more intrinsically luminous SNe have slower light-curve decline rates. This is consistent with the correlation claimed by Graur et al., which was based on a comparison of pseudo-bolometric light curves. By conducting a direct comparison of the late-time light curves in the same filters, I remove any systematic uncertainties introduced by the assumptions that go into constructing the pseudo-bolometric light curves, thus strengthening the Graur et al. claim.

Late-time observations of Type Ia supernovae (SNe Ia), >900 days after explosion, have shown that this type of SN does not suffer an "IR catastrophe" at 500 days as previously predicted. Instead, several groups have observed a slow-down in the optical light curves of these SNe. A few reasons have been suggested for this slow-down, from a changing fraction of positrons reprocessed by the expanding ejecta, through a boost of energy from slow radioactive decay chains such as 57Co --> 57Fe, to atomic "freeze-out." Discovering which of these (or some other) heating mechanisms is behind the slow-down will directly impact studies of SN Ia progenitors, explosion models, and nebular-stage physics. Recently, Graur et al. (2018) suggested a possible correlation between the shape of the late-time light curves of four SNe Ia and their stretch values, which are proxies for their intrinsic luminosities. Here, we present Hubble Space Telescope observations of the SN Ia ASASSN-14lp at ~850-960 days past maximum light. With a stretch of s = 1.15 +/- 0.05, it is the most luminous normal SN Ia observed so far at these late times. We rule out contamination by light echoes and show that the late-time, optical light curve of ASASSN-14lp is flatter than that of previous SNe Ia observed at late times. This result is in line with-and strengthens-the Graur et al. (2018) correlation, but additional SNe are needed to verify it.

Many spiral galaxies have extended outer H~I disks and display low levels of star formation, inferred from the far-ultraviolet emission detected by {\it GALEX}, well beyond the optical radius. Here, we investigate the supernova (SN) rate in the outskirts of galaxies, using the largest and most homogeneous set of nearby supernovae (SNe) from the Lick Observatory Supernova Search (LOSS). While SN rates have been measured with respect to various galaxy properties, such as stellar mass and metallicity, their relative frequency in the outskirts versus the inner regions has not yet been studied. Understanding the SN rate as a function of intragalactic environment has many ramifications, including the interpretation of LIGO observations, the formation of massive stars, and the puzzlingly high velocity dispersion of the outer H~I disk. Using data from the LOSS survey, we find that the rate beyond the optical radius of spiral galaxies is $2.5 \pm 0.5$ SNe per millennium, while dwarf galaxies host $ 4.0 \pm 2.2$ SNe per millennium. The rates of core-collapse SNe (that may collapse to form the massive black holes detected by LIGO/Virgo) in the outer disks of spirals is $1.5 \pm 0.15$ SNe per millennium and in dwarf galaxies is $2.6 \pm 1.5$ SNe per millennium. Core-collapse SNe in spiral outskirts contribute $7600 \pm 1700$\,SNe\,~Gpc$^{-3}$\,yr$^{-1}$ to the volumetric rate, and dwarf galaxies have a rate of $31,000 \pm 18,000$\,SNe\,Gpc$^{-3}$\,yr$^{-1}$. The relative ratio of core-collapse to Type Ia SNe is comparable in the inner and outer parts of spirals, and in dwarf galaxies.

The late-time light curves of Type Ia supernovae (SNe Ia), observed $>900$ days after explosion, present the possibility of a new diagnostic for SN Ia progenitor and explosion models. First, however, we must discover what physical process (or combination of processes) leads to the slow-down of the late-time light curve relative to a pure $^{56}$Co decay, as observed in SNe 2011fe, 2012cg, and 2014J. We present Hubble Space Telescope observations of SN 2015F, taken $\approx 600-1040$ days past maximum light. Unlike those of the three other SNe Ia, the light curve of SN 2015F remains consistent with being powered solely by the radioactive decay of $^{56}$Co. We fit the light curves of these four SNe Ia in a consistent manner and measure possible correlations between the light curve stretch - a proxy for the intrinsic luminosity of the SN - and the parameters of the physical model used in the fit (e.g., the mass ratio of $^{56}$Co and $^{57}$Co produced in the explosion, or the time at which freeze-out sets in). We propose a new, late-time Phillips-like correlation between the stretch of the SNe and the shape of their late-time light curves, which we parametrize as the difference between their pseudo-bolometric luminosities at 600 and 900 days: $\Delta L_{900} = {\rm log}(L_{600}/L_{900})$. This model-independent correlation provides a new way to test which physical process lies behind the slow-down of SN Ia light curves $>900$ days after explosion, and, ultimately, fresh constraints on the various SN Ia progenitor and explosion models.

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.

The rate of tidal disruption events (TDEs), $R_\text{TDE}$, is predicted to depend on stellar conditions near the super-massive black hole (SMBH), which are on difficult-to-measure sub-parsec scales. We test whether $R_\text{TDE}$ depends on kpc-scale global galaxy properties, which are observable. We concentrate on stellar surface mass density, $\Sigma_{M_\star}$, and velocity dispersion, $\sigma_v$, which correlate with the stellar density and velocity dispersion of the stars around the SMBH. We consider 35 TDE candidates, with and without known X-ray emission. The hosts range from star-forming to quiescent to quiescent with strong Balmer absorption lines. The last (often with post-starburst spectra) are overrepresented in our sample by a factor of $35^{+21}_{-17}$ or $18^{+8}_{-7}$, depending on the strength of the H$\delta$ absorption line. For a subsample of hosts with homogeneous measurements, $\Sigma_{M_\star}=10^9$-$10^{10}~{\rm M_\odot / kpc^2}$, higher on average than for a volume-weighted control sample of Sloan Digital Sky Survey galaxies with similar redshifts and stellar masses. This is because: (1) most of the TDE hosts here are quiescent galaxies, which tend to have higher $\Sigma_{M_\star}$ than the star-forming galaxies that dominate the control, and (2) the star-forming hosts have higher average $\Sigma_{M_\star}$ than the star-forming control. There is also a weak suggestion that TDE hosts have lower $\sigma_v$ than for the quiescent control. Assuming that $R_{\rm TDE}\propto \Sigma_{M_\star}^\alpha \times \sigma_v^\beta$, and applying a statistical model to the TDE hosts and control sample, we estimate $\hat{\alpha}=0.9 \pm 0.2$ and $\hat{\beta}=-1.0 \pm 0.6$. This is broadly consistent with $R_\text{TDE}$ being tied to the dynamical relaxation of stars surrounding the SMBH.

Type Ia supernovae (SNe Ia) exhibit a wide diversity of peak luminosities and light curve shapes: the faintest SNe Ia are 10 times less luminous and evolve more rapidly than the brightest SNe Ia. Their differing characteristics also extend to their stellar age distributions, with fainter SNe Ia preferentially occurring in old stellar populations and vice versa. In this Letter, we quantify this SN Ia luminosity - stellar age connection using data from the Lick Observatory Supernova Search (LOSS). Our binary population synthesis calculations agree qualitatively with the observed trend in the >1 Gyr-old populations probed by LOSS if the majority of SNe Ia arise from prompt detonations of sub-Chandrasekhar mass white dwarfs (WDs) in double WD systems. Under appropriate assumptions, we show that double WD systems with less massive primaries, which yield fainter SNe Ia, interact and explode at older ages than those with more massive primaries. We find that prompt detonations in double WD systems are capable of reproducing the observed evolution of the SN Ia luminosity function, a constraint that any SN Ia progenitor scenario must confront.

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.

Recent versions of the observed cosmic star-formation history (SFH) have resolved an inconsistency with the stellar mass density history. We show that the revised SFH also scales up the delay-time distribution (DTD) of Type Ia supernovae (SNe Ia), as determined from the observed volumetric SN Ia rate history, aligning it with other field-galaxy SN Ia DTD measurements. The revised-SFH-based DTD has a $t^{-1.1 \pm 0.1}$ form and a Hubble-time-integrated production efficiency of $N/M_\star=1.3\pm0.1$ SNe Ia per $1000~{\rm M_\odot}$ of formed stellar mass. Using these revised histories and updated empirical iron yields of the various SN types, we re-derive the cosmic iron accumulation history. Core-collapse SNe and SNe Ia have contributed about equally to the total mass of iron in the Universe today. We find the track of the average cosmic gas element in the [$\alpha$/Fe] vs. [Fe/H] abundance-ratio plane. The track is broadly similar to the observed main locus of Galactic stars in this plane, indicating a Milky Way (MW) SFH similar in form to the cosmic one. We easily find a simple MW SFH that makes the track closely match this stellar locus. Galaxy clusters appear to have a higher-normalization DTD. This cluster DTD, combined with a short-burst MW SFH peaked at $z=3$, produces a track that matches remarkably well the observed "high-$\alpha$" locus of MW stars, suggesting the halo/thick-disk population has had a galaxy-cluster-like formation mode. Thus, a simple two-component SFH, combined with empirical DTDs and SN iron yields, suffices to closely reproduce the MW's stellar abundance patterns.

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 present accurate photometric redshifts for galaxies observed by the Cluster Lensing and Supernova survey with Hubble (CLASH). CLASH observed 25 massive galaxy cluster cores with the Hubble Space Telescope in 16 filters spanning 0.2 - 1.7 $\mu$m. Photometry in such crowded fields is challenging. Compared to our previously released catalogs, we make several improvements to the photometry, including smaller apertures, ICL subtraction, PSF matching, and empirically measured uncertainties. We further improve the Bayesian Photometric Redshift (BPZ) estimates by adding a redder elliptical template and by inflating the photometric uncertainties of the brightest galaxies. The resulting photometric redshift accuracies are dz/(1+z) $\sim$ 0.8\%, 1.0\%, and 2.0\% for galaxies with I-band F814W AB magnitudes $<$ 18, 20, and 23, respectively. These results are consistent with our expectations. They improve on our previously reported accuracies by a factor of 4 at the bright end and a factor of 2 at the faint end. Our new catalog includes 1257 spectroscopic redshifts, including 382 confirmed cluster members. We also provide stellar mass estimates. Finally, we include lensing magnification estimates of background galaxies based on our public lens models. Our new catalog of all 25 CLASH clusters is available via MAST. The analysis techniques developed here will be useful in other surveys of crowded fields, including the Frontier Fields and surveys carried out with J-PAS and JWST.

Most types of supernovae (SNe) have yet to be connected with their progenitor stellar systems. Here, we reanalyze the 10-year SN sample collected during 1998-2008 by the Lick Observatory Supernova Search (LOSS) in order to constrain the progenitors of SNe Ia and stripped-envelope SNe (SE SNe, i.e., SNe~IIb, Ib, Ic, and broad-lined Ic). We matched the LOSS galaxy sample with spectroscopy from the Sloan Digital Sky Survey and measured SN rates as a function of galaxy stellar mass, specific star formation rate, and oxygen abundance (metallicity). We find significant correlations between the SN rates and all three galaxy properties. The SN Ia correlations are consistent with other measurements, as well as with our previous explanation of these measurements in the form of a combination of the SN Ia delay-time distribution and the correlation between galaxy mass and age. The ratio between the SE SN and SN II rates declines significantly in low-mass galaxies. This rules out single stars as SE~SN progenitors, and is consistent with predictions from binary-system progenitor models. Using well-known galaxy scaling relations, any correlation between the rates and one of the galaxy properties examined here can be expressed as a correlation with the other two. These redundant correlations preclude us from establishing causality - that is, from ascertaining which of the galaxy properties (or their combination) is the physical driver for the difference between the SE SN and SN II rates. We outline several methods that have the potential to overcome this problem in future works.

In Paper I of this series, we showed that the ratio between stripped-envelope (SE) supernova (SN) and Type II SN rates reveals a significant SE SN deficiency in galaxies with stellar masses $\lesssim 10^{10}~{\rm M}_\odot$. Here, we test this result by splitting the volume-limited subsample of the Lick Observatory Supernova Search (LOSS) SN sample into low- and high-mass galaxies and comparing the relative rates of various SN types found in them. The LOSS volume-limited sample contains 180 SNe and SN impostors and is complete for SNe Ia out to 80 Mpc and core-collapse SNe out to 60 Mpc. All of these transients were recently reclassified by us in Shivvers et al. (2017). We find that the relative rates of some types of SNe differ between low- and high-mass galaxies: SNe Ib and Ic are underrepresented by a factor of ~3 in low-mass galaxies. These galaxies also contain the only examples of SN 1987A-like SNe in the sample and host about 9 times as many SN impostors. Normal SNe Ia seem to be ~30\% more common in low-mass galaxies, making these galaxies better sources for homogeneous SN Ia cosmology samples. The relative rates of SNe IIb are consistent in both low- and high-mass galaxies. The same is true for broad-line SNe Ic, although our sample includes only two such objects. The results presented here are in tension with a similar analysis from the Palomar Transient Factory, especially as regards SNe IIb.

We consider a novel mechanism to account for the observed distance-redshift relation. This is done by presenting a toy model for the large-scale matter distribution in a static Universe. Our model mainly concerns particles with masses far below those in the Standard Model of Particle Physics. The model is founded on three main assumptions: (1) a mass spectrum dN$_{\rm i}$/dm$_{\rm i}$ $=$ $\beta$m$_{\rm i}^{-\alpha}$ (where $\alpha$ and $\beta$ are both positive constants) for low-mass particles with m$_{\rm i}$ $\ll$ 10$^{-22}$ eV $\ll$ M$_{\rm P}$, where M$_{\rm P}$ is the Planck mass; (2) a particle mass-wavelength relation of the form $\lambda_{\rm i} =$ $\hbar$/$\delta_{\rm i}$m$_{\rm i}$c, where $\delta_{\rm i} =$ $\eta$m$_{\rm i}^{\gamma}$ and $\eta$ and $\gamma$ are both constants; and (3) For such low-mass particles, locality can only be defined on large spatial scales, comparable to or exceeding the particle wavelengths. We use our model to derive the cosmological redshift characteristic of the Standard Model of Cosmology, which becomes a gravitational redshift in our model. We compare the results of our model to empirical data and show that, in order to reproduce the sub-linear form of the observed distance-redshift relation, our model requires $\alpha$ + $\gamma$ $<$ 0. We further place our toy model in the context of the Friedmann Universe via a superposition of Einstein Universes, each with its own scale factor a$_{\rm i}$. Given the overwhelming evidence supporting an expanding Universe, we then address possible modifications to our base model that would be required to account for the available empirical constraints, including the addition of some initial expansion. Finally, we consider potentially observable distinctions between the cosmological redshift and our proposed mechanism to account for the observed distance-redshift relation.

We re-examine the classifications of supernovae (SNe) presented in the Lick Observatory Supernova Search (LOSS) volume-limited sample with a focus on the stripped-envelope SNe. The LOSS volume-limited sample, presented by Leaman et al. (2011) and Li et al. (2011b), was calibrated to provide meaningful measurements of SN rates in the local universe; the results presented therein continue to be used for comparisons to theoretical and modeling efforts. Many of the objects from the LOSS sample were originally classified based upon only a small subset of the data now available, however, and recent studies have both updated some subtype distinctions and improved our ability to perform robust classications, especially for stripped-envelope SNe. We re-examine the spectroscopic classifications of all events in the LOSS volume-limited sample (180 SNe and SN impostors) and update them if necessary. We discuss the populations of rare objects in our sample including broad-lined Type Ic SNe, Ca-rich SNe, SN 1987A-like events (we identify SN 2005io as SN 1987A-like here for the first time), and peculiar subtypes. The relative fractions of Type Ia SNe, Type II SNe, and stripped-envelope SNe in the local universe are not affected, but those of some subtypes are. Most significantly, after discussing the often unclear boundary between SNe Ib and Ic when only noisy spectra are available, we find a higher SN Ib fraction and a lower SN Ic fraction than calculated by Li et al. (2011b): spectroscopically normal SNe Ib occur in the local universe $1.7 \pm 0.9$ times more often than do normal SNe Ic.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

We use time-domain optical spectroscopy to distinguish between broad emission lines powered by accreting black holes (BHs) or stellar processes (i.e., supernovae) for 16 galaxies identified as AGN candidates by Reines \etal (2013). Our study is primarily focused on those objects with narrow emission-line ratios dominated by star formation. Based on follow-up spectra taken with the Magellan Echellette Spectrograph (MagE), the Dual Imaging Spectrograph, and the Ohio State Multi-Object Spectrograph, we find that the broad H$\alpha$ emission has faded or was ambiguous for all of the star-forming objects (14/16) over baselines ranging from 5 to 14 years. For the two objects in our follow-up sample with narrow-line AGN signatures (RGG 9 and RGG 119), we find persistent broad H$\alpha$ emission consistent with an AGN origin. Additionally, we use our MagE observations to measure stellar velocity dispersions for 15 objects in the Reines et al. (2013) sample, all with narrow-line ratios indicating the presence of an AGN. Stellar masses range from $\sim5\times10^{8}$ to $3\times10^{9}$~\msun, and we measure $\sigma_{\ast}$ ranging from $28-71~{\rm km~s^{-1}}$. These $\sigma_{\ast}$ correspond to some of the lowest-mass galaxies with optical signatures of AGN activity. We show that RGG 119, the one object which has both a measured $\sigma_{\ast}$ and persistent broad H$\alpha$ emission, falls near the extrapolation of the $\rm M_{BH}-\sigma_{\star}$ relation to the low-mass end.

Using the largest spectroscopic dataset of stripped-envelope core-collapse supernovae (stripped SNe), we present a systematic investigation of spectral properties of Type IIb SNe (SNe IIb), Type Ib SNe (SNe Ib), and Type Ic SNe (SNe Ic). Prior studies have been based on individual objects or small samples. Here, we analyze 227 spectra of 14 SNe IIb, 258 spectra of 21 SNe Ib, and 207 spectra of 17 SNe Ic based on the stripped SN dataset of Modjaz et al. (2014) and other published spectra of individual SNe. Each SN in our sample has a secure spectroscopic ID, a date of $V$-band maximum light, and multiple spectra at different phases. We analyze these spectra as a function of subtype and phase in order to improve the SN identification scheme and constrain the progenitors of different kinds of stripped SNe. By comparing spectra of SNe IIb with those of SNe Ib, we find that the strength of H$\alpha$ can be used to quantitatively differentiate between these two subtypes at all epochs. Moreover, we find a continuum in observational properties between SNe IIb and Ib. We address the question of hidden He in SNe Ic by comparing our observations with predictions from various models that either include hidden He or in which He has been burnt. Our results favor the He-free progenitor models for SNe Ic. Finally, we construct continuum-divided average spectra as a function of subtype and phase to quantify the spectral diversity of the different types of stripped SNe.

We present the first systematic investigation of spectral properties of 17 Type Ic Supernovae (SNe Ic), 10 broad-lined SNe Ic (SNe Ic-bl) without observed Gamma-Ray Bursts (GRBs) and 11 SNe Ic-bl with GRBs (SN-GRBs) as a function of time in order to probe their explosion conditions and progenitors. We analyze a total of 407 spectra, which were drawn from published spectra of individual SNe as well as from the densely time-sampled spectra data of Modjaz et al. (2014). In order to quantify the diversity of the SN spectra as a function of SN subtype, we construct average spectra of SNe Ic, SNe Ic-bl without GRBs and SNe Ic-bl with GRBs. We find that SN 1994I is not a typical SN Ic, in contrast to common belief, while the spectra of SN 1998bw/GRB 980425 are representative of mean spectra of SNe Ic-bl. We measure the ejecta absorption and width velocities using a new method described here and find that SNe Ic-bl with GRBs, on average, have quantifiably higher absorption velocities, as well as broader line widths than SNe without observed GRBs. In addition, we search for correlations between SN-GRB spectral properties and the energies of their accompanying GRBs. Finally, we show that the absence of clear He lines in optical spectra of SNe Ic-bl, and in particular of SN-GRBs, is not due to them being too smeared out due to the high velocities present in the ejecta. This implies that the progenitor stars of SN-GRBs are probably He-free, in addition to being H-free, which puts strong constraints on the stellar evolutionary paths needed to produce such SN-GRB progenitors at the observed low metallicities.

We present the open-source Python code pyMCZ that determines oxygen abundance and its distribution from strong emission lines in the standard metallicity calibrators, based on the original IDL code of Kewley & Dopita (2002) with updates from Kewley & Ellison (2008), and expanded to include more recently developed calibrators. The standard strong-line diagnostics have been used to estimate the oxygen abundance in the interstellar medium through various emission line ratios in many areas of astrophysics, including galaxy evolution and supernova host galaxy studies. We introduce a Python implementation of these methods that, through Monte Carlo sampling, better characterizes the statistical oxygen abundance confidence region including the effect due to the propagation of observational uncertainties. These uncertainties are likely to dominate the error budget in the case of distant galaxies, hosts of cosmic explosions. Given line flux measurements and their uncertainties, our code produces synthetic distributions for the oxygen abundance in up to 15 metallicity calibrators simultaneously, as well as for E(B-V), and estimates their median values and their 68% confidence regions. We test our code on emission line measurements from a sample of nearby supernova host galaxies (z < 0.15) and compare our metallicity results with those from previous methods. Our metallicity estimates are consistent with previous methods but yield smaller statistical uncertainties. Systematic uncertainties are not taken into account. We offer visualization tools to assess the spread of the oxygen abundance in the different calibrators, as well as the shape of the estimated oxygen abundance distribution in each calibrator, and develop robust metrics for determining the appropriate Monte Carlo sample size. The code is open access and open source and can be found at https://github.com/nyusngroup/pyMCZ (Abridged)

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.

Seitenzahl et al. (2009) have predicted that roughly three years after its explosion, the light we receive from a Type Ia supernova (SN Ia) will come mostly from reprocessing of electrons and X-rays emitted by the radioactive decay chain $^{57}{\rm Co}~\to~^{57}{\rm Fe}$, instead of positrons from the decay chain $^{56}{\rm Co}~\to~^{56}{\rm Fe}$ that dominates the SN light at earlier times. Using the {\it Hubble Space Telescope}, we followed the light curve of the SN Ia SN 2012cg out to $1055$ days after maximum light. Our measurements are consistent with the light curves predicted by the contribution of energy from the reprocessing of electrons and X-rays emitted by the decay of $^{57}$Co, offering evidence that $^{57}$Co is produced in SN Ia explosions. However, the data are also consistent with a light echo $\sim14$ mag fainter than SN 2012cg at peak. Assuming no light-echo contamination, the mass ratio of $^{57}$Ni and $^{56}$Ni produced by the explosion, a strong constraint on any SN Ia explosion model, is $0.043^{+0.012}_{-0.011}$, roughly twice Solar. In the context of current explosion models, this value favors a progenitor white dwarf with a mass near the Chandrasekhar limit.

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.

The progenitors of stripped-envelope supernovae (SNe Ibc) remain to be conclsuively identified, but correlations between SN rates and host-galaxy properties can constrain progenitor models. Here, we present one result from a re-analysis of the rates from the Lick Observatory Supernova Search. Galaxies with stellar masses less than $\sim 10^{10}~{\rm M_\odot}$ are less efficient at producing SNe Ibc than more massive galaxies. Any progenitor scenario must seek to explain this new observation.

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.

Using a method to discover and classify supernovae (SNe) in galaxy spectra, we detect 91 Type Ia SNe (SNe Ia) and 16 Type II SNe (SNe II) among 740,000 galaxies of all types and 215,000 star-forming galaxies without active galactic nuclei, respectively, in Data Release 9 of the Sloan Digital Sky Survey. Of these SNe, 15 SNe Ia and 8 SNe II are new discoveries reported here for the first time. We use our SN samples to measure SN rates per unit mass as a function of galaxy stellar mass, star-formation rate (SFR), and specific SFR (sSFR), as derived by the MPA-JHU Galspec pipeline. We show that correlations between SN Ia and SN II rates per unit mass and galaxy stellar mass, SFR, and sSFR can be explained by a combination of the respective SN delay-time distributions (the distributions of times that elapse between the formation of a stellar population and all ensuing SNe), the ages of the surveyed galaxies, the redshifts at which they are observed, and their star-formation histories. This model was first suggested by Kistler et al. for the SN Ia rate-mass correlation, but is expanded here to SNe II and to correlations with galaxy SFR and sSFR. Finally, we measure a volumetric SN II rate at redshift 0.075 of $R_{\rm II,V} = 0.621^{+0.197}_{-0.154}({\rm stat})^{+0.024}_{-0.063}({\rm sys}) \times 10^{-4}$ yr$^{-1}$ Mpc$^{-3}$. Assuming that SNe IIP and IIL account for 60 per cent of all core-collapse (CC) SNe, the CC SN rate is $R_{\rm CC,V} = 1.04^{+0.33}_{-0.26}({\rm stat})^{+0.04}_{-0.11}({\rm sys}) \times 10^{-4}$ yr$^{-1}$ Mpc$^{-3}$.

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.

Context: The Subaru Deep Field (SDF) Supernova Survey discovered 10 Type Ia supernovae (SNe Ia) in the redshift range 1.5<z<2.0, as determined solely from photometric redshifts of the host galaxies. However, photometric redshifts might be biased, and the SN sample could be contaminated by active galactic nuclei (AGNs). Aims: We aim to obtain the first robust redshift measurement and classification of a z > 1.5 SDF SN Ia host galaxy candidate Methods: We use the X-shooter (U-to-K-band) spectrograph on the Very Large Telescope to allow the detection of different emission lines in a wide spectral range. Results: We measure a spectroscopic redshift of 1.54563 +/- 0.00027 of hSDF0705.25, consistent with its photometric redshift of 1.552 +/- 0.018. From the strong emission-line spectrum we rule out AGN activity, thereby confirming the optical transient as a SN. The host galaxy follows the fundamental metallicity relation defined in Mannucci et al. (2010, 2011) showing that the properties of this high-redshift SN Ia host galaxy is similar to other field galaxies. Conclusions: Spectroscopic confirmation of additional SDF SN hosts would be required to confirm the cosmic SN rate evolution measured in the SDF.

We present Hubble Space Telescope (HST) imaging observations of the site of the Type-Ia supernova SN2011fe in the nearby galaxy M101, obtained about one year prior to the event, in a narrow band centred on the HeII 4686 \AA{} emission line. In a "single-degenerate" progenitor scenario, the hard photon flux from an accreting white dwarf (WD), burning hydrogen on its surface over $\sim1$ Myr should, in principle, create a HeIII Str\"{o}mgren sphere or shell surrounding the WD. Depending on the WD luminosity, the interstellar density, and the velocity of an outflow from the WD, the HeIII region could appear unresolved, extended, or as a ring, with a range of possible surface brightnesses. We find no trace of HeII 4686 \AA{} line emission in the HST data. Using simulations, we set $2\sigma$ upper limits on the HeII 4686 \AA{} luminosity of $L_{\rm HeII} < 3.4 \times 10^{34}$ erg s$^{-1}$ for a point source, corresponding to an emission region of radius $r < 1.8$ pc. The upper limit for an extended source is $L_{\rm HeII} < 1.7 \times 10^{35}$ erg s$^{-1}$, corresponding to an extended region with $r\sim11$ pc. The largest detectable shell, given an interstellar-medium density of 1 cm$^{-3}$, has a radius of $\sim6$ pc. Our results argue against the presence, within the $\sim10^5$ yr prior to the explosion, of a supersoft X-ray source of luminosity $L_{\rm bol} \ge 3 \times 10^{37}$ erg s$^{-1}$, or of a super-Eddington accreting WD that produces an outflowing wind capable of producing cavities with radii of 2-6 pc.

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.

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.

Using a method to discover and classify supernovae (SNe) in galaxy spectra, we find 90 Type Ia SNe (SNe Ia) and 10 Type II SNe among the ~700,000 galaxy spectra in the Sloan Digital Sky Survey Data Release 7 that have VESPA-derived star-formation histories (SFHs). We use the SN Ia sample to measure SN Ia rates per unit stellar mass. We confirm, at the median redshift of the sample, z = 0.1, the inverse dependence on galaxy mass of the SN Ia rate per unit mass, previously reported by Li et al. (2011b) for a local sample. We further confirm, following Kistler et al. (2011), that this relation can be explained by the combination of galaxy "downsizing" and a power-law delay-time distribution (DTD; the distribution of times that elapse between a hypothetical burst of star formation and the subsequent SN Ia explosions) with an index of -1, inherent to the double-degenerate progenitor scenario. We use the method of Maoz et al. (2011) to recover the DTD by comparing the number of SNe Ia hosted by each galaxy in our sample with the VESPA-derived SFH of the stellar population within the spectral aperture. In this galaxy sample, which is dominated by old and massive galaxies, we recover a "delayed" component to the DTD of 4.5 +/- 0.6 (statistical) +0.3 -0.5 (systematic) X 10^-14 SNe Msun^-1 yr^-1 for delays in the range > 2.4 Gyr. The mass-normalized SN Ia rate, averaged over all masses and redshifts in our galaxy sample, is R(Ia,M,z=0.1) = 0.10 +/- 0.01 (statistical) +/- 0.01 (systematic) SNuM, and the volumetric rate is R(Ia,V,z=0.1) = 0.247 +0.029 -0.026 (statistical) +0.016 -0.031 (systematic) X 10^-4 SNe yr^-1 Mpc^-3. This rate is consistent with the rates and rate evolution from other recent SN Ia surveys, which together also indicate a ~t^-1 DTD.

We examine the inner mass distribution of the relaxed galaxy cluster Abell 383 in deep 16-band HST/ACS+WFC3 imaging taken as part of the CLASH multi-cycle treasury program. Our program is designed to study the dark matter distribution in 25 massive clusters, and balances depth with a wide wavelength coverage to better identify lensed systems and generate precise photometric redshifts. This information together with the predictive strength of our strong-lensing analysis method identifies 13 new multiply-lensed images and candidates, so that a total of 27 multiple-images of 9 systems are used to tightly constrain the inner mass profile, $d\log \Sigma/d\log r\simeq -0.6\pm 0.1$ (r<160 kpc). We find consistency with the standard distance-redshift relation for the full range spanned by the lensed images, 1.01<z<6.03, with the higher redshift sources deflected through larger angles as expected. The inner mass profile derived here is consistent with the results of our independent weak-lensing analysis of wide-field Subaru images, with good agreement in the region of overlap. The overall mass profile is well fitted by an NFW profile with M_{vir}=(5.37^{+0.70}_{-0.63}\pm 0.26) x 10^{14}M_{\odot}/h and a relatively high concentration, c_{vir}=8.77^{+0.44}_{-0.42}\pm 0.23, which lies above the standard c-M relation similar to other well-studied clusters. The critical radius of Abell 383 is modest by the standards of other lensing clusters, r_{E}\simeq16\pm2\arcsec (for z_s=2.55), so the relatively large number of lensed images uncovered here with precise photometric redshifts validates our imaging strategy for the CLASH survey. In total we aim to provide similarly high-quality lensing data for 25 clusters, 20 of which are X-ray selected relaxed clusters, enabling a precise determination of the representative mass profile free from lensing bias. (ABRIDGED)

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.

We precisely constrain the inner mass profile of Abell 2261 (z=0.225) for the first time and determine this cluster is not "over-concentrated" as found previously, implying a formation time in agreement with {\Lambda}CDM expectations. These results are based on strong lensing analyses of new 16-band HST imaging obtained as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). Combining this with revised weak lensing analyses of Subaru wide field imaging with 5-band Subaru + KPNO photometry, we place tight new constraints on the halo virial mass M_vir = 2.2\pm0.2\times10^15 M\odot/h70 (within r \approx 3 Mpc/h70) and concentration c = 6.2 \pm 0.3 when assuming a spherical halo. This agrees broadly with average c(M,z) predictions from recent {\Lambda}CDM simulations which span 5 <~ <c> <~ 8. Our most significant systematic uncertainty is halo elongation along the line of sight. To estimate this, we also derive a mass profile based on archival Chandra X-ray observations and find it to be ~35% lower than our lensing-derived profile at r2500 ~ 600 kpc. Agreement can be achieved by a halo elongated with a ~2:1 axis ratio along our line of sight. For this elongated halo model, we find M_vir = 1.7\pm0.2\times10^15 M\odot/h70 and c_vir = 4.6\pm0.2, placing rough lower limits on these values. The need for halo elongation can be partially obviated by non-thermal pressure support and, perhaps entirely, by systematic errors in the X-ray mass measurements. We estimate the effect of background structures based on MMT/Hectospec spectroscopic redshifts and find these tend to lower Mvir further by ~7% and increase cvir by ~5%.

The Cluster Lensing And Supernova survey with Hubble (CLASH) is a 524-orbit multi-cycle treasury program to use the gravitational lensing properties of 25 galaxy clusters to accurately constrain their mass distributions. The survey, described in detail in this paper, will definitively establish the degree of concentration of dark matter in the cluster cores, a key prediction of CDM. The CLASH cluster sample is larger and less biased than current samples of space-based imaging studies of clusters to similar depth, as we have minimized lensing-based selection that favors systems with overly dense cores. Specifically, twenty CLASH clusters are solely X-ray selected. The X-ray selected clusters are massive (kT > 5 keV; 5 - 30 x 10^14 M_solar) and, in most cases, dynamically relaxed. Five additional clusters are included for their lensing strength (Einstein radii > 35 arcsec at z_source = 2) to further quantify the lensing bias on concentration, to yield high resolution dark matter maps, and to optimize the likelihood of finding highly magnified high-redshift (z > 7) galaxies. The high magnification, in some cases, provides angular resolutions unobtainable with any current UVOIR facility and can yield z > 7 candidates bright enough for spectroscopic follow-up. A total of 16 broadband filters, spanning the near-UV to near-IR, are employed for each 20-orbit campaign on each cluster. These data are used to measure precise (sigma_phz < 0.02(1+z)) photometric redshifts for dozens of newly discovered multiply-lensed images per cluster. Observations of each cluster are spread over 8 epochs to enable a search, primarily in the parallel fields, for Type Ia supernovae at z > 1 to improve constraints on the time dependence of the dark energy equation of state and the evolution of such supernovae in an epoch when the universe is matter dominated.

We present a strong-lensing analysis of the galaxy cluster MACS J1206.2-0847 ($z$=0.44) using UV, Optical, and IR, HST/ACS/WFC3 data taken as part of the CLASH multi-cycle treasury program, with VLT/VIMOS spectroscopy for some of the multiply-lensed arcs. The CLASH observations, combined with our mass-model, allow us to identify 47 new multiply-lensed images of 12 distant sources. These images, along with the previously known arc, span the redshift range $1\la z\la5.5$, and thus enable us to derive a detailed mass distribution and to accurately constrain, for the first time, the inner mass-profile of this cluster. We find an inner profile slope of $d\log \Sigma/d\log \theta\simeq -0.55\pm 0.1$ (in the range [1\arcsec, 53\arcsec], or $5\la r \la300$ kpc), as commonly found for relaxed and well-concentrated clusters. Using the many systems uncovered here we derive credible critical curves and Einstein radii for different source redshifts. For a source at $z_{s}\simeq2.5$, the critical curve encloses a large area with an effective Einstein radius of $\theta_{E}=28\pm3\arcsec$, and a projected mass of $1.34\pm0.15\times10^{14} M_{\odot}$. From the current understanding of structure formation in concordance cosmology, these values are relatively high for clusters at $z\sim0.5$, so that detailed studies of the inner mass distribution of clusters such as MACS J1206.2-0847 can provide stringent tests of the $\Lambda$CDM paradigm.

The type Ia supernova (SN Ia) rate, when compared to the cosmic star formation history (SFH), can be used to derive the delay-time distribution (DTD) of SNe Ia, which can distinguish among progenitor models. We present the results of a SN survey in the Subaru Deep Field (SDF). Over a period of 3 years, we have observed the SDF on 4 independent epochs with Suprime-Cam on the Subaru 8.2-m telescope, with 2 nights of exposure per epoch, in the R, i', and z' bands. We have discovered 150 SNe out to redshift z~2. Using 11 photometric bands from the observer-frame far-ultraviolet to the near-infrared, we derive photometric redshifts for the SN host galaxies (for 24 we also have spectroscopic redshifts). This information is combined with the SN photometry to determine the type and redshift distribution of the SN sample. Our final sample includes 28 SNe Ia in the range 1.0<z<1.5 and 10 in the range 1.5<z<2.0. Our SN Ia rate measurements are consistent with those derived from the HST/GOODS sample, but the overall uncertainty of our 1.5<z<2.0 measurement is a factor of 2 smaller, of 35-50%. We find that the SN Ia rate evolution levels off at 1.0<z<2.0, but shows no sign of declining. Combining our SN Ia rate measurements and those from the literature, and comparing to a wide range of possible SFHs, the best-fit DTD is a power law of the form Psi(t) ~ t^beta, with beta = -1.1 +/- 0.1 (statistical) +/- 0.17 (systematic). This result is consistent with other recent DTD measurements at various redshifts and environments, and is in agreement with a generic prediction of the double-degenerate progenitor scenario for SNe Ia. By combining the contribution to iron production from core-collapse SNe, based on the wide range of SFHs, with that from SNe Ia, calculated with the best-fit DTD, we predict that the mean present-day cosmic iron abundance is in the range Z(Fe) = (0.09-0.37) Z(Fe,Sun) (abridged).

We discuss the properties of an object in the Subaru Deep Field (SDF) classified as a galaxy in on-line data bases and revealed on the Subaru images as a genuine polar-ring galaxy (PRG) candidate. We analyse available photometric data and conclude that this object consists of a >5 Gyr old early-type central body surrounded by a faint, narrow inner ring tilted at a ~25 deg angle relative to the polar axis of the host galaxy. The halo surrounding the main stellar body exhibits a diversity of spatially extended stellar features of low surface brightness, including a faint asymmetric stellar cloud and two prominent loops. These faint features, together with the unperturbed morphology of the central host, are clear signs of a recent coalescence of two highly unequal mass galaxies, most likely a pre-existing early-type galaxy and a close-by gas-rich dwarf galaxy. The presumed stellar remnants observed near the edges of the ring, including possibly the surviving captured companion itself, indicate that the merger is still taking place.