While microlensing is very rare, occurring on average once per million stars observed, current and near-future surveys are coming online with the capability of providing photometry of almost the entire visible sky to depths up to R ~ 22 mag or fainter every few days, which will contribute to the detection of black holes and exoplanets through follow-up observations of microlensing events. Based on galactic models, we can expect microlensing events across a vastly wider region of the galaxy, although the cadence of these surveys (2-3 per day ) is lower than traditional microlensing surveys, making efficient detection a challenge. Rapid advances are being made in the utility of time-series data to detect and classify transient events in real-time using very high data-rate surveys, but limited work has been published regarding the detection of microlensing events, particularly for when the data streams are of relatively low-cadence. In this research, we explore the utility of a Random Forest algorithm for identifying microlensing signals using time-series data, with the goal of creating an efficient machine learning classifier that can be applied to search for microlensing in wide-field surveys even with low-cadence data. We have applied and optimized our classifier using the OGLE-II microlensing dataset, in addition to testing with PTF/iPTF survey data and the currently operating ZTF, which applies the same data handling infrastructure that is envisioned for the upcoming LSST.

With the advent of the Large Synoptic Survey Telescope (LSST), time-domain astronomy will be faced with an unprecedented volume and rate of data. Real-time processing of variables and transients detected by such large-scale surveys is critical to identifying the more unusual events and allocating scarce follow-up resources efficiently. We develop an algorithm to identify these novel events within a given population of variable sources. We determine the distributions of magnitude changes (dm) over time intervals (dt) for a given passband f, pf(dm|dt), and use these distributions to compute the likelihood of a test source being consistent with the population, or an outlier. We demonstrate our algorithm by applying it to the DECam multi-band time-series data of more than 2000 variable stars identified by Saha et al. (2019) in the Galactic Bulge that are largely dominated by long-period variables and pulsating stars. Our algorithm discovers 18 outlier sources in the sample, including a microlensing event, a dwarf nova, and two chromospherically active RS CVn stars, as well as sources in the Blue Horizontal Branch region of the color-magnitude diagram without any known counterparts. We compare the performance of our algorithm for novelty detection with multivariate KDE and Isolation Forest on the simulated PLAsTiCC dataset. We find that our algorithm yields comparable results despite its simplicity. Our method provides an efficient way for flagging the most unusual events in a real-time alert-broker system.

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.

A Kavli foundation sponsored workshop on the theme \emph{Petabytes to Science} was held 12$^{th}$ to 14$^{th}$ of February 2019 in Las Vegas. The aim of the this workshop was to discuss important trends and technologies which may support astronomy. We also tackled how to better shape the workforce for the new trends and how we should approach education and public outreach. This document was coauthored during the workshop and edited in the weeks after. It comprises the discussions and highlights many recommendations which came out of the workshop. We shall distill parts of this document and formulate potential white papers for the decadal survey.

The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our technical and collaborative frameworks to promote efficient algorithmic development and take advantage of opportunities for scientific discovery in the petabyte era. We discuss challenges for discovery in large and complex data sets; challenges and requirements for the next stage of development of statistical methodologies and algorithmic tool sets; how we might change our paradigms of collaboration and education; and the ethical implications of scientists' contributions to widely applicable algorithms and computational modeling. We start with six distinct recommendations that are supported by the commentary following them. This white paper is related to a larger corpus of effort that has taken place within and around the Petabytes to Science Workshops (https://petabytestoscience.github.io/).

Fundamental coding and software development skills are increasingly necessary for success in nearly every aspect of astronomical and astrophysical research as large surveys and high resolution simulations become the norm. However, professional training in these skills is inaccessible or impractical for many members of our community. Students and professionals alike have been expected to acquire these skills on their own, apart from formal classroom curriculum or on-the-job training. Despite the recognized importance of these skills, there is little opportunity to develop them - even for interested researchers. To ensure a workforce capable of taking advantage of the computational resources and the large volumes of data coming in the next decade, we must identify and support ways to make software development training widely accessible to community members, regardless of affiliation or career level. To develop and sustain a technology capable astronomical and astrophysical workforce, we recommend that agencies make funding and other resources available in order to encourage, support and, in some cases, require progress on necessary training, infrastructure and policies. In this white paper, we focus on recommendations for how funding agencies can lead in the promotion of activities to support the astronomy and astrophysical workforce in the 2020s.

We present RAPID (Real-time Automated Photometric IDentification), a novel time-series classification tool capable of automatically identifying transients from within a day of the initial alert, to the full lifetime of a light curve. Using a deep recurrent neural network with Gated Recurrent Units (GRUs), we present the first method specifically designed to provide early classifications of astronomical time-series data, typing 12 different transient classes. Our classifier can process light curves with any phase coverage, and it does not rely on deriving computationally expensive features from the data, making RAPID well-suited for processing the millions of alerts that ongoing and upcoming wide-field surveys such as the Zwicky Transient Facility (ZTF), and the Large Synoptic Survey Telescope (LSST) will produce. The classification accuracy improves over the lifetime of the transient as more photometric data becomes available, and across the 12 transient classes, we obtain an average area under the receiver operating characteristic curve of 0.95 and 0.98 at early and late epochs, respectively. We demonstrate RAPID's ability to effectively provide early classifications of observed transients from the ZTF data stream. We have made RAPID available as an open-source software package (https://astrorapid.readthedocs.io) for machine learning-based alert-brokers to use for the autonomous and quick classification of several thousand light curves within a few seconds.

Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era of astronomical science, citing software and measuring its impact enables academia to retain and reward researchers that make significant software contributions. These highly skilled researchers must be retained to maximize the scientific return from petabyte-scale datasets. Evolving beyond the one-dimensional credit model requires overcoming several key challenges, including the current scholarly ecosystem and scientific culture issues. This white paper will present these challenges and suggest practical solutions for elevating the role of software as a product of the research enterprise.

We describe the simulated data sample for the "Photometric LSST Astronomical Time Series Classification Challenge" (PLAsTiCC), a publicly available challenge to classify transient and variable events that will be observed by the Large Synoptic Survey Telescope (LSST), a new facility expected to start in the early 2020s. The challenge was hosted by Kaggle, ran from 2018 September 28 to 2018 December 17, and included 1,094 teams competing for prizes. Here we provide details of the 18 transient and variable source models, which were not revealed until after the challenge, and release the model libraries at https://doi.org/10.5281/zenodo.2612896. We describe the LSST Operations Simulator used to predict realistic observing conditions, and we describe the publicly available SNANA simulation code used to transform the models into observed fluxes and uncertainties in the LSST passbands (ugrizy). Although PLAsTiCC has finished, the publicly available models and simulation tools are being used within the astronomy community to further improve classification, and to study contamination in photometrically identified samples of type Ia supernova used to measure properties of dark energy. Our simulation framework will continue serving as a platform to improve the PLAsTiCC models, and to develop new models.

Measurements of the dark energy equation-of-state parameter, $w$, have been limited by uncertainty in the selection effects and photometric calibration of $z<0.1$ Type Ia supernovae (SNe Ia). The Foundation Supernova Survey is designed to lower these uncertainties by creating a new sample of $z<0.1$ SNe Ia observed on the Pan-STARRS system. Here, we combine the Foundation sample with SNe from the Pan-STARRS Medium Deep Survey and measure cosmological parameters with 1,338 SNe from a single telescope and a single, well-calibrated photometric system. For the first time, both the low-$z$ and high-$z$ data are predominantly discovered by surveys that do not target pre-selected galaxies, reducing selection bias uncertainties. The $z>0.1$ data include 875 SNe without spectroscopic classifications and we show that we can robustly marginalize over CC SN contamination. We measure Foundation Hubble residuals to be fainter than the pre-existing low-$z$ Hubble residuals by $0.046 \pm 0.027$ mag (stat+sys). By combining the SN Ia data with cosmic microwave background constraints, we find $w=-0.938 \pm 0.053$, consistent with $\Lambda$CDM. With 463 spectroscopically classified SNe Ia alone, we measure $w=-0.933\pm0.061$. Using the more homogeneous and better-characterized Foundation sample gives a 55% reduction in the systematic uncertainty attributed to SN Ia sample selection biases. Although use of just a single photometric system at low and high redshift increases the impact of photometric calibration uncertainties in this analysis, previous low-$z$ samples may have correlated calibration uncertainties that were neglected in past studies. The full Foundation sample will observe up to 800 SNe to anchor the LSST and WFIRST Hubble diagrams.

We identify minimal observing cadence requirements that enable photometric astronomical surveys to detect and recognize fast and explosive transients and fast transient features. Observations in two different filters within a short time window (e.g., g-and-i, or r-and-z, within < 0.5 hr) and a repeat of one of those filters with a longer time window (e.g., > 1.5 hr) are desirable for this purpose. Such an observing strategy delivers both the color and light curve evolution of transients on the same night. This allows the identification and initial characterization of fast transient -- or fast features of longer timescale transients -- such as rapidly declining supernovae, kilonovae, and the signatures of SN ejecta interacting with binary companion stars or circumstellar material. Some of these extragalactic transients are intrinsically rare and generally all hard to find, thus upcoming surveys like the Large Synoptic Survey Telescope (LSST) could dramatically improve our understanding of their origin and properties. We colloquially refer to such a strategy implementation for the LSST as the Presto-Color strategy (rapid-color). This cadence's minimal requirements allow for overall optimization of a survey for other science goals.

The Large Synoptic Survey Telescope (LSST) will open a discovery frontier for faint and fast transients with its ability to detect variable flux components down to $\sim$24.5 mag in a $\sim$30 second exposure. Spectroscopic follow-up of such phenomena - which are necessary for understanding the physics of stellar explosions - can require a rapid response and several hours with a 8-10m telescope, making it both expensive and difficult to acquire. The future Extremely Large Telescopes (ELTs) would be able to provide not only spectroscopy but capabilities such as spectropolarimetry and high-resolution diffraction-limited imaging that would contribute to future advances in our physical understanding of stellar explosions. In this white paper we focus on several specific scientific impacts in the field of explosive transient astrophysics that will be generated by the combination of LSST's discovery abilities and ELTs' follow-up capacities. First, we map the uncharted frontier of discovery phase-space in terms of intrinsic luminosity and timescales for explosive transients, where we expect the unexpected. We then focus on six areas with open science questions for known transients: the progenitors of thermonuclear supernovae (SNe), mass loss prior to core collapse, asymmetry in stellar explosions, light echoes, high-$z$ transients, and strongly lensed SNe. We conclude with a brief discussion of the practical aspects of ELT & LSST synergy.

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

The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

We have obtained repeated images of 6 fields towards the Galactic bulge in 5 passbands (u, g, r, i, z) with the DECam imager on the Blanco 4m telescope at CTIO. From over 1.6 billion individual photometric measurements in the field centered on Baade's window, we have detected 4877 putative variable stars. 474 of these have been confirmed as fundamental mode RR Lyrae stars, whose colors at minimum light yield line-of-sight reddening determinations as well as a reddening law towards the Galactic Bulge which differs significantly from the standard R_V = 3.1 formulation. Assuming that the stellar mix is invariant over the 3 square-degree field, we are able to derive a line-of-sight reddening map with sub-arcminute resolution, enabling us to obtain de-reddened and extinction corrected color-magnitude diagrams (CMD's) of this bulge field using up to 2.5 million well-measured stars. The corrected CMD's show unprecedented detail and expose sparsely populated sequences: e.g., delineation of the very wide red giant branch, structure within the red giant clump, the full extent of the horizontal branch, and a surprising bright feature which is likely due to stars with ages younger than 1 Gyr. We use the RR Lyrae stars to trace the spatial structure of the ancient stars, and find an exponential decline in density with Galactocentric distance. We discuss ways in which our data products can be used to explore the age and metallicity properties of the bulge, and how our larger list of all variables is useful for learning to interpret future LSST alerts.

We have established a network of 19 faint (16.5 mag $< V < $19 mag) northern and equatorial DA white dwarfs as spectrophotometric standards for present and future wide-field observatories. Our analysis infers SED models for the stars that are tied to the three CALSPEC primary standards. Our SED models are consistent with panchromatic Hubble Space Telescope ($HST$) photometry to better than 1%. The excellent agreement between observations and models validates the use of non-local-thermodynamic-equilibrium (NLTE) DA white dwarf atmospheres extinguished by interstellar dust as accurate spectrophotometric references. Our standards are accessible from both hemispheres and suitable for ground and space-based observatories covering the ultraviolet to the near infrared. The high-precision of these faint sources make our network of standards ideally suited for any experiment that has very stringent requirements on flux calibration, such as studies of dark energy using the Large Synoptic Survey Telescope (LSST) and the Wide-Field Infrared Survey Telescope ($WFIRST$).

This paper describes the organization of the database and the catalog data products from the Pan-STARRS1 $3\pi$ Steradian Survey. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI. The database is described in detail, including the construction of the database, the provenance of the data, the schema, and how the database tables are related. Examples of queries for a range of science goals are included. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI.

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

We present precise photometry and spectroscopy for 23 candidate spectrophotometric standard white dwarfs. The selected stars are distributed in the Northern hemisphere and around the celestial equators and are all fainter than r ~ 16.5 mag. This network of stars, when established as standards, together with the three Hubble Space Telescope primary CALSPEC white dwarfs, will provide a set of spectrophotometric standards to directly calibrate data products to better than 1%. These new faint standard white dwarfs will have enough signal-to-noise ratio in future deep photometric surveys and facilities to be measured accurately while still avoiding saturation in such surveys. They will also fall within the dynamic range of large telescopes and their instruments for the foreseeable future. This paper discusses the provenance of the observational data for our candidate standard stars. The comparison with models, reconciliation with reddening, and the consequent derivation of the full spectral energy density distributions for each of them is reserved for a subsequent paper.

Type Iax supernovae (SNe~Iax) are the most common class of peculiar SNe. While they are thought to be thermonuclear white-dwarf (WD) SNe, SNe~Iax are observationally similar to, but distinct from SNe~Ia. Unlike SNe~Ia, where roughly 30\% occur in early-type galaxies, only one SN~Iax has been discovered in an early-type galaxy, suggesting a relatively short delay time and a distinct progenitor system. Furthermore, one SN~Iax progenitor system has been detected in pre-explosion images with its properties consistent with either of two models: a short-lived (<100 Myr) progenitor system consisting of a WD primary and a He-star companion, or a singular Wolf-Rayet progenitor star. Using deep \textit{Hubble Space Telescope} images of nine nearby SN~Iax host galaxies, we measure the properties of stars within 200 pc of the SN position. The ages of local stars, some of which formed with the SN progenitor system, can constrain the time between star formation and SN, known as the delay time. We compare the local stellar properties to synthetic photometry of single-stellar populations, fitting to a range of possible delay times for each SN. With this sample, we uniquely constrain the delay-time distribution for SNe~Iax, with a median and $1-\sigma$ confidence interval delay time of $63_{- 15}^{+ 58} \times 10^{6}$ years. The measured delay-time distribution provides an excellent constraint on the progenitor system for the class, indicating a preference for a WD progenitor system over a Wolf-Rayet progenitor star.

The discovery of the electromagnetic counterparts to the binary neutron star merger GW170817 has opened the era of GW+EM multi-messenger astronomy. Exploiting this breakthrough requires increasing samples to explore the diversity of kilonova behaviour and provide more stringent constraints on the Hubble constant, and tests of fundamental physics. LSST can play a key role in this field in the 2020s, when the gravitational wave detector network is expected to detect higher rates of merger events involving neutron stars ($\sim$10s per year) out to distances of several hundred Mpc. Here we propose comprehensive target-of-opportunity (ToOs) strategies for follow-up of gravitational-wave sources that will make LSST the premiere machine for discovery and early characterization for neutron star mergers and other gravitational-wave sources.

Supernova (SN) 2018oh (ASASSN-18bt) is the first spectroscopically-confirmed type Ia supernova (SN Ia) observed in the $Kepler$ field. The $Kepler$ data revealed an excess emission in its early light curve, allowing to place interesting constraints on its progenitor system (Dimitriadis et al. 2018, Shappee et al. 2018b). Here, we present extensive optical, ultraviolet, and near-infrared photometry, as well as dense sampling of optical spectra, for this object. SN 2018oh is relatively normal in its photometric evolution, with a rise time of 18.3$\pm$0.3 days and $\Delta$m$_{15}(B)=0.96\pm$0.03 mag, but it seems to have bluer $B - V$ colors. We construct the "uvoir" bolometric light curve having peak luminosity as 1.49$\times$10$^{43}$erg s$^{-1}$, from which we derive a nickel mass as 0.55$\pm$0.04M$_{\odot}$ by fitting radiation diffusion models powered by centrally located $^{56}$Ni. Note that the moment when nickel-powered luminosity starts to emerge is +3.85 days after the first light in the Kepler data, suggesting other origins of the early-time emission, e.g., mixing of $^{56}$Ni to outer layers of the ejecta or interaction between the ejecta and nearby circumstellar material or a non-degenerate companion star. The spectral evolution of SN 2018oh is similar to that of a normal SN Ia, but is characterized by prominent and persistent carbon absorption features. The C II features can be detected from the early phases to about 3 weeks after the maximum light, representing the latest detection of carbon ever recorded in a SN Ia. This indicates that a considerable amount of unburned carbon exists in the ejecta of SN 2018oh and may mix into deeper layers.

We present an exquisite, 30-min cadence Kepler (K2) light curve of the Type Ia supernova (SN Ia) 2018oh (ASASSN-18bt), starting weeks before explosion, covering the moment of explosion and the subsequent rise, and continuing past peak brightness. These data are supplemented by multi-color Pan-STARRS1 and CTIO 4-m DECam observations obtained within hours of explosion. The K2 light curve has an unusual two-component shape, where the flux rises with a steep linear gradient for the first few days, followed by a quadratic rise as seen for typical SNe Ia. This "flux excess" relative to canonical SN Ia behavior is confirmed in our $i$-band light curve, and furthermore, SN 2018oh is especially blue during the early epochs. The flux excess peaks 2.14$\pm0.04$ days after explosion, has a FWHM of 3.12$\pm0.04$ days, a blackbody temperature of $T=17,500^{+11,500}_{-9,000}$ K, a peak luminosity of $4.3\pm0.2\times10^{37}\,{\rm erg\,s^{-1}}$, and a total integrated energy of $1.27\pm0.01\times10^{43}\,{\rm erg}$. We compare SN 2018oh to several models that may provide additional heating at early times, including collision with a companion and a shallow concentration of radioactive nickel. While all of these models generally reproduce the early K2 light curve shape, we slightly favor a companion interaction, at a distance of $\sim$$2\times10^{12}\,{\rm cm}$ based on our early color measurements, although the exact distance depends on the uncertain viewing angle. Additional confirmation of a companion interaction in future modeling and observations of SN 2018oh would provide strong support for a single-degenerate progenitor system.

On 2018 Feb. 4.41, the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered ASASSN-18bt in the K2 Campaign 16 field. With a redshift of z=0.01098 and a peak apparent magnitude of B_{max}=14.31, ASASSN-18bt is the nearest and brightest SNe Ia yet observed by the Kepler spacecraft. Here we present the discovery of ASASSN-18bt, the K2 light curve, and pre-discovery data from ASAS-SN and the Asteroid Terrestrial-impact Last Alert System (ATLAS). The K2 early-time light curve has an unprecedented 30-minute cadence and photometric precision for an SN~Ia light curve, and it unambiguously shows a ~4 day nearly linear phase followed by a steeper rise. Thus, ASASSN-18bt joins a growing list of SNe Ia whose early light curves are not well described by a single power law. We show that a double-power-law model fits the data reasonably well, hinting that two physical processes must be responsible for the observed rise. However, we find that current models of the interaction with a non-degenerate companion predict an abrupt rise and cannot adequately explain the initial, slower linear phase. Instead, we find that existing, published models with shallow 56Ni are able to span the observed behavior and, with tuning, may be able to reproduce the ASASSN-18bt light curve. Regardless, more theoretical work is needed to satisfactorily model this and other early-time SNe~Ia light curves. Finally, we use Swift X-ray non-detections to constrain the presence of circumstellar material (CSM) at much larger distances and lower densities than possible with the optical light curve. For a constant density CSM these non-detections constrain rho<4.5 * 10^5 cm^-3 at a radius of 4 *10^15 cm from the progenitor star. Assuming a wind-like environment, we place mass-loss limits of Mdot< 8 * 10^-6 M_sun yr^-1 for v_w=100 km s^-1, ruling out some symbiotic progenitor systems.

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

Classification of transient and variable light curves is an essential step in using astronomical observations to develop an understanding of their underlying physical processes. However, upcoming deep photometric surveys, including the Large Synoptic Survey Telescope (LSST), will produce a deluge of low signal-to-noise data for which traditional labeling procedures are inappropriate. Probabilistic classification is more appropriate for the data but are incompatible with the traditional metrics used on deterministic classifications. Furthermore, large survey collaborations intend to use these classification probabilities for diverse science objectives, indicating a need for a metric that balances a variety of goals. We describe the process used to develop an optimal performance metric for an open classification challenge that seeks probabilistic classifications and must serve many scientific interests. The Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC) is an open competition aiming to identify promising techniques for obtaining classification probabilities of transient and variable objects by engaging a broader community both within and outside astronomy. Using mock classification probability submissions emulating archetypes of those anticipated of PLAsTiCC, we compare the sensitivity of metrics of classification probabilities under various weighting schemes, finding that they yield qualitatively consistent results. We choose as a metric for PLAsTiCC a weighted modification of the cross-entropy because it can be meaningfully interpreted. Finally, we propose extensions of our methodology to ever more complex challenge goals and suggest some guiding principles for approaching the choice of a metric of probabilistic classifications.

The Photometric LSST Astronomical Time Series Classification Challenge (PLAsTiCC) is an open data challenge to classify simulated astronomical time-series data in preparation for observations from the Large Synoptic Survey Telescope (LSST), which will achieve first light in 2019 and commence its 10-year main survey in 2022. LSST will revolutionize our understanding of the changing sky, discovering and measuring millions of time-varying objects. In this challenge, we pose the question: how well can we classify objects in the sky that vary in brightness from simulated LSST time-series data, with all its challenges of non-representativity? In this note we explain the need for a data challenge to help classify such astronomical sources and describe the PLAsTiCC data set and Kaggle data challenge, noting that while the references are provided for context, they are not needed to participate in the challenge.

The past year has witnessed discovery of the first identified counterparts to a gravitational wave transient (GW 170817A) and a very high-energy neutrino (IceCube-170922A). These source identifications, and ensuing detailed studies, have realized longstanding dreams of astronomers and physicists to routinely carry out observations of cosmic sources by other than electromagnetic means, and inaugurated the era of "multi-messenger" astronomy. While this new era promises extraordinary physical insights into the universe, it brings with it new challenges, including: highly heterogeneous, high-volume, high-velocity datasets; globe-spanning cross-disciplinary teams of researchers, regularly brought together into transient collaborations; an extraordinary breadth and depth of domain-specific knowledge and computing resources required to anticipate, model, and interpret observations; and the routine need for adaptive, distributed, rapid-response observing campaigns to fully exploit the scientific potential of each source. We argue, therefore, that the time is ripe for the community to conceive and propose an Institute for Multi-Messenger Astrophysics that would coordinate its resources in a sustained and strategic fashion to efficiently address these challenges, while simultaneously serving as a center for education and key supporting activities. In this fashion, we can prepare now to realize the bright future that we see, beyond, through these newly opened windows onto the universe.

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.

Much of the progress made in time-domain astronomy is accomplished by relating observational multi-wavelength time series data to models derived from our understanding of physical laws. This goal is typically accomplished by dividing the task in two: collecting data (observing), and constructing models to represent that data (theorizing). Owing to the natural tendency for specialization, a disconnect can develop between the best available theories and the best available data, potentially delaying advances in our understanding new classes of transients. We introduce MOSFiT: the Modular Open-Source Fitter for Transients, a Python-based package that downloads transient datasets from open online catalogs (e.g., the Open Supernova Catalog), generates Monte Carlo ensembles of semi-analytical light curve fits to those datasets and their associated Bayesian parameter posteriors, and optionally delivers the fitting results back to those same catalogs to make them available to the rest of the community. MOSFiT is designed to help bridge the gap between observations and theory in time-domain astronomy; in addition to making the application of existing models and creation of new models as simple as possible, MOSFiT yields statistically robust predictions for transient characteristics, with a standard output format that includes all the setup information necessary to reproduce a given result. As large-scale surveys such as LSST discover entirely new classes of transients, tools such as MOSFiT will be critical for enabling rapid comparison of models against data in statistically consistent, reproducible, and scientifically beneficial ways.

The unprecedented volume and rate of transient events that will be discovered by the Large Synoptic Survey Telescope (LSST) demands that the astronomical community update its followup paradigm. Alert-brokers -- automated software system to sift through, characterize, annotate and prioritize events for followup -- will be critical tools for managing alert streams in the LSST era. The Arizona-NOAO Temporal Analysis and Response to Events System (ANTARES) is one such broker. In this work, we develop a machine learning pipeline to characterize and classify variable and transient sources only using the available multiband optical photometry. We describe three illustrative stages of the pipeline, serving the three goals of early, intermediate and retrospective classification of alerts. The first takes the form of variable vs transient categorization, the second, a multi-class typing of the combined variable and transient dataset, and the third, a purity-driven subtyping of a transient class. While several similar algorithms have proven themselves in simulations, we validate their performance on real observations for the first time. We quantitatively evaluate our pipeline on sparse, unevenly sampled, heteroskedastic data from various existing observational campaigns, and demonstrate very competitive classification performance. We describe our progress towards adapting the pipeline developed in this work into a real-time broker working on live alert streams from time-domain surveys.

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.

We characterize the absolute magnitudes and colors of RR Lyrae stars in the globular cluster M5 in the ugriz filter system of the Dark Energy Camera (DECam). We provide empirical Period-Luminosity (P-L) relationships in all 5 bands based on 47 RR Lyrae stars of the type ab and 14 stars of the type c. The P-L relationships were found to be better constrained for the fundamental mode RR Lyrae stars in the riz passbands, with dispersion of 0.03, 0.02 and 0.02 magnitudes, respectively. The dispersion of the color at minimum light was found to be small, supporting the use of this parameter as a means to obtain accurate interstellar extinctions along the line of sight up to the distance of the RR Lyrae star. We found a trend of color at minimum light with pulsational period that, if taken into account, brings the dispersion in color at minimum light to < 0.016 magnitudes for the (r-i), (i-z), and (r-z) colors. These calibrations will be very useful for using RR Lyrae stars from DECam observations as both standard candles for distance determinations and color standards for reddening measurements.

The Arizona-NOAO Temporal Analysis and Response to Events System (ANTARES) is a joint effort of NOAO and the Department of Computer Science at the University of Arizona to build prototype software to process alerts from time-domain surveys, especially LSST, to identify those alerts that must be followed up immediately. Value is added by annotating incoming alerts with existing information from previous surveys and compilations across the electromagnetic spectrum and from the history of past alerts. Comparison against a knowledge repository of properties and features of known or predicted kinds of variable phenomena is used for categorization. The architecture and algorithms being employed are described.

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 analyze the wavelength-dependent variability of a sample of spectroscopically confirmed active galactic nuclei (AGN) selected from near-UV ($NUV$) variable sources in the GALEX Time Domain Survey that have a large amplitude of optical variability (difference-flux S/N $>$ 3) in the Pan-STARRS1 Medium Deep Survey (PS1 MDS). By matching GALEX and PS1 epochs in 5 bands ($NUV$, $g_{P1}$, $r_{P1}$, $i_{P1}$, $z_{P1}$) in time, and taking their flux difference, we create co-temporal difference-flux spectral energy distributions ($\Delta f$SEDs) using two chosen epochs for each of the 23 objects in our sample on timescales of about a year. We confirm the "bluer-when-brighter" trend reported in previous studies, and measure a median spectral index of the $\Delta f$SEDs of $\alpha_{\lambda}$ = 2.1 that is consistent with an accretion disk spectrum. We further fit the $\Delta f$SEDs of each source with a standard accretion disk model in which the accretion rate changes from one epoch to the other. In our sample, 17 out of 23 ($\sim$74 %) sources are well described by this variable accretion-rate disk model, with a median average characteristic disk temperature $\bar{T}^*$ of $1.2\times 10^5$~K that is consistent with the temperatures expected given the distribution of accretion rates and black hole masses inferred for the sample. Our analysis also shows that the variable accretion rate model is a better fit to the $\Delta f$SEDs than a simple power law.

We present initial results from a program aimed at establishing a network of hot DA white dwarfs to serve as spectrophotometric standards for present and future wide-field surveys. These stars span the equatorial zone and are faint enough to be conveniently observed throughout the year with large-aperture telescopes. Spectra of these white dwarfs are analyzed to generate a non-local-thermodynamic-equilibrium (NLTE) model atmosphere normalized to HST colors, including adjustments for wavelength-dependent interstellar extinction. Once established, this standard star network will serve ground-based observatories in both hemispheres as well as space-based instrumentation from the UV to the near IR. We demonstrate the effectiveness of this concept and show how two different approaches to the problem using somewhat different assumptions produce equivalent results. We discuss lessons learned and the resulting corrective actions applied to our program.

The ESSENCE survey discovered 213 Type Ia supernovae at redshifts 0.1 < z < 0.81 between 2002 and 2008. We present their R and I-band photometry, measured from images obtained using the MOSAIC II camera at the CTIO 4 m Blanco telescope, along with rapid-response spectroscopy for each object. We use our spectroscopic follow-up observations to determine an accurate, quantitative classification and a precise redshift. Through an extensive calibration program we have improved the precision of the CTIO Blanco natural photometric system. We use several empirical metrics to measure our internal photometric consistency and our absolute calibration of the survey. We assess the effect of various potential sources of systematic bias on our measured fluxes, and we estimate that the dominant term in the systematic error budget from the photometric calibration on our absolute fluxes is ~1%.

We present late-time observations of the tidal disruption event candidate PS1-10jh. UV and optical imaging with HST/WFC3 localize the transient to be coincident with the host galaxy nucleus to an accuracy of 0.023 arcsec, corresponding to 66 pc. The UV flux in the F225W filter, measured 3.35 rest-frame years after the peak of the nuclear flare, is consistent with a decline that continues to follow a $t^{-5/3}$ power-law with no spectral evolution. Late epochs of optical spectroscopy obtained with MMT ~ 2 and 4 years after the peak, enable a clean subtraction of the host galaxy from the early spectra, revealing broad helium emission lines on top of a hot continuum, and placing stringent upper limits on the presence of hydrogen line emission. We do not measure Balmer H\delta absorption in the host galaxy strong enough to be indicative of a rare, post-starburst "E+A" galaxy as reported by Arcavi et al. (2014). The light curve of PS1-10jh over a baseline of 3.5 yr is best modeled by fallback accretion of a tidally disrupted star. Its strong broad helium emission relative to hydrogen (He II \lambda 4686/H\alpha > 5) could be indicative of either the hydrogen-poor chemical composition of the disrupted star, or certain conditions in the tidal debris of a solar-composition star in the presence of an optically-thick, extended reprocessing envelope.

CfAIR2 is a large homogeneously reduced set of near-infrared (NIR) light curves for Type Ia supernovae (SN Ia) obtained with the 1.3m Peters Automated InfraRed Imaging TELescope (PAIRITEL). This data set includes 4607 measurements of 94 SN Ia and 4 additional SN Iax observed from 2005-2011 at the Fred Lawrence Whipple Observatory on Mount Hopkins, Arizona. CfAIR2 includes JHKs photometric measurements for 88 normal and 6 spectroscopically peculiar SN Ia in the nearby universe, with a median redshift of z~0.021 for the normal SN Ia. CfAIR2 data span the range from -13 days to +127 days from B-band maximum. More than half of the light curves begin before the time of maximum and the coverage typically contains ~13-18 epochs of observation, depending on the filter. We present extensive tests that verify the fidelity of the CfAIR2 data pipeline, including comparison to the excellent data of the Carnegie Supernova Project. CfAIR2 contributes to a firm local anchor for supernova cosmology studies in the NIR. Because SN Ia are more nearly standard candles in the NIR and are less vulnerable to the vexing problems of extinction by dust, CfAIR2 will help the supernova cosmology community develop more precise and accurate extragalactic distance probes to improve our knowledge of cosmological parameters, including dark energy and its potential time variation.

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.

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.

We present a novel method for the light-curve characterization of Pan-STARRS1 Medium Deep Survey (PS1 MDS) extragalactic sources into stochastic variables (SV) and burst-like (BL) transients, using multi-band image-differencing time-series data. We select detections in difference images associated with galaxy hosts using a star/galaxy catalog extracted from the deep PS1 MDS stacked images, and adopt a maximum a posteriori formulation to model their difference-flux time-series in four Pan-STARRS1 photometric bands g,r,i, and z. We use three deterministic light-curve models to fit burst-like transients and one stochastic light curve model, the Ornstein-Uhlenbeck process, in order to fit variability that is characteristic of active galactic nuclei (AGN). We assess the quality of fit of the models band-wise source-wise, using their estimated leave-out-one cross-validation likelihoods and corrected Akaike information criteria. We then apply a K-means clustering algorithm on these statistics, to determine the source classification in each band. The final source classification is derived as a combination of the individual filter classifications. We use our clustering method to characterize 4361 extragalactic image difference detected sources in the first 2.5 years of the PS1 MDS, into 1529 BL, and 2262 SV, with a purity of 95.00% for AGN, and 90.97% for SN based on our verification sets. We combine our light-curve classifications with their nuclear or off-nuclear host galaxy offsets, to define a robust photometric sample of 1233 active galactic nuclei and 812 supernovae. We use these samples to identify simple photometric priors that would enable their real-time identification in future wide-field synoptic surveys.

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.

The Arizona-NOAO Temporal Analysis and Response to Events System (ANTARES) is a joint project of the National Optical Astronomy Observatory and the Department of Computer Science at the University of Arizona. The goal is to build the software infrastructure necessary to process and filter alerts produced by time-domain surveys, with the ultimate source of such alerts being the Large Synoptic Survey Telescope (LSST). The ANTARES broker will add value to alerts by annotating them with information from external sources such as previous surveys from across the electromagnetic spectrum. In addition, the temporal history of annotated alerts will provide further annotation for analysis. These alerts will go through a cascade of filters to select interesting candidates. For the prototype, `interesting' is defined as the rarest or most unusual alert, but future systems will accommodate multiple filtering goals. The system is designed to be flexible, allowing users to access the stream at multiple points throughout the process, and to insert custom filters where necessary. We describe the basic architecture of ANTARES and the principles that will guide development and implementation.

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).

UV and optical photometry of Type Ia supernovae (SNe Ia) at low redshift have revealed the existence of two distinct color groups, NUV-red and NUV-blue events. The color curves differ primarily by an offset, with the NUV-blue u- color curves bluer than the NUV-red curves by 0.4 mag. For a sample of 23 low-z SNe~Ia observed with Swift, the NUV-red group dominates by a ratio of 2:1. We compare rest-frame UV/optical spectrophotometry of intermediate and high-z SNe Ia with UVOT photometry and HST spectrophotometry of low-z SNe Ia, finding that the same two color groups exist at higher-z, but with the NUV-blue events as the dominant group. Within each red/blue group, we do not detect any offset in color for different redshifts, providing insight into how SN~Ia UV emission evolves with redshift. Through spectral comparisons of SNe~Ia with similar peak widths and phase, we explore the wavelength range that produces the UV/OPT color differences. We show that the ejecta velocity of NUV-red SNe is larger than that of NUV-blue objects by roughly 12% on average. This velocity difference can explain some of the UV/optical color difference, but differences in the strengths of spectral features seen in meanspectra require additional explanation. Because of the different b-v colors for these groups, NUV-red SNe will have their extinction underestimated using common techniques. This, in turn, leads to under-estimation of the optical luminosity of the NUV-blue SNe~Ia, in particular, for the high-redshift cosmological sample. Not accounting for this effect should thus produce a distance bias that increases with redshift and could significantly bias measurements of cosmological parameters.

Type Iax supernovae (SNe Iax) are thermonuclear explosions that are related to SNe Ia, but are physically distinct. The most important differences are that SNe Iax have significantly lower luminosity (1% - 50% that of typical SNe Ia), lower ejecta mass (~0.1 - 0.5 M_sun), and may leave a bound remnant. The most extreme SN Iax is SN 2008ha, which peaked at M_V = -14.2 mag, about 5 mag below that of typical SNe Ia. Here, we present Hubble Space Telescope (HST) images of UGC 12682, the host galaxy of SN 2008ha, taken 4.1 years after the peak brightness of SN 2008ha. In these deep, high-resolution images, we detect a source coincident (0.86 HST pixels; 0.043"; 1.1 sigma) with the position of SN 2008ha with M_F814W = -5.4 mag. We determine that this source is unlikely to be a chance coincidence, but that scenario cannot be completely ruled out. If this source is directly related to SN 2008ha, it is either the luminous bound remnant of the progenitor white dwarf or its companion star. The source is consistent with being an evolved >3 M_sun initial mass star, and is significantly redder than the SN Iax 2012Z progenitor system, the first detected progenitor system for a thermonuclear SN. If this source is the companion star for SN 2008ha, there is a diversity in SN Iax progenitor systems, perhaps related to the diversity in SN Iax explosions. If the source is the bound remnant of the white dwarf, it must have expanded significantly. Regardless of the nature of this source, we constrain the progenitor system of SN 2008ha to have an age of <80 Myr.

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 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.

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 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 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 present multi-band optical photometry of 94 spectroscopically-confirmed Type Ia supernovae (SN Ia) in the redshift range 0.0055 to 0.073, obtained between 2006 and 2011. There are a total of 5522 light curve points. We show that our natural system SN photometry has a precision of roughly 0.03 mag or better in BVr'i', 0.06 mag in u', and 0.07 mag in U for points brighter than 17.5 mag and estimate that it has a systematic uncertainty of 0.014, 0.010, 0.012, 0.014, 0.046, and 0.073 mag in BVr'i'u'U, respectively. Comparisons of our standard system photometry with published SN Ia light curves and comparison stars reveal mean agreement across samples in the range of ~0.00-0.03 mag. We discuss the recent measurements of our telescope-plus-detector throughput by direct monochromatic illumination by Cramer et al (in prep.). This technique measures the whole optical path through the telescope, auxiliary optics, filters, and detector under the same conditions used to make SN measurements. Extremely well-characterized natural-system passbands (both in wavelength and over time) are crucial for the next generation of SN Ia photometry to reach the 0.01 mag accuracy level. The current sample of low-z SN Ia is now sufficiently large to remove most of the statistical sampling error from the dark energy error budget. But pursuing the dark-energy systematic errors by determining highly-accurate detector passbands, combining optical and near-infrared (NIR) photometry and spectra, using the nearby sample to illuminate the population properties of SN Ia, and measuring the local departures from the Hubble flow will benefit from larger, carefully measured nearby samples.

The flare of radiation from the tidal disruption and accretion of a star can be used as a marker for supermassive black holes that otherwise lie dormant and undetected in the centres of distant galaxies. Previous candidate flares have had declining light curves in good agreement with expectations, but with poor constraints on the time of disruption and the type of star disrupted, because the rising emission was not observed. Recently, two `relativistic' candidate tidal disruption events were discovered, each of whose extreme X-ray luminosity and synchrotron radio emission were interpreted as the onset of emission from a relativistic jet. Here we report the discovery of a luminous ultraviolet-optical flare from the nuclear region of an inactive galaxy at a redshift of 0.1696. The observed continuum is cooler than expected for a simple accreting debris disk, but the well-sampled rise and decline of its light curve follows the predicted mass accretion rate, and can be modelled to determine the time of disruption to an accuracy of two days. The black hole has a mass of about 2 million solar masses, modulo a factor dependent on the mass and radius of the star disrupted. On the basis of the spectroscopic signature of ionized helium from the unbound debris, we determine that the disrupted star was a helium-rich stellar core.

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.

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.

We have constructed a comprehensive statistical model for Type Ia supernova (SN Ia) light curves spanning optical through near infrared (NIR) data. A hierarchical framework coherently models multiple random and uncertain effects, including intrinsic supernova light curve covariances, dust extinction and reddening, and distances. An improved BayeSN MCMC code computes probabilistic inferences for the hierarchical model by sampling the global probability density of parameters describing individual supernovae and the population. We have applied this hierarchical model to optical and NIR data of 127 SN Ia from PAIRITEL, CfA3, CSP, and the literature. We find an apparent population correlation between the host galaxy extinction A_V and the the ratio of total-to-selective dust absorption R_V. For SN with low dust extinction, A_V < 0.4, we find R_V = 2.5 - 2.9, while at high extinctions, A_V > 1, low values of R_V < 2 are favored. The NIR luminosities are excellent standard candles and are less sensitive to dust extinction. They exhibit low correlation with optical peak luminosities, and thus provide independent information on distances. The combination of NIR and optical data constrains the dust extinction and improves the predictive precision of individual SN Ia distances by about 60%. Using cross-validation, we estimate an rms distance modulus prediction error of 0.11 mag for SN with optical and NIR data versus 0.15 mag for SN with optical data alone. Continued study of SN Ia in the NIR is important for improving their utility as precise and accurate cosmological distance indicators.

Small-scale magnetic features are present everywhere in the solar photosphere. Theoretical models, numerical calculations, and simulations describing the formation of these features have existed for a few decades, but there are only a few observational studies in direct support of the simulations. In this study we present the evolution of small-scale magnetic features with a spatial resolution close to 0.15 arcsecond and compare these observations with those predicted by numerical simulations and also with previous observational work of a similar nature. We analyze a 40 min time sequence of full Stokes spectropolarimetric 630.25 nm data from a plage region near the Sun center. We use line-of-sight velocities and magnetic field measurements obtained using Milne-Eddington inversion techniques with and without stray-light compensation along with measured continuum and line minimum intensities. We discuss the results in relation to earlier observations and simulations. We present eight cases involving strong downflows and magnetic field intensification. All cases studied are associated with the formation of a bright point in the continuum. In three out of the eight cases we find the presence of weak opposite polarity field in close proximity to the downflow. Our data are consistent with earlier simulations describing flux tube collapse, but the transition to a state with stronger field appears transient and short-lived, rather than resulting in a permanent field intensification. Three cases of weak opposite polarity field found adjacent to the downflows do not appear related to reconnection but may be related to overturning convection pulling down some field lines and leading to up/down "serpentine" field, as seen in some simulations.

The light echo systems of historical supernovae in the Milky Way and local group galaxies provide an unprecedented opportunity to reveal the effects of asymmetry on observables, particularly optical spectra. Scattering dust at different locations on the light echo ellipsoid witnesses the supernova from different perspectives and the light consequently scattered towards Earth preserves the shape of line profile variations introduced by asymmetries in the supernova photosphere. However, the interpretation of supernova light echo spectra to date has not involved a detailed consideration of the effects of outburst duration and geometrical scattering modifications due to finite scattering dust filament dimension, inclination, and image point-spread function and spectrograph slit width. In this paper, we explore the implications of these factors and present a framework for future resolved supernova light echo spectra interpretation, and test it against Cas A and SN 1987A light echo spectra. We conclude that the full modeling of the dimensions and orientation of the scattering dust using the observed light echoes at two or more epochs is critical for the correct interpretation of light echo spectra. Indeed, without doing so one might falsely conclude that differences exist when none are actually present.

We report the first detection of asymmetry in a supernova (SN) photosphere based on SN light echo (LE) spectra of Cas A from the different perspectives of dust concentrations on its LE ellipsoid. New LEs are reported based on difference images, and optical spectra of these LEs are analyzed and compared. After properly accounting for the effects of finite dust-filament extent and inclination, we find one field where the He I and H alpha features are blueshifted by an additional ~4000 km/s relative to other spectra and to the spectra of the Type IIb SN 1993J. That same direction does not show any shift relative to other Cas A LE spectra in the Ca II near-infrared triplet feature. We compare the perspectives of the Cas A LE dust concentrations with recent three-dimensional modeling of the SN remnant (SNR) and note that the location having the blueshifted He I and H alpha features is roughly in the direction of an Fe-rich outflow and in the opposite direction of the motion of the compact object at the center of the SNR. We conclude that Cas A was an intrinsically asymmetric SN. Future LE spectroscopy of this object, and of other historical SNe, will provide additional insight into the connection of explosion mechanism to SN to SNR, as well as give crucial observational evidence regarding how stars explode.

We report the discovery of a supernova (SN) with the highest apparent energy output to date and conclude that it represents an extreme example of the Type IIn subclass. The SN, which was discovered behind the Large Magellanic Cloud at z = 0.289 by the SuperMACHO microlensing survey, peaked at M_R = -21.5 mag and only declined by 2.9 mag over 4.7 years after the peak. Over this period, SN 2003ma had an integrated bolometric luminosity of 4 x 10^51 ergs, more than any other SN to date. The radiated energy is close to the limit allowed by conventional core-collapse explosions. Optical spectra reveal that SN 2003ma has persistent single-peaked intermediate-width hydrogen lines, a signature of interaction between the SN and a dense circumstellar medium. The light curves show further evidence for circumstellar interaction, including a long plateau with a shape very similar to the classic SN IIn 1988Z -- however, SN 2003ma is ten times more luminous at all epochs. The fast velocity measured for the intermediate-width H_alpha component (~6000 km/s) points towards an extremely energetic explosion (> 10^52 ergs), which imparts a faster blast-wave speed to the post-shock material and a higher luminosity from the interaction than is observed in typical SNe IIn. Mid-infrared observations of SN 2003ma suggest an infrared light echo is produced by normal interstellar dust at a distance ~0.5 pc from the SN.

SN2009ku, discovered by Pan-STARRS-1, is a Type Ia supernova (SNIa), and a member of the distinct SN2002cx-like class of SNeIa. Its light curves are similar to the prototypical SN2002cx, but are slightly broader and have a later rise to maximum in g. SN2009ku is brighter (~0.6 mag) than other SN2002cx-like objects, peaking at M_V = -18.4 mag - which is still significantly fainter than typical SNeIa. SN2009ku, which had an ejecta velocity of ~2000 kms^-1 at 18 days after maximum brightness is spectroscopically most similar to SN2008ha, which also had extremely low-velocity ejecta. However, SN2008ha had an exceedingly low luminosity, peaking at M_V = -14.2 mag, ~4 mag fainter than SN2009ku. The contrast of high luminosity and low ejecta velocity for SN2009ku is contrary to an emerging trend seen for the SN2002cx class. SN2009ku is a counter-example of a previously held belief that the class was more homogeneous than typical SNeIa, indicating that the class has a diverse progenitor population and/or complicated explosion physics. As the first example of a member of this class of objects from the new generation of transient surveys, SN2009ku is an indication of the potential for these surveys to find rare and interesting objects.

We present the early UV and optical light curve of Type IIP supernova (SN) 2010aq at z=0.0862, and compare it to analytical models for thermal emission following SN shock breakout in a red supergiant star. SN 2010aq was discovered in joint monitoring between the Galaxy Evolution Explorer (GALEX) Time Domain Survey (TDS) in the NUV and the Pan-STARRS1 Medium Deep Survey (PS1 MDS) in the g, r, i, and z bands. The GALEX and Pan-STARRS1 observations detect the SN less than 1 day after shock breakout, measure a diluted blackbody temperature of 31,000 +/- 6,000 K 1 day later, and follow the rise in the UV/optical light curve over the next 2 days caused by the expansion and cooling of the SN ejecta. The high signal-to-noise ratio of the simultaneous UV and optical photometry allows us to fit for a progenitor star radius of 700 +/- 200 R_sun, the size of a red supergiant star. An excess in UV emission two weeks after shock breakout compared to SNe well fitted by model atmosphere-code synthetic spectra with solar metallicity, is best explained by suppressed line blanketing due to a lower metallicity progenitor star in SN 2010aq. Continued monitoring of PS1 MDS fields by the GALEX TDS will increase the sample of early UV detections of Type II SNe by an order of magnitude, and probe the diversity of SN progenitor star properties.

In this work, we study and quantify properties of strong-field small-scale convection and compare observed properties with those predicted by numerical simulations. We analyze spectropolarimetric 630.25 nm data from a unipolar ephemeral region near sun center. We use line-of-sight velocities and magnetic field measurements obtained with Milne-Eddington inversion techniques along with measured continuum intensities and Stokes V amplitude asymmetry at a spatial resolution of 0.15 arcseconds to establish statistical relations between the measured quantities. We also study these properties for different types of distinct magnetic features, such as micropores, bright points, ribbons, flowers and strings. We present the first direct observations of a small-scale granular magneto-convection pattern within extended regions of strong (more than 600 G average) magnetic field. Along the boundaries of the flux concentrations we see mostly downflows and asymmetric Stokes V profiles, consistent with synthetic line profiles calculated from MHD simulations. We note the frequent occurrence of bright downflows along these boundaries. In the interior of the flux concentrations, we observe an up/down flow pattern that we identify as small-scale magnetoconvection, appearing similar to that of field-free granulation but with scales 4 times smaller. Measured RMS velocities are 70% of those of nearby field-free granulation, even though the average radiative flux is not reduced. The interiors of these flux concentrations are dominated by upflows.

We present photometric and spectroscopic observations of a luminous type IIP Supernova 2009kf discovered by the Pan-STARRS 1 (PS1) survey and detected also by GALEX. The SN shows a plateau in its optical and bolometric light curves, lasting approximately 70 days in the rest frame, with absolute magnitude of M_V = -18.4 mag. The P-Cygni profiles of hydrogen indicate expansion velocities of 9000km/s at 61 days after discovery which is extremely high for a type IIP SN. SN 2009kf is also remarkably bright in the near-ultraviolet (NUV) and shows a slow evolution 10-20 days after optical discovery. The NUV and optical luminosity at these epochs can be modelled with a black-body with a hot effective temperature (T ~16,000 K) and a large radius (R ~1x10^{15} cm). The bright bolometric and NUV luminosity, the lightcurve peak and plateau duration, the high velocities and temperatures suggest that 2009kf is a type IIP SN powered by a larger than normal explosion energy. Recently discovered high-z SNe (0.7 < z < 2.3) have been assumed to be IIn SNe, with the bright UV luminosities due to the interaction of SN ejecta with a dense circumstellar medium (CSM). UV bright SNe similar to SN 2009kf could also account for these high-z events, and its absolute magnitude M_NUV = -21.5 +/- 0.5 mag suggests such SNe could be discovered out to z ~2.5 in the PS1 survey.

We have used a precision calibrated photodiode as the fundamental metrology reference in order to determine the relative throughput of the PanSTARRS telescope and the Gigapixel imager, from 400 nm to 1050 nm. Our technique uses a tunable laser as a source of illumination on a transmissive flat-field screen. We determine the full-aperture system throughput as a function of wavelength, including (in a single integral measurement) the mirror reflectivity, the transmission functions of the filters and the corrector optics, and the detector quantum efficiency, by comparing the light seen by each pixel in the CCD array to that measured by a precision-calibrated silicon photodiode. This method allows us to determine the relative throughput of the entire system as a function of wavelength, for each pixel in the instrument, without observations of celestial standards. We present promising initial results from this characterization of the PanSTARRS system, and we use synthetic photometry to assess the photometric perturbations due to throughput variation across the field of view.

SN 2006bt displays characteristics unlike those of any other known Type Ia supernova (SN Ia). We present optical light curves and spectra of SN 2006bt which demonstrate the peculiar nature of this object. SN 2006bt has broad, slowly declining light curves indicative of a hot, high-luminosity SN, but lacks a prominent second maximum in the i band as do low-luminosity SNe Ia. Its spectra are similar to those of low-luminosity SNe Ia, containing features that are only present in cool SN photospheres. Light-curve fitting methods suggest that SN 2006bt is reddened by a significant amount of dust; however, it occurred in the outskirts of its early-type host galaxy and has no strong Na D absorption in any of its spectra, suggesting a negligible amount of host-galaxy dust absorption. C II is possibly detected in our pre-maximum spectra, but at a much lower velocity than other elements. The progenitor was likely very old, being a member of the halo population of a galaxy that shows no signs of recent star formation. SNe Ia have been very successfully modeled as a one-parameter family, and this is fundamental to their use as cosmological distance indicators. SN 2006bt is a challenge to that picture, yet its relatively normal light curves allowed SN 2006bt to be included in cosmological analyses. We generate mock SN Ia datasets which indicate that contamination by similar objects will both increase the scatter of a SN Ia Hubble diagram and systematically bias measurements of cosmological parameters. However, spectra and rest-frame i-band light curves should provide a definitive way to identify and eliminate such objects.

We discuss penumbral fine structure in a small part of a pore, observed with the CRISP imaging spectropolarimeter at the Swedish 1-m Solar Telescope (SST), close to its diffraction limit of 0.16 arcsec. Milne-Eddington inversions applied to these Stokes data reveal large variations of field strength and inclination angle over dark-cored penumbral intrusions and a dark-cored light bridge. The mid-outer part of this penumbra structure shows 0.3 arcsec wide spines, separated by 1.6 arcsec (1200 km) and associated with 30 deg inclination variations. Between these spines, there are no small-scale magnetic structures that easily can be be identified with individual flux tubes. A structure with nearly 10 deg more vertical and weaker magnetic field is seen midways between two spines. This structure is co-spatial with the brightest penumbral filament, possibly indicating the location of a convective upflow from below.

We report the discovery and orbit determination of 14 trans-Neptunian objects (TNOs) from the ESSENCE Supernova Survey difference imaging dataset. Two additional objects discovered in a similar search of the SDSS-II Supernova Survey database were recovered in this effort. ESSENCE repeatedly observed fields far from the Solar System ecliptic (-21 deg < beta < -5 deg), reaching limiting magnitudes per observation of I~23.1 and R~23.7. We examine several of the newly detected objects in detail, including 2003 UC_414 which orbits entirely between Uranus and Neptune and lies very close to a dynamical region that would make it stable for the lifetime of the Solar System. 2003 SS_422 and 2007 TA_418 have high eccentricities and large perihelia, making them candidate members of an outer class of trans-Neptunian objects. We also report a new member of the ''extended'' or ''detached'' scattered disk, 2004 VN_112, and verify the stability of its orbit using numerical simulations. This object would have been visible to ESSENCE for only ~2% of its orbit, suggesting a vast number of similar objects across the sky. We emphasize that off-ecliptic surveys are optimal for uncovering the diversity of such objects, which in turn will constrain the history of gravitational influences that shaped our early Solar System.

We present multiepoch spectra of 13 high-redshift Type Ia supernovae (SNe Ia) drawn from the literature, the ESSENCE and SNLS projects, and our own separate dedicated program on the ESO Very Large Telescope. We use the Supernova Identification (SNID) code of Blondin & Tonry to determine the spectral ages in the supernova rest frame. Comparison with the observed elapsed time yields an apparent aging rate consistent with the 1/(1+z) factor (where z is the redshift) expected in a homogeneous, isotropic, expanding universe. These measurements thus confirm the expansion hypothesis, while unambiguously excluding models that predict no time dilation, such as Zwicky's "tired light" hypothesis. We also test for power-law dependencies of the aging rate on redshift. The best-fit exponent for these models is consistent with the expected 1/(1+z) factor.

In this paper we advance the simple analytic photometric redshift estimator for Type Ia supernovae (SNe Ia) proposed by Wang (2007), and use it to study simulated SN Ia data. We find that better than 0.5% accuracy in z_phot (with RMS[(z_phot-z_spec)/(1+z_spec)]<0.005) is possible for SNe Ia with well sampled lightcurves in three observed passbands (riz) with a signal-to-noise ratio of 25 at peak brightness, if the extinction by dust is negligible. The corresponding bias in z_phot (the mean of (z_phot-z_spec)) is 5.4\times 10^{-4}. If dust extinction is taken into consideration in the riz observer-frame lightcurves, the accuracy in z_phot deteriorates to 4.4%, with a bias in z_phot of 8.0\times 10^{-3}. Adding the g band lightcurve improves the accuracy in z_phot to 2.5%, and reduces the bias in z_phot to -1.5\times 10^{-3}. Our results have significant implications for the design of future photometric surveys of SNe Ia from both ground and space telescopes. Accurate and precise photometric redshifts boost the cosmological utility of such surveys.

We describe the implementation and optimization of the ESSENCE supernova survey, which we have undertaken to measure the equation of state parameter of the dark energy. We present a method for optimizing the survey exposure times and cadence to maximize our sensitivity to the dark energy equation of state parameter w=P/rho c^2 for a given fixed amount of telescope time. For our survey on the CTIO 4m telescope, measuring the luminosity distances and redshifts for supernovae at modest redshifts (z~0.5 +- 0.2) is optimal for determining w. We describe the data analysis pipeline based on using reliable and robust image subtraction to find supernovae automatically and in near real-time. Since making cosmological inferences with supernovae relies crucially on accurate measurement of their brightnesses, we describe our efforts to establish a thorough calibration of the CTIO 4m natural photometric system. In its first four years, ESSENCE has discovered and spectroscopically confirmed 102 type Ia SNe, at redshifts from 0.10 to 0.78, identified through an impartial, effective methodology for spectroscopic classification and redshift determination. We present the resulting light curves for the all type Ia supernovae found by ESSENCE and used in our measurement of w, presented in Wood-Vasey et al, 2007.

We present constraints on the dark energy equation-of-state parameter, w=P/(rho c^2), using 60 Type Ia supernovae (SNe Ia) from the ESSENCE supernova survey. We derive a set of constraints on the nature of the dark energy assuming a flat Universe. By including constraints on (Omega_M, w) from baryon acoustic oscillations, we obtain a value for a static equation-of-state parameter w=-1.05^{+0.13}_{-0.12} (stat; 1 sigma) +- 0.11 (sys) and Omega_M=0.274^{+0.033}_{-0.020} (stat; 1 sigma) with a best-fit chi^2/DoF of 0.96. These results are consistent with those reported by the SuperNova Legacy Survey in a similar program measuring supernova distances and redshifts. We evaluate sources of systematic error that afflict supernova observations and present Monte Carlo simulations that explore these effects. Currently, the largest systematic currently with the potential to affect our measurements is the treatment of extinction due to dust in the supernova host galaxies. Combining our set of ESSENCE SNe Ia with the SuperNova Legacy Survey SNe Ia, we obtain a joint constraint of w=-1.07^{+0.09}_{-0.09} (stat; 1 sigma) +- 0.12 (sys), Omega_M=0.267^{+0.028}_{-0.018} (stat; 1 sigma) with a best-fit chi^2/DoF of 0.91. The current SN Ia data are fully consistent with a cosmological constant.