Photometric classification of supernovae (SNe) is imperative as recent and upcoming optical time-domain surveys, such as the Large Synoptic Survey Telescope (LSST), overwhelm the available resources for spectrosopic follow-up. Here we develop a range of light curve classification pipelines, trained on 518 spectroscopically-classified SNe from the Pan-STARRS1 Medium-Deep Survey (PS1-MDS): 357 Type Ia, 93 Type II, 25 Type IIn, 21 Type Ibc, and 17 Type I SLSNe. We present a new parametric analytical model that can accommodate a broad range of SN light curve morphologies, including those with a plateau, and fit this model to data in four PS1 filters (griz). We test a number of feature extraction methods, data augmentation strategies, and machine learning algorithms to predict the class of each SN. Our best pipelines result in 90% average accuracy, 70% average purity, and 80% average completeness for all SN classes, with the highest success rates for Type Ia SNe and SLSNe and the lowest for Type Ibc SNe. Despite the greater complexity of our classification scheme, the purity of our Type Ia SN classification, 95%, is on par with methods developed specifically for Type Ia versus non-Type Ia binary classification. As the first of its kind, this study serves as a guide to developing and training classification algorithms for a wide range of SN types with a purely empirical training set, particularly one that is similar in its characteristics to the expected LSST main survey strategy. Future work will implement this classification pipeline on ~3000 PS1/MDS light curves that lack spectroscopic classification.

We present optical light curves, redshifts, and classifications for 365 spectroscopically confirmed Type Ia supernovae (SNe Ia) discovered by the Pan-STARRS1 (PS1) Medium Deep Survey. We detail improvements to the PS1 SN photometry, astrometry and calibration that reduce the systematic uncertainties in the PS1 SN Ia distances. We combine the subset of 279 PS1 SN Ia ($0.03 < z < 0.68$) with useful distance estimates of SN Ia from SDSS, SNLS, various low-z and HST samples to form the largest combined sample of SN Ia consisting of a total of 1048 SN Ia ranging from $0.01 < z < 2.3$, which we call the `Pantheon Sample'. When combining Planck 2015 CMB measurements with the Pantheon SN sample, we find $\Omega_m=0.307\pm0.012$ and $w = -1.026\pm0.041$ for the wCDM model. When the SN and CMB constraints are combined with constraints from BAO and local H0 measurements, the analysis yields the most precise measurement of dark energy to date: $w0 = -1.007\pm 0.089$ and $wa = -0.222 \pm0.407$ for the w0waCDM model. Tension with a cosmological constant previously seen in an analysis of PS1 and low-z SNe has diminished after an increase of $2\times$ in the statistics of the PS1 sample, improved calibration and photometry, and stricter light-curve quality cuts. We find the systematic uncertainties in our measurements of dark energy are almost as large as the statistical uncertainties, primarily due to limitations of modeling the low-redshift sample. This must be addressed for future progress in using SN Ia to measure dark energy.

We present light curves and classification spectra of 17 hydrogen-poor superluminous supernovae (SLSNe) from the Pan-STARRS1 Medium Deep Survey (PS1 MDS). Our sample contains all objects from the PS1 MDS sample with spectroscopic classification that are similar to either of the prototypes SN2005ap or SN2007bi, without an explicit limit on luminosity. With a redshift range $0.3 < z < 1.6$, PS1MDS is the first SLSN sample primarily probing the high-redshift population; our multi-filter PS1 light curves probe the rest-frame UV emission, and hence the peak of the spectral energy distribution. We measure the temperature evolution and construct bolometric light curves, and find peak luminosities of $(0.5-5) \times 10^{44}$ erg s$^{-1}$ and lower limits on the total radiated energies of $(0.3-2) \times 10^{51}$ erg. The light curve shapes are diverse, with both rise- and decline times spanning a factor of $\sim 5$, and several examples of double-peaked light curves. When correcting for the flux-limited nature of our survey, we find a median peak luminosity at 4000 {\AA} of $M_{\rm 4000} = -21.1$ mag, and a spread of $\sigma = 0.7$ mag.

The progenitor systems of the class of "Ca-rich transients" is a key open issue in time domain astrophysics. These intriguing objects exhibit unusually strong calcium line emissions months after explosion, fall within an intermediate luminosity range, are often found at large projected distances from their host galaxies, and may play a vital role in enriching galaxies and the intergalactic medium. Here we present multi-wavelength observations of iPTF15eqv in NGC 3430, which exhibits a unique combination of properties that bridge those observed in Ca-rich transients and Type Ib/c supernovae. iPTF15eqv has among the highest [Ca II]/[O I] emission line ratios observed to date, yet is more luminous and decays more slowly than other Ca-rich transients. Optical and near-infrared photometry and spectroscopy reveal signatures consistent with the supernova explosion of a < 10 solar mass star that was stripped of its H-rich envelope via binary interaction. Distinct chemical abundances and ejecta kinematics suggest that the core collapse occurred through electron capture processes. Deep limits on possible radio emission made with the Jansky Very Large Array imply a clean environment ($n <$ 0.1 cm$^{-3}$) within a radius of $\sim 10^{17}$ cm. Chandra X-ray Observatory observations rule out alternative scenarios involving tidal disruption of a white dwarf by a black hole, for masses > 100 solar masses). Our results challenge the notion that spectroscopically classified Ca-rich transients only originate from white dwarf progenitor systems, complicate the view that they are all associated with large ejection velocities, and indicate that their chemical abundances may vary widely between events.

This paper presents a large amount of observations for the $\delta$ Scuti star AN Lyncis carried out in 2001--2012. The extensive observations include two tri-continent campaigns coordinated in 2002 and 2011, respectively, and several single-site contributions throughout the period. The data in total have more than 104100 raw CCD frames and photoelectric records, which consist of 165 nights (about 816 hours) spanning over 3778 days. The final reduced light curves have more than 26500 data points (including those 3462 unpublished BYU data), from which we determined 306 new timings of maximum light. A time-dependent behaviour study based on all available data indicates cyclic amplitude variability as well as period change [for the main periodicity]. Orbital sinusoid fittings to $(O-C)$ residuals and pulsation amplitudes may account for their variations being caused by the light-time effect of AN Lyn in a binary system. The orbital period is about 26--30 years. Current results support the binarity of AN Lyn, first suspected by Zhou (2002). We further show the detailed time evolution structure of the pulsation of AN Lyn as function of both time and period through wavelet analyses.

Diffuse interstellar bands (DIBs) are absorption features commonly observed in optical/near-infrared spectra of stars and thought to be associated with polyatomic molecules that comprise a significant reservoir of organic material in the universe. However, the central wavelengths of almost all DIBs do not correspond with electronic transitions of known atomic or molecular species and the specific physical nature of their carriers remains inconclusive despite decades of observational, theoretical, and experimental research. It is well established that DIB carriers are located in the interstellar medium, but the recent discovery of time-varying DIBs in the spectra of the extragalactic supernova SN2012ap suggests that some may be created in massive star environments. Here we report evidence of short time-scale (10-60 d) changes in DIB absorption line substructure toward 3 of 17 massive stars observed as part of a pathfinder survey of variable DIBs conducted with the 1.5-m Tillinghast telescope and Tillinghast Reflector Echelle Spectrograph (TRES) at Fred L. Whipple Observatory. The detections are made in high-resolution optical spectra (R = 44000) having signal-to-noise ratios of 5-15 around the 5797 and 6614 angstrom features, and are considered significant but requiring further investigation. We find that these changes are potentially consistent with interactions between stellar winds and DIB carriers in close proximity. Our findings motivate a larger survey to further characterize these variations and may establish a powerful new method for probing the poorly understood physical characteristics of DIB carriers.

A primary goal of many undergraduate- and graduate-level courses in the physical sciences is to prepare students to engage in scientific research, or to prepare students for careers that leverage skillsets similar to those used by research scientists. Even for students who may not intend to pursue a career with these characteristics, exposure to the context of applications in modern research can be a valuable tool for teaching and learning. However, a persistent barrier to student participation in research is familiarity with the technical language, format, and context that academic researchers use to communicate research methods and findings with each other: the literature of the field. Astrobites, an online web resource authored by graduate students, has published brief and accessible summaries of more than 1,300 articles from the astrophysical literature since its founding in 2010. This article presents three methods for introducing students at all levels within the formal higher education setting to approaches and results from modern research. For each method, we provide a sample lesson plan that integrates content and principles from Astrobites, including step-by-step instructions for instructors, suggestions for adapting the lesson to different class levels across the undergraduate and graduate spectrum, sample student handouts, and a grading rubric.

Super-luminous supernovae that radiate more than 10^44 ergs per second at their peak luminosity have recently been discovered in faint galaxies at redshifts of 0.1-4. Some evolve slowly, resembling models of 'pair-instability' supernovae. Such models involve stars with original masses 140-260 times that of the Sun that now have carbon-oxygen cores of 65-30 solar masses. In these stars, the photons that prevent gravitational collapse are converted to electron-positron pairs, causing rapid contraction and thermonuclear explosions. Many solar masses of 56Ni are synthesized; this isotope decays to 56Fe via 56Co, powering bright light curves. Such massive progenitors are expected to have formed from metal-poor gas in the early Universe. Recently, supernova 2007bi in a galaxy at redshift 0.127 (about 12 billion years after the Big Bang) with a metallicity one-third that of the Sun was observed to look like a fading pair-instability supernova. Here we report observations of two slow-to-fade super-luminous supernovae that show relatively fast rise times and blue colours, which are incompatible with pair-instability models. Their late-time light-curve and spectral similarities to supernova 2007bi call the nature of that event into question. Our early spectra closely resemble typical fast-declining super-luminous supernovae, which are not powered by radioactivity. Modelling our observations with 10-16 solar masses of magnetar-energized ejecta demonstrates the possibility of a common explosion mechanism. The lack of unambiguous nearby pair-instability events suggests that their local rate of occurrence is less than 6x10^-6 times that of the core-collapse rate.

SN 2012dn is a super-Chandrasekhar mass candidate in a purportedly normal spiral (SAcd) galaxy, and poses a challenge for theories of type Ia supernova diversity. Here we utilize the fast and highly parameterized spectrum synthesis tool, SYNAPPS, to estimate relative expansion velocities of species inferred from optical spectra obtained with six facilities. As with previous studies of normal SN Ia, we find that both unburned carbon and intermediate mass elements are spatially coincident within the ejecta near and below 14,000 km/s. Although the upper limit on SN 2012dn's peak luminosity is comparable to some of the most luminous normal SN Ia, we find a progenitor mass exceeding ~1.6 Msun is not strongly favored by leading merger models since these models do not accurately predict spectroscopic observations of SN 2012dn and more normal events. In addition, a comparison of light curves and host-galaxy masses for a sample of literature and Palomar Transient Factory SN Ia reveals a diverse distribution of SN Ia subtypes where carbon-rich material remains unburned in some instances. Such events include SN 1991T, 1997br, and 1999aa where trace signatures of C III at optical wavelengths are presumably detected.

We present optical observations of supernova SN 2014C, which underwent an unprecedented slow metamorphosis from H-poor type Ib to H-rich type IIn over the course of one year. The observed spectroscopic evolution is consistent with the supernova having exploded in a cavity before encountering a massive shell of the progenitor star's stripped hydrogen envelope. Possible origins for the circumstellar shell include a brief Wolf-Rayet fast wind phase that overtook a slower red supergiant wind, eruptive ejection, or confinement of circumstellar material by external influences of neighboring stars. An extended high velocity Halpha absorption feature seen in near-maximum light spectra implies that the progenitor star was not completely stripped of hydrogen at the time of core collapse. Archival pre-explosion Subaru Telescope Suprime-Cam and Hubble Space Telescope Wide Field Planetary Camera 2 images of the region obtained in 2009 show a coincident source that is most likely a compact massive star cluster in NGC 7331 that hosted the progenitor system. By comparing the emission properties of the source with stellar population models that incorporate interacting binary stars we estimate the age of the host cluster to be 30 - 300 Myr, and favor ages closer to 30 Myr in light of relatively strong Halpha emission. SN 2014C is the best-observed member of a class of core-collapse supernovae that fill the gap between events that interact strongly with dense, nearby environments immediately after explosion and those that never show signs of interaction. Better understanding of the frequency and nature of this intermediate population can contribute valuable information about the poorly understood final stages of stellar evolution.

We present Hubble Space Telescope (HST) rest-frame ultraviolet imaging of the host galaxies of 16 hydrogen-poor superluminous supernovae (SLSNe), including 11 events from the Pan-STARRS Medium Deep Survey. Taking advantage of the superb angular resolution of HST, we characterize the galaxies' morphological properties, sizes and star formation rate (SFR) densities. We determine the supernova (SN) locations within the host galaxies through precise astrometric matching, and measure physical and host-normalized offsets, as well as the SN positions within the cumulative distribution of UV light pixel brightness. We find that the host galaxies of H-poor SLSNe are irregular, compact dwarf galaxies, with a median half-light radius of just 0.9 kpc. The UV-derived SFR densities are high (<Sigma_SFR> ~ 0.1 M_sun/yr/kpc^2), suggesting that SLSNe form in overdense environments. Their locations trace the UV light of their host galaxies, with a distribution intermediate between that of long-duration gamma-ray bursts (LGRBs) (which are strongly clustered on the brightest regions of their hosts) and a uniform distribution (characteristic of normal core-collapse SNe), though cannot be statistically distinguished from either with the current sample size. Taken together, this strengthens the picture that SLSN progenitors require different conditions than those of ordinary core-collapse SNe to form, and that they explode in broadly similar galaxies as do LGRBs. If the tendency for SLSNe to be less clustered on the brightest regions than are LGRBs is confirmed by a larger sample, this would indicate a different, potentially lower-mass progenitor for SLSNe than LRGBs.

We present the GALEX detection of a UV burst at the time of explosion of an optically normal Type II-P supernova (PS1-13arp) from the Pan-STARRS1 survey at z=0.1665. The temperature and luminosity of the UV burst match the theoretical predictions for shock breakout in a red supergiant, but with a duration a factor of ~50 longer than expected. We compare the $NUV$ light curve of PS1-13arp to previous GALEX detections of Type IIP SNe, and find clear distinctions that indicate that the UV emission is powered by shock breakout, and not by the subsequent cooling envelope emission previously detected in these systems. We interpret the ~ 1 d duration of the UV signal with a shock breakout in the wind of a red supergiant with a pre-explosion mass-loss rate of ~ 10^-3 Msun yr^-1. This mass-loss rate is enough to prolong the duration of the shock breakout signal, but not enough to produce an excess in the optical plateau light curve or narrow emission lines powered by circumstellar interaction. This detection of non-standard, potentially episodic high mass-loss in a RSG SN progenitor has favorable consequences for the prospects of future wide-field UV surveys to detect shock breakout directly in these systems, and provide a sensitive probe of the pre-explosion conditions of SN progenitors.

In recent years, wide-field sky surveys providing deep multi-band imaging have presented a new path for indirectly characterizing the progenitor populations of core-collapse supernovae (SN): systematic light curve studies. We assemble a set of 76 grizy-band Type IIP SN light curves from Pan-STARRS1, obtained over a constant survey program of 4 years and classified using both spectroscopy and machine learning-based photometric techniques. We develop and apply a new Bayesian model for the full multi-band evolution of each light curve in the sample. We find no evidence of a sub-population of fast-declining explosions (historically referred to as "Type IIL" SNe). However, we identify a highly significant relation between the plateau phase decay rate and peak luminosity among our SNe IIP. These results argue in favor of a single parameter, likely determined by initial stellar mass, predominantly controlling the explosions of red supergiants. This relation could also be applied for supernova cosmology, offering a standardizable candle good to an intrinsic scatter of 0.2 mag. We compare each light curve to physical models from hydrodynamic simulations to estimate progenitor initial masses and other properties of the Pan-STARRS1 Type IIP SN sample. We show that correction of systematic discrepancies between modeled and observed SN IIP light curve properties and an expanded grid of progenitor properties, are needed to enable robust progenitor inferences from multi-band light curve samples of this kind. This work will serve as a pathfinder for photometric studies of core-collapse SNe to be conducted through future wide field transient searches.

Deep late-time X-ray observations of the relativistic, engine-driven, type Ic SN2012ap allow us to probe the nearby environment of the explosion and reveal the unique properties of relativistic SNe. We find that on a local scale of ~0.01 pc the environment was shaped directly by the evolution of the progenitor star with a pre-explosion mass-loss rate <5x10^-6 Msun yr-1 in line with GRBs and the other relativistic SN2009bb. Like sub-energetic GRBs, SN2012ap is characterized by a bright radio emission and evidence for mildly relativistic ejecta. However, its late time (t~20 days) X-ray emission is ~100 times fainter than the faintest sub-energetic GRB at the same epoch, with no evidence for late-time central engine activity. These results support theoretical proposals that link relativistic SNe like 2009bb and 2012ap with the weakest observed engine-driven explosions, where the jet barely fails to breakout. Furthermore, our observations demonstrate that the difference between relativistic SNe and sub-energetic GRBs is intrinsic and not due to line-of-sight effects. This phenomenology can either be due to an intrinsically shorter-lived engine or to a more extended progenitor in relativistic SNe.

We probe the systematic uncertainties from 113 Type Ia supernovae (SNIa) in the Pan-STARRS1 (PS1) sample along with 197 SN Ia from a combination of low-redshift surveys. The companion paper by Rest et al. (2013) describes the photometric measurements and cosmological inferences from the PS1 sample. The largest systematic uncertainty stems from the photometric calibration of the PS1 and low-z samples. We increase the sample of observed Calspec standards from 7 to 10 used to define the PS1 calibration system. The PS1 and SDSS-II calibration systems are compared and discrepancies up to ~0.02 mag are recovered. We find uncertainties in the proper way to treat intrinsic colors and reddening produce differences in the recovered value of w up to 3%. We estimate masses of host galaxies of PS1 supernovae and detect an insignificant difference in distance residuals of the full sample of 0.037\pm0.031 mag for host galaxies with high and low masses. Assuming flatness in our analysis of only SNe measurements, we find $w = {-1.120^{+0.360}_{-0.206}\textrm{(Stat)} ^{+0.269}_{-0.291}\textrm{(Sys)}}$. With additional constraints from BAO, CMB(Planck) and H0 measurements, we find $w = -1.166^{+0.072}_{-0.069}$ and $\Omega_M=0.280^{+0.013}_{-0.012}$ (statistical and systematic errors added in quadrature). Significance of the inconsistency with $w=-1$ depends on whether we use Planck or WMAP measurements of the CMB: $w_{\textrm{BAO+H0+SN+WMAP}}=-1.124^{+0.083}_{-0.065}$.

We present griz light curves of 146 spectroscopically confirmed Type Ia Supernovae ($0.03 < z <0.65$) discovered during the first 1.5 years of the Pan-STARRS1 Medium Deep Survey. The Pan-STARRS1 natural photometric system is determined by a combination of on-site measurements of the instrument response function and observations of spectrophotometric standard stars. We find that the systematic uncertainties in the photometric system are currently 1.2\% without accounting for the uncertainty in the HST Calspec definition of the AB system. A Hubble diagram is constructed with a subset of 113 out of 146 SNe Ia that pass our light curve quality cuts. The cosmological fit to 310 SNe Ia (113 PS1 SNe Ia + 222 light curves from 197 low-z SNe Ia), using only SNe and assuming a constant dark energy equation of state and flatness, yields $w=-1.120^{+0.360}_{-0.206}\textrm{(Stat)} ^{+0.269}_{-0.291}\textrm{(Sys)}$. When combined with BAO+CMB(Planck)+$H_0$, the analysis yields $\Omega_{\rm M}=0.280^{+0.013}_{-0.012}$ and $w=-1.166^{+0.072}_{-0.069}$ including all identified systematics (see also Scolnic et al. 2014). The value of $w$ is inconsistent with the cosmological constant value of $-1$ at the 2.3$\sigma$ level. Tension endures after removing either the BAO or the $H_0$ constraint, though it is strongest when including the $H_0$ constraint. If we include WMAP9 CMB constraints instead of those from Planck, we find $w=-1.124^{+0.083}_{-0.065}$, which diminishes the discord to $<2\sigma$. We cannot conclude whether the tension with flat $\Lambda$CDM is a feature of dark energy, new physics, or a combination of chance and systematic errors. The full Pan-STARRS1 supernova sample with $\sim\!\!$3 times as many SNe should provide more conclusive results.

We present ultraviolet, optical, and near-infrared observations of SN 2012ap, a broad-lined Type Ic supernova in the galaxy NGC 1729 that produced a relativistic and rapidly decelerating outflow without a gamma-ray burst signature. Photometry and spectroscopy follow the flux evolution from -13 to +272 days past the B-band maximum of -17.4 +/- 0.5 mag. The spectra are dominated by Fe II, O I, and Ca II absorption lines at ejecta velocities of 20,000 km/s that change slowly over time. Other spectral absorption lines are consistent with contributions from photospheric He I, and hydrogen may also be present at higher velocities (> 27,000 km/s). We use these observations to estimate explosion properties and derive a total ejecta mass of 2.7 Msolar, a kinetic energy of 1.0x10^{52} erg, and a 56Ni mass of 0.1-0.2 Msolar. Nebular spectra (t > 200d) exhibit an asymmetric double-peaked [OI] 6300,6364 emission profile that we associate with absorption in the supernova interior, although toroidal ejecta geometry is an alternative explanation. SN 2012ap joins SN 2009bb as another exceptional supernova that shows evidence for a central engine (e.g., black-hole accretion or magnetar) capable of launching a non-negligible portion of ejecta to relativistic velocities without a coincident gamma-ray burst detection. Defining attributes of their progenitor systems may be related to notable properties including above-average environmental metallicities of Z > Zsolar, moderate to high levels of host-galaxy extinction (E(B-V) > 0.4 mag), detection of high-velocity helium at early epochs, and a high relative flux ratio of [Ca II]/[O I] > 1 at nebular epochs. These events support the notion that jet activity at various energy scales may be present in a wide range of supernovae.

In the past decade, several rapidly-evolving transients have been discovered whose timescales and luminosities are not easily explained by traditional supernovae (SN) models. The sample size of these objects has remained small due, at least in part, to the challenge of detecting short timescale transients with traditional survey cadences. Here we present the results from a search within the Pan-STARRS1 Medium Deep Survey (PS1-MDS) for rapidly-evolving and luminous transients. We identify 10 new transients with a time above half-maximum of less than 12 days and -16.5 > M > -20 mag. This increases the number of known events in this region of SN phase space by roughly a factor of three. The median redshift of the PS1-MDS sample is z=0.275 and they all exploded in star forming galaxies. In general, the transients possess faster rise than decline timescale and blue colors at maximum light (g - r < -0.2). Best fit blackbodies reveal photospheric temperatures/radii that expand/cool with time and explosion spectra taken near maximum light are dominated by a blue continuum, consistent with a hot, optically thick, ejecta. We find it difficult to reconcile the short timescale, high peak luminosity (L > 10^43 erg/s), and lack of UV line blanketing observed in many of these transients with an explosion powered mainly by the radioactive decay of Ni-56. Rather, we find that many are consistent with either (1) cooling envelope emission from the explosion of a star with a low-mass extended envelope which ejected very little (<0.03 M_sun) radioactive material, or (2) a shock breakout within a dense, optically thick, wind surrounding the progenitor star. After calculating the detection efficiency for objects with rapid timescales in the PS1-MDS we find a volumetric rate of 4800 - 8000 events/yr/Gpc^3 (4-7% of the core-collapse SN rate at z=0.2).

We present optical spectroscopy and optical/near-IR photometry of 31 host galaxies of hydrogen-poor superluminous supernovae (SLSNe), including 15 events from the Pan-STARRS1 Medium Deep Survey. Our sample spans the redshift range 0.1 < z < 1.6 and is the first comprehensive host galaxy study of this specific subclass of cosmic explosions. Combining the multi-band photometry and emission-line measurements, we determine the luminosities, stellar masses, star formation rates and metallicities. We find that as a whole, the hosts of SLSNe are a low-luminosity (<M_B> ~ -17.3 mag), low stellar mass (<M_*> ~ 2 x 10^8 M_sun) population, with a high median specific star formation rate (<sSFR> ~ 2 Gyr^-1). The median metallicity of our spectroscopic sample is low, 12 + log(O/H}) ~ 8.35 ~ 0.45 Z_sun, although at least one host galaxy has solar metallicity. The host galaxies of H-poor SLSNe are statistically distinct from the hosts of GOODS core-collapse SNe (which cover a similar redshift range), but resemble the host galaxies of long-duration gamma-ray bursts (LGRBs) in terms of stellar mass, SFR, sSFR and metallicity. This result indicates that the environmental causes leading to massive stars forming either SLSNe or LGRBs are similar, and in particular that SLSNe are more effectively formed in low metallicity environments. We speculate that the key ingredient is large core angular momentum, leading to a rapidly-spinning magnetar in SLSNe and an accreting black hole in LGRBs.

Historically, light curve studies of supernovae (SNe) and other transient classes have focused on individual objects with copious and high signal-to-noise observations. In the nascent era of wide field transient searches, objects with detailed observations are decreasing as a fraction of the overall known SN population, and this strategy sacrifices the majority of the information contained in the data about the underlying population of transients. A population level modeling approach, simultaneously fitting all available observations of objects in a transient sub-class of interest, fully mines the data to infer the properties of the population and avoids certain systematic biases. We present a novel hierarchical Bayesian statistical model for population level modeling of transient light curves, and discuss its implementation using an efficient Hamiltonian Monte Carlo technique. As a test case, we apply this model to the Type IIP SN sample from the Pan-STARRS1 Medium Deep Survey, consisting of 18,837 photometric observations of 76 SNe, corresponding to a joint posterior distribution with 9,176 parameters under our model. Our hierarchical model fits provide improved constraints on light curve parameters relevant to the physical properties of their progenitor stars relative to modeling individual light curves alone. Moreover, we directly evaluate the probability for occurrence rates of unseen light curve characteristics from the model hyperparameters, addressing observational biases in survey methodology. We view this modeling framework as an unsupervised machine learning technique with the ability to maximize scientific returns from data to be collected by future wide field transient searches like LSST. \smallskip

The diffuse interstellar bands (DIBs) are absorption features observed in optical and near-infrared spectra that are thought to be associated with carbon-rich polyatomic molecules in interstellar gas. However, because the central wavelengths of these bands do not correspond with electronic transitions of any known atomic or molecular species, their nature has remained uncertain since their discovery almost a century ago. Here we report on unusually strong DIBs in optical spectra of the broad-lined Type Ic supernova SN 2012ap that exhibit changes in equivalent width over short (~30 days) timescales. The 4428 and 6283 Angstrom DIB features get weaker with time, whereas the 5780 Angstrom feature shows a marginal increase. These nonuniform changes suggest that the supernova is interacting with a nearby source of the DIBs and that the DIB carriers possess high ionization potentials, such as small cations or charged fullerenes. We conclude that moderate-resolution spectra of supernovae with DIB absorptions obtained within weeks of outburst could reveal unique information about the mass-loss environment of their progenitor systems and provide new constraints on the properties of DIB carriers.

Observations spanning a large wavelength range, from X-ray to radio, of the Type IIb supernova 2011hs are presented, covering its evolution during the first year after explosion. The optical light curve presents a narrower shape and a fainter luminosity at peak than previously observed for Type IIb SNe. High expansion velocities are measured from the broad absorption H I and He I lines. From the comparison of the bolometric light curve and the time evolution of the photospheric velocities with hydrodynamical models, we found that SN 2011hs is consistent with the explosion of a 3-4 Msun He-core progenitor star, corresponding to a main sequence mass of 12-15 Msun, that ejected a mass of 56Ni of about 0.04 Msun, with an energy of E= 8.5 x 10^50 erg. Such a low-mass progenitor scenario is in full agreement with the modelling of the nebular spectrum taken at $\sim$215 days from maximum. From the modelling of the adiabatic cooling phase, we infer a progenitor radius of $\approx$500-600 Rsun, clearly pointing to an extended progenitor star. The radio light curve of SN 2011hs yields a peak luminosity similar to that of SN 1993J, but with a higher mass loss rate and a wind density possibly more similar to that of SN 2001ig. Although no significant deviations from a smooth decline have been found in the radio light curves, we cannot rule out the presence of a binary companion star.

Educational technology has attained significant importance as a mechanism for supporting experiential learning of science concepts. However, the growth of this mechanism is limited by the significant time and technical expertise needed to develop such products, particularly in specialized fields of science. We sought to test whether interactive, educational, online software modules can be developed effectively by students as a curriculum component of an advanced science course. We discuss a set of fifteen such modules developed by Harvard University graduate students to demonstrate various concepts related to astronomy and physics. Their successful development of these modules demonstrates that online software tools for education and outreach on specialized topics can be produced while simultaneously fulfilling project-based learning objectives. We describe a set of technologies suitable for module development and present in detail four examples of modules developed by the students. We offer recommendations for incorporating educational software development within a graduate curriculum and conclude by discussing the relevance of this novel approach to new online learning environments like edX

We present extensive radio and X-ray observations of SN\, 2012au, the energetic radio luminous supernova of type Ib that may be a link between subsets of hydrogen-poor superluminous and normal core-collapse supernovae. The observations closely follow models of synchrotron emission from shock heated circum-burst medium that has a wind density profile ($\rho \propto r^{-2}$). We infer a sub-relativistic velocity for the shock wave $v \approx 0.2\,c$ and a radius of $r \approx 1.4 \times 10^{16} \rm cm$ at 25 days after the estimated date of explosion. For a constant wind velocity of 1000 km/s we determine the constant mass loss rate of the progenitor to be $\dot{M} = 3.6 \times 10^{-6} \rm M_{\odot} yr^{-1}$, consistent with the estimates from X-ray observations. We estimate the total internal energy of the radio emitting material to be $E \approx 10^{47} \rm erg$, which is intermediate to SN\,1998bw and SN\,2002ap. Evolution of the radio light curves of SN\,2012au is consistent with interaction with a smoothly distributed circum-burst medium and absence of stellar shells ejected from previous outbursts out to $r \approx 10^{17} \rm cm$ from the supernova site. Based on this we conclude that the evolution of the SN\,2012au progenitor star was relatively quiet during the final years preceding explosion. We find that the bright radio emission from SN2012au was not dissimilar from other core collapse supernovae despite it's extraordinary optical properties. We speculate that it was the nature of the explosion that led to the unusual demise of the SN2012au progenitor star.

We present the Pan-STARRS1 discovery of the long-lived and blue transient PS1-11af, which was also detected by GALEX with coordinated observations in the near-ultraviolet (NUV) band. PS1-11af is associated with the nucleus of an early-type galaxy at redshift z=0.4046 that exhibits no evidence for star formation or AGN activity. Four epochs of spectroscopy reveal a pair of transient broad absorption features in the UV on otherwise featureless spectra. Despite the superficial similarity of these features to P-Cygni absorptions of supernovae (SNe), we conclude that PS1-11af is not consistent with the properties of known types of SNe. Blackbody fits to the spectral energy distribution are inconsistent with the cooling, expanding ejecta of a SN, and the velocities of the absorption features are too high to represent material in homologous expansion near a SN photosphere. However, the constant blue colors and slow evolution of the luminosity are similar to previous optically-selected tidal disruption events (TDEs). The shape of the optical light curve is consistent with models for TDEs, but the minimum accreted mass necessary to power the observed luminosity is only ~0.002M_sun, which points to a partial disruption model. A full disruption model predicts higher bolometric luminosities, which would require most of the radiation to be emitted in a separate component at high energies where we lack observations. In addition, the observed temperature is lower than that predicted by pure accretion disk models for TDEs and requires reprocessing to a constant, lower temperature. Three deep non-detections in the radio with the VLA over the first two years after the event set strict limits on the production of any relativistic outflow comparable to Swift J1644+57, even if off-axis.

There is a well known correlation between the mass and metallicity of star-forming galaxies. Because mass is correlated with luminosity, this relation is often exploited, when spectroscopy is not available, to estimate galaxy metallicities based on single band photometry. However, we show that galaxy color is typically more effective than luminosity as a predictor of metallicity. This is a consequence of the correlation between color and the galaxy mass-to-light ratio and the recently discovered correlation between star formation rate (SFR) and residuals from the mass-metallicity relation. Using Sloan Digital Sky Survey spectroscopy of 180,000 nearby galaxies, we derive "LZC relations," empirical relations between metallicity (in seven common strong line diagnostics), luminosity, and color (in ten filter pairs and four methods of photometry). We show that these relations allow photometric metallicity estimates, based on luminosity and a single optical color, that are 50% more precise than those made based on luminosity alone; galaxy metallicity can be estimated to within 0.05 - 0.1 dex of the spectroscopically-derived value depending on the metallicity diagnostic used. Including color information in metallicity estimates also reduces systematic biases for populations skewed toward high or low SFR environments, as we illustrate using the host galaxy of the supernova SN 2010ay. This new tool will lend more statistical power to studies of galaxy populations, such as supernova and gamma-ray burst (GRB) host environments, in ongoing and future wide field imaging surveys.

We present extensive multi-wavelength observations of the extremely rapidly declining Type Ic supernova, SN 2005ek. Reaching a peak magnitude of M_R = -17.3 and decaying by ~3 mag in the first 15 days post-maximum, SN 2005ek is among the fastest Type I supernovae observed to date. The spectra of SN 2005ek closely resemble those of normal SN Ic, but with an accelerated evolution. There is evidence for the onset of nebular features at only nine days post-maximum. Spectroscopic modeling reveals an ejecta mass of ~0.3 Msun that is dominated by oxygen (~80%), while the pseudo-bolometric light curve is consistent with an explosion powered by ~0.03 Msun of radioactive Ni-56. Although previous rapidly evolving events (e.g., SN 1885A, SN 1939B, SN 2002bj, SN 2010X) were hypothesized to be produced by the detonation of a helium shell on a white dwarf, oxygen-dominated ejecta are difficult to reconcile with this proposed mechanism. We find that the properties of SN 2005ek are consistent with either the edge-lit double detonation of a low-mass white dwarf or the iron-core collapse of a massive star, stripped by binary interaction. However, if we assume that the strong spectroscopic similarity of SN 2005ek to other SN Ic is an indication of a similar progenitor channel, then a white-dwarf progenitor becomes very improbable. SN 2005ek may be one of the lowest mass stripped-envelope core-collapse explosions ever observed. We find that the rate of such rapidly declining Type I events is at least 1-3% of the normal SN Ia rate.

We present observations and analysis of PS1-10bzj, a superluminous supernova (SLSN) discovered in the Pan-STARRS Medium Deep Survey at a redshift z = 0.650. Spectroscopically, PS1-10bzj is similar to the hydrogen-poor SLSNe 2005ap and SCP 06F6, though with a steeper rise and lower peak luminosity (M_bol = -21.4 mag) than previous events. We construct a bolometric light curve, and show that while PS1-10bzj's energetics were less extreme than previous events, its luminosity still cannot be explained by radioactive nickel decay alone. We explore both a magnetar spin-down and circumstellar interaction scenario and find that either can fit the data. PS1-10bzj is located in the Extended Chandra Deep Field South and the host galaxy is imaged in a number of surveys, including with the Hubble Space Telescope. The host is a compact dwarf galaxy (M_B ~ -18 mag, diameter < 800 pc), with a low stellar mass (M_* ~ 2.4 * 10^7 M_sun), young stellar population (\tau_* ~ 5 Myr), and a star formation rate of ~ 2-3 M_sun/yr. The specific star formation rate is the highest seen in a SLSN host so far (~ 100 Gyr^{-1}). We detect the [O III]\lambda 4363 line, and find a low metallicity: 12+(O/H) = 7.8 +/- 0.2 (~ 0.1 Z_sun). Together, this indicates that at least some of the progenitors of SLSNe come from young, low-metallicity populations.

We present optical and near-infrared observations of SN 2012au, a slow-evolving supernova (SN) with properties that suggest a link between subsets of energetic and H-poor SNe and superluminous SNe. SN 2012au exhibited conspicuous SN Ib-like He I lines and other absorption features at velocities reaching 2 x 10^4 km/s in its early spectra, and a broad light curve that peaked at M_B = -18.1 mag. Models of these data indicate a large explosion kinetic energy of 10^{52} erg and 56Ni mass ejection of 0.3 Msolar on par with SN 1998bw. SN 2012au's spectra almost one year after explosion show a blend of persistent Fe II P-Cyg absorptions and nebular emissions originating from two distinct velocity regions. These late-time emissions include strong [Fe II], [Ca II], [O I], Mg I], and Na I lines at velocities > 4500 km/s, as well as O I and Mg I lines at noticeably smaller velocities of 2000 km/s. Many of the late-time properties of SN 2012au are similar to the slow-evolving hypernovae SN 1997dq and SN 1997ef, and the superluminous SN 2007bi. Our observations suggest that a single explosion mechanism may unify all of these events that span -21 < M_B < -17 mag. The aspherical and possibly jetted explosion was most likely initiated by the core collapse of a massive progenitor star and created substantial high-density, low-velocity Ni-rich material.

The 2012 explosion of SN2009ip raises questions about our understanding of the late stages of massive star evolution. Here we present a comprehensive study of SN2009ip during its remarkable re-brightening(s). High-cadence photometric and spectroscopic observations from the GeV to the radio band obtained from a variety of ground-based and space facilities (including the VLA, Swift, Fermi, HST and XMM) constrain SN2009ip to be a low energy (E~ 10^50 erg for an ejecta mass ~ 0.5 Msun) and likely asymmetric explosion in a complex medium shaped by multiple eruptions of the restless progenitor star. Most of the energy is radiated as a result of the shock breaking out through a dense shell of material located at 5x10^14 cm with M~0.1 Msun, ejected by the precursor outburst ~40 days before the major explosion. We interpret the NIR excess of emission as signature of dust vaporization of material located further out (R>4x 10^15 cm), the origin of which has to be connected with documented mass loss episodes in the previous years. Our modeling predicts bright neutrino emission associated with the shock break-out if the cosmic ray energy is comparable to the radiated energy. We connect this phenomenology with the explosive ejection of the outer layers of the massive progenitor star, that later interacted with material deposited in the surroundings by previous eruptions. Future observations will reveal if the luminous blue variable (LBV) progenitor star survived. Irrespective of whether the explosion was terminal, SN2009ip brought to light the existence of new channels for sustained episodic mass-loss, the physical origin of which has yet to be identified.

We report on our discovery and observations of the Pan-STARRS1 supernova (SN) PS1-12sk, a transient with properties that indicate atypical star formation in its host galaxy cluster or pose a challenge to popular progenitor system models for this class of explosion. The optical spectra of PS1-12sk classify it as a Type Ibn SN (c.f. SN 2006jc), dominated by intermediate-width (3x10^3 km/s) and time variable He I emission. Our multi-wavelength monitoring establishes the rise time dt = 9-23 days and shows an NUV-NIR SED with temperature > 17x10^3 K and a peak rise magnitude of Mz = -18.9 mag. SN Ibn spectroscopic properties are commonly interpreted as the signature of a massive star (17 - 100 M_sun) explosion within a He-enriched circumstellar medium. However, unlike previous Type Ibn supernovae, PS1-12sk is associated with an elliptical brightest cluster galaxy, CGCG 208-042 (z = 0.054) in cluster RXC J0844.9+4258. The expected probability of an event like PS1-12sk in such environments is low given the measured infrequency of core-collapse SNe in red sequence galaxies compounded by the low volumetric rate of SN Ibn. Furthermore, we find no evidence of star formation at the explosion site to sensitive limits (Sigma Halpha < 2x10^-3 M_sun/yr/kpc^2). We therefore discuss white dwarf binary systems as a possible progenitor channel for SNe Ibn. We conclude that PS1-12sk represents either a fortuitous and statistically unlikely discovery, evidence for a top-heavy IMF in galaxy cluster cooling flow filaments, or the first clue suggesting an alternate progenitor channel for Type Ibn SNe.

We present the Pan-STARRS1 discovery of PS1-10afx, a unique hydrogen-deficient superluminous supernova (SLSN) at z=1.388. The light curve peaked at z_P1=21.7 mag, making PS1-10afx comparable to the most luminous known SNe, with M_u = -22.3 mag. Our extensive optical and NIR observations indicate that the bolometric light curve of PS1-10afx rose on the unusually fast timescale of ~12 d to the extraordinary peak luminosity of 4.1e44 erg/s (M_bol = -22.8 mag) and subsequently faded rapidly. Equally important, the SED is unusually red for a SLSN, with a color temperature of 6800 K near maximum light, in contrast to previous H-poor SLSNe, which are bright in the UV. The spectra more closely resemble those of a normal SN Ic than any known SLSN, with a photospheric velocity of 11,000 km/s and evidence for line blanketing in the rest-frame UV. Despite the fast rise, these parameters imply a very large emitting radius (>5e15 cm). We demonstrate that no existing theoretical model can satisfactorily explain this combination of properties: (i) a nickel-powered light curve cannot match the combination of high peak luminosity with the fast timescale; (ii) models powered by the spindown energy of a rapidly-rotating magnetar predict significantly hotter and faster ejecta; and (iii) models invoking shock breakout through a dense circumstellar medium cannot explain the observed spectra or color evolution. The host galaxy is well detected in pre-explosion imaging with a luminosity near L*, a star formation rate of 15 M_sun/yr, and is fairly massive (2e10 M_sun), with a stellar population age of 1e8 yr, also in contrast to the dwarf hosts of known H-poor SLSNe. PS1-10afx is distinct from known examples of SLSNe in its spectra, colors, light-curve shape, and host galaxy properties, suggesting that it resulted from a different channel than other hydrogen-poor SLSNe.

We present X-ray, UV/optical, and radio observations of the stripped-envelope, core-collapse supernova (SN) 2011ei, one of the least luminous SNe IIb or Ib observed to date. Our observations begin with a discovery within 1 day of explosion and span several months afterward. Early optical spectra exhibit broad, Type II-like hydrogen Balmer profiles that subside rapidly and are replaced by Type Ib-like He-rich features on the timescale of one week. High-cadence monitoring of this transition suggests that absorption attributable to a high velocity (> 12,000 km/s) H-rich shell is not rare in Type Ib events. Radio observations imply a shock velocity of v = 0.13c and a progenitor star mass-loss rate of 1.4 x 10^{-5} Msun yr^{-1} (assuming wind velocity v_w=10^3 km/s). This is consistent with independent constraints from deep X-ray observations with Swift-XRT and Chandra. Overall, the multi-wavelength properties of SN 2011ei are consistent with the explosion of a lower-mass (3-4 Msun), compact (R* <= 1x10^{11} cm), He core star. The star retained a thin hydrogen envelope at the time of explosion, and was embedded in an inhomogeneous circumstellar wind suggestive of modest episodic mass-loss. We conclude that SN 2011ei's rapid spectral metamorphosis is indicative of time-dependent classifications that bias estimates of explosion rates for Type IIb and Ib objects, and that important information about a progenitor star's evolutionary state and mass-loss immediately prior to SN explosion can be inferred from timely multi-wavelength observations.

We present the largest spectroscopic study of the host environments of Type Ibc supernovae (SN Ibc) discovered exclusively by untargeted SN searches. Past studies of SN Ibc host environments have been biased towards high-mass, high-metallicity galaxies by focusing on SNe discovered in galaxy-targeted SN searches. Our new observations more than double the total number of spectroscopic stellar population age and metallicity measurements published for untargeted SN Ibc host environments, and extend to a median redshift about twice as large as previous statistical studies (z = 0.04). For the 12 SNe Ib and 21 SNe Ic in our metallicity sample, we find median metallicities of log(O/H)+12 = 8.48 and 8.61, respectively, but determine that the discrepancy in the full distribution of metallicities is not statistically significant. This median difference would correspond to only a small difference in the mass loss via metal-line driven winds (<30%), suggesting this does not play the dominant role in distinguishing SN Ib and Ic progenitors. However, the median metallicity of the 7 broad-lined SN Ic (SN Ic-BL) in our sample is significantly lower, log(O/H)+12 = 8.34. The age of the young stellar population of SN Ic-BL host environments also seems to be lower than for SN Ib and Ic, but our age sample is small. A synthesis of SN Ibc host environment spectroscopy to date does not reveal a significant difference in SN Ib and Ic metallicities, but reinforces the significance of the lower metallicities for SN Ic-BL. This combined sample demonstrates that galaxy-targeted SN searches introduce a significant bias for studies seeking to infer the metallicity distribution of SN progenitors, and we identify and discuss other systematic effects that play smaller roles. We discuss the path forward for making progress on SN Ibc progenitor studies in the LSST era.

The oxygen abundance gradients among nebular emission line regions in spiral galaxies have been used as important constraints for models of chemical evolution. We present the largest ever full-wavelength optical spectroscopic sample of emission line nebulae in a spiral galaxy (M31). We have collected spectra of 253 HII regions and 407 planetary nebulae with the Hectospec multi-fiber spectrograph of the MMT. We measure the line-of-sight extinction for 199 HII regions and 333 PNe; we derive oxygen abundance directly, based on the electron temperature, for 51 PNe; and we use strong line methods to estimate oxygen abundance for 192 HII regions and nitrogen abundance for 52 HII regions. The relatively shallow oxygen abundance gradient of the more extended HII regions in our sample is generally in agreement with the result of Zaritsky et al. (1994), based on only 19 M31 HII regions, but varies with the strong-line diagnostic employed. Our large sample size demonstrates that there is significant intrinsic scatter around this abundance gradient, as much as 3 times the systematic uncertainty in the strong line diagnostics. The intrinsic scatter is similar in the nitrogen abundances, although the gradient is significantly steeper. On small scales (deprojected distance < 0.5 kpc), HII regions exhibit local variations in oxygen abundance that are larger than 0.3 dex in 33% of neighboring pairs. We do not identify a significant oxygen abundance gradient among PNe, but we do find a significant gradient in the [N II] ratio that varies systematically with surface brightness. Our results underscore the complex and inhomogeneous nature of the ISM of M31, and our dataset illustrates systematic effects relevant to future studies of the metallicity gradients in nearby spiral galaxies.

Because undergraduate participation in research is a longstanding and increasingly important aspect of the career path for future scientists, students can benefit from additional resources to introduce them to the culture and process of research. We suggest the adoption of the web resource Astrobites as a classroom tool to increase the preparation of undergraduate physics and astronomy students for careers in research. We describe the content and development of the website, discuss previous university courses that have made use of Astrobites, and suggest additional strategies for using Astrobites in the classroom.

We present multi-wavelength observations of the afterglow of the short GRB111117A, and follow-up observations of its host galaxy. From rapid optical and radio observations we place limits of r \gtrsim 25.5 mag at \deltat \approx 0.55 d and F_nu(5.8 GHz) < 18 \muJy at \deltat \approx 0.50 d, respectively. However, using a Chandra observation at t~3.0 d we locate the absolute position of the X-ray afterglow to an accuracy of 0.22" (1 sigma), a factor of about 6 times better than the Swift-XRT position. This allows us to robustly identify the host galaxy and to locate the burst at a projected offset of 1.25 +/- 0.20" from the host centroid. Using optical and near-IR observations of the host galaxy we determine a photometric redshift of z=1.3 (+0.3,-0.2), one of the highest for any short GRB, and leading to a projected physical offset for the burst of 10.5 +/- 1.7 kpc, typical of previous short GRBs. At this redshift, the isotropic gamma-ray energy is E_{gamma,iso} \approx 3\times10^51 erg (rest-frame 23-2300 keV) with a peak energy of E_{pk} \approx 850-2300 keV (rest-frame). In conjunction with the isotropic X-ray energy, GRB111117A appears to follow our recently-reported E_x,iso-E_gamma,iso-E_pk universal scaling. Using the X-ray data along with the optical and radio non-detections we find that for a blastwave kinetic energy of E_{K,iso} \approx E_{gamma,iso}, the circumburst density is n_0 \sim 3x10^(-4)-1 cm^-3 (for a range of epsilon_B=0.001-0.1). Similarly, from the non-detection of a break in the X-ray light curve at t<3 d, we infer a minimum opening angle for the outflow of theta_j> 3-10 degrees (depending on the circumburst density). We conclude that Chandra observations of short GRBs are effective at determining precise positions and robust host galaxy associations in the absence of optical and radio detections.

We present the Pan-STARRS1 discovery and light curves, and follow-up MMT and Gemini spectroscopy of an ultra-luminous supernova (ULSN; dubbed PS1-11bam) at a redshift of z=1.566 with a peak brightness of M_UV=-22.3 mag. PS1-11bam is one of the highest redshift spectroscopically-confirmed SNe known to date. The spectrum is characterized by broad absorption features typical of previous ULSNe (e.g., CII, SiIII), and by strong and narrow MgII and FeII absorption lines from the interstellar medium (ISM) of the host galaxy, confirmed by an [OII]3727 emission line at the same redshift. The equivalent widths of the FeII2600 and MgII2803 lines are in the top quartile of the quasar intervening absorption system distribution, but are weaker than those of gamma-ray burst intrinsic absorbers (i.e., GRB host galaxies). We also detect the host galaxy in pre-explosion Pan-STARRS1 data and find that its UV spectral energy distribution is best fit with a young stellar population age of tau~15-45 Myr and a stellar mass of M \sim (1.1-2.6)x10^9 M_sun (for Z=0.05-1 Z_sun). The star formation rate inferred from the UV continuum and [OII]3727 emission line is ~10 M_sun/yr, higher than in any previous ULSN host. PS1-11bam provides the first direct demonstration that ULSNe can serve as probes of the interstellar medium in distant galaxies. At the present, the depth and red sensitivity of PS1 are uniquely suited to finding such events at cosmologically interesting redshifts (z~1-2); the future combination of LSST and 30-m class telescopes promises to extend this technique to z~4.

We report on our serendipitous pre-discovery detection and detailed follow-up of the broad-lined Type Ic supernova (SN) 2010ay at z = 0.067 imaged by the Pan-STARRS1 3pi survey just ~4 days after explosion. The SN had a peak luminosity, M_R ~ -20.2 mag, significantly more luminous than known GRB-SNe and one of the most luminous SNe Ib/c ever discovered. The absorption velocity of SN 2010ay is v_Si ~ 19,000 km/s at ~40 days after explosion, 2-5 times higher than other broad-lined SNe and similar to the GRB-SN 2010bh at comparable epochs. Moreover, the velocity declines ~2 times slower than other SNe Ic-BL and GRB-SNe. Assuming that the optical emission is powered by radioactive decay, the peak magnitude implies the synthesis of an unusually large mass of 56 Ni, M_Ni = 0.9 M_solar. Modeling of the light-curve points to a total ejecta mass, M_ej ~ 4.7 M_sol, and total kinetic energy, E_K ~ 11x10^51 ergs. The ratio of M_Ni to M_ej is ~2 times as large for SN 2010ay as typical GRB-SNe and may suggest an additional energy reservoir. The metallicity (log(O/H)_PP04 + 12 = 8.19) of the explosion site within the host galaxy places SN 2010ay in the low-metallicity regime populated by GRB-SNe, and ~0.5(0.2) dex lower than that typically measured for the host environments of normal (broad-lined) Ic supernovae. We constrain any gamma-ray emission with E_gamma < 6x10^{48} erg (25-150 keV) and our deep radio follow-up observations with the Expanded Very Large Array rule out relativistic ejecta with energy, E > 10^48 erg. We therefore rule out the association of a relativistic outflow like those which accompanied SN 1998bw and traditional long-duration GRBs, but place less-stringent constraints on a weak afterglow like that seen from XRF 060218. These observations challenge the importance of progenitor metallicity for the production of a GRB, and suggest that other parameters also play a key role.

To understand how best to use observations of Type Ia supernovae (SNe Ia) to obtain precise and accurate distances, we investigate the relations between spectra of SNe Ia and their intrinsic colors. Using a sample of 1630 optical spectra of 255 SNe, based primarily on data from the CfA Supernova Program, we examine how the velocity evolution and line strengths of Si II 6355 and Ca II H&K are related to the B-V color at peak brightness. We find that the maximum-light velocity of Si II 6355 and Ca II H&K and the maximum-light pseudo-equivalent width of Si II 6355 are correlated with intrinsic color, with intrinsic color having a linear relation with the Si II 6355 measurements. Ca II H&K does not have a linear relation with intrinsic color, but lower-velocity SNe tend to be intrinsically bluer. Combining the spectroscopic measurements does not improve intrinsic color inference. The intrinsic color scatter is larger for higher-velocity SNe Ia --- even after removing a linear trend with velocity --- indicating that lower-velocity SNe Ia are more "standard crayons." Employing information derived from SN Ia spectra has the potential to improve the measurements of extragalactic distances and the cosmological properties inferred from them.

We present the discovery of two ultra-luminous supernovae (SNe) at z ~ 0.9 with the Pan-STARRS1 Medium-Deep Survey. These SNe, PS1-10ky and PS1-10awh, are amongst the most luminous SNe ever discovered, comparable to the unusual transients SN 2005ap and SCP 06F6. Like SN 2005ap and SCP 06F6, they show characteristic high luminosities (M_bol ~ -22.5 mag), blue spectra with a few broad absorption lines, and no evidence for H or He. We have constructed a full multi-color light curve sensitive to the peak of the spectral energy distribution in the rest-frame ultraviolet, and we have obtained time-series spectroscopy for these SNe. Given the similarities between the SNe, we combine their light curves to estimate a total radiated energy over the course of explosion of (0.9-1.4) x 10^51 erg. We find photospheric velocities of 12,000-19,000 km/s with no evidence for deceleration measured across ~3 rest-frame weeks around light-curve peak, consistent with the expansion of an optically-thick massive shell of material. We show that, consistent with findings for other ultra-luminous SNe in this class, radioactive decay is not sufficient to power PS1-10ky, and we discuss two plausible origins for these events: the initial spin-down of a newborn magnetar in a core-collapse SN, or SN shock breakout from the dense circumstellar wind surrounding a Wolf-Rayet star.