• ### The Pan-STARRS1 Surveys(1612.05560)

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.
• ### The Pan-STARRS1 Database and Data Products(1612.05243)

Jan. 29, 2019 astro-ph.IM
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-STARRS Pixel Processing: Detrending, Warping, Stacking(1612.05245)

Jan. 28, 2019 astro-ph.IM
The Pan-STARRS1 Science Consortium have carried out a set of imaging surveys using the 1.4 giga-pixel GPC1 camera on the PS1 telescope. As this camera is composed of many individual electronic readouts, and covers a very large field of view, great care was taken to ensure that the many instrumental effects were corrected to produce the most uniform detector response possible. We present the image detrending steps used as part of the processing of the data contained within the public release of the Pan-STARRS1 Data Release 1 (DR1). In addition to the single image processing, the methods used to transform the 375,573 individual exposures into a common sky-oriented grid are discussed, as well as those used to produce both the image stack and difference combination products.
• ### The lowest metallicity type II supernova from the highest mass red-supergiant progenitor(1805.04434)

May 11, 2018 astro-ph.SR, astro-ph.HE
Red supergiants have been confirmed as the progenitor stars of the majority of hydrogen-rich type II supernovae. However, while such stars are observed with masses >25M$_\odot$, detections of >18M$_\odot$ progenitors remain elusive. Red supergiants are also expected to form at all metallicities, but discoveries of explosions from low-metallicity progenitors are scarce. Here, we report observations of the type II supernova, SN 2015bs, for which we infer a progenitor metallicity of $\leq$0.1Z$_\odot$ from comparison to photospheric-phase spectral models, and a Zero Age Main-Sequence mass of 17-25M$_\odot$ through comparison to nebular-phase spectral models. SN 2015bs displays a normal 'plateau' light-curve morphology, and typical spectral properties, implying a red supergiant progenitor. This is the first example of such a high mass progenitor for a 'normal' type II supernova, suggesting a link between high mass red supergiant explosions and low-metallicity progenitors.
• ### The Pan-STARRS1 Proper-motion Survey for Young Brown Dwarfs in Nearby Star-forming Regions. I. Taurus Discoveries and a Reddening-free Classification Method for Ultracool Dwarfs(1804.01533)

We are conducting a proper-motion survey for young brown dwarfs in the Taurus-Auriga molecular cloud based on the Pan-STARRS1 3$\pi$ Survey. Our search uses multi-band photometry and astrometry to select candidates, and is wider (370 deg$^{2}$) and deeper (down to $\approx$3 M$_{\rm Jup}$) than previous searches. We present here our search methods and spectroscopic follow-up of our high-priority candidates. Since extinction complicates spectral classification, we have developed a new approach using low-resolution ($R \approx 100$) near-infrared spectra to quantify reddening-free spectral types, extinctions, and gravity classifications for mid-M to late-L ultracool dwarfs ($\approx 100-3$ M$_{\rm Jup}$ in Taurus). We have discovered 25 low-gravity (VL-G) and the first 11 intermediate-gravity (INT-G) substellar (M6-L1) members of Taurus, constituting the largest single increase of Taurus brown dwarfs to date. We have also discovered 1 new Pleiades member and 13 new members of the Perseus OB2 association, including a candidate very wide separation (58 kAU) binary. We homogeneously reclassify the spectral types and extinctions of all previously known Taurus brown dwarfs. Altogether our discoveries have thus far increased the substellar census in Taurus by $\approx 40\%$ and added three more L-type members ($\approx 5-10$ M$_{\rm Jup}$). Most notably, our discoveries reveal an older ($>$10 Myr) low-mass population in Taurus, in accord with recent studies of the higher-mass stellar members. The mass function appears to differ between the younger and older Taurus populations, possibly due to incompleteness of the older stellar members or different star formation processes.
• ### The Complete Light-curve Sample of Spectroscopically Confirmed Type Ia Supernovae from Pan-STARRS1 and Cosmological Constraints from The Combined Pantheon Sample(1710.00845)

March 26, 2018 astro-ph.CO
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.
• ### Measuring Dark Energy Properties with Photometrically Classified Pan-STARRS Supernovae. II. Cosmological Parameters(1710.00846)

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

Jan. 12, 2018 astro-ph.HE
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.
• ### Photometry and Proper Motions of M, L, and T Dwarfs from the Pan-STARRS1 3$\pi$ Survey(1701.00490)

Nov. 9, 2017 astro-ph.SR
We present a catalog of 9888 M, L and T dwarfs detected in the Pan-STARRS1 3$\pi$ Survey (PS1), covering three-quarters of the sky. Our catalog contains nearly all known objects of spectral types L0-T2 in the PS1 field, with objects as early as M0 and as late as T9, and includes PS1, 2MASS, AllWISE, and Gaia DR1 photometry. We analyze the different types of photometry reported by PS1, and use two types in our catalog to maximize both depth and accuracy. Using parallaxes from the literature, we construct empirical SEDs for field ultracool dwarfs spanning 0.5-12 $\mu$m. We determine typical colors of M0-T9 dwarfs, and we highlight the distinctive colors of subdwarfs and young objects. Our catalog includes 492 L dwarfs detected in $r_{\rm P1}$, the largest sample of L dwarfs detected at such blue wavelengths. We combine astrometry from PS1, 2MASS, and Gaia DR1 to calculate new proper motions for our catalog. We achieve a median precision of 2.9 mas yr$^{-1}$, a factor of $\approx$3-10 improvement over previous large catalogs. Our catalog contains proper motions for 2405 M6-T9 dwarfs and includes the largest set of homogeneous proper motions for L and T dwarfs published to date, 406 objects for which there were no previous measurements, and 1176 objects for which we improve upon previous literature values. We analyze the kinematics of ultracool dwarfs in our catalog and find evidence that bluer but otherwise generic late-M and L field dwarfs (i.e., not subdwarfs) tend to have higher tangential velocities compared to typical field objects. With the public release of the PS1 data, this survey will continue to be an essential tool for characterizing the ultracool dwarf population.
• ### The Foundation Supernova Survey: Motivation, Design, Implementation, and First Data Release(1711.02474)

Nov. 6, 2017 astro-ph.CO, astro-ph.HE
The Foundation Supernova Survey aims to provide a large, high-fidelity, homogeneous, and precisely-calibrated low-redshift Type Ia supernova (SN Ia) sample for cosmology. The calibration of the current low-redshift SN sample is the largest component of systematic uncertainties for SN cosmology, and new data are necessary to make progress. We present the motivation, survey design, observation strategy, implementation, and first results for the Foundation Supernova Survey. We are using the Pan-STARRS telescope to obtain photometry for up to 800 SNe Ia at z < 0.1. This strategy has several unique advantages: (1) the Pan-STARRS system is a superbly calibrated telescopic system, (2) Pan-STARRS has observed 3/4 of the sky in grizy making future template observations unnecessary, (3) we have a well-tested data-reduction pipeline, and (4) we have observed ~3000 high-redshift SNe Ia on this system. Here we present our initial sample of 225 SN Ia griz light curves, of which 180 pass all criteria for inclusion in a cosmological sample. The Foundation Supernova Survey already contains more cosmologically useful SNe Ia than all other published low-redshift SN Ia samples combined. We expect that the systematic uncertainties for the Foundation Supernova Sample will be 2-3 times smaller than other low-redshift samples. We find that our cosmologically useful sample has an intrinsic scatter of 0.111 mag, smaller than other low-redshift samples. We perform detailed simulations showing that simply replacing the current low-redshift SN Ia sample with an equally sized Foundation sample will improve the precision on the dark energy equation-of-state parameter by 35%, and the dark energy figure-of-merit by 72%.
• ### A kilonova as the electromagnetic counterpart to a gravitational-wave source(1710.05841)

Oct. 17, 2017 astro-ph.HE
Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements.
• ### Physical properties of 15 quasars at $z\gtrsim 6.5$(1710.01251)

Oct. 3, 2017 astro-ph.GA
Quasars are galaxies hosting accreting supermassive black holes; due to their brightness, they are unique probes of the early universe. To date, only few quasars have been reported at $z > 6.5$ ($<$800 Myr after the Big Bang). In this work, we present six additional $z \gtrsim 6.5$ quasars discovered using the Pan-STARRS1 survey. We use a sample of 15 $z \gtrsim 6.5$ quasars to perform a homogeneous and comprehensive analysis of this highest-redshift quasar population. We report four main results: (1) the majority of $z\gtrsim$6.5 quasars show large blueshifts of the broad CIV 1549\AA$\,$emission line compared to the systemic redshift of the quasars, with a median value $\sim$3$\times$ higher than a quasar sample at $z\sim$1; (2) we estimate the quasars' black hole masses (M$\rm_{BH}\sim$0.3$-$5 $\times$ 10$^{9}$ M$_{\odot}$) via modeling of the MgII 2798\AA$\,$emission line and rest-frame UV continuum; we find that quasars at high redshift accrete their material (with $\langle (L_{\mathrm{bol}}/L_{\mathrm{Edd}}) \rangle = 0.39$) at a rate comparable to a luminosity-matched sample at lower$-$redshift, albeit with significant scatter ($0.4$ dex); (3) we recover no evolution of the FeII/MgII abundance ratio with cosmic time; (4) we derive near zone sizes; together with measurements for $z\sim6$ quasars from recent work, we confirm a shallow evolution of the decreasing quasar near zone sizes with redshift. Finally, we present new millimeter observations of the [CII] 158 $\mu$m emission line and underlying dust continuum from NOEMA for four quasars, and provide new accurate redshifts and [CII]/infrared luminosities estimates. The analysis presented here shows the large range of properties of the most distant quasars.
• ### A transient search using combined human and machine classifications(1707.05223)

July 17, 2017 astro-ph.IM
Large modern surveys require efficient review of data in order to find transient sources such as supernovae, and to distinguish such sources from artefacts and noise. Much effort has been put into the development of automatic algorithms, but surveys still rely on human review of targets. This paper presents an integrated system for the identification of supernovae in data from Pan-STARRS1, combining classifications from volunteers participating in a citizen science project with those from a convolutional neural network. The unique aspect of this work is the deployment, in combination, of both human and machine classifications for near real-time discovery in an astronomical project. We show that the combination of the two methods outperforms either one used individually. This result has important implications for the future development of transient searches, especially in the era of LSST and other large-throughput surveys.
• ### Supernovae 2016bdu and 2005gl, and their link with SN 2009ip-like transients: another piece of the puzzle(1707.00611)

July 3, 2017 astro-ph.SR, astro-ph.HE
Supernova (SN) 2016bdu is an unusual transient resembling SN 2009ip. SN 2009ip-like events are characterized by a long-lasting phase of erratic variability which ends with two luminous outbursts a few weeks apart. The second outburst is significantly more luminous (about 3 mag) than the first. In the case of SN 2016bdu, the first outburst (Event A) reached an absolute magnitude M(r) ~ -15.3 mag, while the second one (Event B) occurred over one month later and reached M(r) ~ -18 mag. By inspecting archival data, a faint source at the position of SN 2016bdu is detectable several times in the past few years. We interpret these detections as signatures of a phase of erratic variability, similar to that experienced by SN 2009ip between 2008 and mid-2012, and resembling the currently observed variability of the luminous blue variable SN 2000ch in NGC 3432. Spectroscopic monitoring of SN 2016bdu during the second peak initially shows features typical of a SN IIn. One month after the Event B maximum, the spectra develop broad Balmer lines with P Cygni profiles and broad metal features. At these late phases, the spectra resemble those of a typical Type II SN. All members of this SN 2009ip-like group are remarkably similar to the Type IIn SN 2005gl. For this object, the claim of a terminal SN explosion is supported by the disappearance of the progenitor star. The similarity with SN 2005gl suggests that all members of this family may finally explode as genuine SNe, although the unequivocal detection of nucleosynthesised elements in their nebular spectra is still missing.
• ### Measuring the Properties of Dark Energy with Photometrically Classified Pan-STARRS Supernovae. I. Systematic Uncertainty from Core-Collapse Supernova Contamination(1611.07042)

June 19, 2017 astro-ph.CO
The Pan-STARRS (PS1) Medium Deep Survey discovered over 5,000 likely supernovae (SNe) but obtained spectral classifications for just 10% of its SN candidates. We measured spectroscopic host galaxy redshifts for 3,147 of these likely SNe and estimate that $\sim$1,000 are Type Ia SNe (SNe Ia) with light-curve quality sufficient for a cosmological analysis. We use these data with simulations to determine the impact of core-collapse SN (CC SN) contamination on measurements of the dark energy equation of state parameter, $w$. Using the method of Bayesian Estimation Applied to Multiple Species (BEAMS), distances to SNe Ia and the contaminating CC SN distribution are simultaneously determined. We test light-curve based SN classification priors for BEAMS as well as a new classification method that relies upon host galaxy spectra and the association of SN type with host type. By testing several SN classification methods and CC SN parameterizations on large SN simulations, we estimate that CC SN contamination gives a systematic error on $w$ ($\sigma_w^{CC}$) of 0.014, 29% of the statistical uncertainty. Our best method gives $\sigma_w^{CC} = 0.004$, just 8% of the statistical uncertainty, but could be affected by incomplete knowledge of the CC SN distribution. This method determines the SALT2 color and shape coefficients, $\alpha$ and $\beta$, with $\sim$3% bias. However, we find that some variants require $\alpha$ and $\beta$ to be fixed to known values for BEAMS to yield accurate measurements of $w$. Finally, the inferred abundance of bright CC SNe in our sample is greater than expected based on measured CC SN rates and luminosity functions.
• ### Observations of the GRB afterglow ATLAS17aeu and its possible association with GW170104(1706.00175)

June 1, 2017 astro-ph.HE
We report the discovery and multi-wavelength data analysis of the peculiar optical transient, ATLAS17aeu. This transient was identified in the skymap of the LIGO gravitational wave event GW170104 by our ATLAS and Pan-STARRS coverage. ATLAS17aeu was discovered 23.1hrs after GW170104 and rapidly faded over the next 3 nights, with a spectrum revealing a blue featureless continuum. The transient was also detected as a fading x-ray source by Swift and in the radio at 6 and 15\,GHz. A gamma ray burst GRB170105 was detected by 3 satellites 19.04\,hrs after GW170104 and 4.10\,hrs before our first optical detection. We analyse the multi-wavelength fluxes in the context of the known GRB population and discuss the observed sky rates of GRBs and their afterglows. We find it statistically likely that ATLAS17aeu is an afterglow associated with GRB170105, with a chance coincidence ruled out at the 99\% confidence or 2.6$\sigma$. A long, soft GRB within a redshift range of $1 \lesssim z \lesssim 2.9$ would be consistent with all the observed multi-wavelength data. The Poisson probability of a chance occurrence of GW170104 and ATLAS17aeu is $p=0.04$. This is the probability of a chance coincidence in 2D sky location and in time. These observations indicate that ATLAS17aeu is plausibly a normal GRB afterglow at significantly higher redshift than the distance constraint for GW170104 and therefore a chance coincidence. However if a redshift of the faint host were to place it within the GW170104 distance range, then physical association with GW170104 should be considered.
• ### A Search for L/T Transition Dwarfs With Pan-STARRS1 and WISE. III. Young L Dwarf Discoveries and Proper Motion Catalogs in Taurus and Scorpius-Centaurus(1702.00789)

Feb. 8, 2017 astro-ph.SR
We present the discovery of eight young M7-L2 dwarfs in the Taurus star-forming region and the Scorpius-Centaurus OB Association, serendipitously found during a wide-field search for L/T transition dwarfs using Pan-STARRS1 (optical) and WISE (mid-infrared) photometry. We identify PSO J060.3200+25.9644 (near-infrared spectral type L1) and PSO J077.1033+24.3809 (L2) as new members of Taurus based on their VL-G gravity classifications, the consistency of their photometry and proper motions with previously known Taurus objects, and the low probability of contamination by field objects. PSO J077.1033+24.3809 is the coolest substellar member of Taurus found to date. Both Taurus objects are among the lowest mass free-floating objects ever discovered, with estimated masses $\approx$6 M$_{\rm Jup}$, and provide further evidence that isolated planetary-mass objects can form as part of normal star-formation processes. PSO J060.3200+25.9644 (a.k.a. DANCe J040116.80+255752.2) was previously identified as a likely member of the Pleiades (age $\approx125$ Myr) based on photometry and astrometry, but its VL-G gravity classification and near-infrared photometry imply a much younger age and thus point to Taurus membership. We have also discovered six M7-L1 dwarfs in outlying regions of Scorpius-Centaurus with photometry, proper motions, and low-gravity spectral signatures consistent with membership. These objects have estimated masses $\approx$15-36 M$_{\rm Jup}$. The M7 dwarf, PSO J237.1470-23.1489, shows excess mid-infrared flux implying the presence of a circumstellar disk. Finally, we present catalogs of Pan-STARRS1 proper motions for low-mass members of Taurus and Upper Scorpius with median precisions of $\approx$3 mas yr$^{-1}$, including 67 objects with no previous proper motion and 359 measurements that improve on literature values.
• ### A Search for L/T Transition Dwarfs With Pan-STARRS1 and WISE. II. L/T Transition Atmospheres and Young Discoveries(1612.02824)

Dec. 8, 2016 astro-ph.SR
The evolution of brown dwarfs from L to T spectral types is one of the least understood aspects of the ultracool population, partly for lack of a large, well-defined, and well-characterized sample in the L/T transition. To improve the existing census, we have searched $\approx$28,000 deg$^2$ using the Pan-STARRS1 and WISE surveys for L/T transition dwarfs within 25 pc. We present 130 ultracool dwarf discoveries with estimated distances $\approx9-130$ pc, including 21 that were independently discovered by other authors and 3 that were previously identified as photometric candidates. Seventy-nine of our objects have near-IR spectral types of L6-T4.5, the most L/T transition dwarfs from any search to date, and we have increased the census of L9-T1.5 objects within 25 pc by over 50%. The color distribution of our discoveries provides further evidence for the "L/T gap," a deficit of objects with $(J-K)_{\rm MKO}\approx0.0-0.5$ mag in the L/T transition, and thus reinforces the idea that the transition from cloudy to clear photospheres occurs rapidly. Among our discoveries are 31 candidate binaries based on their low-resolution spectral features. Two of these candidates are common proper motion companions to nearby main sequence stars; if confirmed as binaries, these would be rare benchmark systems with the potential to stringently test ultracool evolutionary models. Our search also serendipitously identified 23 late-M and L dwarfs with spectroscopic signs of low gravity implying youth. Finally, we identify 10 candidate members of nearby young moving groups (YMG) with spectral types L7-T4.5, including three showing spectroscopic signs of low gravity. If confirmed, any of these would be among the coolest known YMG members and would help to determine the effective temperature at which young brown dwarfs cross the L/T transition. (Abridged)
• ### PS1-14bj: A Hydrogen-Poor Superluminous Supernova With a Long Rise and Slow Decay(1605.05235)

Oct. 10, 2016 astro-ph.HE
We present photometry and spectroscopy of PS1-14bj, a hydrogen-poor superluminous supernova (SLSN) at redshift $z=0.5215$ discovered in the last months of the Pan-STARRS1 Medium Deep Survey. PS1-14bj stands out by its extremely slow evolution, with an observed rise of $\gtrsim 125$ rest-frame days, and exponential decline out to $\sim 250$ days past peak at a measured rate of $0.01~{\rm mag~day}^{-1}$, consistent with fully-trapped $^{56}$Co decay. This is the longest rise time measured in a SLSN to date, and the first SLSN to show a rise time consistent with pair-instability supernova (PISN) models. Compared to other slowly-evolving SLSNe, it is spectroscopically similar to the prototype SN2007bi at maximum light, though lower in luminosity ($L_{\rm peak} \simeq 4.6 \times 10^{43} {\rm erg s}^{-1}$) and with a flatter peak than previous events. PS1-14bj shows a number of peculiar properties, including a near-constant color temperature for $>200$ days past peak, and strong emission lines from [O III] $\lambda$5007 and [O III] $\lambda$4363 with a velocity width of $\sim$3400 km/s, in its late-time spectra. These both suggest there is a sustained source of heating over very long timescales, and are incompatible with a simple $^{56}$Ni-powered/PISN interpretation. A modified magnetar model including emission leakage at late times can reproduce the light curve, in which case the blue continuum and [O III] features are interpreted as material heated and ionized by the inner pulsar wind nebula becoming visible at late times. Alternatively, the late-time heating could be due to interaction with a shell of H-poor circumstellar material.
• ### Dust in three dimensions in the Galactic Plane(1609.03177)

Sept. 17, 2016 astro-ph.GA
We present three dimensional maps in monochromatic extinction $A_{\rm 0}$ and the extinction parameter $R_0$ within a few degrees of the Galactic plane. These are inferred using photometry from the Pan-STARRS1 and Spitzer Glimpse surveys of nearly $20$ million stars located in the region $l = 0-250^{\circ}$ and from $b = -4.5^{\circ}$ to $b=4.5^{\circ}$. Given the available stellar number density, we use an angular resolution of $7' \times 7'$ and steps of $1{\rm mag}$ in distance modulus. We simultaneously estimate distance modulus and effective temperature $T_{\rm eff}$ alongside the other parameters for stars individually using the method of \citet{Hanson2014} before combining these estimates to a complete map. The full maps are available via the MNRAS website.
• ### Superluminous supernova 2015bn in the nebular phase: evidence for the engine-powered explosion of a stripped massive star(1608.02995)

Sept. 8, 2016 astro-ph.SR, astro-ph.HE
We present nebular-phase imaging and spectroscopy for the hydrogen-poor superluminous supernova SN 2015bn, at redshift z=0.1136, spanning +250-400 d after maximum light. The light curve exhibits a steepening in the decline rate from 1.4 mag/(100 d) to 1.7 mag/(100 d), suggestive of a significant decrease in the opacity. This change is accompanied by a transition from a blue continuum superposed with photospheric absorption lines to a nebular spectrum dominated by emission lines of oxygen, calcium and magnesium. There are no obvious signatures of circumstellar interaction or large nickel mass. We show that the spectrum at +400 d is virtually identical to a number of energetic Type Ic supernovae such as SN 1997dq, SN 2012au, and SN 1998bw, indicating similar core conditions and strengthening the link between hypernovae'/long gamma-ray bursts and superluminous supernovae. A single explosion mechanism may unify these events that span absolute magnitudes of -22 < M_B < -17. Both the light curve and spectrum of SN 2015bn are consistent with an engine-driven explosion ejecting 7-30 M$_\odot$ of oxygen-dominated ejecta (for reasonable choices in temperature and opacity). A strong and relatively narrow O II $\lambda$7774 line, seen in a number of these energetic events but not in normal supernovae, may point to an inner shell that is the signature of a central engine.
• ### A search for an optical counterpart to the gravitational wave event GW151226(1606.04795)

Aug. 8, 2016 astro-ph.SR, astro-ph.HE
We present a search for an electromagnetic counterpart of the gravitational wave source GW151226. Using the Pan-STARRS1 telescope we mapped out 290 square degrees in the optical i_ps filter starting 11.5hr after the LIGO information release and lasting for a further 28 days. The first observations started 49.5hr after the time of the GW151226 detection. We typically reached sensitivity limits of i_ps = 20.3-20.8 and covered 26.5% of the LIGO probability skymap. We supplemented this with ATLAS survey data, reaching 31% of the probability region to shallower depths of m~19. We found 49 extragalactic transients (that are not obviously AGN), including a faint transient in a galaxy at 7Mpc (a luminous blue variable outburst) plus a rapidly decaying M-dwarf flare. Spectral classification of 20 other transient events showed them all to be supernovae. We found an unusual transient, PS15dpn, with an explosion date temporally coincident with GW151226 which evolved into a type Ibn supernova. The redshift of the transient is secure at z=0.1747 +/- 0.0001 and we find it unlikely to be linked, since the luminosity distance has a negligible probability of being consistent with that of GW151226. In the 290 square degrees surveyed we therefore do not find a likely counterpart. However we show that our survey strategy would be sensitive to NS-NS mergers producing kilonovae at D < 100 Mpc which is promising for future LIGO/Virgo searches.
• ### Pan-STARRS and PESSTO search for an optical counterpart to the LIGO gravitational wave source GW150914(1602.04156)

July 28, 2016 astro-ph.SR, astro-ph.HE
We searched for an optical counterpart to the first gravitational wave source discovered by LIGO (GW150914), using a combination of the Pan-STARRS1 wide-field telescope and the PESSTO spectroscopic follow-up programme. As the final LIGO sky maps changed during analysis, the total probability of the source being spatially coincident with our fields was finally only 4.2 per cent. Therefore we discuss our results primarily as a demonstration of the survey capability of Pan-STARRS and spectroscopic capability of PESSTO. We mapped out 442 square degrees of the northern sky region of the initial map. We discovered 56 astrophysical transients over a period of 41 days from the discovery of the source. Of these, 19 were spectroscopically classified and a further 13 have host galaxy redshifts. All transients appear to be fairly normal supernovae and AGN variability and none is obviously linked with GW150914. We illustrate the sensitivity of our survey by defining parameterised lightcurves with timescales of 4, 20 and 40 days and use the sensitivity of the Pan-STARRS1 images to set limits on the luminosities of possible sources. The Pan-STARRS1 images reach limiting magnitudes of i = 19.2, 20.0 and 20.8 respectively for the three timescales. For long timescale parameterised lightcurves (with FWHM=~40d) we set upper limits of M_i <= -17.2 -0.9/+1.4 if the distance to GW150914 is D = 400 +/- 200Mpc. The number of type Ia SN we find in the survey is similar to that expected from the cosmic SN rate, indicating a reasonably complete efficiency in recovering supernova like transients out to D = 400 +/- 200 Mpc.
• ### SN 2015bn: a detailed multi-wavelength view of a nearby superluminous supernova(1603.04748)

We present observations of SN 2015bn (= PS15ae = CSS141223-113342+004332 = MLS150211-113342+004333), a Type I superluminous supernova (SLSN) at redshift $z=0.1136$. As well as being one of the closest SLSNe I yet discovered, it is intrinsically brighter ($M_U\approx-23.1$) and in a fainter galaxy ($M_B\approx-16.0$) than other SLSNe at $z\sim0.1$. We used this opportunity to collect the most extensive dataset for any SLSN I to date, including densely-sampled spectroscopy and photometry, from the UV to the NIR, spanning $-$50 to +250 days from optical maximum. SN 2015bn fades slowly, but exhibits surprising undulations in the light curve on a timescale of 30-50 days, especially in the UV. The spectrum shows extraordinarily slow evolution except for a rapid transformation between +7 and +20-30 days. No narrow emission lines from slow-moving material are observed at any phase. We derive physical properties including the bolometric luminosity, and find slow velocity evolution and non-monotonic temperature and radial evolution. A deep radio limit rules out a healthy off-axis gamma-ray burst, and places constraints on the pre-explosion mass loss. The data can be consistently explained by a $\gtrsim10\,{\rm M}_\odot$ stripped progenitor exploding with $\sim 10^{51}\,$erg kinetic energy, forming a magnetar with a spin-down timescale of $\sim20$ days (thus avoiding a gamma-ray burst) that reheats the ejecta and drives ionization fronts. The most likely alternative scenario -- interaction with $\sim20\,{\rm M}_\odot$ of dense, inhomogeneous circumstellar material -- can be tested with continuing radio follow-up.
• This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.