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 use 1169 Pan-STARRS supernovae (SNe) and 195 low-$z$ ($z < 0.1$) SNe Ia to measure cosmological parameters. Though most Pan-STARRS SNe lack spectroscopic classifications, in a previous paper (I) we demonstrated that photometrically classified SNe can be used to infer unbiased cosmological parameters by using a Bayesian methodology that marginalizes over core-collapse (CC) SN contamination. Our sample contains nearly twice as many SNe as the largest previous SN Ia compilation. Combining SNe with Cosmic Microwave Background (CMB) constraints from Planck, we measure the dark energy equation of state parameter $w$ to be -0.989$\pm$0.057 (stat$+$sys). If $w$ evolves with redshift as $w(a) = w_0 + w_a(1-a)$, we find $w_0 = -0.912 \pm 0.149$ and $w_a =$ -0.513$\pm$0.826. These results are consistent with cosmological parameters from the Joint Lightcurve Analysis and the Pantheon sample. We try four different photometric classification priors for Pan-STARRS SNe and two alternate ways of modeling CC SN contamination, finding that no variant gives a $w$ differing by more than 2% from the baseline measurement. The systematic uncertainty on $w$ due to marginalizing over CC SN contamination, $\sigma_w^{\textrm{CC}} = 0.012$, is the third-smallest source of systematic uncertainty in this work. We find limited (1.6$\sigma$) evidence for evolution of the SN color-luminosity relation with redshift, a possible systematic that could constitute a significant uncertainty in future high-$z$ analyses. Our data provide one of the best current constraints on $w$, demonstrating that samples with $\sim$5% CC SN contamination can give competitive cosmological constraints when the contaminating distribution is marginalized over in a Bayesian framework.

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.

We present the complete set of ultra-violet, optical and near-infrared photometry and spectroscopy for SN 2012ca, covering the period from 6~days prior to maximum light, until 531 days after maximum. The spectroscopic time series for SN 2012ca is essentially unchanged over 1.5 years, and appear to be dominated at all epochs by signatures of interaction with a dense circumstellar medium rather than the underlying supernova (SN). At late phases, we see a near infrared excess in flux which is possibly associated with dust formation, although without any signs of accompanying line shifts. SN 2012ca is a member of the set of type of the ambiguous IIn/Ia-CSM SNe, the nature of which have been debated extensively in the literature. The two leading scenarios are either a type Ia SN exploding within a dense CSM from a non-degenerate, evolved companion, or a core-collapse SN from a massive star. While some members of the population have been unequivocally associated with type Ia SNe, in other cases the association is less certain. While it is possible that SN 2012ca does arise from a thermonuclear SN, this would require a relatively high (between 20 and 70 per cent) efficiency in converting kinetic energy to optical luminosity, and a massive (~2.3-2.6 Msun) circumstellar medium. On the basis of energetics, and the results of simple modelling, we suggest that SN 2012ca is more likely associated with a core-collapse SN. This would imply that the observed set of similar SNe to SN 2012ca is in fact originated by two populations, and while these are drawn from physically distinct channels, they can have observationally similar properties.

The Public European Southern Observatory Spectroscopic Survey of Transient Objects (PESSTO) began as a public spectroscopic survey in April 2012. We describe the data reduction strategy and data products which are publicly available through the ESO archive as the Spectroscopic Survey Data Release 1 (SSDR1). PESSTO uses the New Technology Telescope with EFOSC2 and SOFI to provide optical and NIR spectroscopy and imaging. We target supernovae and optical transients brighter than 20.5mag for classification. Science targets are then selected for follow-up based on the PESSTO science goal of extending knowledge of the extremes of the supernova population. The EFOSC2 spectra cover 3345-9995A (at resolutions of 13-18 Angs) and SOFI spectra cover 0.935-2.53 micron (resolutions 23-33 Angs) along with JHK imaging. This data release contains spectra from the first year (April 2012 - 2013), consisting of all 814 EFOSC2 spectra and 95 SOFI spectra (covering 298 distinct objects), in standard ESO Phase 3 format. We estimate the accuracy of the absolute flux calibrations for EFOSC2 to be typically 15%, and the relative flux calibration accuracy to be about 5%. The PESSTO standard NIR reduction process does not yet produce high accuracy absolute spectrophotometry but the SOFI JHK imaging will improve this. Future data releases will focus on improving the automated flux calibration of the data products.

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.

The Pan-STARRS1 (PS1) survey has obtained imaging in 5 bands (grizy_P1) over 10 Medium Deep Survey (MDS) fields covering a total of 70 square degrees. This paper describes the search for apparently hostless supernovae (SNe) within the first year of PS1 MDS data with an aim of discovering new superluminous supernovae (SLSNe). A total of 249 hostless transients were discovered down to a limiting magnitude of M_AB ~ 23.5, of which 76 were classified as Type Ia SNe. There were 57 SNe with complete light curves that are likely core-collapse SNe (CCSNe) or SLSNe and 12 of these have had spectra taken. Of these 12 hostless, non-Type Ia SNe, 7 were SLSNe of Type Ic at redshifts between 0.5-1.4. This illustrates that the discovery rate of Type Ic SLSNe can be maximised by concentrating on hostless transients and removing normal SNe Ia. We present data for two new possible SLSNe; PS1-10pm (z = 1.206) and PS1-10ahf (z = 1.1), and estimate the rate of SLSNe-Ic to be between 3^{+3}_{-2} * 10^{-5} and 8^{+2}_{-1} * 10^{-5} of the CCSNe rate within 0.3 <= z <= 1.4 by applying a Monte-Carlo technique. The rate of slowly evolving, SN2007bi-like explosions is estimated as a factor of 10 lower than this range.

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 optical spectra and light curves for three hydrogen-poor super-luminous supernovae followed by the Public ESO Spectroscopic Survey of Transient Objects (PESSTO). Time series spectroscopy from a few days after maximum light to 100 days later shows them to be fairly typical of this class, with spectra dominated by Ca II, Mg II, Fe II and Si II, which evolve slowly over most of the post-peak photospheric phase. We determine bolometric light curves and apply simple fitting tools, based on the diffusion of energy input by magnetar spin-down, 56Ni decay, and collision of the ejecta with an opaque circumstellar shell. We investigate how the heterogeneous light curves of our sample (combined with others from the literature) can help to constrain the possible mechanisms behind these events. We have followed these events to beyond 100-200 days after peak, to disentangle host galaxy light from fading supernova flux and to differentiate between the models, which predict diverse behaviour at this phase. Models powered by radioactivity require unrealistic parameters to reproduce the observed light curves, as found by previous studies. Both magnetar heating and circumstellar interaction still appear to be viable candidates. A large diversity is emerging in observed tail-phase luminosities, with magnetar models failing in some cases to predict the rapid drop in flux. This would suggest either that magnetars are not responsible, or that the X-ray flux from the magnetar wind is not fully trapped. The light curve of one object shows a distinct re-brightening at around 100d after maximum light. We argue that this could result either from multiple shells of circumstellar material, or from a magnetar ionisation front breaking out of the ejecta.

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.

We present optical photometric and spectroscopic coverage of the superluminous supernova (SLSN) PS1-11ap, discovered with the Pan-STARRS1 Medium Deep Survey at z = 0.524. This intrinsically blue transient rose slowly to reach a peak magnitude of M_u = -21.4 mag and bolometric luminosity of 8 x 10^43 ergs^-1 before settling onto a relatively shallow gradient of decline. The observed decline is significantly slower than those of the superluminous type Ic SNe which have been the focus of much recent attention. Spectroscopic similarities with the lower redshift SN2007bi and a decline rate similar to 56Co decay timescale initially indicated that this transient could be a candidate for a pair instability supernova (PISN) explosion. Overall the transient appears quite similar to SN2007bi and the lower redshift object PTF12dam. The extensive data set, from 30 days before peak to 230 days after, allows a detailed and quantitative comparison with published models of PISN explosions. We find that the PS1-11ap data do not match these model explosion parameters well, supporting the recent claim that these SNe are not pair instability explosions. We show that PS1-11ap has many features in common with the faster declining superluminous Ic supernovae and the lightcurve evolution can also be quantitatively explained by the magnetar spin down model. At a redshift of z = 0.524 the observer frame optical coverage provides comprehensive restframe UV data and allows us to compare it with the superluminous SNe recently found at high redshifts between z = 2-4. While these high-z explosions are still plausible PISN candidates, they match the photometric evolution of PS1-11ap and hence could be counterparts to this lower redshift transient.

We present optical and infrared monitoring data of SN 2012hn collected by the Public ESO Spectroscopic Survey for Transient Objects (PESSTO). We show that SN 2012hn has a faint peak magnitude (MR ~ -15.7) and shows no hydrogen and no clear evidence for helium in its spectral evolution. Instead, we detect prominent Ca II lines at all epochs, which relates this transient to previously described 'Ca-rich' or 'gap' transients. However, the photospheric spectra (from -3 to +32 d with respect to peak) of SN 2012hn show a series of absorption lines which are unique, and a red continuum that is likely intrinsic rather than due to extinction. Lines of Ti II and Cr II are visible. This may be a temperature effect, which could also explain the red photospheric colour. A nebular spectrum at +150d shows prominent CaII, OI, CI and possibly MgI lines which appear similar in strength to those displayed by core-collapse SNe. To add to the puzzle, SN 2012hn is located at a projected distance of 6 kpc from an E/S0 host and is not close to any obvious starforming region. Overall SN 2012hn resembles a group of faint H-poor SNe that have been discovered recently and for which a convincing and consistent physical explanation is still missing. They all appear to explode preferentially in remote locations offset from a massive host galaxy with deep limits on any dwarf host galaxies, favouring old progenitor systems. SN 2012hn adds heterogeneity to this sample of objects. We discuss potential explosion channels including He-shell detonations and double detonations of white dwarfs as well as peculiar core-collapse SNe.

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 report the results of a 3 year-long dedicated monitoring campaign of a restless Luminous Blue Variable (LBV) in NGC 7259. The object, named SN 2009ip, was observed photometrically and spectroscopically in the optical and near-infrared domains. We monitored a number of erupting episodes in the past few years, and increased the density of our observations during eruptive episodes. In this paper we present the full historical data set from 2009-2012 with multi-wavelength dense coverage of the two high luminosity events between August - September 2012. We construct bolometric light curves and measure the total luminosities of these eruptive or explosive events. We label them the 2012a event (lasting ~50 days) with a peak of 3x10^41 erg/s, and the 2012b event (14 day rise time, still ongoing) with a peak of 8x10^42 erg/s. The latter event reached an absolute R-band magnitude of about -18, comparable to that of a core-collapse supernova (SN). Our historical monitoring has detected high-velocity spectral features (~13000 km/s) in September 2011, one year before the current SN-like event. This implies that the detection of such high velocity outflows cannot, conclusively, point to a core-collapse SN origin. We suggest that the initial peak in the 2012a event was unlikely to be due to a faint core-collapse SN. We propose that the high intrinsic luminosity of the latest peak, the variability history of SN 2009ip, and the detection of broad spectral lines indicative of high-velocity ejecta are consistent with a pulsational pair-instability event, and that the star may have survived the last outburst. The question of the survival of the LBV progenitor star and its future fate remain open issues, only to be answered with future monitoring of this historically unique explosion.

Super-luminous supernovae have a tendency to occur in faint host galaxies which are likely to have low mass and low metallicity. While these extremely luminous explosions have been observed from z=0.1 to 1.55, the closest explosions allow more detailed investigations of their host galaxies. We present a detailed analysis of the host galaxy of SN 2010gx (z=0.23), one of the best studied super-luminous type Ic supernovae. The host is a dwarf galaxy (M_g=-17.42+/-0.17) with a high specific star formation rate. It has a remarkably low metallicity of 12+log(O/H)=7.5+/-0.1 dex as determined from the detection of the [OIII] 4363 Angs line. This is the first reliable metallicity determination of a super-luminous stripped-envelope supernova host. We collected deep multi-epoch imaging with Gemini + GMOS between 240-560 days after explosion to search for any sign of radioactive nickel-56, which might provide further insights on the explosion mechanism and the progenitor's nature. We reach griz magnitudes of m_AB~26, but do not detect SN 2010gx at these epochs. The limit implies that any nickel-56 production was similar to or below that of SN 1998bw (a luminous type Ic SN that produced around 0.4 M_sun of nickel-56). The low volumetric rates of these supernovae (~10^-4 of the core-collapse population) could be qualitatively matched if the explosion mechanism requires a combination of low-metallicity (below 0.2 Z_sun), high progenitor mass (>60 M_sun) and high rotation rate (fastest 10% of rotators).

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.