We present the result of a search for orphan Gamma-Ray Burst (GRB) afterglows in the Pan-STARRS1 (PS1) COSMOS survey. There is extensive theoretical and observational evidence suggesting that GRBs are collimated jets; the direct observation of orphan GRB afterglows would further support this model. An optimal survey strategy is designed by coupling the PS1 survey with the Subaru/Hyper-Suprime-Cam (HSC) survey. The PS1 COSMOS survey, one of the survey fields in the PS1 Medium Deep Survey (PS1/MDS), searches a field of 7 deg$^2$ from December 2011 to January 2014, reaching a limiting magnitude R $\sim$ 23. The dense cadence of PS1/MDS is crucial for identifying transients, and the deep magnitude reached by the HSC survey (R $\sim$ 26) is important for evaluating potential GRB hosts. A transient classification method is employed to select potential orphan GRB afterglow candidates. After a thorough analysis of the transient and host galaxy properties, we conclude that there are no candidates in this survey field. The null result implies that the consideration of jet structures is essential for further orphan GRB afterglow surveys.

The intra-night trajectories of asteroids can be approximated by straight lines, and so are their intra-night detections. Therefore, the Hough transform, a line detecting algorithm, can be used to connect the line-up detections to find asteroids. We applied this algorithm to a high-cadence Pan-STARRS 1 (PS1) observation, which was originally designed to collect asteroid light curves for rotation period measurements (Chang et al., 2019). The algorithm recovered most of the known asteroids in the observing fields and, moreover, discovered 3574 new asteroids with magnitude mainly of 21.5 < w_p1 < 22.5 mag. This magnitude range is equivalent to sub-kilometer main-belt asteroids (MBAs), which usually lack of rotation period measurements due to their faintness. Using the light curves of the 3574 new asteroids, we obtained 122 reliable rotation periods, of which 13 are super-fast rotators (SRFs; i.e., rotation period of < 2 hr). The required cohesion to survive these SFRs range from tens to thousands of Pa, a value consistent with the known SFRs and the regolith on the Moon and Mars. The higher chance of discovering SFRs here suggests that sub-kilometer MBAs probably harbor more SFRs.

We present a systematic search for periodically varying quasar and supermassive black hole binary (SMBHB) candidates in the Pan-STARRS1 Medium Deep Survey. From $\sim9,000$ color-selected quasars in a $\sim50$ deg$^{2}$ sky area, we initially identify $26$ candidates with more than $1.5$ cycles of variation. We extend the baseline of observations via our imaging campaign with the Discovery Channel Telescope and the Las Cumbres Observatory network and reevaluate the candidates using a more rigorous, maximum likelihood method. Using a range of statistical criteria and assuming the Damped Random Walk model for normal quasar variability, we identify one statistically significant periodic candidate. We also investigate the capabilities of detecting SMBHBs by the Large Synoptic Survey Telescope using our study with MDS as a benchmark and explore any complementary, multiwavelength evidence for SMBHBs in our sample.

We present the results of the monitoring campaign of the Type IIn supernova (SN) 2018cnf (aka ASASSN-18mr). It was discovered about 10 days before the maximum light (on MJD = 58293.4+-5.7 in the V band, with MV = -18.13+-0.15 mag). The multiband light curves show an immediate post-peak decline with some minor luminosity fluctuations, followed by a flattening starting about 40 days after maximum. The early spectra are relatively blue and show narrow Balmer lines with P Cygni profiles. Additionally, Fe II, O I, He I and Ca II are detected. The spectra show little evolution with time, with intermediate-width features becoming progressively more prominent, indicating stronger interaction of the SN ejecta with the circumstellar medium. The inspection of archival images from the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) survey has revealed a variable source at the SN position, with a brightest detection in December 2015 at Mr = -14.66+-0.17 mag. This was likely an eruptive phase from the massive progenitor star started from at least mid-2011, and that produced the circumstellar environment within which the star exploded as a Type IIn SN. The overall properties of SN 2018cnf closely resemble those of transients such as SN 2009ip. This similarity favours a massive hypergiant, perhaps a luminous blue variable, as progenitor for SN 2018cnf.

We report on the catastrophic disintegration of P/2016 G1 (PANSTARRS), an active asteroid, in April 2016. Deep images over three months show that object is constituted by a central concentration of fragments surrounded by an elongated coma, and presents previously unreported sharp arc-like and a narrow linear features. The morphology and evolution of these characteristics independently point toward a brief event on 2016 March 6. The arc and the linear feature can be reproduced by large particles on a ring, moving at ~2.5 m/s. The expansion of the ring defines a cone with a ~40deg half-opening. We propose that G1 was hit by a small object which caused its (partial or total) disruption, and that the ring corresponds to large fragments ejected during the final stages of the crater formation.}

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 introduce a new technique to estimate the comet nuclear size frequency distribution (SFD) that combines a cometary activity model with a survey simulation and apply it to 150 long period comets (LPC) detected by the Pan-STARRS1 near-Earth object survey. The debiased LPC size-frequency distribution is in agreement with previous estimates for large comets with nuclear diameter $>\sim 1$~km but we measure a significant drop in the SFD slope for small objects with diameters $<1$~km and approaching only $100$~m diameter. Large objects have a slope $\alpha_{big} = 0.72 \pm 0.09 (stat.) \pm 0.15 (sys.)$ while small objects behave as $\alpha_{small} = 0.07 \pm 0.03 (stat.) \pm 0.09 (sys.)$ where the SFD is $\propto 10^{\alpha H_N}$ and $H_N$ represents the cometary nuclear absolute magnitude. The total number of LPCs that are $>1$~km diameter and have perihelia $q<10$~au is $0.46 \pm 0.15 \times 10^9$ while there are only $2.4 \pm 0.5 (stat.) \pm 2 (sys.) \times 10^9$ objects with diameters $>100$~m due to the shallow slope of the SFD for diameters $<1$~m. We estimate that the total number of `potentially active' objects with diameters $\ge 1$~km in the Oort cloud, objects that would be defined as LPCs if their perihelia evolved to $<10$~au, is $(1.5\pm1)\times10^{12}$ with a combined mass of $1.3\pm0.9 \, M_{Earth}$. The debiased LPC orbit distribution is broadly in agreement with expectations from contemporary dynamical models but there are discrepancies that could point towards a future ability to disentangle the relative importance of stellar perturbations and galactic tides in producing the LPC population.

In order to identify the sources of the observed diffuse high-energy neutrino flux, it is crucial to discover their electromagnetic counterparts. IceCube began releasing alerts for single high-energy ($E > 60$ TeV) neutrino detections with sky localisation regions of order 1 deg radius in 2016. We used Pan-STARRS1 to follow-up five of these alerts during 2016-2017 to search for any optical transients that may be related to the neutrinos. Typically 10-20 faint ($m < 22.5$ mag) extragalactic transients are found within the Pan-STARRS1 footprints and are generally consistent with being unrelated field supernovae (SNe) and AGN. We looked for unusual properties of the detected transients, such as temporal coincidence of explosion epoch with the IceCube timestamp. We found only one transient that had properties worthy of a specific follow-up. In the Pan-STARRS1 imaging for IceCube-160427A (probability to be of astrophysical origin of $\sim$50 %), we found a SN PS16cgx, located at 10.0' from the nominal IceCube direction. Spectroscopic observations of PS16cgx showed that it was an H-poor SN at z = 0.2895. The spectra and light curve resemble some high-energy Type Ic SNe, raising the possibility of a jet driven SN with an explosion epoch temporally coincident with the neutrino detection. However, distinguishing Type Ia and Type Ic SNe at this redshift is notoriously difficult. Based on all available data we conclude that the transient is more likely to be a Type Ia with relatively weak SiII absorption and a fairly normal rest-frame r-band light curve. If, as predicted, there is no high-energy neutrino emission from Type Ia SNe, then PS16cgx must be a random coincidence, and unrelated to the IceCube-160427A. We find no other plausible optical transient for any of the five IceCube events observed down to a 5$\sigma$ limiting magnitude of $m \sim 22$ mag, between 1 day and 25 days after detection.

On 2019 January 5 a streamer associated with the 4--10 km main-belt asteroid (6478)~Gault was detected by the ATLAS sky survey, a rare discovery of activity around a main-belt asteroid. Archival data from ATLAS and Pan-STARRS1 show the trail in early December 2018, but not between 2010 and January 2018. The feature has significantly changed over one month, perfectly matching predictions of pure dust dynamical evolution and changes in observing geometry for a short release of dust around 2018 October 28. Follow-up observations with HST show a second narrow trail corresponding to a brief release of dust on 2018 December 30. Both releases occurred with negligible velocity. We find the dust grains to be fairly large, with power-law size distributions in the $10^{-5} - 10^{-3}$~m range and power-law indices of $\sim -1.5$. Three runs of ground-based data find a signature of $\sim 2\,\rm h$ rotation, close to the rotational limit, suggesting that the activity is the result of landslides or reconfigurations after YORP spin-up.

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.

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.

A class of asteroids, called large super-fast rotators (large SFRs), have rotation periods shorter than 2 hours and diameters larger than ~ 0.3 km. They pose challenges to the usual interior rubble-pile structure unless a relatively high bulk density is assumed. So far, only six large SFRs have been found. Therefore, we present a survey of asteroid rotation period using Pan-STARRS 1 telescopes during 2016 October 26 to 31 to search more large SFRs and study their properties. A total of 876 reliable rotation periods are measured, among which seven are large SFRs, thereby increasing the inventory of known large SFRs. These seven newly discovered large SFRs have diverse colors and locations in the main asteroid belt, suggesting that the taxonomic tendency and the location preference in the inner main belt of the six perviously known large SFRs could be a bias due to various observational limits. Interestingly, five out of the seven newly discovered large SFRs are mid main-belt asteroids. Considering the rare discovery rates of large SFR in the previously similar surveys (Chang et al., 2015, 2016) and the survey condition in this work, the chance of detecting a large SFR in the inner main belt seems to be relatively low. This probably suggests that the inner main belt harbors less large SFRs than the mid main belt. From our survey, we also found the drop in number appearing at f > 5 rev/day on the spin-rate distribution for the outer main-belt asteroids of D < 3 km, which was reported for the inner and mid main belt by (Chang et al., 2015, 2016).

We present new spatial models and distance estimates for globular clusters (GC) and dwarf spheroidals (dSphs) orbiting our Galaxy based on RR Lyrae (RRab) stars in the Pan-STARRS1 (PS1) 3$\pi$ survey. Using the PS1 sample of RRab stars from Sesar et al. (2017) in 16 globular clusters and 5 dwarf galaxies, we fit structural models in $(l,b,D)$ space; for 13 globular clusters and 6 dwarf galaxies, we give only their mean heliocentric distance $D$. We verify the accuracy of the period-luminosity (PL) relations used in Sesar et al. (2017) to constrain the distance to those stars, and compare them to period-luminosity-metallicity (PLZ) relations using metallicities from Carretta et al. (2009). We compare our Sesar et al. (2017) distances to the parallax-based \textit{Gaia} DR2 distance estimates from Bailer-Jones et al. (2018), and find our distances to be consistent and considerably more precise.

We present a nuclear transient event, PS1-13cbe, that was first discovered in the Pan-STARRS1 survey in 2013. The outburst occurred in the nucleus of the galaxy SDSS J222153.87+003054.2 at $z = 0.12355$, which was classified as a Seyfert 2 in a pre-outburst archival Sloan Digital Sky Survey (SDSS) spectrum. PS1-13cbe showed the appearance of strong broad H$\alpha$ and H$\beta$ emission lines and a non-stellar continuum in a Magellan spectrum taken 57 days after the peak of the outburst that resembled the characteristics of a Seyfert 1. These broad lines were not present in the SDSS spectrum taken a decade earlier and faded away within two years, as observed in several late-time MDM spectra. We argue that the dramatic appearance and disappearance of the broad lines and factor of $\sim 8$ increase in the optical continuum is most likely caused by variability in the pre-existing accretion disk than a tidal disruption event, supernova, or variable obscuration. The timescale for the turn-on of the optical emission of $\sim 70$ days observed in this transient is among the shortest observed in a "changing look" active galactic nucleus.

A large sample of white dwarfs is selected by both proper motion and colours from the Pan-STARRS 1 3{\pi} Steradian Survey Processing Version 2 to construct the White Dwarf Luminosity Functions of the discs and halo in the solar neighbourhood. Four-parameter astrometric solutions were recomputed from the epoch data. The generalised maximum volume method is then used to calculate the density of the populations. After removal of crowded areas near the Galactic plane and centre, the final sky area used by this work is 7.833 sr, which is 83% of the 3{\pi} sky and 62% of the whole sky. By dividing the sky using Voronoi tessellation, photometric and astrometric uncertainties are recomputed at each step of the integration to improve the accuracy of the maximum volume. Interstellar reddening is considered throughout the work. We find a disc-to-halo white dwarf ratio of about 100.

Super-luminous supernovae (SLSNe) are rare events defined as being significantly more luminous than normal terminal stellar explosions. The source of the extra powering needed to achieve such luminosities is still unclear. Discoveries in the local Universe (i.e. $z<0.1$) are scarce, but afford dense multi-wavelength observations. Additional low-redshift objects are therefore extremely valuable. We present early-time observations of the type I SLSN ASASSN-18km/SN~2018bsz. These data are used to characterise the event and compare to literature SLSNe and spectral models. Host galaxy properties are also analysed. Optical and near-IR photometry and spectroscopy were analysed. Early-time ATLAS photometry was used to constrain the rising light curve. We identified a number of spectral features in optical-wavelength spectra and tracked their time evolution. Finally, we used archival host galaxy photometry together with HII region spectra to constrain the host environment. ASASSN-18km/SN~2018bsz is found to be a type I SLSN in a galaxy at a redshift of 0.0267 (111 Mpc), making it the lowest-redshift event discovered to date. Strong CII lines are identified in the spectra. Spectral models produced by exploding a Wolf-Rayet progenitor and injecting a magnetar power source are shown to be qualitatively similar to ASASSN-18km/SN~2018bsz, contrary to most SLSNe-I that display weak/non-existent CII lines. ASASSN-18km/SN~2018bsz displays a long, slowly rising, red 'plateau' of $>$26 days, before a steeper, faster rise to maximum. The host has an absolute magnitude of --19.8 mag ($r$), a mass of M$_{*}$ = 1.5$^{+0.08}_{-0.33}$ $\times$10$^{9}$ M$_{\odot}$ , and a star formation rate of = 0.50$^{+2.22}_{-0.19}$ M$_{\odot}$ yr$^{-1}$. A nearby HII region has an oxygen abundance (O3N2) of 8.31$\pm$0.01 dex.

We present the largest Cepheid sample in M31 based on the complete Pan-STARRS1 survey of Andromeda (PAndromeda) in the $r_{\mathrm{P1}}$ , $i_{\mathrm{P1}}$ and $g_{\mathrm{P1}}$ bands. We find 2686 Cepheids with 1662 fundamental mode Cepheids, 307 first-overtone Cepheids, 278 type II Cepheids and 439 Cepheids with undetermined Cepheid type. Using the method developed by Kodric et al. (2013) we identify Cepheids by using a three dimensional parameter space of Fourier parameters of the Cepheid light curves combined with a color cut and other selection criteria. This is an unbiased approach to identify Cepheids and results in a homogeneous Cepheid sample. The Period-Luminosity relations obtained for our sample have smaller dispersions than in our previous work. We find a broken slope that we previously observed with HST data in Kodric et al. (2015), albeit with a lower significance.

Superluminous supernovae (SLSNe) are at least $\sim$5 times more luminous than common supernovae (SNe). Especially hydrogen-poor SLSN-I are difficult to explain with conventional powering mechanisms. One possible scenario that might explain such luminosities is that SLSNe-I are powered by an internal engine, such as a magnetar or an accreting black hole. Strong magnetic fields or collimated jets can circularly polarize light. In this work, we measured circular polarization of two SLSNe-I with the FOcal Reducer and low dispersion Spectrograph (FORS2) mounted at the ESO's Very Large Telescope (VLT). PS17bek, a fast evolving SLSN-I, was observed around peak, while OGLE16dmu, a slowly evolving SLSN-I, was observed 100 days after maximum. Neither SLSN shows evidence of circularly polarized light, however, these non-detections do not rule out the magnetar scenario as the powering engine for SLSNe-I. We calculate the strength of the magnetic field and the expected circular polarization as a function of distance from the magnetar, which decreases very fast. Additionally, we observed no significant linear polarization for PS17bek at four epochs, suggesting that the photosphere near peak is close to spherical symmetry.

We characterize the spatial density of the Pan-STARRS1 (PS1) sample of RR Lyrae stars, to study the properties of the old Galactic stellar halo as traced by RRab stars. This sample of 44,403 sources spans Galactocentric radii of $0.55 \; \mathrm{kpc} \leq R_{\mathrm{gc}} \leq 141 \; \mathrm{kpc}$ with a distance precision of 3\% and thus is able to trace the halo out to larger distances than most previous studies. After excising stars that are attributed to dense regions such as stellar streams, the Galactic disc and bulge as well as halo globular clusters, the sample contains ${\sim}11,000$ sources within $20 \; \mathrm{kpc} \leq R_{\mathrm{gc}} \leq 131 \; \mathrm{kpc}$. We then apply forward modeling using ellipsoidal stellar density models $\rho(l,b,R_{\mathrm{gc}})$ both with a constant and a radius-dependent halo flattening $q(R_{\mathrm{gc}})$. Assuming constant flattening $q$, the distribution of the sources is reasonably well fit from $20 \; \mathrm{kpc}$ to $131 \; \mathrm{kpc}$ by a single power law with $n=4.40^{+0.05}_{-0.04}$ and $q=0.918^{+0.016}_{-0.014}$. The distance distribution is fit comparably well by an Einasto profile with $n=9.53^{+0.27}_{-0.28}$, an effective radius $r_{\mathrm{eff}}=1.07 \pm 0.10 \; \mathrm{kpc}$ and a halo flattening of $q=0.923 \pm 0.007$. If we allow for a radius-dependent flattening $q(R_{\mathrm{gc}})$, we find evidence for a distinct flattening of $q{\sim}0.8$ of the inner halo at ${\sim} 25 \; \mathrm{kpc}$. Additionally, we find that the south Galactic hemisphere is more flattened than the north Galactic hemisphere. The results of our work are largely consistent with many earlier results, e.g. \cite{Watkins2009}, \cite{Iorio2017}. We find that the stellar halo, as traced in RR Lyrae stars, exhibits a substantial number of further significant over- and underdensities, even after all known overdensities have been masked.

We show that `Oumuamua's excited spin could be in a high energy LAM state, which implies that its shape could be far from the highly elongated shape found in previous studies. CLEAN and ANOVA algorithms are used to analyze `Oumuamua's lightcurve using 818 observations over 29.3~days. Two fundamental periodicities are found at frequencies (2.77$\pm$0.11) and (6.42$\pm$0.18)~cycles/day, corresponding to (8.67$\pm$0.34)~h and (3.74$\pm$0.11)~h, respectively. The phased data show that the lightcurve does not repeat in a simple manner, but approximately shows a double minimum at 2.77~cycles/day and a single minimum at 6.42~cycles/day. This is characteristic of an excited spin state. `Oumuamua could be spinning in either the long (LAM) or short (SAM) axis mode. For both, the long axis precesses around the total angular momentum vector with an average period of (8.67$\pm$0.34)~h. For the three LAMs we have found, the possible rotation periods around the long axis are 6.58, 13.15, or 54.48~h, with 54.48~h being the most likely. `Oumuamua may also be nutating with respective periods of half of these values. We have also found two possible SAM states where `Oumuamua oscillates around the long axis with possible periods at 13.15 and 54.48~h, the latter as the most likely. In this case any nutation will occur with the same periods. Determination of the spin state, the amplitude of the nutation, the direction of the TAMV, and the average total spin period may be possible with a direct model fit to the lightcurve. We find that `Oumuamua is "cigar-shaped"', if close to its lowest rotational energy, and an extremely oblate spheroid if close to its highest energy state for its total angular momentum.

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.

We present a simple but effective technique for measuring angular variation in $R_V$ across the sky. We divide stars from the Pan-STARRS1 catalog into Healpix pixels and determine the posterior distribution of reddening and $R_V$ for each pixel using two independent Monte Carlo methods. We find the two methods to be self-consistent in the limits where they are expected to perform similarly. We also find some agreement with high-precision photometric studies of $R_V$ in Perseus and Ophiuchus, as well as with a map of reddening near the Galactic plane based on stellar spectra from APOGEE. While current studies of $R_V$ are mostly limited to isolated clouds, we have developed a systematic method for comparing $R_V$ values for the majority of observable dust. This is a proof of concept for a more rigorous Galactic reddening map.x

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.

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 present a new 3D map of interstellar dust reddening, covering three quarters of the sky (declinations greater than -30 degrees) out to a distance of several kiloparsecs. The map is based on high-quality stellar photometry of 800 million stars from Pan-STARRS 1 and 2MASS. We divide the sky into sightlines containing a few hundred stars each, and then infer stellar distances and types, along with the line-of-sight dust distribution. Our new map incorporates a more accurate average extinction law and an additional 1.5 years of Pan-STARRS 1 data, tracing dust to greater extinctions and at higher angular resolutions than our previous map. Out of the plane of the Galaxy, our map agrees well with 2D reddening maps derived from far-infrared dust emission. After accounting for a 15% difference in scale, we find a mean scatter of 10% between our map and the Planck far-infrared emission-based dust map, out to a depth of 0.8 mag in E(r-z), with the level of agreement varying over the sky. Our map can be downloaded at http://argonaut.skymaps.info, or by its DOI: 10.7910/DVN/LCYHJG.

Recent all-sky surveys have led to the discovery of new types of transients. These include stars disrupted by the central supermassive black hole, and supernovae that are 10-100 times more energetic than typical ones. However, the nature of even more energetic transients that apparently occur in the innermost regions of their host galaxies is hotly debated. Here we report the discovery of the most energetic of these to date: PS1-10adi, with a total radiated energy of ~ 2.3 x 10^52 erg. The slow evolution of its light curve and persistently narrow spectral lines over ~3 yr are inconsistent with known types of recurring black hole variability. The observed properties imply powering by shock interaction between expanding material and large quantities of surrounding dense matter. Plausible sources of this expanding material are a star that has been tidally disrupted by the central black hole, or a supernova. Both could satisfy the energy budget. For the former, we would be forced to invoke a new and hitherto unseen variant of a tidally disrupted star, while a supernova origin relies principally on environmental effects resulting from its nuclear location. Remarkably, we also discovered that PS1-10adi is not an isolated case. We therefore surmise that this new population of transients has previously been overlooked due to incorrect association with underlying central black hole activity.

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.

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 GRB170105A was detected by 3 satellites 19.04hrs after GW170104 and 4.10hrs 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 GRB170105A, 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.

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

We present photometry and spectroscopy of nine Type II-P/L supernovae (SNe) with redshifts in the 0.045 < z < 0.335 range, with a view to re-examining their utility as distance indicators. Specifically, we apply the expanding photosphere method (EPM) and the standardized candle method (SCM) to each target, and find that both methods yield distances that are in reasonable agreement with each other. The current record-holder for the highest-redshift spectroscopically confirmed SN II-P is PS1-13bni (z = 0.335 +0.009 -0.012), and illustrates the promise of Type II SNe as cosmological tools. We updated existing EPM and SCM Hubble diagrams by adding our sample to those previously published. Within the context of Type II SN distance measuring techniques, we investigated two related questions. First, we explored the possibility of utilising spectral lines other than the traditionally used Fe II 5169 to infer the photospheric velocity of SN ejecta. Using local well-observed objects, we derive an epoch-dependent relation between the strong Balmer line and Fe II 5169 velocities that is applicable 30 to 40 days post-explosion. Motivated in part by the continuum of key observables such as rise time and decline rates exhibited from II-P to II-L SNe, we assessed the possibility of using Hubble-flow Type II-L SNe as distance indicators. These yield similar distances as the Type II-P SNe. Although these initial results are encouraging, a significantly larger sample of SNe II-L would be required to draw definitive conclusions.

We present a comprehensive and precise description of the Sagittarius (Sgr) stellar stream's 3D geometry as traced by its old stellar population. This analysis draws on the sample of ${\sim}44,000$ RR Lyrae (RRab) stars from the Pan-STARRS1 (PS1) 3$\pi$ survey (Hernitschek et al. 2016,Sesar et al. 2017b), which is ${\sim}80\%$ complete and ${\sim}90\%$ pure within 80~kpc, and extends to ${\gtrsim} 120$~kpc with a distance precision of ${\sim} 3\%$. A projection of RR Lyrae stars within $|\tilde{B}|_{\odot}<9^\circ$ of the Sgr stream's orbital plane reveals the morphology of both the leading and the trailing arms at very high contrast, across much of the sky. In particular, the map traces the stream near-contiguously through the distant apocenters. We fit a simple model for the mean distance and line-of-sight depth of the Sgr stream as a function of the orbital plane angle $\tilde{\Lambda}_{\odot}$, along with a power-law background-model for the field stars. This modeling results in estimates of the mean stream distance precise to ${\sim}1\%$ and it resolves the stream's line-of-sight depth. These improved geometric constraints can serve as new constraints for dynamical stream models.

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.

Comet C/2017 K2 (PANSTARRS) was discovered by the Pan-STARRS1 (PS1) Survey on 2017 May 21 at a distance 16.09 au from the Sun, the second most distant discovery of an active comet. Pre-discovery images in the PS1 archive back to 2014 and additional deep CFHT images between 2013 May 10-13 showed the comet to be active at 23.75 au. We derive an upper limit to the nucleus radius of $R_N$=80 km, assuming a 4\% albedo. The spectral reflectivity of the comet surface is similar to "fresh" regions seen on comet 67P/Churyumov-Gerasimenko using the $Rosetta$ OSIRIS camera. Pre-discovery photometry combined with new data obtained with Megacam on the CFHT show that the activity is consistent with CO-ice sublimation and inconsistent with CO$_2$-ice sublimation. The ice sublimation models were run out to perihelion in 2022 at 1.8 au to predict the CO production rates, assuming that the outgassing area does not change. Assuming a canonical 4\% active surface area for water-ice sublimation, we present production rate ratios, $Q_{\rm CO}$/$Q_{\rm H2O}$, for a range of nucleus sizes. Comparing these results with other CO-rich comets we derive a lower limit to the nucleus radius of $\sim$14 km. We present predictions for $Q_{\rm CO}$ at a range of distances that will be useful for planning observations with JWST and large ground-based facilities.

Thick back-illuminated deep-depletion CCDs have superior quantum efficiency over previous generations of thinned and traditional thick CCDs. As a result, they are being used for wide-field imaging cameras in several major projects. We use observations from the Pan-STARRS $3\pi$ survey to characterize the behavior of the deep-depletion devices used in the Pan-STARRS1 Gigapixel Camera. We have identified systematic spatial variations in the photometric measurements and stellar profiles which are similar in pattern to the so-called "tree rings" identified in devices used by other wide-field cameras (e.g., DECam and Hypersuprime Camera). The tree-ring features identified in these other cameras result from lateral electric fields which displace the electrons as they are transported in the silicon to the pixel location. In contrast, we show that the photometric and morphological modifications observed in the GPC1 detectors are caused by variations in the vertical charge transportation rate and resulting charge diffusion variations.

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.

Context. A lot of photometric data is produced by surveys such as Pan-STARRS, LONEOS, WISE or Catalina. These data are a rich source of information about the physical properties of asteroids. There are several possible approaches for utilizing these data. Lightcurve inversion is a typical method that works with individual asteroids. Our approach in this paper is statistical when we focused on large groups of asteroids like dynamical families and taxonomic classes, and the data were not sufficient for individual models. Aims. Our aim was to study the distributions of shape elongation $b/a$ and the spin axis latitude $\beta$ for various subpopulations of asteroids and to compare our results, based on Pan-STARRS1 survey, with statistics previously done using different photometric data (Lowell database, WISE data). Methods. We use the LEADER algorithm to compare the $b/a$ and $\beta$ distributions for different subpopulations of asteroids. The algorithm creates a cumulative distributive function (CDF) of observed brightness variations, and computes the $b/a$ and $\beta$ distributions using analytical basis functions that yield the observed CDF. A variant of LEADER is used to solve the joint distributions for synthetic populations to test the validity of the method. Results. When comparing distributions of shape elongation for groups of asteroids with different diameters $D$, we found that there are no differences for $D < 25$ km. We also constructed distributions for asteroids with different rotation periods and revealed that the fastest rotators with $P = 0 - 4$ h are more spheroidal than the population with $P = 4 - 8$ h.

We report the discovery of a $H_r = 3.4\pm0.1$ dwarf planet candidate by the Pan-STARRS Outer Solar System Survey. 2010 JO$_{179}$ is red with $(g-r)=0.88 \pm 0.21$, roughly round, and slowly rotating, with a period of $30.6$ hr. Estimates of its albedo imply a diameter of 600--900~km. Observations sampling the span between 2005--2016 provide an exceptionally well-determined orbit for 2010 JO$_{179}$, with a semi-major axis of $78.307\pm0.009$ au, distant orbits known to this precision are rare. We find that 2010 JO$_{179}$ librates securely within the 21:5 mean-motion resonance with Neptune on hundred-megayear time scales, joining the small but growing set of known distant dwarf planets on metastable resonant orbits. These imply a substantial trans-Neptunian population that shifts between stability in high-order resonances, the detached population, and the eroding population of the scattering disk.

We make use of a catalog of 1600 Pan-STARRS1 groups produced by the probability friends-of-friends algorithm to explore how the galaxy properties, i.e. the specific star formation rate (SSFR) and quiescent fraction, depend on stellar mass and group-centric radius. The work is the extension of Lin et al. (2014). In this work, powered by a stacking technique plus a background subtraction for contamination removal, a finer correction and more precise results are obtained than in our previous work. We find that while the quiescent fraction increases with decreasing group-centric radius the median SSFRs of star-forming galaxies in groups at fixed stellar mass drop slightly from the field toward the group center. This suggests that the major quenching process in groups is likely a fast mechanism. On the other hand, a reduction in SSFRs by ~0.2 dex is seen inside clusters as opposed to the field galaxies. If the reduction is attributed to the slow quenching effect, the slow quenching process acts dominantly in clusters. In addition, we also examine the density-color relation, where the density is defined by using a sixth-nearest neighbor approach. Comparing the quiescent fractions contributed from the density and radial effect, we find that the density effect dominates over the massive group or cluster galaxies, and the radial effect becomes more effective in less massive galaxies. The results support mergers and/or starvation as the main quenching mechanisms in the group environment, while harassment and/or starvation dominate in clusters.

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.

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.

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.

Using shape measurement techniques developed for weak lensing surveys we have identified three new ultracool binaries in the Pan-STARRS1 survey. Binary companions which are not completely resolved can still alter the shapes of stellar images. These shape distortions can be measured if PSF anisotropy caused by the telescope is properly accounted for. We show using both a sample of known binary stars and simulated binaries that we can reliably recover binaries wider than around 0.3" and with flux ratios greater than around 0.1. We then applied our method to a sample of ultracool dwarfs within 30pc with 293 objects having sufficient Pan-STARRS1 data for our method. In total we recovered all but one of the 11 binaries wider than 0.3" in this sample. Our one failure was a true binary detected with a significant but erroneously high ellipticity which led it to be rejected in our analysis. We identify three new binaries, one a simultaneous discovery, with primary spectral types M6.5, L1 and T0.5. These latter two were confirmed with Keck/NIRC2 follow-up imaging. This technique will be useful for identifying large numbers of stellar and substellar binaries in the upcoming LSST and DES sky surveys.

We present deep imaging observations, orbital dynamics, and dust tail model analyses of the double-component asteroid P/2016 J1 (J1-A and J1-B). The observations were acquired at the Gran Telescopio Canarias (GTC) and the Canada-France-Hawaii Telescope (CFHT) from mid March to late July, 2016. A statistical analysis of backward-in-time integrations of the orbits of a large sample of clone objects of P/2016 J1-A and J1-B shows that the minimum separation between them occurred most likely $\sim$2300 days prior to the current perihelion passage, i.e., during the previous orbit near perihelion. This closest approach was probably linked to a fragmentation event of their parent body. Monte Carlo dust tail models show that those two components became active simultaneously $\sim$250 days before the current perihelion, with comparable maximum loss rates of $\sim$0.7 kg s$^{-1}$ and $\sim$0.5 kg s$^{-1}$, and total ejected masses of 8$\times$10$^{6}$ kg and 6$\times$10$^{6}$ kg for fragments J1-A and J1-B, respectively. In consequence, the fragmentation event and the present dust activity are unrelated. The simultaneous activation times of the two components and the fact that the activity lasted 6 to 9 months or longer, strongly indicate ice sublimation as the most likely mechanism involved in the dust emission process.

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.

We present the discovery of a 360 AU separation T3 companion to the tight (3.1 AU) M4.5+M6.5 binary 2MASS J02132062+3648506. This companion was identified using Pan-STARRS1 data and, despite its relative proximity to the Sun (22.2$_{-4.0}^{+6.4}$ pc; Pan-STARRS1 parallax) and brightness ($J$=15.3), appears to have been missed by previous studies due to its position near a diffraction spike in 2MASS. The close M~dwarf binary has active X-ray and H$\alpha$ emission and shows evidence for UV flares. The binary's weak {\it GALEX} UV emission and strong Na I 8200\AA Na absorption leads us to an age range of $\sim$1-10Gyr. Applying this age range to evolutionary models implies the wide companion has a mass of 0.063$\pm$0.009\,$M_{\odot}$. 2MASS J0213+3648 C provides a relatively old benchmark close to the L/T transition and acts as a key, older comparison to the much younger early-T companions HN~Peg~B and GU~Psc~b.

We study the time lags between the continuum emission of quasars at different wavelengths, based on more than four years of multi-band ($g$, $r$, $i$, $z$) light-curves in the Pan-STARRS Medium Deep Fields. As photons from different bands emerge from different radial ranges in the accretion disk, the lags constrain the sizes of the accretion disks. We select 240 quasars with redshifts $z \approx 1$ or $z \approx 0.3$ that are relatively emission line free. The light curves are sampled from day to month timescales, which makes it possible to detect lags on the scale of the light crossing time of the accretion disks. With the code JAVELIN, we detect typical lags of several days in the rest frame between the $g$ band and the $riz$ bands. The detected lags are $\sim 2-3$ times larger than the light crossing time estimated from the standard thin disk model, consistent with the recently measured lag in NGC5548 and micro-lensing measurements of quasars. The lags in our sample are found to increase with increasing luminosity. Furthermore, the increase in lags going from $g-r$ to $g-i$ and then to $g-z$ is slower than predicted in the thin disk model, particularly for high luminosity quasars. The radial temperature profile in the disk must be different from what is assumed. We also find evidence that the lags decrease with increasing line ratios between ultraviolet FeII lines and MgII, which may point to changes in the accretion disk structure at higher metallicity.

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)

A key challenge in understanding the feedback mechanism of AGN in Brightest Cluster Galaxies (BCGs) is the inherent rarity of catching an AGN during its strong outburst phase. This is exacerbated by the ambiguity of differentiating between AGN and clusters in X-ray observations. If there is evidence for an AGN then the X-ray emission is commonly assumed to be dominated by the AGN emission, introducing a selection effect against the detection of AGN in BCGs. In order to recover these 'missing' clusters, we systematically investigate the colour-magnitude relation around some ~3500 ROSAT All Sky Survey selected AGN, looking for signs of a cluster red sequence. Amongst our 22 candidate systems, we independently rediscover several confirmed systems, where a strong AGN resides in a central galaxy. We compare the X-ray luminosity to red sequence richness distribution of our AGN candidate systems with that of a similarly selected comparison sample of ~1000 confirmed clusters and identify seven 'best' candidates (all of which are BL Lac objects), where the X-ray flux is likely to be a comparable mix between cluster and AGN emission. We confirm that the colours of the red sequence are consistent with the redshift of the AGN, that the colours of the AGN host galaxy are consistent with AGN, and, by comparing their luminosities with those from our comparison clusters, confirm that the AGN hosts are consistent with BCGs.

We determine the absolute magnitude (H) distribution (or size-frequency distribution, SFD; $N(H) \propto 10^{\alpha H}$ where $\alpha$ is the slope of the distribution) for near-Earth objects (NEO) with $13<H<30$ and Asteroid Retrieval Mission (ARM) targets with $27<H<31$ that were detected by the 1\st\ telescope of the Panoramic Survey Telescope and Rapid Response System - Pan-STARRS1 (e.g. Kaiser et al. 2002, Kaiser 2004, Hodapp et al. 2004). The NEO and ARM target detection efficiencies were calculated using the Greenstreet et al. (2012) NEO orbit distribution. The debiased Pan-STARRS1 NEO absolute magnitude distribution is more complex than a single slope power law - it shows two transitions - at H$\sim$16 from steep to shallow slope, and in the $21<H<23$ interval from a shallow to steep slope, which is consistent with other recent works (e.g. Mainzer et al. 2011c, Brown et al. 2013, Harris and D`Abramo 2015). We fit $\alpha = 0.48\pm0.02$ for NEOs with $13<H<16$, $\alpha = 0.33\pm0.01$ for NEOs with $16<H<22$, and $\alpha = 0.62\pm0.03$ for the smaller objects with $H>22$. There is also another change in slope from steep to shallow around H=27. The three ARM target candidates detected by Pan-STARRS1 in one year of surveying have a corrected SFD with slope $\alpha = 0.40^{+0.33}_{-0.45}$. We also show that the window for follow up observations of small (H$\gtrsim$22) NEOs with the NASA IRTF telescope and Arecibo and Goldstone radars are extremely short - on order of days, and procedures for fast response must be implemented in order to measure physical characteristics of small Earth-approaching objects. CFHT's MegaCam and Pan-STARRS1 have longer observing windows and are capable of following-up more NEOs due to their deeper limiting magnitudes and wider fields of view.

We present a systematic search for periodically varying quasars and supermassive black hole binary (SMBHB) candidates in the Pan-STARRS1 (PS1) Medium Deep Survey's MD09 field. From a color-selected sample of 670 quasars extracted from a multi-band deep-stack catalog of point sources, we locally select variable quasars and look for coherent periods with the Lomb-Scargle periodogram. Three candidates from our sample demonstrate strong variability for more than ~3 cycles, and their PS1 light curves are well fitted to sinusoidal functions. We test the persistence of the candidates' apparent periodic variations detected during the 4.2 years of the PS1 survey with archival photometric data from the SDSS Stripe 82 survey or new monitoring with the Large Monolithic Imager at the Discovery Channel Telescope. None of the three periodic candidates (including PSO J334.2028+1.4075) remain persistent over the extended baseline of 7 - 14 years, corresponding to a detection rate of < 1 in 670 quasars in a search area of 5 deg^2. Even though SMBHBs should be a common product of the hierarchal growth of galaxies, and periodic variability in SMBHBs has been theoretically predicted, a systematic search for such signatures in a large optical survey is strongly limited by its temporal baseline and the "red noise" associated with normal quasar variability. We show that follow-up long-term monitoring (>5 cycles) is crucial to our search for these systems.

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

In this work we report the detection of seven Neptune Trojans (NTs) in the Pan-STARRS 1 (PS1) survey. Five of these are new discoveries, consisting of four L4 Trojans and one L5 Trojan. Our orbital simulations show that the L5 Trojan stably librates for only several million years. This suggests that the L5 Trojan must be of recent capture origin. On the other hand, all four new L4 Trojans stably occupy the 1:1 resonance with Neptune for more than 1 Gyr. They can, therefore, be of primordial origin. Our survey simulation results show that the inclination width of the Neptune Trojan population should be between $7^{\circ}$ and $27^{\circ}$ at $>$ 95% confidence, and most likely $\sim 11^{\circ}$. In this paper, we describe the PS1 survey, the Outer Solar System pipeline, the confirming observations, and the orbital/physical properties of the new Neptune Trojans.

We present initial time-resolved observations of the split comet 332P/Ikeya-Murakami taken using the Hubble Space Telescope. Our images reveal a dust-bathed cluster of fragments receding from their parent nucleus at projected speeds in the range 0.06 to 3.5 m s$^{-1}$ from which we estimate ejection times from October to December 2015. The number of fragments with effective radii $\gtrsim$20 m follows a differential power law with index $\gamma$ = -3.6$\pm$0.6, while smaller fragments are less abundant than expected from an extrapolation of this power-law. We argue that, in addition to losses due to observational selection, torques from anisotropic outgassing are capable of destroying the small fragments by driving them quickly to rotational instability. Specifically, the spin-up times of fragments $\lesssim$20 m in radius are shorter than the time elapsed since ejection from the parent nucleus. The effective radius of the parent nucleus is $r_e \le$ 275 m (geometric albedo 0.04 assumed). This is about seven times smaller than previous estimates and results in a nucleus mass at least 300 times smaller than previously thought. The mass in solid pieces, $2\times10^9$ kg, is about 4% of the mass of the parent nucleus. As a result of its small size, the parent nucleus also has a short spin-up time. Brightness variations in time-resolved nucleus photometry are consistent with rotational instability playing a role in the release of fragments.

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.

We describe 2016 January to April observations of the fragments of 332P/Ikeya-Murakami, a comet earlier observed in a 2010 October outburst (Ishiguro et al 2014). We present photometry of the fragments, and perform simulations to infer the time of breakup. We argue that the eastern-most rapidly brightening fragment ($F4$) best corresponds to the original nucleus, rather than the initial bright fragment $F1$. We compute radial and tangential non-gravitational parameters, $A_1 = (1.5 \pm 0.4) \times 10^{-8}$ AU day$^{-2}$ and $(7.2 \pm 1.9) \times 10^{-9}$ AU day$^{-2}$; both are consistent with zero at the $4\sigma$ level. Monte Carlo simulations indicate that the fragments were emitted on the outbound journey well after the 2010 outburst, with bright fragment $F1$ splitting in mid--2013 and the fainter fragments within months of the 2016 January recovery. Western fragment $F7$ is the oldest, dating from 2011. We suggest that the delayed onset of the splitting is consistent with a self-propagating crystallization of water ice.

We present a panoramic map of the entire Milky Way halo north of dec~-30 degrees (~30,000 deg^2), constructed by applying the matched-filter technique to the Pan-STARRS1 3Pi Survey dataset. Using single-epoch photometry reaching to g~22, we are sensitive to stellar substructures with heliocentric distances between 3.5 and ~35 kpc. We recover almost all previously-reported streams in this volume and demonstrate that several of these are significantly more extended than earlier datasets have indicated. In addition, we also report five new candidate stellar streams. One of these features appears significantly broader and more luminous than the others and is likely the remnant of a dwarf galaxy. The other four streams are consistent with a globular cluster origin, and three of these are rather short in projection (<10 degrees), suggesting that streams like Ophiuchus may not be that rare. Finally, a significant number of more marginal substructures are also revealed by our analysis; many of these features can also be discerned in matched-filter maps produced by other authors from SDSS data, and hence they are very likely to be genuine. However, the extant 3Pi data is currently too shallow to determine their properties or produce convincing CMDs. The global view of the Milky Way provided by Pan-STARRS1 provides further evidence for the important role of both globular cluster disruption and dwarf galaxy accretion in building the Milky Way's stellar halo.

Luminous quasars at z>5.6 can be studied in detail with the current generation of telescopes and provide us with unique information on the first gigayear of the universe. Thus far these studies have been statistically limited by the number of quasars known at these redshifts. Such quasars are rare and therefore wide-field surveys are required to identify them and multiwavelength data are needed to separate them efficiently from their main contaminants, the far more numerous cool dwarfs. In this paper, we update and extend the selection for z~6 quasars presented in Banados et al. (2014) using the Pan-STARRS1 (PS1) survey. We present the PS1 distant quasar sample, which currently consists of 124 quasars in the redshift range 5.6<z<6.7 that satisfy our selection criteria. Seventy-seven of these quasars have been discovered with PS1, and 63 of them are newly identified in this paper. We present composite spectra of the PS1 distant quasar sample. This sample spans a factor of ~20 in luminosity and shows a variety of emission line properties. The number of quasars at z>5.6 presented in this work almost double the quasars previously known at these redshifts, marking a transition phase from studies of individual sources to statistical studies of the high-redshift quasar population, which was impossible with earlier, smaller samples.

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.

Although the majority of Centaurs are thought to have originated in the scattered disk, with the high-inclination members coming from the Oort cloud, the origin of the high inclination component of trans-Neptunian objects (TNOs) remains uncertain. We report the discovery of a retrograde TNO, which we nickname "Niku", detected by the Pan-STARRS 1 Outer Solar System Survey. Our numerical integrations show that the orbital dynamics of Niku are very similar to that of 2008 KV$_{42}$ (Drac), with a half-life of $\sim 500$ Myr. Comparing similar high inclination TNOs and Centaurs ($q > 10$ AU, $a < 100$ AU and $i > 60^\circ$), we find that these objects exhibit a surprising clustering of ascending node, and occupy a common orbital plane. This orbital configuration has high statistical significance: 3.8-$\sigma$. An unknown mechanism is required to explain the observed clustering. This discovery may provide a pathway to investigate a possible reservoir of high-inclination objects.

We discuss 76 large amplitude transients (Delta-m>1.5) occurring in the nuclei of galaxies, nearly all with no previously known Active Galactic Nucleus (AGN). They have been discovered as part of the Pan-STARRS1 (PS1) 3pi survey, by comparison with SDSS photometry a decade earlier, and then monitored with the Liverpool Telescope, and studied spectroscopically with the William Herschel Telescope (WHT). Based on colours, light curve shape, and spectra, these transients fall into four groups. A few are misclassified stars or objects of unknown type. Some are red/fast transients and are known or likely nuclear supernovae. A few are either radio sources or erratic variables and so likely blazars. However the majority (~66%) are blue and evolve slowly, on a timescale of years. Spectroscopy shows them to be AGN at z~ 0.3 - 1.4, which must have brightened since the SDSS photometry by around an order of magnitude. It is likely that these objects were in fact AGN a decade ago, but too weak to be recognised by SDSS; they could then be classed as "hypervariable" AGN. By searching the SDSS Stripe 82 quasar database, we find 15 similar objects. We discuss several possible explanations for these slow blue hypervariables - (i) unusually luminous tidal disruption events; (ii) extinction events; (iii) changes in accretion state; and (iv) large amplitude microlensing by stars in foreground galaxies. A mixture of explanations (iii) and (iv) seems most likely. Both hold promise of considerable new insight into the AGN phenomenon.

We report on an ongoing project to investigate activity in the M dwarf stellar population observed by the Pan-STARRS 1 Medium Deep Survey (PS1-MDS). Using a custom-built pipeline, we refine an initial sample of $\approx$ 4 million sources in PS1-MDS to a sample of 184,148 candidate cool stars using color cuts. Motivated by the well-known relationship between rotation and stellar activity, we use a multi-band periodogram analysis and visual vetting to identify 271 sources that are likely rotating M dwarfs. We derive a new set of polynomials relating M dwarf PS1 colors to fundamental stellar parameters and use them to estimate the masses, distances, effective temperatures, and bolometric luminosities of our sample. We present a catalog containing these values, our measured rotation periods, and cross-matches to other surveys. Our final sample spans periods of $\lesssim$1-130 days in stars with estimated effective temperatures of $\approx$ 2700-4000 K. Twenty-two of our sources have X-ray cross-matches, and they are found to be relatively X-ray bright as would be expected from selection effects. Our data set provides evidence that Kepler-based searches have not been sensitive to very slowly-rotating stars ($P_{\rm rot} \gtrsim 70$ d), implying that the observed emergence of very slow rotators in studies of low-mass stars may be a systematic effect. We also see a lack of low-amplitude ($<$2%) variability in objects with intermediate (10-40 d) rotation periods, which, considered in conjunction with other observational results, may be a signpost of a loss of magnetic complexity associated with a phase of rapid spin-down in intermediate-age M dwarfs. This work represents just a first step in exploring stellar variability in data from the PS1-MDS and, in the farther future, LSST.

We present a method to identify distant solar system objects in long-term wide-field asteroid survey data, and conduct a search for them in the Pan-STARRS1 (PS1) image data acquired from 2010 to mid-2015. We demonstrate that our method is able to find multi-opposition orbital links, and we present the resulting orbital distributions which consist of 154 Centaurs, 255 classical Trans-Neptunian Objects (TNOs), 121 resonant TNOs, 89 Scattered Disc Objects (SDOs) and 10 comets. Our results show more than half of these are new discoveries, including a newly discovered 19th magnitude TNO. Our identified objects do not show clustering in their argument of perihelia, which if present, might support the existence of a large unknown planetary-sized object in the outer solar system.

We describe development and application of a Global Astrometric Solution (GAS) to the problem of Pan-STARRS1 (PS1) astrometry. Current PS1 astrometry is based on differential astrometric measurements using 2MASS reference stars, thus PS1 astrometry inherits the errors of the 2MASS catalog. The GAS, based on a single, least squares adjustment to approximately 750k grid stars using over 3000 extragalactic objects as reference objects, avoids this catalog-to-catalog propagation of errors to a great extent. The GAS uses a relatively small number of Quasi-Stellar Objects (QSOs, or distant AGN) with very accurate (<1 mas) radio positions, referenced to the ICRF2. These QSOs provide a hard constraint in the global least squares adjustment. Solving such a system provides absolute astrometry for all the stars simultaneously. The concept is much cleaner than conventional astrometry but is not easy to perform for large catalogs. In this paper we describe our method and its application to Pan-STARRS1 data. We show that large-scale systematic errors are easily corrected but our solution residuals for position (~60 mas) are still larger than expected based on simulations (~10 mas). We provide a likely explanation for the reason the small-scale residual errors are not corrected in our solution as would be expected.

Despite their reputation as being "red and dead", the unique environment inhabited by Brightest Cluster Galaxies (BCGs) can often lead to a self-regulated feedback cycle between radiatively cooling intracluster gas and star formation and AGN activity in the BCG. However the prevalence of "active" BCGs, and details of the feedback involved, are still uncertain. We have performed an optical, UV and Mid-IR photometric analysis of the BCGs in 981 clusters at 0.03 < z < 0.5, selected from the ROSAT All Sky Survey. Using Pan-STARRS PS1 3pi, GALEX and WISE survey data we look for BCGs with photometric colours which deviate from that of the bulk population of passive BCGs - indicative of AGN and/or star formation activity within the BCG. We find that whilst the majority of BCGs are consistent with being passive, at least 14% of our BCGs show a significant colour offset from passivity in at least one colour index. And, where available, supplementary spectroscopy reveals the majority of these particular BCGs show strong optical emission lines. On comparing BCG "activity" with the X-ray luminosity of the host cluster, we find that BCGs showing a colour offset are preferentially found in the more X-ray luminous clusters, indicative of the connection between BCG "activity" and the intracluster medium.

We study the properties of 975 active galactic nuclei (AGN) selected by variability in the Pan-STARRS1 Medium-Deep Survey. Using complementary multi wavelength data from the ultraviolet to the far-infrared, we use SED fitting to determine the AGN and host properties at $z<1$, and compare to a well-matched control sample. We confirm the trend previously observed that the variability amplitude decreases with AGN luminosity, but on the other hand, we observe that the slope of this relation steepens with wavelength resulting in a "redder when brighter" trend at low luminosities. Our results show that AGN are hosted by more massive hosts than control sample galaxies, while the restframe, dust-corrected $NUV-r$ color distribution of AGN hosts is similar to control galaxies. We find a positive correlation between the AGN luminosity and star formation rate (SFR), independent of redshift. AGN hosts populate the whole range of SFRs within and outside the Main Sequence of star forming galaxies. Comparing the distribution of AGN hosts and control galaxies, we show that AGN hosts are less likely to be hosted by quiescent galaxies, but more likely to be hosted by Main Sequence or starburst galaxies.

Using the Pan-STARRS1 survey, we derive limiting magnitude, spatial completeness and density maps that we use to probe the three dimensional structure and estimate the stellar mass of the so-called Monoceros Ring. The Monoceros Ring is an enormous and complex stellar sub-structure in the outer Milky Way disk. It is most visible across the large Galactic Anticenter region, 120 < l < 240 degrees, -30 < b < +40 degrees. We estimate its stellar mass density profile along every line of sight in 2 X 2 degree pixels over the entire 30,000 square degree Pan-STARRS1 survey using the previously developed MATCH software. By parsing this distribution into a radially smooth component and the Monoceros Ring, we obtain its mass and distance from the Sun along each relevant line of sight. The Monoceros Ring is significantly closer to us in the South (6 kpc) than in the North (9 kpc). We also create 2D cross sections parallel to the Galactic plane that show 135 degrees of the Monoceros Ring in the South and 170 degrees of the Monoceros Ring in the North. We show that the Northern and Southern structures are also roughly concentric circles, suggesting that they may be a wave rippling from a common origin. Excluding the Galactic plane, we observe an excess stellar mass of 4 million solar masses across 120 < l < 240 degrees. If we interpolate across the Galactic plane, we estimate that this region contains 8 million solar masses. If we assume (somewhat boldly) that the Monoceros Ring is a set of two Galactocentric rings, its total stellar mass is 60 million solar masses. Finally, if we assume that it is a set of two circles centered at a point 4 kpc from the Galactic center in the anti-central direction, as our data suggests, we estimate its stellar mass to be 40 million solar masses.

The rotational state of asteroids is controlled by various physical mechanisms including collisions, internal damping and the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect. We have analysed the changes in magnitude between consecutive detections of approximately 60,000 asteroids measured by the PanSTARRS 1 survey during its first 18 months of operations. We have attempted to explain the derived brightness changes physically and through the application of a simple model. We have found a tendency toward smaller magnitude variations with decreasing diameter for objects of 1 < D < 8 km. Assuming the shape distribution of objects in this size range to be independent of size and composition our model suggests a population with average axial ratios 1 : 0.85 \pm 0.13 : 0.71 \pm 0.13, with larger objects more likely to have spin axes perpendicular to the orbital plane.

Substellar members of young ($\lesssim$150 Myr) moving groups are valuable benchmarks to empirically define brown dwarf evolution with age and to study the low-mass end of the initial mass function. We have combined Pan-STARRS1 (PS1) proper motions with optical$-$IR photometry from PS1, 2MASS and $\textit{WISE}$ to search for substellar members of the AB Dor Moving Group within $\approx$50 pc and with spectral types of late-M to early-L, corresponding to masses down to $\approx$30 M$_{Jup}$ at the age of the group ($\approx$125 Myr). Including both photometry and proper motions allows us to better select candidates by excluding field dwarfs whose colors are similar to young AB~Dor Moving Group members. Our near-IR spectroscopy has identified six ultracool dwarfs (M6$-$L4; $\approx$30$-$100 M$_{Jup}$) with intermediate surface gravities (INT-G) as candidate members of the AB Dor Moving Group. We find another two candidate members with spectra showing hints of youth but consistent with field gravities. We also find four field brown dwarfs unassociated with the AB Dor Moving Group, three of which have INT-G gravity classification. While signatures of youth are present in the spectra of our $\approx$125 Myr objects, neither their $J-K$ nor $W1-W2$ colors are significantly redder than field dwarfs with the same spectral types, unlike younger ultracool dwarfs. We also determined PS1 parallaxes for eight of our candidates and one previously identified AB Dor Moving Group candidate. Although radial velocities (and parallaxes, for some) are still needed to fully assess membership, these new objects provide valuable insight into the spectral characteristics and evolution of young brown dwarfs.

The dust extinction curve is a critical component of many observational programs and an important diagnostic of the physics of the interstellar medium. Here we present new measurements of the dust extinction curve and its variation towards tens of thousands of stars, a hundred-fold larger sample than in existing detailed studies. We use data from the APOGEE spectroscopic survey in combination with ten-band photometry from Pan-STARRS1, 2MASS, and WISE. We find that the extinction curve in the optical through infrared is well characterized by a one-parameter family of curves described by R(V). The extinction curve is more uniform than suggested in past works, with sigma(R(V)) = 0.18, and with less than one percent of sight lines having R(V) > 4. Our data and analysis have revealed two new aspects of Galactic extinction: first, we find significant, wide-area variations in R(V) throughout the Galactic plane. These variations are on scales much larger than individual molecular clouds, indicating that R(V) variations must trace much more than just grain growth in dense molecular environments. Indeed, we find no correlation between R(V) and dust column density up to E(B-V) ~ 2. Second, we discover a strong relationship between R(V) and the far-infrared dust emissivity.

We present the first spectroscopic analysis of the faint and compact stellar system Draco II (Dra II, M_V=-2.9+/-0.8, r_h=19^{+8}_{-6} pc), recently discovered in the Pan-STARRS1 3{\pi} survey. The observations, conducted with DEIMOS on the Keck II telescope, establish some of its basic characteristics: the velocity data reveal a narrow peak with 9 member stars at a systemic heliocentric velocity <v_r>=-347.6^{+1.7}_{-1.8} km/s, thereby confirming Dra II is a satellite of the Milky Way; we infer a velocity dispersion with \sigma_{vr}=2.9+/-2.1 km/s (<8.4 km/s at the 95% confidence level), which implies log_{10}(M_{1/2})=5.5^{+0.4}_{-0.6} and log_{10}((M/L)_{1/2})=2.7^{+0.5}_{-0.8}, in Solar units; furthermore, very weak Calcium triplet lines in the spectra of the high signal-to-noise member stars imply [Fe/H]<-2.1, whilst variations in the line strengths of two stars with similar colours and magnitudes suggest a metallicity spread in Dra II. These new data cannot clearly discriminate whether Draco II is a star cluster or amongst the faintest, most compact, and closest dwarf galaxies. However, the sum of the three --- individually inconclusive --- pieces of evidence presented here, seems to favour the dwarf galaxy interpretation.

We present a systematic search for changing-look quasars based on repeat photometry from SDSS and Pan-STARRS1, along with repeat spectra from SDSS and SDSS-III BOSS. Objects with large, |\Delta g|>1 mag photometric variations in their light curves are selected as candidates to look for changes in broad emission line (BEL) features. Out of a sample of 1011 objects that satisfy our selection criteria and have more than one epoch of spectroscopy, we find 10 examples of quasars that have variable and/or "changing-look" BEL features. Four of our objects have emerging BELs; five have disappearing BELs, and one object shows tentative evidence for having both emerging and disappearing BELs. With redshifts in the range 0.20 < z < 0.63, this sample includes the highest-redshift changing-look quasars discovered to date. We highlight the quasar J102152.34+464515.6 at z = 0.204. Here, not only have the Balmer emission lines strongly diminished in prominence, including H$\beta$ all but disappearing, but the blue continuum $f_{\nu} \propto \nu^{1/3}$ typical of an AGN is also significantly diminished in the second epoch of spectroscopy. Using our selection criteria, we estimate that >15% of strongly variable luminous quasars display changing-look BEL features on rest-frame timescales of 8 to 10 years. Plausible timescales for variable dust extinction are factors of 2-10 too long to explain the dimming and brightening in these sources, and simple dust reddening models cannot reproduce the BEL changes. On the other hand, an advancement such as disk reprocessing is needed if the observed variations are due to accretion rate changes.

The Pan-Planets survey observed an area of 42 sq deg. in the galactic disk for about 165 hours. The main scientific goal of the project is the detection of transiting planets around M dwarfs. We establish an efficient procedure for determining the stellar parameters $T_{eff}$ and log$g$ of all sources using a method based on SED fitting, utilizing a three-dimensional dust map and proper motion information. In this way we identify more than 60000 M dwarfs, which is by far the largest sample of low-mass stars observed in a transit survey to date. We present several planet candidates around M dwarfs and hotter stars that are currently being followed up. Using Monte-Carlo simulations we calculate the detection efficiency of the Pan-Planets survey for different stellar and planetary populations. We expect to find $3.0^{+3.3}_{-1.6}$ hot Jupiters around F, G, and K dwarfs with periods lower than 10 days based on the planet occurrence rates derived in previous surveys. For M dwarfs, the percentage of stars with a hot Jupiter is under debate. Theoretical models expect a lower occurrence rate than for larger main sequence stars. However, radial velocity surveys find upper limits of about 1\% due to their small sample, while the Kepler survey finds a occurrence rate that we estimate to be at least $0.17(^{+0.67}_{-0.04})$%, making it even higher than the determined fraction from OGLE-III for F, G and K stellar types, $0.14(^{+0.15}_{-0.076})\%$. With the large sample size of Pan-Planets, we are able to determine an occurrence rate of $0.11(^{+0.37}_{-0.02})$% in case one of our candidates turns out to be a real detection. If, however, none of our candidates turn out to be true planets, we are able to put an upper limit of 0.34% with a 95% confidence on the hot Jupiter occurrence rate of M dwarfs. Therefore we cannot yet confirm the theoretical prediction of a lower occurrence rate for cool stars.

The Ophiuchus stellar stream presents a dynamical puzzle: its old stellar populations ($\sim 12$ Gyr) cannot be reconciled with (1) its orbit in a simple model for the Milky Way potential and (2) its short angular extent, both of which imply that the observed stream formed within the last $<1$ Gyr. Recent theoretical work has shown that streams on chaotic orbits may abruptly fan out near their apparent ends; stars in these fans are dispersed in both position and velocity and may be difficult to associate with the stream. Here we present the first evidence of such stream-fanning in the Ophiuchus stream, traced by four blue horizontal branch (BHB) stars beyond the apparent end of the stream. These stars stand out from the background by their high velocities ($v_{\rm los} > 230$ km s$^{-1}$) against $\sim 40$ other stars: their velocities are comparable to those of the stream, but would be exceptional if they were unrelated halo stars. Their positions and velocities are, however, inconsistent with simple extrapolation of the observed cold, high-density portion of the stream. These observations suggest that stream-fanning may be a real, observable effect and, therefore, that Ophiuchus may be on a chaotic orbit. They also show that the Ophiuchus stream is more extended and hence dynamically older than previously thought, easing the stellar population vs. dynamical age tension.

The NASA Kepler mission has revolutionised time-domain astronomy and has massively expanded the number of known extrasolar planets. However, the effect of wide multiplicity on exoplanet occurrence has not been tested with this dataset. We present a sample of 401 wide multiple systems containing at least one Kepler target star. Our method uses Pan-STARRS1 and archival data to produce an accurate proper motion catalogue of the Kepler field. Combined with Pan-STARRS1 SED fits and archival proper motions for bright stars, we use a newly developed probabilistic algorithm to identify likely wide binary pairs which are not chance associations. As by-products of this we present stellar SED templates in the Pan-STARRS1 photometric system and conversions from this system to Kepler magnitudes. We find that Kepler target stars in our binary sample with separations above 6 arcseconds are no more or less likely to be identified as confirmed or candidate planet hosts than a weighted comparison sample of Kepler targets of similar brightness and spectral type. Therefore we find no evidence that binaries with projected separations greater than 3,000AU affect the occurrence rate of planets with P<300days around FGK stars.

In area and depth, the Pan-STARRS1 (PS1) 3$\pi$ survey is unique among many-epoch, multi-band surveys and has enormous potential for all-sky identification of variable sources. PS1 has observed the sky typically seven times in each of its five bands ($grizy$) over 3.5 years, but unlike SDSS not simultaneously across the bands. Here we develop a new approach for quantifying statistical properties of non-simultaneous, sparse, multi-color lightcurves through light-curve structure functions, effectively turning PS1 into a $\sim 35$-epoch survey. We use this approach to estimate variability amplitudes and timescales $(\omega_r, \tau)$ for all point-sources brighter than $r_{\mathrm{P1}}=21.5$ mag in the survey. With PS1 data on SDSS Stripe 82 as ``ground truth", we use a Random Forest Classifier to identify QSOs and RR Lyrae based on their variability and their mean PS1 and WISE colors. We find that, aside from the Galactic plane, QSO and RR Lyrae samples of purity $\sim$75\% and completeness $\sim$92\% can be selected. On this basis we have identified a sample of $\sim 1,000,000$ QSO candidates, as well as an unprecedentedly large and deep sample of $\sim$150,000 RR Lyrae candidates with distances from $\sim$10 kpc to $\sim$120 kpc. Within the Draco dwarf spheroidal, we demonstrate a distance precision of 6\% for RR Lyrae candidates. We provide a catalog of all likely variable point sources and likely QSOs in PS1, a total of $25.8\times 10^6$ sources.

We present the discovery of three new Milky Way satellites from our search for compact stellar overdensities in the photometric catalog of the Panoramic Survey Telescope and Rapid Response System 1 (Pan-STARRS 1, or PS1) 3pi survey. The first satellite, Laevens 3, is located at a heliocentric distance of d=67+/-3 kpc. With a total magnitude of Mv=-4.4+/-0.3 and a half-light radius rh=7+/-2 pc, its properties resemble those of outer halo globular clusters. The second system, Draco II/Laevens 4 (Dra II), is a closer and fainter satellite (d~20 kpc, Mv =-2.9+/-0.8), whose uncertain size (rh = 19 +8/-6 pc) renders its classification difficult without kinematic information; it could either be a faint and extended globular cluster or a faint and compact dwarf galaxy. The third satellite, Sagittarius II/Laevens 5 (Sgr II), has an ambiguous nature as it is either the most compact dwarf galaxy or the most extended globular cluster in its luminosity range (rh = 37 +9/-8 pc and Mv=-5.2+/-0.4). At a heliocentric distance of 67+/-5 kpc, this satellite lies intriguingly close to the expected location of the trailing arm of the Sagittarius stellar stream behind the Sagittarius dwarf spheroidal galaxy (Sgr dSph). If confirmed through spectroscopic follow up, this connection would locate this part of the trailing arm of the Sagittarius stellar stream that has so far gone undetected. It would further suggest that Sgr II was brought into the Milky Way halo as a satellite of the Sgr dSph.

[Abbreviated] Upcoming large area sky surveys like EUCLID and eROSITA crucially depend on accurate photometric redshifts (photo-z). The identification of variable sources, such as AGNs, and the achievable redshift accuracy for varying objects are important in view of the science goals of the EUCLID and eROSITA missions. We probe AGN optical variability for a large sample of X-ray-selected AGNs in the XMM-COSMOS field, using the light curves provided by the Pan-STARRS1 (PS1) 3pi and MDF04 surveys. Utilizing two different variability parameters, we defined a sample of varying AGNs for every PS1 band. We investigated the influence of variability on the calculation of photo-z by applying three different input photometry sets for our fitting procedure. For each of the five PS1 bands, we chose either the epochs minimizing the interval in observing time, the median magnitude values, or randomly drawn light curve points to compute the redshift. In addition, we derived photo-z using PS1 photometry extended by GALEX/IRAC bands. We find that the photometry produced by the 3pi survey is sufficient to reliably detect variable sources provided that the fractional variability amplitude is at least 3%. Considering the photo-z of variable AGNs, we observe that minimizing the time spacing of the chosen points yields superior photo-z in terms of the percentage of outliers (33%) and accuracy (0.07), outperforming the other two approaches. Drawing random points from the light curve gives rise to typically 57% of outliers and an accuracy of 0.4. Adding GALEX/IRAC bands for the redshift determination weakens the influence of variability. Although the redshift quality generally improves when adding these bands, we still obtain not less than 26% of outliers and an accuracy of 0.05 at best, therefore variable sources should receive a flag stating that their photo-z may be low quality.

The fortuitous occurrence of a type II-Plateau (IIP) supernova, SN 2014bc, in a galaxy for which distance estimates from a number of primary distance indicators are available provides a means with which to cross-calibrate the standardised candle method (SCM) for type IIP SNe. By applying calibrations from the literature we find distance estimates in line with the most precise measurement to NGC 4258 based on the Keplerian motion of masers (7.6$\pm$0.23\,Mpc), albeit with significant scatter. We provide an alternative local SCM calibration by only considering type IIP SNe that have occurred in galaxies for which a Cepheid distance estimate is available. We find a considerable reduction in scatter ($\sigma_I = 0.16$\, mag.), but note that the current sample size is limited. Applying this calibration, we estimate a distance to NGC 4258 of $7.08\pm0.86$ Mpc.

We present the results of a Monte Carlo technique to calculate the absolute magnitudes (H) and slope parameters (G) of about 240000 asteroids observed by the Pan-STARRS1 telescope during the first 15 months of its 3-year all-sky survey mission. The system's exquisite photometry with photometric errors < 0.04 mag, and well-defined filter and photometric system, allowed us to derive accurate H and G even with a limited number of observations and restricted range in phase angles. Our Monte Carlo method simulates each asteroid's rotation period, amplitude and color to derive the most-likely H and G, but its major advantage is in estimating realistic statistical+systematic uncertainties and errors on each parameter. The method was tested by comparison with the well-established and accurate results for about 500 asteroids provided by Pravec et al. (2012) and then applied to determining H and G for the Pan-STARRS1 asteroids using both the Muinonen et al. (2010) and Bowell et al. (1989) phase functions. Our results confirm the bias in MPC photometry discovered by (Juric et al. 2002).

We present a three-dimensional map of interstellar dust reddening, covering three-quarters of the sky out to a distance of several kiloparsecs, based on Pan-STARRS 1 and 2MASS photometry. The map reveals a wealth of detailed structure, from filaments to large cloud complexes. The map has a hybrid angular resolution, with most of the map at an angular resolution of 3.4' to 13.7', and a maximum distance resolution of ~25%. The three-dimensional distribution of dust is determined in a fully probabilistic framework, yielding the uncertainty in the reddening distribution along each line of sight, as well as stellar distances, reddenings and classifications for 800 million stars detected by Pan-STARRS 1. We demonstrate the consistency of our reddening estimates with those of two-dimensional emission-based maps of dust reddening. In particular, we find agreement with the Planck 353 GHz optical depth-based reddening map to within 0.05 mag in E(B-V) to a depth of 0.5 mag, and explore systematics at reddenings less than E(B-V) ~ 0.08 mag. We validate our per-star reddening estimates by comparison with reddening estimates for stars with both SDSS photometry and SEGUE spectral classifications, finding per-star agreement to within 0.1 mag out to a stellar E(B-V) of 1 mag. We compare our map to two existing three-dimensional dust maps, by Marshall et al. (2006) and Lallement et al. (2013), demonstrating our finer angular resolution, and better distance resolution compared to the former within ~3 kpc. The map can be queried or downloaded at http://argonaut.skymaps.info. We expect the three-dimensional reddening map presented here to find a wide range of uses, among them correcting for reddening and extinction for objects embedded in the plane of the Galaxy, studies of Galactic structure, calibration of future emission-based dust maps and determining distances to objects of known reddening.

We report results of a study of Planck Sunyaev-Zel'dovich effect (SZE) selected galaxy cluster candidates using the Panoramic Survey Telescope & Rapid Response System (Pan-STARRS) imaging data. We first examine 150 Planck confirmed galaxy clusters with spectroscopic redshifts to test our algorithm for identifying optical counterparts and measuring their redshifts; our redshifts have a typical accuracy of $\sigma_{z/(1+z)} \sim 0.022$ for this sample. Using 60 random sky locations, we estimate that our chance of contamination through a random superposition is ~ 3 per cent. We then examine an additional 237 Planck galaxy cluster candidates that have no redshift in the source catalogue. Of these 237 unconfirmed cluster candidates we are able to confirm 60 galaxy clusters and measure their redshifts. A further 83 candidates are so heavily contaminated by stars due to their location near the Galactic plane that we do not attempt to identify counterparts. For the remaining 94 candidates we find no optical counterpart but use the depth of the Pan-STARRS1 data to estimate a redshift lower limit $z_{\text{lim}(10^{15})}$ beyond which we would not have expected to detect enough galaxies for confirmation. Scaling from the already published Planck sample, we expect that $\sim$12 of these unconfirmed candidates may be real clusters.

Supermassive black hole binaries (SMBHBs) should be an inevitable consequence of the hierarchical growth of massive galaxies through mergers, and the strongest sirens of gravitational waves (GWs) in the cosmos. And yet, their direct detection has remained elusive due to the compact (sub-parsec) orbital separations of gravitationally bound SMBHBs. Here we exploit a theoretically predicted signature of a SMBHB in the time domain: periodic variability caused by a mass accretion rate that is modulated by the binary's orbital motion. We report our first significant periodically varying quasar detection from the systematic search in the Pan-STARRS1 (PS1) Medium Deep Survey. Our SMBHB candidate, PSO J334.2028+01.4075, is a luminous radio-loud quasar at $z=2.060$, with extended baseline photometry from the Catalina Real-Time Transient Survey, as well as archival spectroscopy from the FIRST Bright Quasar Survey. The observed period ($542 \pm 15$ days) and estimated black hole mass ($\log (M_{\rm BH}/M_\odot) = 9.97 \pm 0.50$), correspond to an orbital separation of $7^{+8}_{-4}$ Schwarzschild radii ($\sim 0.006^{+0.007}_{-0.003}$ pc), assuming the rest-frame period of the quasar variability traces the orbital period of the binary. This SMBHB candidate, discovered at the peak redshift for SMBH mergers, is in a physically stable configuration for a circumbinary accretion disk, and within the regime of GW-driven orbital decay. Our search with PS1 is a benchmark study for the exciting capabilities of LSST, which will have orders of magnitude larger survey power, and will potentially pinpoint the locations of thousands of SMBHBs in the variable night sky.

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

We present the discovery of a faint Milky Way satellite, Laevens 2/Triangulum II, found in the Panoramic Survey Telescope And Rapid Response System (Pan-STARRS 1) 3 pi imaging data and confirmed with follow-up wide-field photometry from the Large Binocular Cameras. The stellar system, with an absolute magnitude of M_V=-1.8 +/-0.5, a heliocentric distance of 30 +2/-2 kpc, and a half-mass radius of 34 +9/-8 pc, shows remarkable similarity to faint, nearby, small satellites such as Willman 1, Segue 1, Segue 2, and Bo\"otes II. The discovery of Laevens 2/Triangulum II further populates the region of parameter space for which the boundary between dwarf galaxies and globular clusters becomes tenuous. Follow-up spectroscopy will ultimately determine the nature of this new satellite, whose spatial location hints at a possible connection with the complex Triangulum-Andromeda stellar structures.

Radio-loud Active Galactic Nuclei at z~2-4 are typically located in dense environments and their host galaxies are among the most massive systems at those redshifts, providing key insights for galaxy evolution. Finding radio-loud quasars at the highest accessible redshifts (z~6) is important to study their properties and environments at even earlier cosmic time. They would also serve as background sources for radio surveys intended to study the intergalactic medium beyond the epoch of reionization in HI 21 cm absorption. Currently, only five radio-loud ($R=f_{\nu,5{\rm GHz}}/f_{\nu,4400\AA}>10$) quasars are known at z~6. In this paper we search for 5.5 < z < 7.2 quasars by cross-matching the optical Pan-STARRS1 and radio FIRST surveys. The radio information allows identification of quasars missed by typical color-based selections. While we find no good 6.4 < z <7.2 quasar candidates at the sensitivities of these surveys, we discover two new radio-loud quasars at z~6. Furthermore, we identify two additional z~6 radio-loud quasars which were not previously known to be radio-loud, nearly doubling the current z~6 sample. We show the importance of having infrared photometry for z>5.5 quasars to robustly classify them as radio-quiet or radio-loud. Based on this, we reclassify the quasar J0203+0012 (z=5.72), previously considered radio-loud, to be radio-quiet. Using the available data in the literature, we constrain the radio-loud fraction of quasars at z~6, using the Kaplan--Meier estimator, to be $8.1^{+5.0}_{-3.2}\%$. This result is consistent with there being no evolution of the radio-loud fraction with redshift, in contrast to what has been suggested by some studies at lower redshifts.

We use the WISE-2MASS infrared galaxy catalog matched with Pan-STARRS1 (PS1) galaxies to search for a supervoid in the direction of the Cosmic Microwave Background Cold Spot. Our imaging catalog has median redshift $z\simeq 0.14$, and we obtain photometric redshifts from PS1 optical colours to create a tomographic map of the galaxy distribution. The radial profile centred on the Cold Spot shows a large low density region, extending over 10's of degrees. Motivated by previous Cosmic Microwave Background results, we test for underdensities within two angular radii, $5^\circ$, and $15^\circ$. The counts in photometric redshift bins show significantly low densities at high detection significance, $\gtrsim 5 \sigma$ and $\gtrsim 6 \sigma$, respectively, for the two fiducial radii. The line-of-sight position of the deepest region of the void is $z\simeq 0.15-0.25$. Our data, combined with an earlier measurement by Granett et al. 2010, are consistent with a large $R_{\rm void}=(220 \pm 50) h^{-1}Mpc $ supervoid with $\delta_{m} \simeq -0.14 \pm 0.04$ centered at $z=0.22\pm0.03$. Such a supervoid, constituting at least a $\simeq 3.3\sigma$ fluctuation in a Gaussian distribution of the $\Lambda CDM$ model, is a plausible cause for the Cold Spot.

Hypervelocity stars (HVS) travel with velocities so high, that they exceed the escape velocity of the Galaxy. Several acceleration mechanisms have been discussed. Only one HVS (US 708, HVS 2) is a compact helium star. Here we present a spectroscopic and kinematic analysis of US\,708. Travelling with a velocity of $\sim1200\,{\rm km\,s^{-1}}$, it is the fastest unbound star in our Galaxy. In reconstructing its trajectory, the Galactic center becomes very unlikely as an origin, which is hardly consistent with the most favored ejection mechanism for the other HVS. Furthermore, we discovered US\,708 to be a fast rotator. According to our binary evolution model it was spun-up by tidal interaction in a close binary and is likely to be the ejected donor remnant of a thermonuclear supernova.

Efficient identification and follow-up of astronomical transients is hindered by the need for humans to manually select promising candidates from data streams that contain many false positives. These artefacts arise in the difference images that are produced by most major ground-based time domain surveys with large format CCD cameras. This dependence on humans to reject bogus detections is unsustainable for next generation all-sky surveys and significant effort is now being invested to solve the problem computationally. In this paper we explore a simple machine learning approach to real-bogus classification by constructing a training set from the image data of ~32000 real astrophysical transients and bogus detections from the Pan-STARRS1 Medium Deep Survey. We derive our feature representation from the pixel intensity values of a 20x20 pixel stamp around the centre of the candidates. This differs from previous work in that it works directly on the pixels rather than catalogued domain knowledge for feature design or selection. Three machine learning algorithms are trained (artificial neural networks, support vector machines and random forests) and their performances are tested on a held-out subset of 25% of the training data. We find the best results from the random forest classifier and demonstrate that by accepting a false positive rate of 1%, the classifier initially suggests a missed detection rate of around 10%. However we also find that a combination of bright star variability, nuclear transients and uncertainty in human labelling means that our best estimate of the missed detection rate is approximately 6%.

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

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

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

We present an observational and dynamical study of newly discovered main-belt comet 313P/Gibbs. We find that the object is clearly active both in observations obtained in 2014 and in precovery observations obtained in 2003 by the Sloan Digital Sky Survey, strongly suggesting that its activity is sublimation-driven. This conclusion is supported by a photometric analysis showing an increase in the total brightness of the comet over the 2014 observing period, and dust modeling results showing that the dust emission persists over at least three months during both active periods, where we find start dates for emission no later than 2003 July 24+/-10 for the 2003 active period and 2014 July 28+/-10 for the 2014 active period. From serendipitous observations by the Subaru Telescope in 2004 when the object was apparently inactive, we estimate that the nucleus has an absolute R-band magnitude of H_R=17.1+/-0.3, corresponding to an effective nucleus radius of r_e~1.00+/-0.15 km. The object's faintness at that time means we cannot rule out the presence of activity, and so this computed radius should be considered an upper limit. We find that 313P's orbit is intrinsically chaotic, having a Lyapunov time of T_l=12000 yr and being located near two 3-body mean-motion resonances with Jupiter and Saturn, 11J-1S-5A and 10J+12S-7A, yet appears stable over >50 Myr in an apparent example of stable chaos. We furthermore find that 313P is the second main-belt comet, after P/2012 T1 (PANSTARRS), to belong to the ~155 Myr old Lixiaohua asteroid family.

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 present the results of an analysis of Pan-STARRS1 Medium Deep Survey multi-band (g,r,i,z,y) images of a sample of 698 low-redshift disk galaxies that span broad ranges in stellar mass, star-formation rate, and bulge/disk ratio. We use population synthesis SED fitting techniques to explore the radial distribution of the light, color, surface mass density, mass/light ratio, and age of the stellar populations. We characterize the structure and stellar content of the galaxy disks out to radii of about twice Petrosian r90, beyond which the halo light becomes significant. We measure normalized radial profiles for sub-samples of galaxies in three bins each of stellar mass and concentration. We also fit radial profiles to each galaxy. The majority of galaxies have down-bending radial surface brightness profiles in the bluer bands with a break radius at roughly r90. However, they typically show single unbroken exponentials in the reddest bands and in the stellar surface mass density. We find that the mass/light ratio and stellar age radial profiles have a characteristic "U-shape". There is a good correlation between the amplitude of the down-bend in the surface brightness profile and the rate of the increase in the M/L ratio in the outer disk. As we move from late- to early-type galaxies, the amplitude of the down-bend and the radial gradient in M/L both decrease. Our results imply a combination of stellar radial migration and suppression of recent star formation can account for the stellar populations of the outer disk.

The Orion Molecular Complex is the nearest site of ongoing high-mass star formation, making it one of the most extensively studied molecular complexes in the Galaxy. We have developed a new technique for mapping the 3D distribution of dust in the Galaxy using Pan-STARRS1 photometry. We isolate the dust at the distance to Orion using this technique, revealing a large (100 pc, 14 degree diameter), previously unrecognized ring of dust, which we term the "Orion dust ring." The ring includes Orion A and B, and is not coincident with current H-alpha features. The circular morphology suggests formation as an ancient bubble in the interstellar medium, though we have not been able to conclusively identify the source of the bubble. This hint at the history of Orion may have important consequences for models of high-mass star formation and triggered star formation.

We analyze 760475 observations of 333026 main-belt objects obtained by the Pan-STARRS1 (PS1) survey telescope between 2012 May 20 and 2013 November 9, a period during which PS1 discovered two main-belt comets, P/2012 T1 (PANSTARRS) and P/2013 R3 (Catalina-PANSTARRS). PS1 comet detection procedures currently consist of the comparison of the point spread functions (PSFs) of moving objects to those of reference stars, and the flagging of objects that show anomalously large radial PSF widths. Based on the number of missed discovery opportunities among comets discovered by other observers, we estimate an upper limit comet discovery efficiency rate of ~70% for PS1. Additional analyses that could improve comet discovery yields in future surveys include linear PSF analysis, modeling of trailed stellar PSFs for comparison to trailed moving object PSFs, searches for azimuthally localized activity, comparison of point-source-optimized photometry to extended-source-optimized photometry, searches for photometric excesses in objects with known absolute magnitudes, and crowd-sourcing. Analysis of PS1 survey statistics indicates an expected fraction of 59 MBCs per 10^6 outer main-belt asteroids, and a 95% confidence upper limit of 96 MBCs per 10^6 outer main-belt asteroids. We note that more sensitive future surveys could detect many more MBCs than estimated here. We find an excess of high eccentricities (0.1 < e < 0.3) among all known MBCs relative to the background asteroid population. Theoretical calculations show that, given these eccentricities, the sublimation rate for a typical MBC is orders of magnitude larger at perihelion than at aphelion, providing a plausible physical explanation for the observed behavior of MBCs peaking in observed activity strength near perihelion.

We present a search for eclipses of $\sim$1700 white dwarfs in the Pan-STARRS1 medium-deep fields. Candidate eclipse events are selected by identifying low outliers in over 4.3 million light curve measurements. We find no short-duration eclipses consistent with being caused by a planetary size companion. This large dataset enables us to place strong constraints on the close-in planet occurrence rates around white dwarfs for planets as small as 2 R$_\oplus$. Our results indicate that gas giant planets orbiting just outside the Roche limit are rare, occurring around less than 0.5% of white dwarfs. Habitable-zone super-Earths and hot super-Earths are less abundant than similar classes of planets around main-sequence stars. These constraints give important insight into the ultimate fate of the large population of exoplanets orbiting main sequence stars.

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

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

We have calculated 90% confidence limits on the steady-state rate of catastrophic disruptions of main belt asteroids in terms of the absolute magnitude at which one catastrophic disruption occurs per year (HCL) as a function of the post-disruption increase in brightness (delta m) and subsequent brightness decay rate (tau). The confidence limits were calculated using the brightest unknown main belt asteroid (V = 18.5) detected with the Pan-STARRS1 (Pan-STARRS1) telescope. We measured the Pan-STARRS1's catastrophic disruption detection efficiency over a 453-day interval using the Pan-STARRS moving object processing system (MOPS) and a simple model for the catastrophic disruption event's photometric behavior in a small aperture centered on the catastrophic disruption event. Our simplistic catastrophic disruption model suggests that delta m = 20 mag and 0.01 mag d-1 < tau < 0.1 mag d-1 which would imply that H0 = 28 -- strongly inconsistent with H0,B2005 = 23.26 +/- 0.02 predicted by Bottke et al. (2005) using purely collisional models. We postulate that the solution to the discrepancy is that > 99% of main belt catastrophic disruptions in the size range to which this study was sensitive (100 m) are not impact-generated, but are instead due to fainter rotational breakups, of which the recent discoveries of disrupted asteroids P/2013 P5 and P/2013 R3 are probable examples. We estimate that current and upcoming asteroid surveys may discover up to 10 catastrophic disruptions/year brighter than V = 18.5.

We present Keck II/DEIMOS spectroscopy of the three distant dwarf galaxies of M31 Lacerta I, Cassiopeia III, and Perseus I, recently discovered within the Pan-STARRS1 3\pi imaging survey. The systemic velocities of the three systems (v_{r,helio} = -198.4 +/- 1.1 km/s, -371.6 +/- 0.7 km/s, and -326 +/- 3 km/s, respectively) confirm that they are satellites of M31. In the case of Lacerta I and Cassiopeia III, the high quality of the data obtained for 126 and 212 member stars, respectively, yields reliable constraints on their global velocity dispersions (\sigma_{vr} = 10.3 +/- 0.9 km/s and 8.4 +/- 0.6 km/s, respectively), leading to dynamical-mass estimates for both of ~4x10^7 Msun within their half-light radius. These translate to V-band mass-to-light ratios of 15^{+12}_{-9} and 8^{+9}_{-5} in solar units. We also use our spectroscopic data to determine the average metallicity of the 3 dwarf galaxies ([Fe/H] = -2.0 +/- 0.1, -1.7 +/- 0.1, and -2.0 +/- 0.2, respectively). All these properties are typical of dwarf galaxy satellites of Andromeda with their luminosity and size.

We present the results of a search for the reactivation of active asteroid 176P/LINEAR during its 2011 perihelion passage using deep optical observations obtained before, during, and after that perihelion passage. Deep composite images of 176P constructed from data obtained between June 2011 and December 2011 show no visible signs of activity, while photometric measurements of the object during this period also show no significant brightness enhancements similar to that observed for 176P between November 2005 and December 2005 when it was previously observed to be active. An azimuthal search for dust emission likewise reveals no evidence for directed emission (i.e., a tail, as was previously observed for 176P), while a one-dimensional surface brightness profile analysis shows no indication of a spherically symmetric coma at any time in 2011. We conclude that 176P did not in fact exhibit activity in 2011, at least not on the level on which it exhibited activity in 2005, and suggest that this could be due to the devolatization or mantling of the active site responsible for its activity in 2005.

We present the discovery of 61 wide (>5 arcsecond) separation, low-mass (stellar and substellar) companions to stars in the solar neighborhood identified from Pan-STARRS\,1 (PS1) data and the spectral classification of 27 previously known companions. Our companions represent a selective subsample of promising candidates and span a range in spectral type of K7-L9 with the addition of one DA white dwarf. These were identified primarily from a dedicated common proper motion search around nearby stars, along with a few as serendipitous discoveries from our Pan-STARRS1 brown dwarf search. Our discoveries include 24 new L dwarf companions and one known L dwarf not previously identified as a companion. The primary stars around which we searched for companions come from a list of bright stars with well-measured parallaxes and large proper motions from the Hipparcos catalog (8583 stars, mostly A-K~dwarfs) and fainter stars from other proper motion catalogues (79170 stars, mostly M~dwarfs). We examine the likelihood that our companions are chance alignments between unrelated stars and conclude that this is unlikely for the majority of the objects that we have followed-up spectroscopically. We also examine the entire population of ultracool (>M7) dwarf companions and conclude that while some are loosely bound, most are unlikely to be disrupted over the course of $\sim$10 Gyr. Our search increases the number of ultracool M dwarf companions wider than 300 AU by 88% and increases the number of L dwarf companions in the same separation range by 96%. Finally, we resolve our new L dwarf companion to HIP 6407 into a tight (0.13 arcsecond, 7.4 AU) L1+T3 binary, making the system a hierarchical triple. Our search for these key benchmarks against which brown dwarf and exoplanet atmosphere models are tested has yielded the largest number of discoveries to date.

We use a WISE-2MASS-Pan-STARRS1 galaxy catalog to search for a supervoid in the direction of the Cosmic Microwave Background Cold Spot. We obtain photometric redshifts using our multicolor data set to create a tomographic map of the galaxy distribution. The radial density profile centred on the Cold Spot shows a large low density region, extending over 10's of degrees. Motivated by previous Cosmic Microwave Background results, we test for underdensities within two angular radii, $5^\circ$, and $15^\circ$. Our data, combined with an earlier measurement by Granett et al 2010, are consistent with a large $R_{\rm void}=(192 \pm 15)h^{-1} Mpc $ $(2\sigma)$ supervoid with $\delta \simeq -0.13 \pm 0.03$ centered at $z=0.22\pm0.01$. Such a supervoid, constituting a $\sim3.5 \sigma$ fluctuation in the $\Lambda CDM$ model, is a plausible cause for the Cold Spot.

We present a panoptic view of the stellar structure in the Galactic disk's outer reaches commonly known as the Monoceros Ring, based on data from Pan-STARRS1. These observations clearly show the large extent of the stellar overdensities on both sides of the Galactic disk, extending between b = -25 and b = +35 degrees and covering over 130 degrees in Galactic longitude. The structure exhibits a complex morphology with both stream-like features and a sharp edge to the structure in both the north and the south. We compare this map to mock observations of two published simulations aimed at explaining such structures in the outer stellar disk, one postulating an origin as a tidal stream and the other demonstrating a scenario where the disk is strongly distorted by the accretion of a satellite. These morphological comparisons of simulations can link formation scenarios to observed structures, such as demonstrating that the distorted-disk model can produce thin density features resembling tidal streams. Although neither model produces perfect agreement with the observations--the tidal stream predicts material at larger distances which is not detected while in the distorted disk model the midplane is warped to an excessive degree--future tuning of the models to accommodate these latest data may yield better agreement.

We present a study on the discoverability of temporarily captured orbiters (TCOs) by present day or near-term anticipated ground-based and space-based facilities. TCOs (Granvik et al. 2012) are potential targets for spacecraft ren- dezvous or human exploration (Chyba et al. 2014) and provide an opportunity to study the population of the smallest asteroids in the solar system. We find that present day ground-based optical surveys such as Pan-STARRS and ATLAS can discover the largest TCOs over years of operation. A targeted survey conducted with the Subaru telescope can discover TCOs in the 0.5 m to 1.0 m diameter size range in about 5 nights of observing. Furthermore, we discuss the application of space-based infrared surveys, such as NEOWISE, and ground-based meteor detection systems such as CAMS, CAMO and ASGARD in discovering TCOs. These systems can detect TCOs but at a uninteresting rate. Finally, we discuss the application of bi-static radar at Arecibo and Green Bank to discover TCOs. Our radar simulations are strongly dependent on the rotation rate distribution of the smallest asteroids but with an optimistic distribution we find that these systems have > 80% chance of detecting a > 10 cm diameter TCO in about 40 h of operation.

Standard inflationary hot big bang cosmology predicts small fluctuations in the Cosmic Microwave Background (CMB) with isotropic Gaussian statistics. All measurements support the standard theory, except for a few anomalies discovered in the Wilkinson Microwave Anisotropy Probe maps and confirmed recently by the Planck satellite. The Cold Spot is one of the most significant of such anomalies, and the leading explanation of it posits a large void that imprints this extremely cold area via the linear Integrated Sachs-Wolfe (ISW) effect due to the decay of gravitational potentials over cosmic time, or via the Rees-Sciama (RS) effect due to late-time non-linear evolution. Despite several observational campaigns targeting the Cold Spot region, to date no suitably large void was found at higher redshifts $z > 0.3$. Here we report the detection of an $R =(192 \pm 15) h^{-1}Mpc$ size supervoid of depth $\delta = -0.13 \pm 0.03$, and centred at redshift $z = 0.22$. This supervoid, possibly the largest ever found, is large enough to significantly affect the CMB via the non-linear RS effect, as shown in our Lemaitre-Tolman-Bondi framework. This discovery presents the first plausible explanation for any of the physical CMB anomalies, and raises the possibility that local large-scale structure could be responsible for other anomalies as well.

We report the discovery of a thin stellar stream found in Pan-STARRS1 photometry near the Galactic bulge in the constellation of Ophiuchus. It appears as a coherent structure in the colour-selected stellar density maps produced to search for tidal debris around nearby globular clusters. The stream is exceptionally short and narrow; it is about 2.5{\deg} long and 6' wide in projection. The colour-magnitude diagram of this object, which harbours a blue horizontal-branch, is consistent with an old and relatively metal-poor population ([Fe/H]~-1.3) located 9.5 +/- 0.9 kpc away at (l,b) ~ (5{\deg},+32{\deg}), and 5.0 +/- 1.0 kpc from the Galactic centre. These properties argue for a globular cluster as progenitor. The finding of such a prominent, nearby stream suggests that many streams could await discovery in the more densely populated regions of our Galaxy.

We present the fiducial sequences of a sample of Galactic star clusters in the five bands of the Pan-STARRS1 (PS1) photometric system (g_P1, r_P1, i_P1, z_P1, y_P1). These empirical sequences -- which include the red giant and sub-giant branches, the main sequence, and the horizontal branch -- were defined from deep colour-magnitude diagrams reaching below the oldest main-sequence turn-offs of 13 globular and 3 old open clusters covering a wide range of metallicities (-2.4 < [Fe/H] < +0.4). We find excellent agreement for the nine clusters in common with previous studies in similar photometric systems when transformed to the PS1 system. Because the photometric and spectroscopic properties of these stellar populations are accurately known, the fiducials provide a solid basis for the interpretation of observations in the PS1 system, as well as valuable constraints to improve the empirical colour--$T_{eff}$ relations.

We have calculated the coherence and detectable lifetimes of synthetic near-Earth object (NEO) families created by catastrophic disruption of a progenitor as it suffers a very close Earth approach. The closest or slowest approaches yield the most violent `s-class' disruption events. We found that the average slope of the absolute magnitude (H) distribution, $N(H)\propto10^{(0.55\pm0.04)\,H}$, for the fragments in the s-class families is steeper than the slope of the NEO population \citep{mainzer2011} in the same size range. The families remain coherent as statistically significant clusters of orbits within the NEO population for an average of $\bar\tau_c = (14.7\pm0.6)\times10^3$ years after disruption. The s-class families are detectable with the techniques developed by \citet{fu2005} and \citet{Schunova2012} for an average duration ($\bar\tau_{det}$) ranging from about 2,000 to about 12,000 years for progenitors in the absolute magnitude ($H_p$) range from 20 to 13 corresponding to diameters in the range from about 0.5 to 10$\km$ respectively. The short detectability lifetime explains why zero NEO families have been discovered to-date. Nonetheless, every tidal disruption event of a progenitor with D$>0.5\km$ is capable of producing several million fragments in the $1\meter$ to $10\meter$ diameter range that can contribute to temporary local density enhancements of small NEOs in Earth's vicinity. We expect that there are about 1,200 objects in the steady state NEO population in this size range due to tidal disruption assuming that one $1\km$ diameter NEO tidally disrupts at Earth every 2,500 years. These objects may be suitable targets for asteroid retrieval missions due to their Earth-like orbits with corresponding low $v_{\infty}$. The fragments from the tidal disruptions at Earth have $\sim5\times$ the collision probability with Earth compared to the background NEO population.

High-redshift quasars are currently the only probes of the growth of supermassive black holes and potential tracers of structure evolution at early cosmic time. Here we present our candidate selection criteria from the Panoramic Survey Telescope & Rapid Response System 1 and follow-up strategy to discover quasars in the redshift range 5.7<z<6.2. With this strategy we discovered eight new 5.7<z<6.0 quasars, increasing the number of known quasars at z>5.7 by more than 10%. We additionally recovered 18 previously known quasars. The eight quasars presented here span a large range of luminosities (-27.3 < M_{1450} < -25.4; 19.6 < z_ps1 < 21.2) and are remarkably heterogeneous in their spectral features: half of them show bright emission lines whereas the other half show a weak or no Ly$\alpha$ emission line (25% with rest-frame equivalent width of the Ly$\alpha$ + Nv line lower than 15{\AA}). We find a larger fraction of weak-line emission quasars than in lower redshift studies. This may imply that the weak-line quasar population at the highest redshifts could be more abundant than previously thought. However, larger samples of quasars are needed to increase the statistical significance of this finding.

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 a map of the dust reddening to 4.5 kpc derived from Pan-STARRS1 stellar photometry. The map covers almost the entire sky north of declination -30 degrees at a resolution of 7' to 14', and is based on the estimated distances and reddenings to more than 500 million stars. The technique is designed to map dust in the Galactic plane, where many other techniques are stymied by the presence of multiple dust clouds at different distances along each line of sight. This reddening-based dust map agrees closely with the Schlegel, Finkbeiner, and Davis (SFD; 1998) far-infrared emission-based dust map away from the Galactic plane, and the most prominent differences between the two maps stem from known limitations of SFD in the plane. We also compare the map with Planck, finding likewise good agreement in general at high latitudes. The use of optical data from Pan-STARRS1 yields reddening uncertainty as low as 25 mmag E(B-V).

We present the detection of 6,371 RR Lyrae (RRL) stars distributed across ~14,000 deg^2 of the sky from the combined data of the Sloan Digital Sky Survey (SDSS), the Panoramic Survey Telescope and Rapid Response System 1 (PS1), and the second photometric catalogue from the Catalina Survey (CSDR2), out of these, ~2,021 RRL stars (~572 RRab and 1,449 RRc) are new discoveries. The RRL stars have heliocentric distances in the 4--28 kpc distance range. RRL-like color cuts from the SDSS and variability cuts from the PS1 are used to cull our candidate list. We then use the CSDR2 multi-epoch data to refine our sample. Periods were measured using the Analysis of Variance technique while the classification process is performed with the Template Fitting Method in addition to the visual inspection of the light curves. A cross-match of our RRL star discoveries with previous published catalogs of RRL stars yield completeness levels of ~50% for both RRab and RRc stars, and an efficiency of ~99% and ~87% for RRab and RRc stars, respectively. We show that our method for selecting RRL stars allows us to recover halo structures. The full lists of all the RRL stars are made publicly available.

We present a new satellite in the outer halo of the Galaxy, the first Milky Way satellite found in the stacked photometric catalog of the Panoramic Survey Telescope and Rapid Response System 1 (Pan-STARRS1) Survey. From follow-up photometry obtained with WFI on the MPG/ESO 2.2m telescope, we argue that the object, located at a heliocentric distance of 145+/-17 kpc, is the most distant Milky Way globular cluster yet known. With a total magnitude of M_V=-4.3+/-0.2 and a half-light radius of 20+/-2 pc, it shares the properties of extended globular clusters found in the outer halo of our Galaxy and the Andromeda galaxy. The discovery of this distant cluster shows that the full spatial extent of the Milky Way globular cluster system has not yet been fully explored.

Dynamically new comet C/2011 L4 (PanSTARRS) is one of the brightest comets since the great comet C/1995 O1 (Hale-Bopp). Here, we present our multi-wavelength observations of C/2011 L4 during its in-bound passage to the inner Solar system. A strong absorption band of water ice at 2.0 $\mu$m was detected in the near infrared spectra, taken with the 8-m Gemini-North and 3-m IRTF telescopes. The companion 1.5 $\mu$m band of water ice, however, was not observed. Spectral modeling show that the absence of the 1.5 $\mu$m feature can be explained by the presence of sub-micron-sized fine ice grains. No gas lines (i.e. CN, HCN or CO) were observed pre-perihelion either in the optical or in the sub-millimeter. 3-$\sigma$ upper limits to the CN and CO production rates were derived. The comet exhibited a very strong continuum in the optical and its slope seemed to become redder as the comet approached the Sun. Our observations suggest that C/2011 L4 is an unusually dust-rich comet with a dust-to-gas mass ratio $>$ 4.

Membership identification is the first step to determine the properties of a star cluster. Low-mass members in particular could be used to trace the dynamical history, such as mass segregation, stellar evaporation, or tidal stripping, of a star cluster in its Galactic environment. We identified member candidates with stellar masses $\sim$0.11--2.4 Msun of the intermediate-age Praesepe cluster (M44), by using Pan-STARRS and 2MASS photometry, and PPMXL proper motions. Within a sky area of 3 deg radius, 1040 candidates are identified, of which 96 are new inclusions. Using the same set of selection criteria on field stars led to an estimate of a false positive rate 16%, suggesting 872 of the candidates being true members. This most complete and reliable membership list allows us to favor the BT-Settl model in comparison with other stellar models. The cluster shows a distinct binary track above the main sequence, with a binary frequency of 20--40%, and a high occurrence rate of similar mass pairs. The mass function is consistent with that of the disk population but shows a deficit of members below 0.3 solar masses. A clear mass segregation is evidenced, with the lowest-mass members in our sample being evaporated from this disintegrating cluster.

We propose a low-cost robotic optical survey aimed at $1-300$ m Near Earth Objects (NEO) based on four state-of-the-art telescopes having extremely wide field of view. The small Near-Earth Asteroids (NEA) represent a potential risk but also easily accessible space resources for future robotic or human space in-situ exploration, or commercial activities. The survey system will be optimized for the detection of fast moving - trailed - asteroids, space debris and will provide real-time alert notifications. The expected cost of the system including 1-year development and 2-year operation is 1,000,000 EUR. The successful demonstration of the system will promote cost-efficient ADAM-WFS (Automatic Detection of Asteroids and Meteoroids - A Wide Field Survey) systems to be built around the world.

P/2013 P5 PANSTARRS was discovered in Aug. 2013, displaying a cometary tail, but with orbital elements typical for a member of the inner asteroid Main Belt. We monitored the object from 2013 Aug. 30 until Oct. 05 using the CFHT, NTT, CA 1.23m, Perkins 1.8m (Lowell), and the 0.6m TRAPPIST telescopes. We measured its nuclear radius to be r < 0.25-0.29km, and its colours g-r = 0.58+/-0.05 and r-i = 0.23+/-0.06, typical for an S-class asteroid. We failed to detect any rotational light curve, with an amplitude < 0.05mag and a double-peaked rotation period < 20h. A detailed Finson-Probstein analysis of deep NTT and CFHT images indicated that the object was active since at least late January 2013 until the time of the latest observations in 2013 September, with at least two peaks of activity around 2013 June 14+/-10d and 2013 July 22+/-3d. The changes of activity level and the activity peaks were extremely sharp and short, shorter than the temporal resolution of our observations (about 1d). The dust distribution was similar during these two events, with dust grains covering at least the 1-1000{\mu}m range. The total mass ejected in grains <1mm was estimated to be 3.0 10$^6$kg and 2.6 10$^7$kg around the two activity peaks. Rotational disruption cannot be ruled out as the cause of the dust ejection. We also propose that the components of a contact binary might gently rub and produce the observed emission. Volatile sublimation might also explain what appears as cometary activity over a period of 8 months. However, while Main Belt comets best explained by ice sublimation are found in the outskirts of the Main Belt, where water ice is believed to be able to survive buried in moderately large objects for the age of the solar system deeply, the presence of volatiles in an object smaller than 300m in radius would be very surprising in the inner asteroid belt.

We present a method to infer reddenings and distances to stars, based only on their broad-band photometry, and show how this method can be used to produce a three-dimensional dust map of the Galaxy. Our method samples from the full probability density function of distance, reddening and stellar type for individual stars, as well as the full uncertainty in reddening as a function of distance in the 3D dust map. We incorporate prior knowledge of the distribution of stars in the Galaxy and the detection limits of the survey. For stars in the Pan-STARRS 1 (PS1) 3 pi survey, we demonstrate that our reddening estimates are unbiased, and accurate to ~0.13 mag in E(B-V) for the typical star. Based on comparisons with mock catalogs, we expect distances for main-sequence stars to be constrained to within ~20% - 60%, although this range can vary, depending on the reddening of the star, the precise stellar type and its position on the sky. A further paper will present a 3D map of dust over the three quarters of the sky surveyed by PS1. Both the individual stellar inferences and the 3D dust map will enable a wealth of Galactic science in the plane. The method we present is not limited to the passbands of the PS1 survey, but may be extended to incorporate photometry from other surveys, such as 2MASS, SDSS (where available), and in the future, LSST and Gaia.

Using a large sample of field and group galaxies drawn from the Pan-STARRS1 Medium-Deep Survey, we present the specific star formation rate (SSFR) - stellar mass (M*) relation, as well as the quiescent fraction versus M* relation in different environments. We confirm that the fraction of quiescent galaxies is strongly dependent on environment at a fixed stellar mass, but that the amplitude and the slope of the star-forming sequence is similar between the field and groups: in other words, the SSFR-density relation at a fixed stellar mass is primarily driven by the change in the star-forming and quiescent fractions between different environments rather than a global suppression in the star formation rate for the star-forming population. However, when we restrict our sample to the cluster-scale environments ($M>10^{14}M_{solar}$), we find a global reduction in the SSFR of the star forming sequence of $17\%$ at 4$\sigma$ confidence as opposed to its field counterpart. After removing the stellar mass dependence of the quiescent fraction seen in field galaxies, the excess in the quiescent fraction due to the environment quenching in groups and clusters is found to increase with stellar mass. We argue that these results are in favor of galaxy mergers to be the primary environment quenching mechanism operating in galaxy groups whereas strangulation is able to reproduce the observed trend in the environment quenching efficiency and stellar mass relation seen in clusters. Our results also suggest that the relative importance between mass quenching and environment quenching depends on stellar mass -- the mass quenching plays a dominant role in producing quiescent galaxies for more massive galaxies, while less massive galaxies are quenched mostly through the environmental effect, with the transition mass around $1-2\times10^{10}M_{solar}$ in the group/cluster environment. (abridged)

Using imaging from the Pan-STARRS1 survey, we identify a precursor outburst at epochs 287 and 170 days prior to the reported explosion of the purported Type IIn supernova (SN) 2011ht. In the Pan-STARRS data, a source coincident with SN 2011ht is detected exclusively in the \zps\ and \yps-bands. An absolute magnitude of M$_z\simeq$-11.8 suggests that this was an outburst of the progenitor star. Unfiltered, archival Catalina Real Time Transient survey images also reveal a coincident source from at least 258 to 138 days before the main event. We suggest that the outburst is likely to be an intrinsically red eruption, although we cannot conclusively exclude a series of erratic outbursts which were observed only in the redder bands by chance. This is only the fourth detection of an outburst prior to a claimed SN, and lends credence to the possibility that many more interacting transients have pre-explosion outbursts, which have been missed by current surveys.

We have used Pan-STARRS1 to discover an extremely red late-L dwarf, which has (J-K)_MKO = 2.84 and (J-K)_2MASS = 2.78, making it the reddest known field dwarf and second only to 2MASS J1207-39b among substellar companions. Near-IR spectroscopy shows a spectral type of L7 and reveals a triangular H-band continuum and weak alkali (K I and Na I) lines, hallmarks of low surface gravity. Near-IR astrometry from the Hawaii Infrared Parallax Program gives a distance of 24.6+/-1.4 pc and indicates a much fainter J-band absolute magnitude than field L dwarfs. The position and kinematics of PSO J318-22 point to membership in the beta Pictoris moving group. Evolutionary models give a temperature of 1160 (-40,+30) K and a mass of 6.5 (-1.0, +1.3) Mjup, making PSO J318-22 one of the lowest mass free-floating objects in the solar neighborhood. This object adds to the growing list of low-gravity field L dwarfs and is the first to be strongly deficient in methane relative to its estimated temperature. Comparing their spectra suggests that young L dwarfs with similar ages and temperatures can have different spectral signatures of youth. For the two objects with well constrained ages (PSO J318-22 and 2MASS J0355+11), we find their temperatures are ~400 K cooler than field objects of similar spectral type but their luminosities are comparable, i.e., these young L dwarfs are very red and unusually cool but not "underluminous." Altogether, PSO J318-22 is the first free-floating object with the colors, magnitudes, spectrum, luminosity, and mass that overlap the young dusty planets around HR 8799 and 2MASS J1207-39.

The Pan-STARRS1 survey is acquiring multi-epoch imaging in 5 bands (grizy) over the entire sky north of declination -30deg (the $3\pi$ survey). In July 2011 a test area of about 70 sq.deg. was observed to the expected final depth of the main survey. In this, the first of a series of papers targetting the galaxy count and clustering properties of the combined multi-epoch test area data, we present a detailed investigation into the depth of the survey and the reliability of the Pan-STARRS1 analysis software. We show that the Pan-STARRS1 reduction software can recover the properties of fake sources, and show good agreement between the magnitudes measured by Pan-STARRS1 and those from Sloan Digital Sky Survey. We also examine the number of false detections apparent in the Pan-STARRS1 data. Our comparisons show that the test area survey is somewhat deeper than the Sloan Digital Sky Survey in all bands, and, in particular, the z band approaches the depth of the stacked Sloan Stripe 82 data.

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

We present the discovery of a new dwarf galaxy, Perseus I/Andromeda XXXIII, found in the vicinity of Andromeda (M31) in stacked imaging data from the Pan-STARRS1 3{\pi} survey. Located 27.9{\deg} away from M31, Perseus I has a heliocentric distance of 785 +/- 65 kpc, compatible with it being a satellite of M31 at 374 +14/-10 kpc from its host. The properties of Perseus I are typical for a reasonably bright dwarf galaxy (M_V = -10.3 +/- 0.7), with an exponential half-light radius of r_h = 1.7 +/- 0.4 arcminutes or r_h = 400 +105/-85 pc at this distance, and a moderate ellipticity (\epsilon = 0.43 +0.15/-0.17). The late discovery of Perseus I is due to its fairly low surface brightness (\mu_0=25.7 +1.0/-0.9 mag/arcsec^2), and to the previous lack of deep, high quality photometric data in this region. If confirmed to be a companion of M31, the location of Perseus I, far east from its host, could place interesting constraints on the bulk motion of the satellite system of M31.

We present a sample of M31 beat Cepheids from the Pan-STARRS 1 PAndromeda campaign. By analyzing three years of PAndromeda data, we identify seventeen beat Cepheids, spreading from a galactocentric distance of 10 to 16 kpc. Since the relation between fundamental mode period and the ratio of fundamental to the first overtone period puts a tight constraint on metallicity we are able to derive the metallicity at the position of the beat Cepheids using the relations from the model of Buchler (2008). Our metallicity estimates show subsolar values within 15 kpc, similar to the metallicities from HII regions (Zurita & Bresolin 2012). We then use the metallicity estimates to calculate the metallicity gradient of the M31 disk, which we find to be closer to the metallicity gradient derived from planetary nebulae (Kwitter et al. 2012) than the metallicity gradient from HII regions (Zurita & Bresolin 2012).

We report photometric observations for comet C/2012 S1 (ISON) obtained during the time period immediately after discovery (r=6.28 AU) until it moved into solar conjunction in mid-2013 June using the UH2.2m, and Gemini North 8-m telescopes on Mauna Kea, the Lowell 1.8m in Flagstaff, the Calar Alto 1.2m telescope in Spain, the VYSOS-5 telescopes on Mauna Loa Hawaii and data from the CARA network. Additional pre-discovery data from the Pan STARRS1 survey extends the light curve back to 2011 September 30 (r=9.4 AU). The images showed a similar tail morphology due to small micron sized particles throughout 2013. Observations at sub-mm wavelengths using the JCMT on 15 nights between 2013 March 9 (r=4.52 AU) and June 16 (r=3.35 AU) were used to search for CO and HCN rotation lines. No gas was detected, with upper limits for CO ranging between (3.5-4.5)E27 molec/s. Combined with published water production rate estimates we have generated ice sublimation models consistent with the photometric light curve. The inbound light curve is likely controlled by sublimation of CO2. At these distances water is not a strong contributor to the outgassing. We also infer that there was a long slow outburst of activity beginning in late 2011 peaking in mid-2013 January (r~5 AU) at which point the activity decreased again through 2013 June. We suggest that this outburst was driven by CO injecting large water ice grains into the coma. Observations as the comet came out of solar conjunction seem to confirm our models.

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 initial results from a wide-field (30,000 deg^2) search for L/T transition brown dwarfs within 25 pc using the Pan-STARRS1 and WISE surveys. Previous large-area searches have been incomplete for L/T transition dwarfs, because these objects are faint in optical bands and have near-infrared colors that are difficult to distinguish from background stars. To overcome these obstacles, we have cross-matched the Pan-STARRS1 (optical) and WISE (mid-IR) catalogs to produce a unique multi-wavelength database for finding ultracool dwarfs. As part of our initial discoveries, we have identified seven brown dwarfs in the L/T transition within 9-15 pc of the Sun. The L9.5 dwarf PSO J140.2308+45.6487 and the T1.5 dwarf PSO J307.6784+07.8263 (both independently discovered by Mace et al. 2013) show possible spectroscopic variability at the Y- and J-bands. Two more objects in our sample show evidence of photometric J-band variability, and two others are candidate unresolved binaries based on their spectra. We expect our full search to yield a well-defined, volume-limited sample of L/T transition dwarfs that will include many new targets for study of this complex regime. PSO J307.6784+07.8263 in particular may be an excellent candidate for in-depth study of variability, given its brightness (J = 14.2 mag) and proximity (11 pc).

The Hyades cluster is an ideal target to study the dynamical evolution of a star cluster over the entire mass range due to its intermediate age and proximity to the Sun. We wanted to extend the Hyades mass function towards lower masses down to 0.1 Msol and to use the full three-dimensional spatial information to characterize the dynamical evolution of the cluster. We performed a kinematic and photometric selection using the PPMXL and Pan-STARRS1 sky surveys, to search for cluster members up to 30 pc from the cluster centre. We determined our detection efficiency and field star contamination rate to derive the cluster luminosity and mass functions down to masses of 0.1 Msol. The thorough astrometric and photometric constraints minimized the contamination. A minimum spanning tree algorithm was used to quantify the mass segregation. We discovered 43 new Hyades member candidates with velocity perpendicular to the Hyades motion up to 2 km/s. They have mass estimates between 0.43 and 0.09 Msol, for a total mass of 10 Msol. This doubles the number of Hyades candidates with masses smaller than 0.15 Msol. We provide an additional list of 11 possible candidates with velocity perpendicular to the Hyades motion up to 4 km/s. The cluster is significantly mass segregated. The extension of the mass function towards lower masses provided an even clearer signature than estimated in the past. We also identified as likely Hyades member an L0 dwarf previously assumed to be a field dwarf. Finally we question the membership of a number of previously published candidates, including a L2.5-type dwarf.

We present initial results from observations and numerical analyses aimed at characterizing main-belt comet P/2012 T1 (PANSTARRS). Optical monitoring observations were made between October 2012 and February 2013 using the University of Hawaii 2.2 m telescope, the Keck I telescope, the Baade and Clay Magellan telescopes, Faulkes Telescope South, the Perkins Telescope at Lowell Observatory, and the Southern Astrophysical Research (SOAR) telescope. The object's intrinsic brightness approximately doubles from the time of its discovery in early October until mid-November and then decreases by ~60% between late December and early February, similar to photometric behavior exhibited by several other main-belt comets and unlike that exhibited by disrupted asteroid (596) Scheila. We also used Keck to conduct spectroscopic searches for CN emission as well as absorption at 0.7 microns that could indicate the presence of hydrated minerals, finding an upper limit CN production rate of QCN<1.5x10^23 mol/s, from which we infer a water production rate of QH2O<5x10^25 mol/s, and no evidence of the presence of hydrated minerals. Numerical simulations indicate that P/2012 T1 is largely dynamically stable for >100 Myr and is unlikely to be a recently implanted interloper from the outer solar system, while a search for potential asteroid family associations reveal that it is dynamically linked to the ~155 Myr-old Lixiaohua asteroid family.

We report the discovery of two new dwarf galaxies, Lacerta I/Andromeda XXXI (Lac I/And XXXI) and Cassiopeia III/Andromeda XXXII (Cas III/And XXXII), in stacked Pan-STARRS1 r_P1- and i_P1-band imaging data. Both are luminous systems (M_V ~ -12) located at projected distances of 20.3{\deg} and 10.5{\deg} from M31. Lac I and Cas III are likely satellites of the Andromeda galaxy with heliocentric distances of 756^{+44}_{-28} kpc and 772^{+61}_{-56} kpc, respectively, and corresponding M31-centric distances of 275+/-7 kpc and 144^{+6}_{-4} kpc . The brightest of recent Local Group member discoveries, these two new dwarf galaxies owe their late discovery to their large sizes (r_h = 4.2^{+0.4}_{-0.5} arcmin or 912^{+124}_{-93} pc for Lac I; r_h = 6.5^{+1.2}_{-1.0} arcmin or 1456+/-267 pc for Cas III), and consequently low surface brightness (\mu_0 ~ 26.0 mag/arcsec^2), as well as to the lack of a systematic survey of regions at large radii from M31, close to the Galactic plane. This latter limitation is now alleviated by the 3{\pi} Pan-STARRS1 survey, which could lead to the discovery of other distant Andromeda satellite dwarf galaxies.

We report extensive observational data for five of the lowest redshift Super-Luminous Type Ic Supernovae (SL-SNe Ic) discovered to date, namely PTF10hgi, SN2011ke, PTF11rks, SN2011kf and SN2012il. Photometric imaging of the transients at +50 to +230 days after peak combined with host galaxy subtraction reveals a luminous tail phase for four of these SL-SNe. A high resolution, optical and near infrared spectrum from xshooter provides detection of a broad He I $\lambda$10830 emission line in the spectrum (+50d) of SN2012il, revealing that at least some SL-SNe Ic are not completely helium free. At first sight, the tail luminosity decline rates that we measure are consistent with the radioactive decay of \co, and would require 1-4M of \ni to produce the luminosity. These \ni masses cannot be made consistent with the short diffusion times at peak, and indeed are insufficient to power the peak luminosity. We instead favour energy deposition by newborn magnetars as the power source for these objects. A semi-analytical diffusion model with energy input from the spin-down of a magnetar reproduces the extensive lightcurve data well. The model predictions of ejecta velocities and temperatures which are required are in reasonable agreement with those determined from our observations. We derive magnetar energies of $0.4\lesssim E$($10^{51}$erg) $\lesssim6.9$ and ejecta masses of $2.3\lesssim M_{ej}$(\M) $\lesssim 8.6$. The sample of five SL-SNe Ic presented here, combined with SN 2010gx - the best sampled SL-SNe Ic so far - point toward an explosion driven by a magnetar as a viable explanation for all SL-SNe Ic.

We describe the Pan-STARRS Moving Object Processing System (MOPS), a modern software package that produces automatic asteroid discoveries and identifications from catalogs of transient detections from next-generation astronomical survey telescopes. MOPS achieves > 99.5% efficiency in producing orbits from a synthetic but realistic population of asteroids whose measurements were simulated for a Pan-STARRS4-class telescope. Additionally, using a non-physical grid population, we demonstrate that MOPS can detect populations of currently unknown objects such as interstellar asteroids. MOPS has been adapted successfully to the prototype Pan-STARRS1 telescope despite differences in expected false detection rates, fill-factor loss and relatively sparse observing cadence compared to a hypothetical Pan-STARRS4 telescope and survey. MOPS remains >99.5% efficient at detecting objects on a single night but drops to 80% efficiency at producing orbits for objects detected on multiple nights. This loss is primarily due to configurable MOPS processing limits that are not yet tuned for the Pan-STARRS1 mission. The core MOPS software package is the product of more than 15 person-years of software development and incorporates countless additional years of effort in third-party software to perform lower-level functions such as spatial searching or orbit determination. We describe the high-level design of MOPS and essential subcomponents, the suitability of MOPS for other survey programs, and suggest a road map for future MOPS development.

Asteroid detections in astronomical images may appear as trails due to a combination of their apparent rate of motion and exposure duration. Nearby asteroids in particular typically have high apparent rates of motion and acceleration. Their recovery, especially on their discovery apparition, depends upon obtaining good astrometry from the trailed detections. We present an analytic function describing a trailed detection under the assumption of a Gaussian point spread function (PSF) and constant rate of motion. We have fit the function to both synthetic and real trailed asteroid detections from the Pan-STARRS1 survey telescope to obtain accurate astrometry and photometry. For short trails our trailing function yields the same astrometric and photometry accuracy as a functionally simpler 2-d Gaussian but the latter underestimates the length of the trail - a parameter that can be important for measuring the object's rate of motion and assessing its cometary activity. For trails longer than about 10 pixels (> 3xPSF) our trail fitting provides 3-times better astrometric accuracy and up to 2 magnitudes improvement in the photometry. The trail fitting algorithm can be implemented at the source detection level for all detections to provide trail length and position angle that can be used to reduce the false tracklet rate.

We present a precise photometric calibration of the first 1.5 years of science imaging from the Pan-STARRS1 survey (PS1), an ongoing optical survey of the entire sky north of declination -30 degrees in five bands. Building on the techniques employed by Padmanabhan et al. (2008) in the Sloan Digital Sky Survey (SDSS), we use repeat PS1 observations of stars to perform the relative calibration of PS1 in each of its five bands, solving simultaneously for the system throughput, the atmospheric transparency, and the large-scale detector flat field. Both internal consistency tests and comparison against the SDSS indicate that we achieve relative precision of <10 mmag in g, r, and i_P1, and ~10 mmag in z and y_P1. The spatial structure of the differences with the SDSS indicates that errors in both the PS1 and SDSS photometric calibration contribute similarly to the differences. The analysis suggests that both the PS1 system and the Haleakala site will enable <1% photometry over much of the sky.

We report the discovery of a wide (approximately 400 AU projected separation), common proper motion companion to the nearby M dwarf LHS 2803 (PSO J207.0300-13.7422). This object was discovered during our census of the local T dwarf population using Pan-STARRS1 and 2MASS data. Using IRTF/SpeX near-infrared spectroscopy, we classify the secondary to be spectral type T5.5. University of Hawai`i 2.2m/SNIFS optical spectroscopy indicates the primary has a spectral type of M4.5, with approximately solar metallicity and no measurable H_alpha emission. We use this lack of activity to set a lower age limit for the system of 3.5 Gyr. Using a comparison with chance alignments of brown dwarfs and nearby stars, we conclude that the two objects are unlikely to be a chance association. The primary's photometric distance of 21 pc and its proper motion implies thin disk kinematics. Based on these kinematics and its metallicity, we set an upper age limit for the system of 10 Gyr. Evolutionary model calculations suggest the secondary has a mass of 72^+4_-7 M_Jup, temperature of 1120+-80 K, and log g=5.4+-0.1 dex. Model atmosphere fitting to the near-IR spectrum gives similar physical parameters of 1100 K and log g=5.0.

We report on our search for genetically related asteroids amongst the near-Earth object (NEO) population - families of NEOs akin to the well known main belt asteroid families. We used the technique proposed by Fu et al. (2005) supplemented with a detailed analysis of the statistical significance of the detected clusters. Their significance was assessed by comparison to identical searches performed on 1,000 'fuzzy-real' NEO orbit distribution models that we developed for this purpose. The family-free 'fuzzy-real' NEO models maintain both the micro and macro distribution of 5 orbital elements (ignoring the mean anomaly). Three clusters were identified that contain four or more NEOs but none of them are statistically significant at \geq 3{\sigma}. The most statistically significant cluster at the \sim 2{\sigma} level contains 4 objects with H < 20 and all members have long observational arcs and concomitant good orbital elements. Despite the low statistical significance we performed several other tests on the cluster to determine if it is likely a genetic family. The tests included examining the cluster's taxonomy, size-frequency distribution, consistency with a family-forming event during tidal disruption in a close approach to Mars, and whether it is detectable in a proper element cluster search. None of these tests exclude the possibility that the cluster is a family but neither do they confirm the hypothesis. We conclude that we have not identified any NEO families.

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

We report the discovery of a wide brown dwarf companion to the mildly metal-poor ([Fe/H]=-0.24), low galactic latitude (b = 1.88 deg) K4V star HIP 38939. The companion was discovered by its common proper motion with the primary and its red optical (Pan-STARRS1) and blue infrared (2MASS) colors. It has a projected separation of 1630 AU and a near-infrared spectral type of T4.5. As such it is one of only three known companions to a main sequence star which have early/mid-T spectral types (the others being HN Peg B and eps Indi B). Using chromospheric activity we estimate an age for the primary of 900{+1900,-600} Myr. This value is also in agreement with the age derived from the star's weak ROSAT detection. Comparison with evolutionary models for this age range indicates that HIP 38939B falls in the mass range 38+/-20 Mjup with an effective temperature range of 1090+/-60 K. Fitting our spectrum with atmospheric models gives a best fitting temperature of 1100 K. We include our object in an analysis of the population of benchmark T dwarfs and find that while older atmospheric models appeared to over-predict the temperature of the coolest objects compared to evolutionary models, more recent atmospheric models provide better agreement.

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.

The discovery of new objects in modern wide-field asteroid and comet surveys can be enhanced by first identifying observations belonging to known solar system objects. The assignation of new observations to a known object is an attribution problem that occurs when a least squares orbit already exists for the object but a separate fit is not possible to just the set of new observations. In this work we explore the strongly asymmetric attribution problem in which the existing least squares orbit is very well constrained and the new data are sparse. We describe an attribution algorithm that introduces new quality control metrics in the presence of strong biases in the astrometric residuals. The main biases arise from the stellar catalogs used in the reduction of asteroid observations and we show that a simple debiasing with measured regional catalog biases significantly improves the results. We tested the attribution algorithm using data from the PS1 survey that used the 2MASS star catalog for the astrometric reduction. We found small but statistically significant biases in the data of up to 0.1 arcsec that are relevant only when the observations reach the level of accuracy made possible by instruments like PS1. The false attribution rate was measured to be < 1/1,000 with a simple additional condition that can reduce it to zero while the attribution efficiency is consistent with 100%.

The Pan-STARRS 1 (PS1) survey of M31 (PAndromeda) is designed to identify gravitational microlensing events, caused by bulge and disk stars (self-lensing) and by compact matter in the halos of M31 and the Milky Way (halo lensing, or lensing by MACHOs). With the 7 deg2 FOV of PS1, the entire disk of M31 can be imaged with one single pointing. Our aim is to monitor M31 with this wide FOV with daily sampling (20 mins/day). In the 2010 season we acquired in total 91 nights towards M31, with 90 nights in the rP1 and 66 nights in the iP1. The total integration time in rP1 and iP1 are 70740s and 36180s, respectively. As a preliminary analysis, we study a 40'\times40' sub-field in the central region of M31, a 20'\times20' sub-field in the disk of M31 and a 20'\times20' sub-field for the investigation of astrometric precision. We demonstrate that the PSF is good enough to detect microlensing events. We present light curves for 6 candidate microlensing events. This is a competitive rate compared to previous M31 microlensing surveys. We finally also present one example light curve for Cepheids, novae and eclipsing binaries in these sub-fields.

The Pan-STARRS1 survey is collecting multi-epoch, multi-color observations of the sky north of declination -30 deg to unprecedented depths. These data are being photometrically and astrometrically calibrated and will serve as a reference for many other purposes. In this paper we present our determination of the Pan-STARRS photometric system: gp1, rp1, ip1, zp1, yp1, and wp1. The Pan-STARRS photometric system is fundamentally based on the HST Calspec spectrophotometric observations, which in turn are fundamentally based on models of white dwarf atmospheres. We define the Pan-STARRS magnitude system, and describe in detail our measurement of the system passbands, including both the instrumental sensitivity and atmospheric transmission functions. Byproducts, including transformations to other photometric systems, galactic extinction, and stellar locus are also provided. We close with a discussion of remaining systematic errors.

We present observations of comet-like main-belt object P/2010 R2 (La Sagra) obtained by Pan-STARRS 1 and the Faulkes Telescope-North on Haleakala in Hawaii, the University of Hawaii 2.2 m, Gemini-North, and Keck I telescopes on Mauna Kea, the Danish 1.54 m telescope at La Silla, and the Isaac Newton Telescope on La Palma. An antisolar dust tail is observed from August 2010 through February 2011, while a dust trail aligned with the object's orbit plane is also observed from December 2010 through August 2011. Assuming typical phase darkening behavior, P/La Sagra is seen to increase in brightness by >1 mag between August 2010 and December 2010, suggesting that dust production is ongoing over this period. These results strongly suggest that the observed activity is cometary in nature (i.e., driven by the sublimation of volatile material), and that P/La Sagra is therefore the most recent main-belt comet to be discovered. We find an approximate absolute magnitude for the nucleus of H_R=17.9+/-0.2 mag, corresponding to a nucleus radius of ~0.7 km, assuming an albedo of p=0.05. Using optical spectroscopy, we find no evidence of sublimation products (i.e., gas emission), finding an upper limit CN production rate of Q_CN<6x10^23 mol/s, from which we infer an H2O production rate of Q_H2O<10^26 mol/s. Numerical simulations indicate that P/La Sagra is dynamically stable for >100 Myr, suggesting that it is likely native to its current location and that its composition is likely representative of other objects in the same region of the main belt, though the relatively close proximity of the 13:6 mean-motion resonance with Jupiter and the (3,-2,-1) three-body mean-motion resonance with Jupiter and Saturn mean that dynamical instability on larger timescales cannot be ruled out.

Main belt asteroid (300163) 2006 VW139 (later designated P/2006 VW139) was discovered to exhibit comet-like activity by the Pan-STARRS1 survey telescope using automated point-spread-function analyses performed by PS1's Moving Object Processing System. Deep follow-up observations show both a short (\sim 10") antisolar dust tail and a longer (\sim 60") dust trail aligned with the object's orbit plane, similar to the morphology observed for another main-belt comet, P/2010 R2 (La Sagra), and other well-established comets, implying the action of a long-lived, sublimation-driven emission event. Photometry showing the brightness of the near-nucleus coma remaining constant over \sim 30 days provides further evidence for this object's cometary nature, suggesting it is in fact a main-belt comet, and not a disrupted asteroid. A spectroscopic search for CN emission was unsuccessful, though we find an upper limit CN production rate of Q_CN < 1.3x10^24 mol/s, from which we infer a water production rate of Q_H2O < 10^26 mol/s. We also find an approximately linear optical spectral slope of 7.2%/1000A, similar to other cometary dust comae. Numerical simulations indicate that P/2006 VW139 is dynamically stable for > 100 Myr, while a search for a potential asteroid family around the object reveals a cluster of 24 asteroids within a cutoff distance of 68 m/s. At 70 m/s, this cluster merges with the Themis family, suggesting that it could be similar to the Beagle family to which another main-belt comet, 133P/Elst-Pizarro, belongs.

The Pan-STARRS1 survey is obtaining multi-epoch imaging in 5 bands (gps rps ips zps yps) over the entire sky North of declination -30deg. We describe here the implementation of the Photometric Classification Server (PCS) for Pan-STARRS1. PCS will allow the automatic classification of objects into star/galaxy/quasar classes based on colors, the measurement of photometric redshifts for extragalactic objects, and constrain stellar parameters for stellar objects, working at the catalog level. We present tests of the system based on high signal-to-noise photometry derived from the Medium Deep Fields of Pan-STARRS1, using available spectroscopic surveys as training and/or verification sets. We show that the Pan-STARRS1 photometry delivers classifications and photometric redshifts as good as the Sloan Digital Sky Survey (SDSS) photometry to the same magnitude limits. In particular, our preliminary results, based on this relatively limited dataset down to the SDSS spectroscopic limits and therefore potentially improvable, show that stars are correctly classified as such in 85% of cases, galaxies in 97% and QSOs in 84%. False positives are less than 1% for galaxies, ~19% for stars and ~28% QSOs. Moreover, photometric redshifts for 1000 luminous red galaxies up to redshift 0.5 are determined to 2.4% precision with just 0.4% catastrophic outliers and small (-0.5%) residual bias. PCS will create a value added catalog with classifications and photometric redshifts for eventually many millions sources.

(abridged) The activities of the Commission have continued to focus on controlling unwanted light and radio emissions at observatory sites, monitoring of conditions at observatory sites, and education and outreach. Commission members have been active in securing new legislation in several locations to further the protection of observatory sites as well as in the international regulation of the use of the radio spectrum and the protection of radio astronomical observations. In 2009, at its XXVII General Assembly in Rio de Janeiro, the IAU passed Resolution B5 in Defence of the Night Sky and the Right to Starlight. The Resolution encourages IAU members to assist in raising public awareness about the contents and objectives of the International Conference in Defence of the Quality of the Night Sky and the Right to Observe Stars [http://www.starlight2007.net/], in particular the importance of preserving access to an unpolluted night sky for all mankind.

The Pan-STARRS1 survey has obtained multi-epoch imaging in five bands (Pan-STARRS1 gps, rps, ips, zps, and yps) on twelve "Medium Deep Fields", each of which spans a 3.3 degree circle. For the period between Apr 2009 and Apr 2011 these fields were observed 50-200 times. Using a reduced proper motion diagram, we have extracted a list of 47 white dwarf (WD) candidates whose Pan-STARRS1 astrometry indicates a non-zero proper motion at the 6-sigma level, with a typical 1-sigma proper motion uncertainty of 10 mas/yr. We also used astrometry from SDSS (when available) and USNO-B to assess our proper motion fits. None of the WD candidates exhibits evidence of statistically significant parallaxes, with a typical 1-sigma uncertainty of 8 mas. Twelve of these candidates are known WDs, including the high proper motion (1.7"/yr) WD LHS 291. We confirm three more objects as WDs through optical spectroscopy. Based on the Pan-STARRS1 colors, ten of the stars are likely to be cool WDs with 4170 K Teff 5000 K and cooling ages <9 Gyr. We classify these objects as likely thick disk WDs based on their kinematics. Our current sample represents only a small fraction of the Pan-STARRS1 data. With continued coverage from the Medium Deep Field Survey and the 3pi survey, Pan-STARRS1 should find many more high proper motion WDs that are part of the old thick disk and halo.

We have searched ~8200 sq. degs for high proper motion (~0.5-2.7"/year) T dwarfs by combining first-epoch data from the Pan-STARRS1 (PS1) 3-Pi Survey, the 2MASS All-Sky Point Source Catalog, and the WISE Preliminary Data Release. We identified two high proper motion objects with the very red (W1-W2) colors characteristic of T dwarfs, one being the known T7.5 dwarf GJ 570D. Near-IR spectroscopy of the other object (PSO J043.5+02 = WISEP J0254+0223) reveals a spectral type of T8, leading to a photometric distance of 7.2+/-0.7 pc. The 2.56"/yr proper motion of PSO J043.5+02 is the second highest among field T dwarfs, corresponding to an tangential velocity of 87+/-8 km/s. According to the Besancon galaxy model, this velocity indicates its galactic membership is probably in the thin disk, with the thick disk an unlikely possibility. Such membership is in accord with the near-IR spectrum, which points to a surface gravity (age) and metallicity typical of the field population. We combine 2MASS, SDSS, WISE, and PS1 astrometry to derive a preliminary parallax of 171+/-45 mas (5.8{+2.0}{-1.2} pc), the first such measurement using PS1 data. The proximity and brightness of PSO J043+02 will facilitate future characterization of its atmosphere, variability, multiplicity, distance, and kinematics. The modest number of candidates from our search suggests that the immediate (~10 pc) solar neighborhood does not contain a large reservoir of undiscovered T dwarfs earlier than about T8.

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.

A complete well-defined sample of ultracool dwarfs is one of the key science programs of the Pan-STARRS 1 optical survey telescope (PS1). Here we combine PS1 commissioning data with 2MASS to conduct a proper motion search (0.1--2.0\arcsec/yr) for nearby T dwarfs, using optical+near-IR colors to select objects for spectroscopic followup. The addition of sensitive far-red optical imaging from PS1 enables discovery of nearby ultracool dwarfs that cannot be identified from 2MASS data alone. We have searched 3700 sq. deg. of PS1 y-band (0.95--1.03 um) data to y$\approx$19.5 mag (AB) and J$\approx$16.5 mag (Vega) and discovered four previously unknown bright T dwarfs. Three of the objects (with spectral types T1.5, T2 and T3.5) have photometric distances within 25 pc and were missed by previous 2MASS searches due to more restrictive color selection criteria. The fourth object (spectral type T4.5) is more distant than 25 pc and is only a single-band detection in 2MASS. We also examine the potential for completing the census of nearby ultracool objects with the PS1 3$\pi$ survey.

Recent searches by unbiased, wide-field surveys have uncovered a group of extremely luminous optical transients. The initial discoveries of SN 2005ap by the Texas Supernova Search and SCP-06F6 in a deep Hubble pencil beam survey were followed by the Palomar Transient Factory confirmation of host redshifts for other similar transients. The transients share the common properties of high optical luminosities (peak magnitudes ~ -21 to -23), blue colors, and a lack of H or He spectral features. The physical mechanism that produces the luminosity is uncertain, with suggestions ranging from jet-driven explosion to pulsational pair-instability. Here we report the most detailed photometric and spectral coverage of an ultra-bright transient (SN 2010gx) detected in the Pan-STARRS 1 sky survey. In common with other transients in this family, early-time spectra show a blue continuum, and prominent broad absorption lines of O II. However, about 25d after discovery, the spectra developed type Ic supernova features, showing the characteristic broad Fe II and Si II absorption lines. Detailed, post-maximum follow-up may show that all SN 2005ap and SCP-06F6 type transients are linked to supernovae Ic. This poses problems in understanding the physics of the explosions: there is no indication from late-time photometry that the luminosity is powered by 56Ni, the broad lightcurves suggest very large ejected masses, and the slow spectral evolution is quite different from typical Ic timescales. The nature of the progenitor stars and the origin of the luminosity are intriguing and open questions.

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.

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 performed a simulation of a next generation sky survey's (Pan-STARRS 1) efficiency for detecting Earth-impacting asteroids. The steady-state sky-plane distribution of the impactors long before impact is concentrated towards small solar elongations (Chesley and Spahr, 2004) but we find that there is interesting and potentially exploitable behavior in the sky-plane distribution in the months leading up to impact. The next generation surveys will find most of the dangerous impactors (>140m diameter) during their decade-long survey missions though there is the potential to miss difficult objects with long synodic periods appearing in the direction of the Sun, as well as objects with long orbital periods that spend much of their time far from the Sun and Earth. A space-based platform that can observe close to the Sun may be needed to identify many of the potential impactors that spend much of their time interior to the Earth's orbit. The next generation surveys have a good chance of imaging a bolide like 2008TC3 before it enters the atmosphere but the difficulty will lie in obtaining enough images in advance of impact to allow an accurate pre-impact orbit to be computed.

(Abridged) We have obtained radial velocities of a sample of 18 ultracool dwarfs (M6.5-T8) using high-resolution, near-infrared spectra obtained with NIRSPEC and the Keck II telescope. We have confirmed that the radial velocity of Gl 570 D is coincident with that of the K-type primary star Gl 570 A, thus providing additional support for their true companionship. The presence of planetary-mass companions around 2MASS J05591914-1404488 (T4.5V) has been analyzed using five NIRSPEC radial velocity measurements obtained over a period of 4.37 yr. We have computed UVW space motions for a total of 21 L and T dwarfs within 20 pc of the Sun. This population shows UVW velocities that nicely overlap the typical kinematics of solar to M-type stars within the same spatial volume. However, the mean Galactic (44.2 km/s) and tangential (36.5 km/s) velocities of the L and T dwarfs appear to be smaller than those of G to M stars. A significant fraction (~40%) of the L and T dwarfs lies near the Hyades moving group (0.4-2 Gyr), which contrasts with the 10-12% found for earlier-type stellar neighbors. Additionally, the distributions of all three UVW components (sigma_{UVW} = 30.2, 16.5, 15.8 km/s) and the distributions of the total Galactic (sigma_{v_tot} = 19.1 km/s) and tangential (sigma_{v_t} = 17.6 km/s) velocities derived for the L and T dwarf sample are narrower than those measured for nearby G, K, and M-type stars, but similar to the dispersions obtained for F stars. This suggests that, in the solar neighborhood, the L- and T-type ultracool dwarfs in our sample (including brown dwarfs) is kinematically younger than solar-type to early M stars with likely ages in the interval 0.5-4 Gyr.

The Panoramic Survey Telescope And Rapid Response System (Pan-STARRS) under development at the University of Hawaii's Institute for Astronomy is creating the first fully automated end-to-end Moving Object Processing System (MOPS) in the world. It will be capable of identifying detections of moving objects in our solar system and linking those detections within and between nights, attributing those detections to known objects, calculating initial and differentially-corrected orbits for linked detections, precovering detections when they exist, and orbit identification. Here we describe new kd-tree and variable-tree algorithms that allow fast, efficient, scalable linking of intra and inter-night detections. Using a pseudo-realistic simulation of the Pan-STARRS survey strategy incorporating weather, astrometric accuracy and false detections we have achieved nearly 100% efficiency and accuracy for intra-night linking and nearly 100% efficiency for inter-night linking within a lunation. At realistic sky-plane densities for both real and false detections the intra-night linking of detections into `tracks' currently has an accuracy of 0.3%. Successful tests of the MOPS on real source detections from the Spacewatch asteroid survey indicate that the MOPS is capable of identifying asteroids in real data.

The nearby brown dwarf LP 944-20 has been monitored for radial velocity variability at optical and near-infrared wavelengths using the VLT/UVES and the Keck/NIRSPEC spectrographs, respectively. The UVES radial velocity data obtained over 14 nights spanning a baseline of 841 days shows significant variability with an amplitude of 3.5 km s$^{-1}$. The periodogram analysis of the UVES data indicates a possible period between 2.5 hours and 3.7 hours, which is likely due to the rotation of the brown dwarf. However, the NIRSPEC data obtained over 6 nights shows an rms dispersion of only 0.36 km s$^{-1}$ and do not follow the periodic trend. These results indicate that the variability seen with UVES is likely to be due to rotationally modulated inhomogeneous surface features. We suggest that future planet searches around very low-mass stars and brown dwarfs using radial velocities will be better conducted in the near-infrared than in the optical.

We have obtained projected rotation velocities (vsini) of a sample of 19 ultracool dwarfs with spectral types in the interval M6.5-T8 using high-resolution, near-infrared spectra obtained with NIRSPEC and the Keck II telescope. Among our targets there are two young brown dwarfs, two likely field stars, and fifteen likely brown dwarfs (30-72 Mjup) of the solar neighborhood. Our results indicate that the T-type dwarfs are fast rotators in marked contrast to M-type stars. We have derived vsini velocities between <15 and 40 km/s for them, and have found no clear evidence for T dwarfs rotating strongly faster than L dwarfs. However, there is a hint for an increasing lower envelope on moving from mid-M to the L spectral types in the vsini-spectral type diagram that was previously reported in the literature; our vsini results extend it to even cooler types. Assuming that field brown dwarfs have a size of 0.08-0.1 Rsol, we can place an upper limit of 12.5 h on the equatorial rotation period of T-type brown dwarfs. In addition, we have compared our vsini measurements to spectroscopic rotational velocities of very young brown dwarfs of similar mass available in the literature. The comparison, although model-dependent, suggests that brown dwarfs lose some angular momentum during their contraction; however, their spin down time seems to be significantly longer than that of solar-type to early-M stars.

The University of Hawaii Wide-Field Imager (UHWFI) is a focal compressor system designed to project the full half-degree field of the UH 2.2 m telescope onto the refurbished UH 8Kx8K CCD camera. The optics use Ohara glasses and are mounted in an oil-filled cell to minimize light losses and ghost images from the large number of internal lens surfaces. The UHWFI is equipped with a six-position filter wheel and a rotating sector blade shutter,both driven by stepper motors. The instrument saw first light in 2004 in an engineering mode. After filling the lens cell with index matching oil, integration of all software components into the user interface, tuning of the CCD performance, and the purchase of the final filter set, UHWFI is now fully commissioned at the UH 2.2 m telescope.

We present a spectrophotometric analysis of S Ori J053825.4-024241, a candidate member close to the substellar boundary of the young (1-8 Myr), nearby (~360 pc) sigma Orionis star cluster. Our optical and near-infrared photometry and low-resolution spectroscopy indicate that S Ori J053825.4-024241 is a likely cluster member with a mass estimated from evolutionary models at 0.06+0.07-0.02 Msol, which makes the object a probable brown dwarf. The radial velocity of S Ori J053825.4-024241 is similar to the cluster systemic velocity. This target, which we have classified as an M 6.0+-1.0 low-gravity object, shows excessemission in the near-infrared and anomalously strong photometric variability for its type (from the blue to the J band), suggesting the presence of a surrounding disc. The optical spectroscopic observations show a continuum excess at short wavelengths and a persistent and resolved Halpha emission (pseudo-equivalent width of ~-250 AA) in addition to the presence of other forbidden and permitted emission lines, which we interpret as indicating accretion from the disc and possibly mass loss. We conclude that despite the low mass of S Ori J053825.4-024241, this object exhibits some of the properties typical of active classical T Tauri stars.

We present photometric and spectroscopic observations of 23 high redshift supernovae spanning a range of z=0.34-1.03, 9 of which are unambiguously classified as Type Ia. These supernovae were discovered during the IfA Deep Survey, which began in September 2001 and observed a total of 2.5 square degrees to a depth of approximately m=25-26 in RIZ over 9-17 visits, typically every 1-3 weeks for nearly 5 months, with additional observations continuing until April 2002. We give a brief description of the survey motivations, observational strategy, and reduction process. This sample of 23 high-redshift supernovae includes 15 at z>0.7, doubling the published number of objects at these redshifts, and indicates that the evidence for acceleration of the universe is not due to a systematic effect proportional to redshift. In combination with the recent compilation of Tonry et al. (2003), we calculate cosmological parameter density contours which are consistent with the flat universe indicated by the CMB (Spergel et al. 2003). Adopting the constraint that Omega_total = 1.0, we obtain best-fit values of (Omega_m, Omega_Lambda)=(0.33, 0.67) using 22 SNe from this survey augmented by the literature compilation. We show that using the empty-beam model for gravitational lensing does not eliminate the need for Omega_Lambda > 0. Experience from this survey indicates great potential for similar large-scale surveys while also revealing the limitations of performing surveys for z>1 SNe from the ground.

We present the discovery of a distant methane dwarf, the first from the Institute for Astronomy (IfA) Deep Survey. The object ("IfA 0230-Z1") was identified from deep optical I and z'-band imaging, being conducted as an IfA-wide collaboration using the prime-focus imager Suprime-Cam on the Subaru 8.2-m Telescope. IfA 0230-Z1 is extremely red in the Iz'J (0.8--1.2 micron) bands but relatively blue in J-H; such colors are uniquely characteristic of T dwarfs. A near-IR spectrum taken with the Keck Telescope shows strong H2O absorption and a continuum break indicative of CH4, confirming the object has a very cool atmosphere. Comparison with nearby T dwarfs gives a spectral type of T3-T4 and a distance of ~45 pc. Simple estimates based on previous T dwarf discoveries suggest that the IfA survey will find a comparable number of T dwarfs as the 2MASS survey, albeit at a much larger average distance. We also discuss the survey's ability to probe the galactic scale height of ultracool (L and T) dwarfs.

One of the main scientific drivers of the Space InfraRed Telescope Facility (SIRTF) is the search for brown dwarfs and extrasolar superplanets. We discuss observational strategies for identification of these objects, and conclude that an optimal strategy is a wide IRAC survey (18 deg$^2$) with a 5 $\sigma$ sensitivity of 3.9 $\mu$Jy in channel 2 (M$\sim19.1^m$). For this sensitivity, we provide estimates of the number of low mass brown dwarfs and isolated planets detected per square degree for power-law mass functions with $\alpha$=1.5, and 1.0. Shallower surveys covering a larger area are inefficient because of large overheads and detector noise. Deeper surveys covering a smaller area become more and more affected by crowding with galaxies. A survey like the one that we propose would determine the field mass function down to a few Jupiter masses through the identification of a large sample of brown dwarfs and isolated planets. The proposed SIRTF survey would also allow the first detection of ultracool substellar objects with temperatures between 700 K and 200 K. The cooling curves of substellar objects with masses less than 20 Jupiters imply that they should spend most of their lifetimes at temperatures below 700 K. Preliminary models indicate that their atmospheres could be dominated by water clouds, which would diminish their optical and near-infrared fluxes. The properties of those objects are still completely unexplored.

We present the results of a wide-field K-selected galaxy survey with complementary optical I- and B-band imaging in six fields with a total coverage of 9.8 square degrees. This survey establishes the bright-end K-band galaxy number counts in the magnitude range 13<K<16 with high precision. We find that our bright-end counts have a significantly steeper slope than the prediction of a no-evolution model, which cannot be accounted for by known observational or theoretical error. We also argue against the likelihood of sufficient evolution at such low redshifts to account for this effect, we describe an alternative picture in which there is a local deficiency of galaxies by a factor of 2 on scale sizes of around 300 Mpc. Taken at face value, this would imply that local measurements of \Omega_0 underestimate the true value of the cosmological mass density by this factor and that local measurements of $H_0$ could be high by as much as 33%.

The Ursa Major Cluster has received remarkably little attention, although it is as near as the Virgo Cluster and contains a comparable number of HI-rich galaxies. In this paper, criteria for group membership are discussed and data are presented for 79 galaxies identified with the group. Of these, all 79 have been imaged at B,R,I bands with CCDs, 70 have been imaged at K' with a HgCdTe array detector, and 70 have been detected in the HI 21cm line. A complete sample of 62 galaxies brighter than M(B)=-16.5 is identified. Images and gradients in surface brightness and color are presented at a common linear scale. As has been seen previously, the galaxies with the reddest global colors are reddest at the centers and get bluer at large radii. However, curiously, among the galaxies with the bluest global colors there are systems with very blue cores that get redder at large radii.