We present the first release of the MaNGA Stellar Library (MaStar), which is a large, well-calibrated, high-quality empirical library covering the wavelength range of 3,622-10,354A at a resolving power of R~1800. The spectra were obtained using the same instrument as used by the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) project, by piggybacking on the SDSS-IV/APOGEE-2N observations. Compared to previous empirical libraries, the MaStar library will have a higher number of stars and a more comprehensive stellar-parameter coverage, especially of cool dwarfs, low-metallicity stars, and stars with different [alpha/Fe], achieved by a sophisticated target selection strategy that takes advantage of stellar-parameter catalogs from the literature. This empirical library will provide a new basis for stellar population synthesis, and is particularly well-suited for stellar-population analysis of MaNGA galaxies. The first version of the library contains 8646 high-quality per-visit spectra for 3321 unique stars. Compared to photometry, the relative flux calibration of the library is accurate to 3.9% in g-r, 2.7% in r-i, and 2.2% in i-z. The data are released as part of Sloan Digital Sky Survey Data Release 15. We expect the final release of the library to contain more than 10,000 stars.

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.

We present a detailed characterization of the 849 broad-line quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Our quasar sample covers a redshift range of 0.1<z<4.5 and is flux-limited to i_PSF<21.7 without any other cuts on quasar properties. The main sample characterization includes: 1) spectral measurements of the continuum and broad emission lines for individual objects from the coadded first-season spectroscopy in 2014; 2) identification of broad and narrow absorption lines in the spectra; 3) optical variability properties for continuum and broad lines from multi-epoch spectroscopy. We provide improved systemic redshift estimates for all quasars, and demonstrate the effects of signal-to-noise ratio on the spectral measurements. We compile measured properties for all 849 quasars along with supplemental multi-wavelength data for subsets of our sample from other surveys. The SDSS-RM sample probes a diverse range in quasar properties, and shows well detected continuum and broad-line variability for many objects from first-season monitoring data. The compiled properties serve as the benchmark for follow-up work based on SDSS-RM data. The spectral fitting tools are made public along with this work.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

Individual highly magnified stars have been recently discovered at lookback times of more than half the age of the Universe, in lensed galaxies that straddle the critical curves of massive galaxy clusters. Having confirmed their detectability, it is now important to carry out systematic searches for them in order to establish their frequency, and in turn learn about the statistical properties of high-redshift stars and of the granularity of matter in the foreground deflector. Here we report the discovery of a highly magnified star at redshift $z = 0.94$ in a strongly lensed arc behind a Hubble Frontier Field galaxy cluster, MACS J0416.1-2403, discovered as part of a systematic archival search. The bright transient (dubbed "Warhol") was discovered in Hubble Space Telescope data taken on 2014 September 15 and 16. This single image faded over a period of two weeks, and observations taken on 2014 September 1 show that the duration of the microlensing event was at most four weeks in total. The light curve may also exhibit slow changes over a period of years consistent with the level of microlensing expected from stars responsible for the intracluster light (ICL) of the cluster. Optical and infrared observations taken near peak brightness can be fit by a stellar spectrum with moderate host-galaxy extinction. A blue supergiant matches the measured spectral energy distribution near peak, implying a temporary magnification of at least several thousand. While the spectrum of an O-type star would also fit the transient's spectral energy distribution, extremely luminous O-type stars are much less common than blue supergiants. The short timescale of the event and the estimated effective temperature indicate that the lensed source is an extremely magnified star.

Active galactic nuclei (AGN) that show strong rest-frame optical/UV variability in their blue continuum and broad line emission are classified as "changing-look" AGN, or at higher luminosities changing look quasars (CLQs). These surprisingly large and sometimes rapid transitions challenge accepted models of quasar physics and duty cycles, offer several new avenues for study of quasar host galaxies, and open a wider interpretation of the cause of differences between broad and narrow line AGN. To better characterize extreme quasar variability, we present follow-up spectroscopy as part of a comprehensive search for CLQs across the full SDSS footprint using spectroscopically confirmed quasars from the SDSS DR7 catalog. Our primary selection requires large-amplitude (|\Delta g|>1 mag, |\Delta r|>0.5 mag) variability over any of the available time baselines probed by the SDSS and Pan-STARRS 1 surveys. We employ photometry from the Catalina Sky Survey to verify variability behavior in CLQ candidates where available, and confirm CLQs using optical spectroscopy from the William Herschel, MMT, Magellan, and Palomar telescopes. For our adopted S/N threshold on variability of broad H\beta emission, we find 17 new CLQs, yielding a confirmation rate of >~ 20%. These candidates are at lower Eddington ratio relative to the overall quasar population which supports a disk-wind model for the broad line region. Based on our sample, the CLQ fraction increases from 10% to roughly half as the continuum flux ratio between repeat spectra at 3420 Angstroms increases from 1.5 to 6. We release a catalog of over 200 highly variable candidates to facilitate future CLQ searches.

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.

Over 3 billion astronomical objects have been detected in the more than 22 million orthogonal transfer CCD images obtained as part of the Pan-STARRS1 $3\pi$ survey. Over 85 billion instances of those objects have been automatically detected and characterized by the Pan-STARRS Image Processing Pipeline photometry software, psphot. This fast, automatic, and reliable software was developed for the Pan-STARRS project, but is easily adaptable to images from other telescopes. We describe the analysis of the astronomical objects by psphot in general as well as for the specific case of the 3rd processing version used for the first two public releases of the Pan-STARRS $3\pi$ survey data, DR1 & DR2.

The Pan-STARRS Data Processing System is responsible for the steps needed to downloaded, archive, and process all images obtained by the Pan-STARRS telescopes, including real-time detection of transient sources such as supernovae and moving objects including potentially hazardous asteroids. With a nightly data volume of up to 4 terabytes and an archive of over 4 petabytes of raw imagery, Pan-STARRS is solidly in the realm of Big Data astronomy. The full data processing system consists of several subsystems covering the wide range of necessary capabilities. This article describes the Image Processing Pipeline and its connections to both the summit data systems and the outward-facing systems downstream. The latter include the Moving Object Processing System (MOPS) & the public database: the Published Science Products Subsystem (PSPS).

Luminous high-redshift quasars can be used to probe of the intergalactic medium (IGM) in the early universe because their UV light is absorbed by the neutral hydrogen along the line of sight. They help us to measure the neutral hydrogen fraction of the high-z universe, shedding light on the end of reionization epoch. In this paper, we present a discovery of a new quasar (PSO J006.1240+39.2219) at redshift $z=6.61\pm0.02$ from Panoramic Survey Telescope & Rapid Response System 1. Including this quasar, there are nine quasars above $z>6.5$ up to date. The estimated continuum brightness is $M_\text{1450}$=$-25.96\pm0.08$. PSO J006.1240+39.2219 has a strong Ly~$\alpha$ emission compared with typical low-redshift quasars, but the measured near-zone region size is $R_\text{NZ}=3.2\pm1.1$ proper megaparsecs, which is consistent with other quasars at z$\sim$6.

Recently Assef et al.(2018) presented two catalogs of AGN candidates over 30,093 deg^2 selected from the Wide-field Infrared Survey Explorer (WISE) observations. From their most reliable sample, Assef et al. (2018) identified 45 AGN candidates with the highest variability levels in the AllWISE catalog, but that are not blazars. Here we present new spectroscopic observations of some of these targets to further constrain their nature. We also study their optical lightcurves using observations from CRTS, and find that only seven show significant optical variability, and that five of those seven are spectroscopically classified as AGN. In one of them, WISEA J094806.56+031801.7 (W0948+0318), we identify a transient event in the CRTS lightcurve. We present a detailed analysis of this transient, and characterize it through its CRTS lightcurve and its multi-wavelength spectral energy distribution obtained from GALEX, Pan-STARRS and WISE observations. We find that the most likely source of the transient is a super-luminous supernova (SLSN) in W0948+0318. We estimate the total radiated energy to be E=1.6\pm 0.3 x 10^52 erg, making it one of the most energetic SLSN observed. Based on the lack of change in mid-IR color throughout and after the transient event, we speculate that the location of the SLSN is within the torus of the AGN. We identify 9 possible analogs to W0948+0318 based on their WISE lightcurves. None show optically detected transients and hence suggest significant dust obscuration. Finally, we estimate a rate of >2x10^-7 yr^-1 per AGN for these transients under the conservative assumption that none of the identified analogs have a common origin with the transient in W0948+0318.

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.

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.

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

Galaxy-cluster gravitational lenses can magnify background galaxies by a total factor of up to ~50. Here we report an image of an individual star at redshift z=1.49 (dubbed "MACS J1149 Lensed Star 1 (LS1)") magnified by >2000. A separate image, detected briefly 0.26 arcseconds from LS1, is likely a counterimage of the first star demagnified for multiple years by a >~3 solar-mass object in the cluster. For reasonable assumptions about the lensing system, microlensing fluctuations in the stars' light curves can yield evidence about the mass function of intracluster stars and compact objects, including binary fractions and specific stellar evolution and supernova models. Dark-matter subhalos or massive compact objects may help to account for the two images' long-term brightness ratio.

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

We present 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

A galaxy cluster acts as a cosmic telescope over background galaxies but also as a cosmic microscope of the lens imperfections. The diverging magnification of lensing caustics enhances the microlensing effect of substructure present within the lensing mass. Fine-scale structure can be accessed as a moving background source brightens and disappears when crossing these caustics. The recent recognition of a distant lensed star near the Einstein radius of the galaxy cluster MACSJ1149.5+2223 (Kelly et al. 2017) allows the rare opportunity to reach subsolar mass microlensing through a super-critical column of cluster matter. Here we compare these observations with high-resolution ray-tracing simulations that include stellar microlensing set by the observed intracluster starlight and also primordial black holes that may be responsible for the recently observed LIGO events. We explore different scenarios with microlenses from the intracluster medium and black holes, including primordial ones, and examine strategies to exploit these unique alignments. We find that the best constraints on the fraction of compact dark matter in the small-mass regime can be obtained in regions of the cluster where the intracluster medium plays a negligible role. This new lensing phenomenon should be widespread and can be detected within modest-redshift lensed galaxies so that the luminosity distance is not prohibitive for detecting individual magnified stars. Continuous {\it Hubble Space Telescope} monitoring of several such optimal arcs will be rewarded by an unprecedented mass spectrum of compact objects that can contribute to uncovering the nature of dark matter.

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.

The recent discovery of fast transient events near critical curves of massive galaxy clusters, which are interpreted as highly magnified individual stars in giant arcs due to caustic crossing, opens up the possibility of using such microlensing events to constrain a range of dark matter models such as primordial black holes and scalar field dark matter. Based on a simple analytic model, we study lensing properties of a point mass lens embedded in a high magnification region, and derive the dependence of the peak brightness, microlensing time scales, and event rates on the mass of the point mass lens as well as the radius of a source star that is magnified. We find that the lens mass and source radius of the first event MACS J1149 Lensed Star 1 (LS1) are constrained, with the lens mass range of $0.1~M_\odot \lesssim M \lesssim 4\times 10^3M_\odot$ and the source radius range of $40~R_\odot \lesssim R \lesssim 260~R_\odot$. In the most plausible case with $M\approx 0.3~M_\odot$ and $R\approx 180~R_\odot$, the source star should have been magnified by a factor of $\approx 4300$ at the peak. The derived lens properties are fully consistent with the interpretation that MACS J1149 LS1 is a microlensing event produced by a star that contributes to the intra-cluster light. We argue that compact dark matter models with high fractional mass densities for the mass range $10^{-5}M_\odot \lesssim M\lesssim 10^2M_\odot$ are inconsistent with the observation of MACS J1149 LS1 because such models predict too low magnifications. Our work demonstrates a potential use of caustic crossing events in giant arcs to constrain compact dark matter.

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.

As astronomers increasingly exploit the information available in the time domain, spectroscopic variability in particular opens broad new channels of investigation. Here we describe the selection algorithms for all targets intended for repeat spectroscopy in the Time Domain Spectroscopic Survey (TDSS), part of the extended Baryon Oscillation Spectroscopic Survey within the Sloan Digital Sky Survey-IV. Also discussed are the scientific rationale and technical constraints leading to these target selections. The TDSS includes a large "Repeat Quasar Spectroscopy" (RQS) program delivering ~13,000 repeat spectra of confirmed SDSS quasars, and several smaller "Few-Epoch Spectroscopy" (FES) programs targeting specific classes of quasars as well as stars. The RQS program aims to provide a large and diverse quasar data set for studying variations in quasar spectra on timescales of years, a comparison sample for the FES quasar programs, and opportunity for discovering rare, serendipitous events. The FES programs cover a wide variety of phenomena in both quasars and stars. Quasar FES programs target broad absorption line quasars, high signal-to-noise ratio normal broad line quasars, quasars with double-peaked or very asymmetric broad emission line profiles, binary supermassive black hole candidates, and the most photometrically variable quasars. Strongly variable stars are also targeted for repeat spectroscopy, encompassing many types of eclipsing binary systems, and classical pulsators like RR Lyrae. Other stellar FES programs allow spectroscopic variability studies of active ultracool dwarf stars, dwarf carbon stars, and white dwarf/M dwarf spectroscopic binaries. We present example TDSS spectra and describe anticipated sample sizes and results.

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

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.

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.

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.

Brueckner theory is used to investigate the in-medium properties of a $\Lambda$-hyperon in nuclear and neutron matter, based on hyperon-nucleon interactions derived within SU(3) chiral effective field theory (EFT). It is shown that the resulting $\Lambda$ single-particle potential $U_\Lambda(p_\Lambda =0,\rho)$ becomes strongly repulsive for densities $\rho$ of two-to-three times that of normal nuclear matter. Adding a density-dependent effective $\Lambda N$-interaction constructed from chiral $\Lambda NN$ three-body forces increases the repulsion further. Consequences of these findings for neutron stars are discussed. It is argued that for hyperon-nuclear interactions with properties such as those deduced from the SU(3) EFT potentials, the onset for hyperon formation in the core of neutron stars is expected to be shifted to extremely high baryon density, thus potentially resolving the so-called hyperon puzzle.

In the conventional framework for cosmological dynamics the scale factor $a(t)$ is assumed to obey the `background' Friedmann equation for a perfectly homogeneous universe while particles move according to equations of motions driven by the gravity of the density fluctuations. It has recently been suggested that the emergence of structure modifies the evolution of $a(t)$ via Newtonian (or `kinematic') backreaction and that this may avoid the need for dark energy. Here we point out that the conventional system of equations is exact in Newtonian gravity and there is no approximation in the use of the homogeneous universe equation for $a(t)$. The recently proposed modification of Racz et al.\ (2017) does not reduce to Newtonian gravity in the limit of low velocities. We discuss the relation of this to the `generalised Friedmann equation' of Buchert and Ehlers. These are quite different things; their formula describes individual regions and is obtained under the restrictive assumption that the matter behaves like a pressure-free fluid whereas our result is exact for collisionless dynamics and is an auxiliary relation appearing in the structure equations. We use the symmetry of the general velocity autocorrelation function to show how Buchert's $\cal Q$ tends very rapidly to zero for large volume and that this does not simply arise `by construction' through the adoption of periodic boundary conditions as has been claimed. We conclude that, to the extent that Newtonian gravity accurately describes the low-$z$ universe, there is no backreaction of structure on $a(t)$ and that the need for dark energy cannot be avoided in this way.

RR Lyrae stars may be the best practical tracers of Galactic halo (sub-)structure and kinematics. The PanSTARRS1 (PS1) $3\pi$ survey offers multi-band, multi-epoch, precise photometry across much of the sky, but a robust identification of RR Lyrae stars in this data set poses a challenge, given PS1's sparse, asynchronous multi-band light curves ($\lesssim 12$ epochs in each of five bands, taken over a 4.5-year period). We present a novel template fitting technique that uses well-defined and physically motivated multi-band light curves of RR Lyrae stars, and demonstrate that we get accurate period estimates, precise to 2~sec in $>80\%$ of cases. We augment these light curve fits with other {\em features} from photometric time-series and provide them to progressively more detailed machine-learned classification models. From these models we are able to select the widest ($3/4$ of the sky) and deepest (reaching 120 kpc) sample of RR Lyrae stars to date. The PS1 sample of $\sim 45,000$ RRab stars is pure (90\%), and complete (80\% at 80 kpc) at high galactic latitudes. It also provides distances precise to 3\%, measured with newly derived period-luminosity relations for optical/near-infrared PS1 bands. With the addition of proper motions from {\em Gaia} and radial velocity measurements from multi-object spectroscopic surveys, we expect the PS1 sample of RR Lyrae stars to become the premier source for studying the structure, kinematics, and the gravitational potential of the Galactic halo. The techniques presented in this study should translate well to other sparse, multi-band data sets, such as those produced by the Dark Energy Survey and the upcoming Large Synoptic Survey Telescope Galactic plane sub-survey.

We analyze four extreme active galactic nuclei (AGN) transients to explore the possibility that they are caused by rare, high-amplitude microlensing events. These previously unknown type-I AGN are located in the redshift range 0.6-1.1 and show changes of > 1.5 mag in the g-band on a time-scale of ~years. Multi-epoch optical spectroscopy, from the William Herschel Telescope, shows clear differential variability in the broad line fluxes with respect to the continuum changes and also evolution in the line profiles. In two cases, a simple point-source, point-lens microlensing model provides an excellent match to the long-term variability seen in these objects. For both models, the parameter constraints are consistent with the microlensing being due to an intervening stellar mass object but as yet there is no confirmation of the presence of an intervening galaxy. The models predict a peak amplification of 10.3/13.5 and an Einstein time-scale of 7.5/10.8 yr, respectively. In one case, the data also allow constraints on the size of the CIII] emitting region, with some simplifying assumptions, to to be ~1.0-6.5 light-days and a lower limit on the size of the MgII emitting region to be > 9 light-days (half-light radii). This CIII] radius is perhaps surprisingly small. In the remaining two objects, there is spectroscopic evidence for an intervening absorber but the extra structure seen in the light curves requires a more complex lensing scenario to adequately explain.

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.

In this study, we present observations of a type Iax supernova, PS1-12bwh, discovered during the Pan-STARRS1 3$\pi$-survey. Our analysis was driven by previously unseen pre-maximum, spectroscopic heterogeneity. While the light curve and post-maximum spectra of PS1-12bwh are virtually identical to those of the well-studied type Iax supernova, SN 2005hk, the $-$2 day spectrum of PS1-12bwh does not resemble SN 2005hk at a comparable epoch; instead, we found it to match a spectrum of SN 2005hk taken over a week earlier ($-$12 day). We are able to rule out the cause as being incorrect phasing, and argue that it is not consistent with orientation effects predicted by existing explosion simulations. To investigate the potential source of this difference, we performed radiative transfer modelling of both supernovae. We found that the pre-maximum spectrum of PS1-12bwh is well matched by a synthetic spectrum generated from a model with a lower density in the high velocity ($\gtrsim$6000 km~s$^{-1}$) ejecta than SN 2005hk. The observed differences between SN 2005hk and PS1-12bwh may therefore be attributed primarily to differences in the high velocity ejecta alone, while comparable densities for the lower velocity ejecta would explain the nearly identical post-maximum spectra. These two supernovae further highlight the diversity within the SNe Iax class, as well as the challenges in spectroscopically identifying and phasing these objects, especially at early epochs.

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 investigate a sample of 2293 ICRF2 extragalactic radio-loud sources with accurate positions determined by VLBI, mostly active galactic nuclei (AGN) and quasars, which are cross-matched with optical sources in the first Gaia release (Gaia DR1). The distribution of offsets between the VLBI sources and their optical counterparts is strongly non-Gaussian, with powerful wings extending beyond 1 arcsecond. Limiting our analysis to only high-confidence difference detections, we find (and publish) a list of 188 objects with normalized variances above 12 and offsets below 1 arcsecond. Pan-STARRS stacked and monochromatic images resolve some of these sources indicating the presence of double sources, confusion sources, or pronounced extended structures. Some 89 high-quality objects, however, do not show any perturbations and appear to be star-like single sources, yet displaced by multiples of the expected error from the radio-loud AGN. We conclude that a fraction of luminous AGN (more than 4%) can be physically dislodged from the optical centers of their parent galaxies.

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.

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

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

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.

We derive the expression for the observed redshift in the weak field limit in the observer's past light cone, including all relativistic terms up to second order in velocity. We then apply it to compute the cluster-galaxy cross-correlation functions (CGCF) using N-body simulations. The CGCF is asymmetric along the line of sight (LOS) owing to the presence of the small second order terms such as the gravitational redshift (GRedshift). We identify two systematics in the modelling of the GRedshift signal in stacked clusters. First, it is affected by the morphology of dark matter haloes and the large-scale cosmic-web. The non-spherical distribution of galaxies around the central halo and the presence of neighbouring clusters systematically reduce the GRedshift signal. This bias is approximately 20% for $M_{\rm min}\simeq 10^{14} {\rm M_{\odot}}/h$, and is more than 50% for haloes with $M_{\rm min}\simeq 2\times 10^{13} {\rm M_{\odot}}/h$ at $r>$4 Mpc/$h$. Second, the best-fit gravitational redshift profiles as well as the profiles of all other relativistic terms are found to be significantly different in velocity space compared to their real space versions. We discuss some subtleties relating to these effects in velocity space. We also find that the S/N of the GRedshift signal increases with decreasing halo mass.

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

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

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

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.

We apply a combination of a Genetic Algorithms (GA) and Support Vector Machines (SVM) machine learning algorithm to solve two important problems faced by the astronomical community: star/galaxy separation, and photometric redshift estimation of galaxies in survey catalogs. We use the GA to select the relevant features in the first step, followed by optimization of SVM parameters in the second step to obtain an optimal set of parameters to classify or regress, in process of which we avoid over-fitting. We apply our method to star/galaxy separation in Pan-STARRS1 data. We show that our method correctly classifies 98% of objects down to i_P1= 24.5, with a completeness (or true positive rate) of 99% for galaxies, and 88% for stars. By combining colors with morphology, our star/classification method yields better results than the new SExtractor classifier spread_model in particular at the faint end (i_P1>22). We also use our method to derive photometric redshifts for galaxies in the COSMOS bright multi-wavelength dataset down to an error in (1+z) of sigma=0.013, which compares well with estimates from SED fitting on the same data (sigma=0.007) while making a significantly smaller number of assumptions.

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.

The Time-Domain Spectroscopic Survey (TDSS) is an SDSS-IV eBOSS subproject primarily aimed at obtaining identification spectra of ~220,000 optically-variable objects systematically selected from SDSS/Pan-STARRS1 multi-epoch imaging. We present a preview of the science enabled by TDSS, based on TDSS spectra taken over ~320 deg^2 of sky as part of the SEQUELS survey in SDSS-III, which is in part a pilot survey for eBOSS in SDSS-IV. Using the 15,746 TDSS-selected single-epoch spectra of photometrically variable objects in SEQUELS, we determine the demographics of our variability-selected sample, and investigate the unique spectral characteristics inherent in samples selected by variability. We show that variability-based selection of quasars complements color-based selection by selecting additional redder quasars, and mitigates redshift biases to produce a smooth quasar redshift distribution over a wide range of redshifts. The resulting quasar sample contains systematically higher fractions of blazars and broad absorption line quasars than from color-selected samples. Similarly, we show that M-dwarfs in the TDSS-selected stellar sample have systematically higher chromospheric active fractions than the underlying M-dwarf population, based on their H-alpha emission. TDSS also contains a large number of RR Lyrae and eclipsing binary stars with main-sequence colors, including a few composite-spectrum binaries. Finally, our visual inspection of TDSS spectra uncovers a significant number of peculiar spectra, and we highlight a few cases of these interesting objects. With a factor of ~15 more spectra, the main TDSS survey in SDSS-IV will leverage the lessons learned from these early results for a variety of time-domain science applications.

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.

We present a recalibration of the Sloan Digital Sky Survey (SDSS) photometry with new flat fields and zero points derived from Pan-STARRS1 (PS1). Using PSF photometry of 60 million stars with $16 < r < 20$, we derive a model of amplifier gain and flat-field corrections with per-run RMS residuals of 3 millimagnitudes (mmag) in $griz$ bands and 15 mmag in $u$ band. The new photometric zero points are adjusted to leave the median in the Galactic North unchanged for compatibility with previous SDSS work. We also identify transient non-photometric periods in SDSS ("contrails") based on photometric deviations co-temporal in SDSS bands. The recalibrated stellar PSF photometry of SDSS and PS1 has an RMS difference of {9,7,7,8} mmag in $griz$, respectively, when averaged over $15'$ regions.

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.

Current cosmological analyses which use Type Ia supernova (SN Ia) observations combine SN samples to expand the redshift range beyond that of a single sample and increase the overall sample size. The inhomogeneous photometric calibration between different SN samples is one of the largest systematic uncertainties of the cosmological parameter estimation. To place these different samples on a single system, analyses currently use observations of a small sample of very bright flux standards on the $HST$ system. We propose a complementary method, called `Supercal', in which we use measurements of secondary standards in each system, compare these to measurements of the same stars in the Pan-STARRS1 (PS1) system, and determine offsets for each system relative to PS1, placing all SN observations on a single, consistent photometric system. PS1 has observed $3\pi$ of the sky and has a relative calibration of better than 5 mmag (for $\sim15<griz<21$ mag), making it an ideal reference system. We use this process to recalibrate optical observations taken by the following SN samples: PS1, SNLS, SDSS, CSP, and CfA1-4. We measure discrepancies on average of 10 mmag, but up to 35 mmag, in various optical passbands. We find that correcting for these differences changes recovered values for the dark energy equation-of-state parameter, $w$, by on average $2.6\%$. This change is roughly half the size of current statistical constraints on $w$. The size of this effect strongly depends on the error in the $B-V$ calibration of the low-$z$ surveys. The Supercal method will allow future analyses to tie past samples to the best calibrated sample.

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

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

The IceCube neutrino observatory pursues a follow-up program selecting interesting neutrino events in real-time and issuing alerts for electromagnetic follow-up observations. In March 2012, the most significant neutrino alert during the first three years of operation was issued by IceCube. In the follow-up observations performed by the Palomar Transient Factory (PTF), a Type IIn supernova (SN) PTF12csy was found $0.2^\circ$ away from the neutrino alert direction, with an error radius of $0.54^\circ$. It has a redshift of $z=0.0684$, corresponding to a luminosity distance of about $300 \, \mathrm{Mpc}$ and the Pan-STARRS1 survey shows that its explosion time was at least 158 days (in host galaxy rest frame) before the neutrino alert, so that a causal connection is unlikely. The a posteriori significance of the chance detection of both the neutrinos and the SN at any epoch is $2.2 \, \sigma$ within IceCube's 2011/12 data acquisition season. Also, a complementary neutrino analysis reveals no long-term signal over the course of one year. Therefore, we consider the SN detection coincidental and the neutrinos uncorrelated to the SN. However, the SN is unusual and interesting by itself: It is luminous and energetic, bearing strong resemblance to the SN IIn 2010jl, and shows signs of interaction of the SN ejecta with a dense circumstellar medium. High-energy neutrino emission is expected in models of diffusive shock acceleration, but at a low, non-detectable level for this specific SN. In this paper, we describe the SN PTF12csy and present both the neutrino and electromagnetic data, as well as their analysis.

The Ophiuchus stream is a recently discovered stellar tidal stream in the Milky Way. We present high-quality spectroscopic data for 14 stream member stars obtained using the Keck and MMT telescopes. We confirm the stream as a fast moving ($v_{los}\sim290$ km s$^{-1}$), kinematically cold group ($\sigma_{v_{los}}\lesssim1$ km s$^{-1}$) of $\alpha$-enhanced and metal-poor stars (${\rm [\alpha/Fe]\sim0.4}$ dex, ${\rm [Fe/H]\sim-2.0}$ dex). Using a probabilistic technique, we model the stream simultaneously in line-of-sight velocity, color-magnitude, coordinate, and proper motion space, and so determine its distribution in 6D phase-space. We find that that the stream extends in distance from 7.5 to 9 kpc from the Sun; it is 50 times longer than wide, merely appearing highly foreshortened in projection. The analysis of the stellar population contained in the stream suggests that it is $\sim12$ Gyr old, and that its initial stellar mass was $\sim2\times10^4$ $M_{\odot}$ (or at least $\gtrsim7\times10^3$ $M_{\odot}$). Assuming a fiducial Milky Way potential, we fit an orbit to the stream which matches the observed phase-space distribution, except for some tension in the proper motions: the stream has an orbital period of $\sim350$ Myr, and is on a fairly eccentric orbit ($e\sim0.66$) with a pericenter of $\sim3.5$ kpc and an apocenter of $\sim17$ kpc. The phase-space structure and stellar population of the stream show that its progenitor must have been a globular cluster that was disrupted only $\sim240$ Myr ago. We do not detect any significant overdensity of stars along the stream that would indicate the presence of a progenitor, and conclude that the stream is all that is left of the progenitor.

We report the discovery and characterisation of a deeply eclipsing AM CVn-system, Gaia14aae (= ASSASN-14cn). Gaia14aae was identified independently by the All-Sky Automated Survey for Supernovae (ASAS-SN; Shappee et al. 2014) and by the Gaia Science Alerts project, during two separate outbursts. A third outburst is seen in archival Pan-STARRS-1 (PS1; Schlafly et al. 2012; Tonry et al. 2012; Magnier et al. 2013) and ASAS-SN data. Spectroscopy reveals a hot, hydrogen-deficient spectrum with clear double-peaked emission lines, consistent with an accreting double degenerate classification. We use follow-up photometry to constrain the orbital parameters of the system. We find an orbital period of 49.71 min, which places Gaia14aae at the long period extremum of the outbursting AM CVn period distribution. Gaia14aae is dominated by the light from its accreting white dwarf. Assuming an orbital inclination of 90 degrees for the binary system, the contact phases of the white dwarf lead to lower limits of 0.78 M solar and 0.015 M solar on the masses of the accretor and donor respectively and a lower limit on the mass ratio of 0.019. Gaia14aae is only the third eclipsing AM CVn star known, and the first in which the WD is totally eclipsed. Using a helium WD model, we estimate the accretor's effective temperature to be 12900+-200 K. The three out-burst events occurred within 4 months of each other, while no other outburst activity is seen in the previous 8 years of Catalina Real-time Transient Survey (CRTS; Drake et al. 2009), Pan-STARRS-1 and ASAS-SN data. This suggests that these events might be rebrightenings of the first outburst rather than individual events.

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 the selection algorithm and anticipated results for the Time Domain Spectroscopic Survey (TDSS). TDSS is an SDSS-IV eBOSS subproject that will provide initial identification spectra of approximately 220,000 luminosity-variable objects (variable stars and AGN) across 7,500 square degrees selected from a combination of SDSS and multi-epoch Pan-STARRS1 photometry. TDSS will be the largest spectroscopic survey to explicitly target variable objects, avoiding pre-selection on the basis of colors or detailed modeling of specific variability characteristics. Kernel Density Estimate (KDE) analysis of our target population performed on SDSS Stripe 82 data suggests our target sample will be 95% pure (meaning 95% of objects we select have genuine luminosity variability of a few magnitudes or more). Our final spectroscopic sample will contain roughly 135,000 quasars and 85,000 stellar variables, approximately 4,000 of which will be RR Lyrae stars which may be used as outer Milky Way probes. The variability-selected quasar population has a smoother redshift distribution than a color-selected sample, and variability measurements similar to those we develop here may be used to make more uniform quasar samples in large surveys. The stellar variable targets are distributed fairly uniformly across color space, indicating that TDSS will obtain spectra for a wide variety of stellar variables including pulsating variables, stars with significant chromospheric activity, cataclysmic variables and eclipsing binaries. TDSS will serve as a pathfinder mission to identify and characterize the multitude of variable objects that will be detected photometrically in even larger variability surveys such as LSST.

A long standing question in cosmology is whether gravitational lensing changes the distance-redshift relation $D(z)$ or the mean flux density of sources. Interest in this has been rekindled by recent studies in non-linear relativistic perturbation theory that find biases in both the area of a surface of constant redshift and in the mean distance to this surface, with a fractional bias in both cases on the order of the mean squared convergence $\langle \kappa^2 \rangle$. Any such area bias could alter CMB cosmology, and the corresponding bias in mean flux density could affect supernova cosmology. Here we show that, in an ensemble averaged sense, the perturbation to the area of a surface of constant redshift is in reality much smaller, being on the order of the cumulative bending angle squared, or roughly a part-in-a-million effect. This validates the arguments of Weinberg (1976) that the mean magnification $\mu$ of sources is unity and of Kibble \& Lieu (2005) that the mean direction-averaged inverse magnification is unity. It also validates the conventional treatment of lensing in analysis of CMB anisotropies. But the existence of a scatter in magnification will cause any non-linear function of these conserved quantities to be statistically biased. The distance $D$, for example, is proportional to $\mu^{-1/2}$ so lensing will bias $\langle D\rangle$ even if $\langle \mu \rangle=1$. The fractional bias in such quantities is generally of order $\langle \kappa^2 \rangle$, which is orders of magnitude larger than the area perturbation. Claims for large bias in area or flux density of sources appear to have resulted from misinterpretation of such effects: they do not represent a new non-Newtonian effect, nor do they invalidate standard cosmological analyses.

We consider some aspects of the perturbation to the luminosity distance $d(z)$ that are of relevance for SN1a cosmology and for future peculiar velocity surveys at non-negligible redshifts. 1) Previous work has shown that the correction to the lowest order perturbation $\delta d / d = -\delta v / c z$ has the peculiar characteristic that it appears to depend on the absolute state of motion of sources, rather than on their motion relative to that of the observer. The resolution of this apparent violation of the equivalence principle is that it is necessary to allow for evolution of the velocities with time, and also, when considering perturbations on the scale of the observer-source separation, to include the gravitational redshift effect. We provide an expression for $\delta d / d$ that provides a physically consistent way to measure peculiar velocities and determine their impact for SN1a cosmology. 2) We then calculate the perturbation to the redshift as a function of source flux density, which has been proposed as an alternative probe of large-scale motions. We show how the inclusion of surface brightness modulation modifies the relation between $\delta z(m)$ and the peculiar velocity, and that, while the noise properties of this method might appear promising, the velocity signal is swamped by the effect of galaxy clustering for most scales of interest. 3) We show how, in linear theory, peculiar velocity measurements are biased downwards by the effect of smaller scale motions or by measurement errors (such as in photometric redshifts). Our results nicely explain the effects seen in simulations by Koda et al.\ 2013. We critically examine the prospects for extending peculiar velocity studies to larger scales with near-term future surveys.

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.

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.

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.

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.

Luminous distant quasars are unique probes of the high redshift intergalactic medium (IGM) and of the growth of massive galaxies and black holes in the early universe. Absorption due to neutral Hydrogen in the IGM makes quasars beyond a redshift of z~6.5 very faint in the optical $z$-band, thus locating quasars at higher redshifts require large surveys that are sensitive above 1 micron. We report the discovery of three new z>6.5 quasars, corresponding to an age of the universe of <850 Myr, selected as z-band dropouts in the Pan-STARRS1 survey. This increases the number of known z>6.5 quasars from 4 to 7. The quasars have redshifts of z=6.50, 6.52, and 6.66, and include the brightest z-dropout quasar reported to date, PSO J036.5078+03.0498 with M_1450=-27.4. We obtained near-infrared spectroscopy for the quasars and from the MgII line we estimate that the central black holes have masses between 5x10^8 and 4x10^9 M_sun, and are accreting close to the Eddington limit (L_Bol/L_Edd=0.13-1.2). We investigate the ionized regions around the quasars and find near zone radii of R_NZ=1.5-5.2 proper Mpc, confirming the trend of decreasing near zone sizes with increasing redshift found for quasars at 5.7<z<6.4. By combining R_NZ of the PS1 quasars with those of 5.7<z<7.1 quasars in the literature, we derive a luminosity corrected redshift evolution of R_NZ,corrected=(7.2+/-0.2)-(6.1+/-0.7)x(z-6) Mpc. However, the large spread in R_NZ in the new quasars implies a wide range in quasar ages and/or a large variation in the neutral Hydrogen fraction along different lines of sight.

Recent calculations using non-linear relativistic cosmological perturbation theory show biases in the mean luminosity distance and distance modulus at low redshift. We show that these effects may be understood very simply as a non-relativistic, and purely kinematic, Malmquist-like bias, and we describe how the effect changes if one averages over sources that are limited by apparent magnitude. This effect is essentially identical to the distance bias from small-scale random velocities that has previously been considered by astronomers, though we find that the standard formula overestimates the homogeneous bias by a factor 2.

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

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

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

The goal of this work is to conduct a photometric study of eclipsing binaries in M31. We apply a modified box-fitting algorithm to search for eclipsing binary candidates and determine their period. We classify these candidates into detached, semi-detached, and contact systems using the Fourier decomposition method. We cross-match the position of our detached candidates with the photometry from Local Group Survey (Massey et al. 2006) and select 13 candidates brighter than 20.5 magnitude in V. The relative physical parameters of these detached candidates are further characterized with Detached Eclipsing Binary Light curve fitter (DEBiL) by Devor (2005). We will followup the detached eclipsing binaries spectroscopically and determine the distance to M31.

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.

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

We present 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 measure quasar variability using the Panoramic Survey Telescope and Rapid Response System 1 Survey (Pan-STARRS1 or PS1) and the Sloan Digital Sky Survey (SDSS) and establish a method of selecting quasars via their variability in 10,000 square degree surveys. We use 100,000 spectroscopically confirmed quasars that have been well measured in both PS1 and SDSS and take advantage of the decadal time scales that separate SDSS measurements and PS1 measurements. A power law model fits the data well over the entire time range tested, 0.01 to 10 years. Variability in the current PS1-SDSS dataset can efficiently distinguish between quasars and non-varying objects. It improves the purity of a griz quasar color cut from 4.1% to 48% while maintaining 67% completeness. Variability will be very effective at finding quasars in datasets with no u band and in redshift ranges where exclusively photometric selection is not efficient. We show that quasars' rest-frame ensemble variability, measured as a root mean squared in delta magnitudes, is consistent with V(z, L, t) = A0 (1+z)^0.37 (L/L0)^-0.16 (t/1yr)^0.246 , where L0 = 10^46 ergs^-1 and A0 = 0.190, 0.162, 0.147 or 0.141 in the gP1 , rP1 , iP1 or zP1 filter, respectively. We also fit across all four filters and obtain median variability as a function of z, L and lambda as V(z, L, lambda, t) = 0.079(1 + z)^0.15 (L/L0 )^-0.2 (lambda/1000 nm)^-0.44 (t/1yr)^0.246 .

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.

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.

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 a method for selecting RR Lyrae (RRL) stars (or other type of variable stars) in the absence of a large number of multi-epoch data and light curve analyses. Our method uses color and variability selection cuts that are defined by applying a Gaussian Mixture Bayesian Generative Method (GMM) on 636 pre-identified RRL stars instead of applying the commonly used rectangular cuts. Specifically, our method selects 8,115 RRL candidates (heliocentric distances < 70 kpc) using GMM color cuts from the Sloan Digital Sky Survey (SDSS) and GMM variability cuts from the Panoramic Survey Telescope and Rapid Response System 1 3pi survey (PS1). Comparing our method with the Stripe 82 catalog of RRL stars shows that the efficiency and completeness levels of our method are ~77% and ~52%, respectively. Most contaminants are either non-variable main-sequence stars or stars in eclipsing systems. The method described here efficiently recovers known stellar halo substructures. It is expected that the current completeness and efficiency levels will further improve with the additional PS1 epochs (~3 epochs per filter) that will be observed before the conclusion of the survey. A comparison between our efficiency and completeness levels using the GMM method to the efficiency and completeness levels using rectangular cuts that are commonly used yielded a significant increase in the efficiency level from ~13% to ~77% and an insignificant change in the completeness levels. Hence, we favor using the GMM technique in future studies. Although we develop it over the SDSS X PS1 footprint, the technique presented here would work well on any multi-band, multi-epoch survey for which the number of epochs is limited.

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 perform a study on the optical and infrared photometric properties of known luminous blue variables (LBVs) in M31 using the sample of LBV candidates from the Local Group Galaxy Survey (Massey et al. 2007). We find that M31 LBV candidates show photometric variability ranging from 0.375 to 1.576 magnitudes in rP1 during a three year time-span observed by the Pan-STARRS 1 Andromeda survey (PAndromeda). Their near-infrared colors also follow the distribution of Galactic LBVs as shown by Oksala et al. (2013). We use these features as selection criteria to search for unknown LBV candidates in M31. We thus devise a method to search for candidate LBVs using both optical color from the Local Group Galaxy Survey and infrared color from Two Micron All Sky Survey, as well as photometric variations observed by PAndromeda. We find four sources exhibiting common properties of known LBVs. These sources also exhibit UV emission as seen from GALEX, which is one of the previously adopted method to search for LBV candidates. The locations of the LBVs are well aligned withM31 spiral arms as seen in the UV light, suggesting they are evolved stars at young age given their high-mass nature. We compare these candidates with the latest Geneva evolutionary tracks, which show that our new M31 LBV candidates are massive evolved stars with an age of 10 to 100 million years.

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.

Distance measurements to molecular clouds are important, but are often made separately for each cloud of interest, employing very different different data and techniques. We present a large, homogeneous catalog of distances to molecular clouds, most of which are of unprecedented accuracy. We determine distances using optical photometry of stars along lines of sight toward these clouds, obtained from PanSTARRS-1. We simultaneously infer the reddenings and distances to these stars, tracking the full probability distribution function using a technique presented in Green et al. (2014). We fit these star-by-star measurements using a simple dust screen model to find the distance to each cloud. We thus estimate the distances to almost all of the clouds in the Magnani et al. (1985) catalog, as well as many other well-studied clouds, including Orion, Perseus, Taurus, Cepheus, Polaris, California, and Monoceros R2, avoiding only the inner Galaxy. Typical statistical uncertainties in the distances are 5%, though the systematic uncertainty stemming from the quality of our stellar models is about 10%. The resulting catalog is the largest catalog of accurate, directly-measured distances to molecular clouds. Our distance estimates are generally consistent with available distance estimates from the literature, though in some cases the literature estimates are off by a factor of more than two.

P/2011 S1 (Gibbs) is an outer solar system comet or active Centaur with a similar orbit to that of the famous 29P/Schwassmann-Wachmann 1. P/2011 S1 (Gibbs) has been observed by the Pan-STARRS 1 (PS1) sky survey from 2010 to 2012. The resulting data allow us to perform multi-color studies of the nucleus and coma of the comet. Analysis of PS1 images reveals that P/2011 S1 (Gibbs) has a small nucleus $< 4$ km radius, with colors $g_{P1}-r_{P1} = 0.5 \pm 0.02$, $r_{P1}-i_{P1} = 0.12 \pm 0.02$ and $i_{P1}-z_{P1} = 0.46 \pm 0.03$. The comet remained active from 2010 to 2012, with a model-dependent mass-loss rate of $\sim100$ kg s$^{-1}$. The mass-loss rate per unit surface area of P/2011 S1 (Gibbs) is as high as that of 29P/Schwassmann-Wachmann 1, making it one of the most active Centaurs. The mass-loss rate also varies with time from $\sim 40$ kg s$^{-1}$ to 150 kg s$^{-1}$. Due to its rather circular orbit, we propose that P/2011 S1 (Gibbs) has 29P/Schwassmann-Wachmann 1-like outbursts that control the outgassing rate. The results indicate that it may have a similar surface composition to that of 29P/Schwassmann-Wachmann 1. Our numerical simulations show that the future orbital evolution of P/2011 S1 (Gibbs) is more similar to that of the main population of Centaurs than to that of 29P/Schwassmann-Wachmann 1. The results also demonstrate that P/2011 S1 (Gibbs) is dynamically unstable and can only remain near its current orbit for roughly a thousand years.

Nova Delphini 2013 was identified on the 14th of August 2013 and eventually rose to be a naked eye object. We sought to study the behaviour of the object in the run-up to outburst and to compare it to the pre-outburst photometric characteristics of other novae. We searched the Pan-STARRS 1 datastore to identify pre-outburst photometry of Nova Del 2013 and identified twenty-four observations in the 1.2 years before outburst. The progenitor of Nova Delphini showed variability of a few tenths of a magnitude but did not brighten significantly in comparison with archival plate photometry. We also found that the object did not vary significantly on the approximately half hour timescale between pairs of Pan-STARRS 1 observations.

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.

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)

In homogeneous cosmological models the wavelength $\lambda$ of a photon exchanged between two fundamental observers changes in proportion to expansion of the space $D$ between them, so $\Delta\log(\lambda / D) = 0$. This is exactly the same as for a pair of observers receding from each other in flat space-time where the effect is purely kinematic. The interpretation of this has been the subject of considerable debate, and it has been suggested that all redshifts are a relative velocity effect, raising the question of whether the wavelength always stretches in proportion to the emitter-receiver separation. Here we show that, for low redshift at least, $\Delta\log(\lambda / D)$ vanishes for a photon exchanged between any two freely-falling observers in a spatially constant tidal field, because such a field stretches wavelengths and the space between the observers identically. But in general there is a non-kinematic, and essentially gravitational, component of the redshift that is given by a weighted average of the gradient of the tidal field along the photon path. While the redshift can always be formally expressed using the Doppler formula, in situations where the gravitational redshift dominates, the `relative velocity' is typically quite different from the rate of change of $D$ and it is misleading to think of the redshift as being a velocity or `kinematic' effect.

We measure the angular clustering of 33 415 extremely red objects (EROs) in the Elais-N1 field covering 5.33 deg$^{2}$, which cover the redshift range $z=0.8$ to $2$. This sample was made by merging the UKIDSS Deep eXtragalactic Survey (DXS) with the optical Subaru and Pan-STARRS PS1 datasets. We confirm the existence of a clear break in the angular correlation function at $\sim 0.02^{\circ}$ corresponding to $1 h^{-1}$ Mpc at $z\sim1$. We find that redder or brighter EROs are more clustered than bluer or fainter ones. Halo Occupation Distribution (HOD) model fits imply that the average mass of dark matter haloes which host EROs is over $10^{13} h^{-1} M_{\odot}$ and that EROs have a bias ranging from 2.7 to 3.5. Compared to EROs at $z\sim1.1$, at $z\sim1.5$ EROs have a higher bias and fewer are expected to be satellite galaxies. Furthermore, EROs reside in similar dark matter haloes to those that host $10^{11.0} M_{\odot}<M_{*}<10^{11.5} M_{\odot}$ galaxies. We compare our new measurement and HOD fits with the predictions of the GALFORM semi-analytical galaxy formation model. Overall, the clustering predicted by GALFORM gives an encouraging match to our results. However, compared to our deductions from the measurements, GALFORM puts EROs into lower mass haloes and predicts that a larger fraction of EROs are satellite galaxies. This suggests that the treatment of gas cooling may need to be revised in the model. Our analysis illustrates the potential of clustering analyses to provide observational constraints on theoretical models of galaxy formation.

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.

The Pan-STARRS1 survey is currently obtaining imaging in 5 bands (grizy) for the $3\pi$ steradian survey, one of the largest optical surveys ever conducted. The finished survey will have spatially varying depth, due to the survey strategy. This paper presents a method to correct galaxy number counts and galaxy clustering for this potential systematic based on a simplified signal to noise measurement. A star and galaxy separation method calibrated using realistic synthetic images is also presented, along with an approach to mask bright stars. By using our techniques on a ~69 sq. degree region of science verification data this paper shows PS1 measurements of the two point angular correlation function as a function of apparent magnitude agree with measurements from deeper, smaller surveys. Clustering measurements appear reliable down to a magnitude limit of rps<22.5. Additionally, stellar contamination and false detection issues are discussed and quantified. This work is the second of two papers which pave the way for the exploitation of the full $3\pi$ survey for studies of large scale structure.

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

Wojtak {\it et al} have stacked 7,800 clusters from the SDSS survey in redshift space. They find a small net blue-shift for the cluster galaxies relative to the brightest cluster galaxies, which agrees quite well with the gravitational redshift from GR. Zhao {\it et al.} have pointed out that, in addition to the gravitational redshift, one would expect to see transverse Doppler (TD) redshifts, and that these two effects are generally of the same order. Here we show that there are other corrections that are also of the same order of magnitude. The fact that we observe galaxies on our past light cone results in a bias such that more of the galaxies observed are moving away from us in the frame of the cluster than are moving towards us. This causes the observed average redshift to be $\langle \delta z \rangle = -\langle \Phi \rangle + \langle \beta^2 \rangle / 2 + \langle \beta_x^2 \rangle$, with $\beta_x$ is the line of sight velocity. That is if we average over galaxies with equal weight. If the galaxies in each cluster are weighted by their fluence, or equivalently if we do not resolve the moving sources, and make an average of the mean redshift giving equal weight per photon, the observed redshift is then opposite to the usual transverse Doppler effect. In the WHH experiment, the weighting is a step-function because of the flux-limit for inclusion in the spectroscopic sample and the result is different again, and depends on the details of the luminosity function and the SEDs of the galaxies. Including these effects substantially modifies the blue-shift profile. We show that in-fall and out-flow have very small effect over the relevant range of impact parameters but become important on larger scales.

We have combined optical and NIR photometry from Pan-STARRS 1 and UKIDSS to search the young (5-10 Myr) star-forming region of Upper Scorpius for wide (~400-4000 AU) substellar companions down to ~5 Mjup. Our search is ~4mag deeper than previous work based on 2MASS. We identified several candidates around known stellar members using a combination of color selection and spectral energy distribution fitting. Our followup spectroscopy has identified two new companions as well as confirmed two companions previously identified from photometry, with spectral types of M7.5-M9 and masses of ~15-60 Mjup, indicating a frequency for such wide substellar companions of ~0.6+/-0.3%. Both USco1610-1913B and USco1612-1800B are more luminous than expected for their spectral type compared with known members of Upper Sco. HIP77900B has an extreme mass ratio (M2/M1~0.005) and an extreme separation of 3200 AU. USco1602-2401B also has a very large separation of 1000 AU. We have also confirmed a low-mass stellar companion, USco1610-2502B (730AU, M5.5). Our substellar companions appear both non-coeval with their primary stars according to evolutionary models and, as a group, are systematically more luminous than the Upper Sco cluster sequence. One possible reason for these luminosity discrepancies could be different formation processes or accretion histories for these objects.

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.

(Abridged) In tandem with observational datasets, we utilize realistic mock catalogs, based on a semi-analytic galaxy formation model, constructed specifically for Pan-STARRS1 Medium Deep Surveys in order to assess the performance of the Probability Friends-of-Friends (PFOF, Liu et al.) group finder, and aim to develop a grouping optimization method applicable to surveys like Pan-STARRS1. Producing mock PFOF group catalogs under a variety of photometric redshift accuracies ({\sigma}{\Delta}z/(1+zs)), we find that catalog purities and completenesses from ``good' {\sigma}{\Delta}z/(1+zs)) ~ 0.01) to ``poor' {\sigma}{\Delta}z/(1+zs)) ~ 0.07) photo-zs gradually degrade respectively from 77% and 70% to 52% and 47%. To avoid model dependency of the mock for use on observational data we apply a ``subset optimization' approach, using spectroscopic-redshift group data from the target field to train the group finder for application to that field, as an alternative method for the grouping optimization. We demonstrate this approach using these spectroscopically identified groups as the training set, i.e. zCOSMOS groups for PFOF searches within PS1 Medium Deep Field04 (PS1MD04) and DEEP2 EGS groups for searches in PS1MD07. We ultimately apply PFOF to four datasets spanning the photo-z uncertainty range from 0.01 to 0.06 in order to quantify the dependence of group recovery performance on photo-z accuracy. We find purities and completenesses calculated from observational datasets broadly agree with their mock analogues. Further tests of the PFOF algorithm are performed via matches to X-ray clusters identified within the PS1MD04 and COSMOS footprints. Across over a decade in group mass, we find PFOF groups match ~85% of X-ray clusters in COSMOS and PS1MD04, but at a lower statistical significance in the latter.

[abridged] We present the results of a pilot study for the extended MACS survey (eMACS), a comprehensive search for distant, X-ray luminous galaxy clusters at z>0.5. Our pilot study applies the eMACS concept to the 71 deg^2 area extended by the ten fields of the Pan-STARRS1 (PS1) Medium Deep Survey (MDS). Candidate clusters are identified by visual inspection of PS1 images in the g,r, i, and z bands in a 5x5 arcmin^2 region around X-ray sources detected in the ROSAT All-Sky Survey (RASS). To test and optimize the eMACS X-ray selection criteria, our pilot study uses the largest possible RASS database, i.e., all RASS sources listed in the Bright and Faint Source Catalogs (BSC and FSC) that fall within the MDS footprint. Scrutiny of PS1/MDS images for 41 BSC and 200 FSC sources combined with dedicated spectroscopic follow-up observations results in a sample of 11 clusters with estimated or spectroscopic redshifts of z>0.3. X-ray follow-up observations will be crucial in order to establish robust cluster luminosities for eMACS clusters. Although the small number of distant X-ray luminous clusters in the MDS does not allow us to make firm predictions for the over 20,000 deg^2 of extragalactic sky covered by eMACS, the identification of two extremely promising eMACS cluster candidates at z>0.6 (both yet to be observed with Chandra) in such a small solid angle is encouraging. Representing a tremendous gain over the presently known two dozen such systems from X-ray, optical, and SZ cluster surveys combined, the sample of over 100 extremely massive clusters at z>0.5 expected from eMACS would be invaluable for the identification of the most powerful gravitational lenses in the Universe, as well as for in-depth and statistical studies of the physical properties of the most massive galaxy clusters out to z~1.

As of 2012 Jan 21, the Pan-STARRS 1 3\pi Survey has observed the 3/4 of the sky visible from Hawaii with a minimum of 2 and mean of 7.6 observations in 5 filters, g_p1, r_p1, i_p1, z_p1, y_p1. Now at the end of the second year of the mission, we are in a position to make an initial public release of a portion of this unprecedented dataset. This article describes the PS1 Photometric Ladder, Release 12.01 This is the first of a series of data releases to be generated as the survey coverage increases and the data analysis improves. The Photometric Ladder has rungs every hour in RA and at 4 intervals in declination. We will release updates with increased area coverage (more rungs) from the latest dataset until the PS1 survey and the final re-reduction are completed. The currently released catalog presents photometry of \approx 1000 objects per square degree in the rungs of the ladder. Saturation occurs at g_p1, r_p1, i_p1 \approx 13.5; z_p1 \approx 13.0; and y_p1 \approx 12.0. Photometry is provided for stars down to g_p1, r_p1, i_p1 \approx 19.1 in the AB system. This data release depends on the rigid `Ubercal' photometric calibration using only the photometric nights, with systematic uncertainties of (8.0, 7.0, 9.0, 10.7, 12.4) millimags in (g_p1, r_p1, i_p1, z_p1, y_p1). Areas covered only with lower quality nights are also included, and have been tied to the Ubercal solution via relative photometry; photometric accuracy of the non-photometric regions is lower and should be used with caution.

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

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

We present the selection and classification of over a thousand ultraviolet (UV) variable sources discovered in ~ 40 deg^2 of GALEX Time Domain Survey (TDS) NUV images observed with a cadence of 2 days and a baseline of observations of ~ 3 years. The GALEX TDS fields were designed to be in spatial and temporal coordination with the Pan-STARRS1 Medium Deep Survey, which provides deep optical imaging and simultaneous optical transient detections via image differencing. We characterize the GALEX photometric errors empirically as a function of mean magnitude, and select sources that vary at the 5\sigma level in at least one epoch. We measure the statistical properties of the UV variability, including the structure function on timescales of days and years. We report classifications for the GALEX TDS sample using a combination of optical host colors and morphology, UV light curve characteristics, and matches to archival X-ray, and spectroscopy catalogs. We classify 62% of the sources as active galaxies (358 quasars and 305 active galactic nuclei), and 10% as variable stars (including 37 RR Lyrae, 53 M dwarf flare stars, and 2 cataclysmic variables). The remaining unclassified sources include UV-bright extragalactic transients, two of which have been spectroscopically confirmed to be a young core-collapse supernova and a flare from the tidal disruption of a star by dormant supermassive black hole. We calculate a surface density for variable sources in the UV with NUV < 23 mag and |\Delta m| > 0.2 mag of ~ 8.0, 7.7, and 1.8 deg^-2 for quasars, AGNs, and RR Lyrae stars, respectively. We also calculate a surface density rate in the UV for transient sources, using the effective survey time at the cadence appropriate to each class, of ~ 15 and 52 deg^-2 yr^-1 for M dwarfs and extragalactic transients, respectively. (Abridged.)

We present a sample of Cepheid variable stars towards M31 based on the first year of regular M31 observations of the PS1 survey in the r_P1 and i_P1 filters. We describe the selection procedure for Cepheid variable stars from the overall variable source sample and develop an automatic classification scheme using Fourier decomposition and the location of the instability strip. We find 1440 fundamental mode (classical \delta) Cep stars, 126 Cepheids in the first overtone mode, and 147 belonging to the Population II types. 296 Cepheids could not be assigned to one of these classes and 354 Cepheids were found in other surveys. These 2009 Cepheids constitute the largest Cepheid sample in M31 known so far and the full catalog is presented in this paper. We briefly describe the properties of our sample in its spatial distribution throughout the M31 galaxy, in its age properties, and we derive an apparent period-luminosity relation (PLR) in our two bands. The Population I Cepheids nicely follow the dust pattern of the M31 disk, whereas the 147 Type II Cepheids are distributed throughout the halo of M31. We outline the time evolution of the star formation in the major ring found previously and find an age gradient. A comparison of our PLR to previous results indicates a curvature term in the PLR.

We use data from the Pan-STARRS1 survey to present a panoramic view of the Sagittarius tidal stream in the southern Galactic hemisphere. As a result of the extensive sky coverage of Pan-STARRS1, the southern stream is visible along more than 60 degrees of its orbit, nearly double the length seen by the SDSS. The recently discovered southern bifurcation of the stream is also apparent, with the fainter branch of the stream visible over at least 30 degrees. Using a combination of fitting both the main sequence turn-off and the red clump, we measure the distance to both arms of the stream in the south. We find that the distances to the bright arm of the stream agree very well with the N-body models of Law & Majewski (2010). We also find that the faint arm lies ~5 kpc closer to the Sun than the bright arm, similar to the behavior seen in the northern hemisphere.

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.

The PS2 telescope is the second in an array of wide-field telescopes that is being built for the Panoramic-Survey Telescope and Rapid Response System (Pan-STARRS) on Haleakala. The PS2 design has evolved incrementally based on lessons learned from PS1, but these changes should result in significant improvements in image quality, tracking performance in windy conditions, and reductions in scattered light. The optics for this telescope are finished save for their coatings and the fabrication for the telescope structure itself is well on the way towards completion and installation on-site late this year (2012). The most significant differences between the two telescopes include the following: secondary mirror support changes, improvements in the optical polishing, changes in the optical coatings to improve throughput and decrease ghosting, removal of heat sources inside the mirror cell, expansion of the primary mirror figure control system, changes in the baffle designs, and an improved cable wrap design. This paper gives a description of each of these design changes and discusses the motivations for making them.

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.

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

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

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.

We describe the design and data sample from the DEEP2 Galaxy Redshift Survey, the densest and largest precision-redshift survey of galaxies at z ~ 1 completed to date. The survey has conducted a comprehensive census of massive galaxies, their properties, environments, and large-scale structure down to absolute magnitude M_B = -20 at z ~ 1 via ~90 nights of observation on the DEIMOS spectrograph at Keck Observatory. DEEP2 covers an area of 2.8 deg^2 divided into four separate fields, observed to a limiting apparent magnitude of R_AB=24.1. Objects with z < 0.7 are rejected based on BRI photometry in three of the four DEEP2 fields, allowing galaxies with z > 0.7 to be targeted ~2.5 times more efficiently than in a purely magnitude-limited sample. Approximately sixty percent of eligible targets are chosen for spectroscopy, yielding nearly 53,000 spectra and more than 38,000 reliable redshift measurements. Most of the targets which fail to yield secure redshifts are blue objects that lie beyond z ~ 1.45. The DEIMOS 1200-line/mm grating used for the survey delivers high spectral resolution (R~6000), accurate and secure redshifts, and unique internal kinematic information. Extensive ancillary data are available in the DEEP2 fields, particularly in the Extended Groth Strip, which has evolved into one of the richest multiwavelength regions on the sky. DEEP2 surpasses other deep precision-redshift surveys at z ~ 1 in terms of galaxy numbers, redshift accuracy, sample number density, and amount of spectral information. We also provide an overview of the scientific highlights of the DEEP2 survey thus far. This paper is intended as a handbook for users of the DEEP2 Data Release 4, which includes all DEEP2 spectra and redshifts, as well as for the publicly-available DEEP2 DEIMOS data reduction pipelines. [Abridged]

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.

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 present the discovery of the first high redshift (z > 5.7) quasar from the Panoramic Survey Telescope and Rapid Response System 1 (Pan-STARRS1 or PS1). This quasar was initially detected as an i dropoutout in PS1, confirmed photometrically with the SAO Widefield InfraRed Camera (SWIRC) at Arizona's Multiple Mirror Telescope (MMT) and the Gamma-Ray Burst Optical/Near-Infrared Detector (GROND) at the MPG 2.2 m telescope in La Silla. The quasar was verified spectroscopically with the the MMT Spectrograph, Red Channel and the Cassegrain Twin Spectrograph (TWIN) at the Calar Alto 3.5 m telescope. It has a redshift of 5.73, an AB z magnitude of 19.4, a luminosity of 3.8 x 10^47 erg/s and a black hole mass of 6.9 x 10^9 solar masses. It is a Broad Absorption Line quasar with a prominent Ly-beta peak and a very blue continuum spectrum. This quasar is the first result from the PS1 high redshift quasar search that is projected to discover more than a hundred i dropout quasars, and could potentially find more than 10 z dropout (z > 6.8) quasars.

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 calculate the one-photon loop radiative corrections to charged pion Compton scattering, $\pi^- \gamma \to \pi^- \gamma $. Ultraviolet and infrared divergencies are both treated in dimensional regularization. Analytical expressions for the ${\cal O}(\alpha)$ corrections to the invariant Compton scattering amplitudes, $A(s,u)$ and $B(s,u)$, are presented for 11 classes of contributing one-loop diagrams. Infrared finiteness of the virtual radiative corrections is achieved (in the standard way) by including soft photon radiation below an energy threshold $\lambda$, and its relation to the experimental detection threshold is discussed. We find that the radiative corrections are maximal in backward directions, reaching e.g. -2.4% for a center-of-mass energy of $\sqrt{s}=4m_\pi$ and $\lambda=5 $MeV. Furthermore, we extend our calculation of the radiative corrections by including the leading pion structure effect (at low energies) in form of its electric and magnetic polarizability difference, $\alpha_\pi - \beta_\pi \simeq 6\cdot 10^{-4} $fm$^3$. We find that this structure effect does not change the relative size and angular dependence of the radiative corrections to pion Compton scattering. Our results are particularly relevant for analyzing the COMPASS experiment at CERN which aims at measuring the pion electric and magnetic polarizabilities with high statistics using the Primakoff effect.

Aims: Massive stars in low-metallicity environments may produce exotic explosions such as long-duration gamma-ray bursts and pair-instability supernovae when they die as core-collapse supernovae (CCSNe). Here we determine the feasibility of searching for these CCSNe in metal-poor dwarf galaxies using various survey strategies. Methods: We determine oxygen abundances and star-formation rates for all spectroscopically typed star-forming galaxies in the Sloan Digital Sky Survey, Data Release 5, within z = 0.04. We then estimate the CCSN rates for sub-samples of galaxies with differing upper-metallicity limits. Using Monte-Carlo simulations we then predict the fraction of these CCSNe that we can expect to detect using different survey strategies. Results: Using a single 2m telescope (with a standard CCD camera) search we predict a detection rate of ~1.3 CCSNe/yr in galaxies with metallicities below 12 + log(O/H) < 8.2 which are within a volume that will allow detailed follow-up with 4m and 8m telescopes (z = 0.04). With a network of seven 2m telescopes we estimate ~9.3 CCSNe/yr could be found, although this would require more than 1,000 hrs of telescope time allocated to the network. Within the same radial distance, a volume-limited search in the future Pan-STARRS PS1 all-sky survey should uncover 12.5 CCSNe/yr in low-metallicity galaxies. Over a period of a few years this would allow a detailed comparison of their properties. We then extend our calculations to determine the total numbers of CCSNe that can potentially be found in magnitude-limited surveys with PS1 (24,000/yr, within z < 0.6), PS4 (69,000/yr, within z < 0.8) and LSST (160,000/yr, within z < 0.9) surveys.

The evolution of the B-band galaxy luminosity function is measured using a sample of more than 11,000 galaxies with spectroscopic redshifts from the DEEP2 Redshift Survey. The rest-frame M_B versus U-B color-magnitude diagram of DEEP2 galaxies shows that the color-magnitude bi-modality seen in galaxies locally is still present at redshifts z > 1. Dividing the sample at the trough of this color bimodality into predominantly red and blue galaxies, we find that the luminosity function of each galaxy color type evolves differently. Blue counts tend to shift to brighter magnitudes at constant number density, while the red counts remain largely constant at a fixed absolute magnitude. Using Schechter functions with fixed faint-end slopes we find that M*_B for blue galaxies brightens by ~ 1.3 magnitudes per unit redshift, with no significant evolution in number density. For red galaxies M*_B brightens somewhat less with redshift, while the formal value of phi* declines. When the population of blue galaxies is subdivided into two halves using the rest-frame color as the criterion, the measured evolution of both blue subpopulations is very similar.

The DEEP2 and COMBO-17 surveys are used to study the evolution of the luminosity function of red and blue galaxies to $z \sim 1$. Schechter function fits show that, since $z = 1$, $M^*_B$ dims by $\sim$ 1.3 mag per unit redshift for both color classes, $\phi^*$ of blue galaxies shows little change, while $\phi^*$ for red galaxies has formally nearly quadrupled. At face value, the number density of blue galaxies has remained roughly constant since $ z = 1$, whereas that of red galaxies has been rising. Luminosity densities support both conclusions, but we note that most red-galaxy evolution occurs between our data and local surveys and in our highest redshift bin, where the data are weakest. We discuss the implications of having most red galaxies emerge after $z = 1$ from precursors among the blue population, taking into account the properties of local and distant E/S0s. We suggest a ``mixed'' scenario in which some blue galaxies have their star-formation quenched in gas-rich mergers, migrate to the red sequence with a variety of masses, and merge further on the red sequence in one or more purely stellar mergers. E/S0s of a given mass today will have formed via different routes, in a manner that may help to explain the fundamental plane and other local scaling laws.

We use the first 25% of the DEEP2 Galaxy Redshift Survey spectroscopic data to identify groups and clusters of galaxies in redshift space. The data set contains 8370 galaxies with confirmed redshifts in the range 0.7 <= z <= 1.4, over one square degree on the sky. Groups are identified using an algorithm (the Voronoi-Delaunay Method) that has been shown to accurately reproduce the statistics of groups in simulated DEEP2-like samples. We optimize this algorithm for the DEEP2 survey by applying it to realistic mock galaxy catalogs and assessing the results using a stringent set of criteria for measuring group-finding success, which we develop and describe in detail here. We find in particular that the group-finder can successfully identify ~78% of real groups and that ~79% of the galaxies that are true members of groups can be identified as such. Conversely, we estimate that ~55% of the groups we find can be definitively identified with real groups and that ~46% of the galaxies we place into groups are interloper field galaxies. Most importantly, we find that it is possible to measure the distribution of groups in redshift and velocity dispersion, n(sigma,z), to an accuracy limited by cosmic variance, for dispersions greater than 350km/s. We anticipate that such measurements will allow strong constraints to be placed on the equation of state of the dark energy in the future. Finally, we present the first DEEP2 group catalog, which assigns 32% of the galaxies to 899 distinct groups with two or more members, 153 of which have velocity dispersions above 350km/s. We provide locations, redshifts and properties for this high-dispersion subsample. This catalog represents the largest sample to date of spectroscopically detected groups at z~1.

We present our image processing system for the reduction of optical imaging data from multi-chip cameras. In the framework of the Garching Bonn Deep Survey (GaBoDS; Schirmer et al. 2003) consisting of about 20 square degrees of high-quality data from WFI@MPG/ESO 2.2m, our group developed an imaging pipeline for the homogeneous and efficient processing of this large data set. Having weak gravitational lensing as the main science driver, our algorithms are optimised to produce deep co-added mosaics from individual exposures obtained from empty field observations. However, the modular design of our pipeline allows an easy adaption to different scientific applications. Our system has already been ported to a large variety of optical instruments and its products have been used in various scientific contexts. In this paper we give a thorough description of the algorithms used and a careful evaluation of the accuracies reached. This concerns the removal of the instrumental signature, the astrometric alignment, photometric calibration and the characterisation of final co-added mosaics. In addition we give a more general overview on the image reduction process and comment on observing strategies where they have significant influence on the data quality.

When applied to deep photometric catalogs, the two-point angular correlation function, w(theta), is a sensitive probe of the evolution of galaxy clustering properties. Here we present measurements of w(theta) as a function of I_AB magnitude and (R-I) color to a depth of I_AB=24, using a sample of ~350,000 galaxies covering 5 degrees^2 in total over 5 separate fields. Using redshifts of 3319 galaxies in early DEEP2 Galaxy Redshift Survey data, we construct robust galaxy redshift distributions as a function of I_AB and R_AB magnitude and (R-I) color for galaxies between 0<z<2. We constrain models of the redshift evolution of galaxy clustering and find that significant growth of clustering has occurred from z>1 to z=0. A model in which the comoving scale-length, x_0, evolves linearly with redshift, x_0(z)=x_0(0)(1-Bz), fits the data better than the 'epsilon' model proposed by Groth and Peebles (1977). The clustering properties depend strongly on observed (R-I) color, with both the reddest and bluest galaxies exhibiting large clustering amplitudes and steeper slopes. Different observed (R-I) color ranges are sensitive to very disparate redshift regimes. Red galaxies with (R-I)~1.5 lie in a narrow redshift range centered at z~0.85 and have a comoving scale length of clustering at z=0.85 of x_0=5.0 +/-0.3 Mpc/h. These galaxies have early-type spectra and are likely progenitors of massive local ellipticals. The bluest galaxies with (R-I)~0 appear to be a mix of star-forming galaxies, both relatively local (z~0.3-0.6) dwarfs and bright z>1.4 galaxies, and broad-line AGN. We find that local blue dwarfs are relatively unclustered, with x_0=1.6 +/-0.2 Mpc/h. The z>1.4 blue galaxies have a larger clustering scale-length, x_0>5 Mpc/h.

We measure the two-point correlation function xi(r) using a sample of 2219 galaxies in an area of 0.32 degrees^2 at z=0.7-1.35 from the first season of the DEEP2 Galaxy Redshift Survey. We find that xi(r) can be approximated by a power-law, xi(r)=(r/r_0)^-gamma, on scales 0.1-20 Mpc/h. In a sample with an effective redshift of z_eff=0.82, for a Lcdm cosmology we find r_0=3.53 +/-0.81 Mpc/h (comoving) and gamma=1.66 +/-0.12, while in a higher-redshift sample with z_eff=1.14 we find r_0=3.14 +/-0.72 Mpc/h and gamma=1.61 +/-0.11. We find that red, absorption-dominated, passively-evolving galaxies have a larger clustering scale length, r_0, and more prominent ``fingers of God'' than blue, emission-line, actively star-forming galaxies. Intrinsically brighter galaxies also cluster more strongly than fainter galaxies at z~1, with a significant luminosity-bias seen for galaxies fainter than M*. Our results are suggestive of evolution in the galaxy clustering within our survey volume and imply that the DEEP2 galaxies, with a median brightness one magnitude fainter than M* have an effective bias b=0.97 +/-0.13 if sigma_{8 DM}=1 today or b=1.20 +/-0.16 if sigma_{8 DM}=0.8 today. Given the strong luminosity-dependence in the bias that we measure at z~1, the galaxy bias at M* may be significantly greater. We note that our star-forming sample at z~1 has very similar selection criteria as the Lyman-break galaxies at z~3 and that our red, absorption-line sample displays a clustering strength comparable to the expected clustering of the Lyman-break galaxy descendants at z~1. Our results demonstrate that the clustering properties in the galaxy distribution seen in the local Universe were largely in place by z~1.

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 a Principal Component Analysis (PCA)-based spectral classification, eta, for the first 5600 galaxies observed in the DEEP2 Redshift Survey. This parameter provides a very pronounced separation between absorption and emission dominated galaxy spectra - corresponding to passively evolving and actively star-forming galaxies in the survey respectively. In addition it is shown that despite the high resolution of the observed spectra, this parameter alone can be used to quite accurately reconstruct any given galaxy spectrum, suggesting there are not many `degrees of freedom' in the observed spectra of this galaxy population. It is argued that this form of classification, eta, will be particularly valuable in making future comparisons between high and low-redshift galaxy surveys for which very large spectroscopic samples are now readily available, particularly when used in conjunction with high-resolution spectral synthesis models which will be made public in the near future. We also discuss the relative advantages of this approach to distant galaxy classification compared to other methods such as colors and morphologies. Finally, we compare the classification derived here with that adopted for the 2dF Galaxy Redshift Survey and in so doing show that the two systems are very similar. This will be particularly useful in subsequent analyses when making comparisons between results from each of these surveys to study evolution in the galaxy populations and large-scale structure.

We use weak lensing shear measurements of six z>0.5 clusters of galaxies to derive the mean lensing redshift of the background galaxies used to measure the shear. Five of these clusters are compared to X-ray mass models and verify a mean lensing redshift for a 23<R<26.3, R-I<0.9 background galaxy population in good agreement with photometric redshift surveys of the HDF-S. The lensing strength of the six clusters is also analyzed as a function of the magnitude of the background galaxies, and an increase in shear with increasing magnitude is detected at moderate significance. The change in the strength of the shear is presumed to be caused by an increase in the mean redshift of the background galaxies with increasing magnitude, and the degree of change detected is also in agreement with those in photometric redshift surveys of the HDF-S.

In the course of a weak gravitational lensing survey of 39 clusters of galaxies,covering a total sky area of ~1 square degree, we have serendipitously discovered mass concentrations in the fields of A1705 and A1722 which are most probably not associated with the main cluster target. By combining weak lensing information with two-color galaxy photometry in fields centered on our sample clusters, we identify a new cluster candidate at z~0.5 in the field of A1705. This cluster candidate also displays strong lensing in the form of a giant luminous arc. The new mass concentration in the field of A1722 also seems to be associated with an optically luminous cluster of galaxies at z~0.5, but in this case there is some evidence for additional structures along the line of sight that may contribute to the lensing signal. A third cluster, A959, has a dark sub-clump which shows interesting morphological evidence in the mass map for being associated with the main cluster. This is the first case where there is any significant evidence for a physical association between a dark sub-clump (discovered from weak lensing) and a normal cluster. Analysis of archival X-ray data shows that the three new mass concentrations are not firmly detected in X-rays and that they are X-ray underluminous.

We introduce the theory of non-linear cosmological perturbations using the correspondence limit of the Schr\"odinger equation. The resulting formalism is equivalent to using the collisionless Boltzman (or Vlasov) equations which remain valid during the whole evolution, even after shell crossing. Other formulations of perturbation theory explicitly break down at shell crossing, e.g. Eulerean perturbation theory, which describes gravitational collapse in the fluid limit. This paper lays the groundwork by introducing the new formalism, calculating the perturbation theory kernels which form the basis of all subsequent calculations. We also establish the connection with conventional perturbation theories, by showing that third order tree level results, such as bispectrum, skewness, cumulant correlators, three-point function are exactly reproduced in the appropriate expansion of our results. We explicitly show that cumulants up to N=5 predicted by Eulerian perturbation theory for the dark matter field $\delta$ are exactly recovered in the corresponding limit. A logarithmic mapping of the field naturally arises in the Schr\"odinger context, which means that tree level perturbation theory translates into (possibly incomplete) loop corrections for the conventional perturbation theory. We show that the first loop correction for the variance is $\sigma^2 = \sigma_L^2+ (-1.14+n)\sigma_L^4$ for a field with spectral index $n$. This yields 1.86 and 0.86 for $n=-3,-2$ respectively, and to be compared with the exact loop order corrections 1.82, and 0.88. Thus our tree-level theory recovers the dominant part of first order loop corrections of the conventional theory, while including (partial) loop corrections to infinite order in terms of $\delta$.

The Galactic Exoplanet Survey Telescope (GEST) will observe a 2 square degree field in the Galactic bulge to search for extra-solar planets using a gravitational lensing technique. This gravitational lensing technique is the only method employing currently available technology that can detect Earth-mass planets at high signal-to-noise, and can measure the frequency of terrestrial planets as a function of Galactic position. GEST's sensitivity extends down to the mass of Mars, and it can detect hundreds of terrestrial planets with semi-major axes ranging from 0.7 AU to infinity. GEST will be the first truly comprehensive survey of the Galaxy for planets like those in our own Solar System.

The DEIMOS spectrograph has now been installed on the Keck-II telescope and commissioning is nearly complete. The DEEP2 Redshift Survey, which will take approximately 120 nights at the Keck Observatory over a three year period and has been designed to utilize the power of DEIMOS, began in the summer of 2002. The multiplexing power and high efficiency of DEIMOS enables us to target 1000 faint galaxies per clear night. Our goal is to gather high-quality spectra of \~60,000 galaxies with z>0.75 in order to study the properties and large scale clustering of galaxies at z ~ 1. The survey will be executed at high spectral resolution, $R=\lambda/\Delta \lambda \approx 5000$, allowing us to work between the bright OH sky emission lines and to infer linewidths for many of the target galaxies (for several thousand objects, we will obtain rotation curves as well). The linewidth data will facilitate the execution of the classical redshift-volume cosmological test, which can provide a precision measurement of the equation of state of the Universe. This talk reviews the project, summarizes our science goals and presents some early DEIMOS data.

We present measurements of the extended dark halo profiles of bright early type galaxies at redshifts 0.1 to 0.9 obtained via galaxy-galaxy lensing analysis of images taken at the CFHT using the UH8K CCD mosaic camera. Six half degree fields were observed for a total of 2 hours each in I and V, resulting in catalogs containing ~20 000 galaxies per field. We used V-I color and I magnitude to select bright early type galaxies as the lens galaxies, yielding a sample of massive lenses with fairly well determined redshifts and absolute magnitudes M ~ M_* pm 1. We paired these with faint galaxies lying at angular distances 20" to 60", corresponding to physical radii of 26 to 77 kpc (z = 0.1) and 105 to 315 kpc (z = 0.9), and computed the mean tangential shear of the faint galaxies. The shear falls off with radius roughly as expected for flat rotation curve halos. The shear values were weighted in proportion to the square root of the luminosity of the lens galaxy. Our results give a value for the average mean rotation velocity of an L_* galaxy halo at r~50-200 kpc of v_* = 238^{+27}_{-30} km per sec for a flat lambda (Omega_m0 = 0.3, Omega_l0 = 0.7) cosmology (v_* = 269^{+34}_{-39} km per sec for Einstein-de Sitter), and with little evidence for evolution with redshift. We find a mass-to-light ratio of \~121\pm28h(r/100 kpc) and these halos constitute Omega ~0.04 pm 0.01 (r/100 kpc) of closure density.

We present measurements of the extended dark halo profiles of bright early type galaxies at redshifts 0.1 to 0.9 obtained via galaxy-galaxy lensing analysis of images taken at the CFHT using the UH8K CCD mosaic camera. Six half degree fields were observed for a total of 2 hours each in I and V, resulting in catalogs containing ~20 000 galaxies per field. We used V-I color and I magnitude to select bright early type galaxies as the lens galaxies, yielding a sample of massive lenses with fairly well determined redshifts and absolute magnitudes M ~ M_* \pm 1. We paired these with faint galaxies lying at angular distances 20" to 60", corresponding to physical radii of 26 to 77 kpc (z = 0.1) and 105 to 315 kpc (z = 0.9), and computed the mean tangential shear of the faint galaxies. The shear falls off with radius roughly as expected for flat rotation curve halos. The shear values were weighted in proportion to the square root of the luminosity of the lens galaxy. Our results give a value for the average mean rotation velocity of an L_* galaxy halo at r~50-200 kpc of v_* = 238^{+27}_{-30} km per sec for a flat lambda (Omega_m0 = 0.3, Omega_l0 = 0.7) cosmology (v_* = 269^{+34}_{-39} km per sec for Einstein-de Sitter), and with little evidence for evolution with redshift. We compare to halo masses measured by other groups/techniques. We find a mass-to-light ratio of ~121\pm28h(r/100 kpc) and these halos constitute Omega ~0.04 \pm 0.01(r/100 kpc) of closure density. (abridged)

We present a weak lensing and photometric study of six half by half degree fields observed at the CFHT using the UH8K CCD mosaic camera. The fields were observed for a total of 2 hours each in I and V, resulting in catalogs containing ~ 20 000 galaxies per passband per field. We use V-I color and I magnitude to select bright early type galaxies at redshifts 0.1 < z < 0.9. We measure the gravitational shear from faint galaxies in the range 21 < m_I < 25 from a composite catalog and find a strong correlation with that predicted from the early types if they trace the mass with mass-to-light ratio 300\pm75 h (in solar units) for a flat (Omega_m0 = 0.3, Omega_l0 = 0.7) lambda cosmology and 400\pm100 h for Einstein-de Sitter. We make two-dimensional reconstructions of the mass surface density. Cross-correlation of the measured mass surface density with that predicted from the early type galaxy distribution shows a strong peak at zero lag (significant at the 5.2-sigma level). We azimuthally average the cross- and auto-correlation functions. We conclude that the profiles are consistent with early type galaxies tracing mass on scales of > 45 arcsec (> 200 kpc at z = 0.5). We sub-divide our bright early type galaxies by redshift and obtain similar conclusions. These mass-to-light ratios imply \Omega_m0 = 0.10\pm0.02 (\Omega_m0 = 0.13\pm0.03 for Einstein-de Sitter) of closure density.

We present estimates of the gravitational lensing shear variance obtained from images taken at the CFHT using the UH8K CCD mosaic camera. Six fields were observed for a total of 1 hour each in V and I, resulting in catalogs containing approximately 20,000$ galaxies per field, with properly calibrated and optimally weighted shear estimates. These were averaged in cells of sizes ranging from 1'.875 to 30' to obtain estimates of the cosmic shear variance $<\bar{\gamma}^2>$, with uncertainty estimated from the scatter among the estimates for the 6 fields. Our most reliable estimator for cosmic shear is provided by the cross-correlation of the shear measured in the two passbands. At scales $\lsim 10'$ the results are in good agreement with those of van Waerbeke et al. 2000, Bacon et al. 2000 and Wittman et al. 2000 and with currently fashionable cosmological models. At larger scales the shear variance falls below the theoretical predictions, and on the largest scales we find a null detection of shear variance averaged in 30' cells of $<\bar{\gamma}^2 > = (0.28 \pm 1.84) \times 10^{-5}$.

We present the results of a weak lensing survey of six high-redshift (z > 0.5), X-ray selected clusters of galaxies. We have obtained ultra-deep R-band images of each cluster with the Keck Telescope, and have measured a weak lensing signal from each cluster. From the background galaxy ellipticities we create two-dimensional maps of the surface mass density of each cluster. We find that the substructure seen in the mass reconstructions typically agree well with substructure in both the cluster galaxy distributions and X-ray images of the clusters. We also measure the one-dimensional radial profiles of the lensing signals and fit these with both isothermal spheres and "universal" CDM profiles. We find that the more massive clusters are less compact and not as well fit by isothermal spheres as the less massive clusters, possibly indicating that they are still in the process of collapse.

In this talk I review some of the key questions that weak lensing observations of clusters can potentially answer, and sketch the progress that has been made to date in extracting quantitative estimates of masses and density profiles. A major difficulty in interpreting these measurements is the diverse methods that have been used to calibrate the shear-polarization relation, none of which are completely satisfactory. I describe some recent developments in the theory of weak lensing calibration. I then present some results from the UH8K camera and discuss their implication for cosmology. Finally, I outline a new strategy for obtaining high resolution optical imaging from the ground using an array of small telescopes, and describe how such an instrument can greatly improve the power of weak lensing observations

We propose a new strategy for obtaining enhanced resolution (FWHM = 0.12 arcsec) deep optical images over a wide field of view. As is well known, this type of image quality can be obtained in principle simply by fast guiding on a small (D = 1.5m) telescope at a good site, but only for target objects which lie within a limited angular distance of a suitably bright guide star. For high altitude turbulence this 'isokinetic angle' is approximately 1 arcminute. With a 1 degree field say one would need to track and correct the motions of thousands of isokinetic patches, yet there are typically too few sufficiently bright guide stars to provide the necessary guiding information. Our proposed solution to these problems has two novel features. The first is to use orthogonal transfer charge-coupled device (OTCCD) technology to effectively implement a wide field 'rubber focal plane' detector composed of an array of cells which can be guided independently. The second is to combine measured motions of a set of guide stars made with an array of telescopes to provide the extra information needed to fully determine the deflection field. We discuss the performance, feasibility and design constraints on a system which would provide the collecting area equivalent to a single 9m telescope, a 1 degree square field and 0.12 arcsec FWHM image quality.

We present results of a weak gravitational lensing survey of six X-ray selected high-redshift clusters of galaxies. We find that the masses of the clusters derived from weak lensing are comparable to those derived from the X-ray observations. We show that many of the clusters have significant substructure not observed in the X-ray observations and that for the more massive clusters a singular isothermal sphere does not provide a good fit to the radial mass profile.

We describe in detail the processing of a set of images of the z = 0.42 supercluster MS0302 taken with the UH8K camera at CFHT. The result of this is a pair of seamless combined V- and I-band images of the field, along with a characterization of the noise properties and of the point spread function (PSF), and catalogs of about 30,000 faint galaxies. The analysis involves the following steps: image preparation; detection of stars and registration to find the transformation from detector to sky coordinates; correction for extinction and/or gain variations; modeling of the PSF; generation of images with a circular PSF; image warping and averaging; modeling of the noise auto-correlation function; faint object detection, aperture photometry, and shape measurement. The shear analysis is described elsewhere.

We present a new shear estimator for weak lensing observations which properly accounts for the effects of a realistic point spread function (PSF). Images of faint galaxies are subject to gravitational shearing followed by smearing with the instrumental and/or atmospheric PSF. We construct a `finite resolution shear operator' which when applied to an observed image has the same effect as a gravitational shear applied prior to smearing. This operator allows one to calibrate essentially any shear estimator. We then specialize to the case of weighted second moment shear estimators. We compute the shear polarizability which gives the response of an individual galaxy's polarization to a gravitational shear. We then compute the response of the population of galaxies, and thereby construct an optimal weighting scheme for combining shear estimates from galaxies of various shapes, luminosities and sizes. We define a figure of merit --- an inverse shear variance per unit solid angle --- which characterizes the quality of image data for shear measurement. The new method is tested with simulated image data. We discuss the correction for anisotropy of the PSF and propose a new technique involving measuring shapes from images which have been convolved with a re-circularizing PSF. We draw attention to a hitherto ignored noise related bias and show how this can be analyzed and corrected for. The analysis here draws heavily on the properties of real PSF's and we include as an appendix a brief review, highlighting those aspects which are relevant for weak lensing.

We describe the 2dF Galaxy Redshift Survey (2dFGRS), and the current status of the observations. In this exploratory paper, we apply a Principal Component Analysis to a preliminary sample of 5869 galaxy spectra and use the two most significant components to split the sample into five spectral classes. These classes are defined by considering visual classifications of a subset of the 2dF spectra, and also by comparing to high quality spectra of local galaxies. We calculate a luminosity function for each of the different classes and find that later-type galaxies have a fainter characteristic magnitude, and a steeper faint-end slope. For the whole sample we find M*=-19.7 (for Omega=1, H_0=100 km/sec/Mpc), alpha=-1.3, phi*=0.017. For class 1 (`early-type') we find M*=-19.6, alpha=-0.7, while for class 5 (`late-type') we find M*=-19.0, alpha=-1.7. The derived 2dF luminosity functions agree well with other recent luminosity function estimates.

In this introductory talk I give an elementary overview of the development, current status, and near term outlook for three probes of large scale structure; galaxy clustering, bulk flows and weak lensing observations.

We perform a weak lensing and photometric study of the z=0.42 supercluster MS0302+17 using deep I and V band images taken with the UH8K CCD mosaic camera at the CFHT. We use archival ROSAT HRI data to estimate fluxes, gas masses and, in one case, the binding mass of the three major clusters. We then use our CCD data to determine the optical richness and luminosities of the clusters and to map out the spatial distribution of the early type galaxies in the supercluster and in other foreground and background structures. We measure the gravitational shear from a sample of \simeq 30,000 faint background galaxies in the range 22 < m_I < 26 and find this correlates strongly with that predicted from the early type galaxies if they trace the mass with M/L_B \simeq 250 h. We make 2-dimensional reconstructions of the mass surface density. These recover all of the major concentrations of galaxies and indicate that most of the supercluster mass, like the early type galaxies, is concentrated in the three X-ray clusters, and we obtain mean mass-to-light ratios for the clusters of M/L_B \simeq 260 h. Cross-correlation of the measured mass surface density with that predicted from the early type galaxy distribution shows a strong peak at zero lag (significant at the \simeq 9-sigma level), and that at separations \gsim 200 h^{-1}$kpc the early galaxies trace the mass very accurately. This conclusion is supported by cross-correlation in Fourier space; we see little evidence for any variation of M/L or `bias' with scale, and from the longest wavelength modes with \lambda = 1.5-6 h^{-1}Mpc we find M/L \simeq (280 \pm 40)h, quite similar to that obtained for the cluster centers. We discuss the implication of these results for the cosmological density parameter.

A tangential distortion of background source galaxies around foreground lens galaxies in the Hubble Deep Field is detected at the 99.3% confidence level. An important element of our analysis is the use of photometric redshifts to determine distances of lens and source galaxies and rest-frame B-band luminosities of the lens galaxies. The lens galaxy halos obey a Tully-Fisher relation between halo circular velocity and luminosity; the typical lens galaxy, at a redshift z = 0.6, has a circular velocity of 210 +/-40 km/s at M_B = -18.5, if q_0 = 0.5. Control tests, in which lens and source positions and source ellipticities are randomized, confirm the significance level of the detection quoted above. Furthermore, a marginal signal is also detected from an independent, fainter sample of source galaxies without photometric redshifts. Potential systematic effects, such as contamination by aligned satellite galaxies, the distortion of source shapes by the light of the foreground galaxies, PSF anisotropies, and contributions from mass distributed on the scale of galaxy groups are shown to be negligible. A comparison of our result with the local Tully-Fisher relation indicates that intermediate-redshift galaxies are fainter than local spirals by 1.0 +/- 0.6 B mag at a fixed circular velocity. This is consistent with some spectroscopic studies of the rotation curves of intermediate-redshift galaxies. This result suggests that the strong increase in the global luminosity density with redshift is dominated by evolution in the galaxy number density.

We report the detection of weak gravitational lensing of distant background galaxies by two rich, X-ray luminous clusters of galaxies, MS 1137+66 at z=0.783 and RXJ 1716+67 at z=0.813. We detect a shear signal in each cluster over the radial range 100<r<700 h^-1 kpc. Assuming the background galaxies lie in a sheet at z=2, we measure a mass of 2.45*10^14 h^-1 Msolar and M/Lv = 270 at 500 h^-1 kpc radius for MS 1137+66, and a mass of 2.6*10^14 h^-1 Msolar and M/Lv = 190 for RXJ 1716+67 at the same radius. Both the optical lightmaps and weak lensing massmaps of RXJ 1716+67 show two spatially distinct sub-clusters, as well as a long filamentary structure extending out of the cluster to the northeast. In contrast, MS 1137+66 is an ultra-compact massive cluster in both massmaps and optical lightmaps, and contains a strongly-lensed arc system in the cluster core. These data add to the growing number of massive clusters at z>~0.8.

We explore the dependence of weak lensing phenomena on the background cosmology. We first generalise the relation between $P_\psi(\omega)$, the angular power spectrum of the distortion, and the power spectrum of density fluctuations to non-flat cosmologies. We then compute $P_\psi$ for various illustrative models. A useful cosmological discriminator is the growth of $P_\psi$ with source redshift which is much stronger in low matter density models, and especially in $\Lambda$-dominated models. With even crude redshift information (say from broad band colours) it should be possible to constrain the cosmological world model. The amplitude of $P_\psi(\omega)$ is also quite sensitive to the cosmology, but requires a reliable external normalisation for the mass fluctuations. If one normalises to galaxy clustering, with $M/L$ fixed by small-scale galaxy dynamics, then low density models predict a much stronger distortion. If, however, one normalises to large-scale bulk-flows, the predicted distortion for sources at redshifts $Z_s \sim 1-3$ is rather insensitive to the background cosmology. The signals predicted here can be detected at a very high level of significance with a photometric survey covering say 10 square degrees, but sparse sampling is needed to avoid large sampling variance and we discuss the factors influencing the design of an optimum survey. Turning to weak lensing by clusters we find that for high lens redshifts ($Z_l\simeq1$) the critical density is substantially reduced in $\Lambda$ models, but that the ratio of the shear or convergence to the velocity dispersions or X-ray temperature of clusters is only very weakly dependent on the cosmology.

We describe the statistical properties of light rays propagating though a random sea of gravity waves and compare with the case for scalar metric perturbations from density inhomogeneities. For scalar fluctuations the deflection angle grows as the square-root of the path length $D$ in the manner of a random walk, and the rms displacement of a ray from the unperturbed trajectory grows as $D^{3/2}$. For gravity waves the situation is very different. The mean square deflection angle remains finite and is dominated by the effect of the metric fluctuations at the ends of the ray, and the mean square displacement grows only as the logarithm of the path length. In terms of power spectra, the displacement for scalar perturbations has $P(k) \propto 1/ k^4$ while for gravity waves the trajectories of photons have $P(k) \propto 1/k$ which is a scale-invariant or `flicker-noise' process, and departures from rectilinear motion are suppressed, relative to the scalar case, by a factor $\sim (\lambda / D)^{3/2}$ where $\lambda$ is the characteristic scale of the metric fluctuations and $D$ is the path length. This result casts doubt on the viability of some recent proposals for detecting or constraining the gravity wave background by astronomical measurements.

In Paper I we studied the theory of gravitational microlensing for a planar distribution of point masses. In this second paper, we extend the analysis to a three-dimensional lens distribution. First we study the lensing properties of three-dimensional lens distributions by considering in detail the critical curves, the caustics, the illumination patterns, and the magnification cross-sections sigma(A) of multiplane configurations with two, three, and four point masses. For N_\ast point masses that are widely separated in Lagrangian space, we find that there are ~ 2^{N_\ast} - 1 critical curves in total, but that only ~ N_\ast of these produce prominent caustic-induced features at the high magnification end of sigma(A). In the case of a low optical depth random distribution of point masses, we show that the multiplane lens equation near a point mass can be reduced to the single plane equation of a point mass perturbed by weak shear. This allows us to calculate the caustic-induced feature in the macroimage magnification distribution P(A) as a weighted sum of the semi-analytic feature derived in Paper I for a planar lens distribution. The resulting semi-analytic caustic-induced feature is similar to the feature in the planar case, but it does not have any simple scaling properties, and it is shifted to higher magnification. The semi-analytic distribution is compared to the results of previous numerical simulations for optical depth tau ~ 0.1, and they are in better agreement than a similar comparison in the planar case. We explain this by estimating the fraction of caustics of individual lenses that merge with those of their neighbors. For tau=0.1, the fraction is approx 20%, much less than the approx 55% for the planar case.

(Abridged) In this paper we refine the theory of microlensing for a planar distribution of point masses. We derive the macroimage magnification distribution P(A) at high magnification (A-1 >> tau^2) for a low optical depth (tau << 1) lens distribution by modeling the illumination pattern as a superposition of the patterns due to individual ``point mass plus weak shear'' lenses. We show that a point mass plus weak shear lens produces an astroid- shaped caustic and that the magnification cross-section obeys a simple scaling property. By convolving this cross-section with the shear distribution, we obtain a caustic-induced feature in P(A) which also exhibits a simple scaling property. This feature results in a 20% enhancement in P(A) at A approx 2/tau. In the low magnification (A-1 << 1) limit, the macroimage consists of a bright primary image and a large number of faint secondary images formed close to each of the point masses. Taking into account the correlations between the primary and secondary images, we derive P(A) for low A. The low-A distribution has a peak of amplitude ~ 1/tau^2 at A-1 ~ tau^2 and matches smoothly to the high-A distribution. We combine the high- and low-A results and obtain a practical semi-analytic expression for P(A). This semi-analytic distribution is in qualitative agreement with previous numerical results, but the latter show stronger caustic-induced features at moderate A for tau as small as 0.1. We resolve this discrepancy by re-examining the criterion for low optical depth. A simple argument shows that the fraction of caustics of individual lenses that merge with those of their neighbors is approx 1-exp(-8 tau). For tau=0.1, the fraction is surprisingly high: approx 55%. For the purpose of computing P(A) in the manner we did, low optical depth corresponds to tau << 1/8.

We report on the detection of dark matter in the cluster of galaxies Abell~2163 using the weak gravitational distortion of background galaxies, and an analysis of the cluster X-ray emission. We find that while the qualitative distributions of the cluster light and the dark matter are similar -- shallow and extended, with significant substructure -- the X-ray morphology shows a more regular overall appearance. We interpret the joint lensing and X-ray observations as a signature of a merger event in the cluster. We present new ROSAT/HRI data and reanalyze ROSAT/PSPC data, accounting for the effect of a varying background to determine the best fit parameters in the $\beta$-model formalism. We combine the surface brightness fits with two determinations of the radial temperature profile to determine the total mass. Although there are slight variations in the total mass determinations introduced by the uncertainties in the $\beta$-fit, the main contributor to the error arises from the uncertainties in the temperature determinations. Even though the morphologies of the dark matter/light and X-ray gas are quite different, we find that the total mass determined from the X-ray and weak lensing estimates are consistent with each other within the $2\sigma$ error bars, with the X-ray inferred mass a factor of $\simeq 2$ larger. However, as the lensing mass estimates are differential (the surface density at any point is determined relative to the mean in a control annulus), the shallow, extended nature of the mass profile biases the lensing inferred mass downwards. We estimate the correction for this effect and find very good agreement between the corrected lensing and X-ray results. We determine the gas mass fraction and find $f_g \simeq 0.07h^{-3/2}$ at all radii and a constant mass-to-light ratio of $M/L_V

We report on the detection of dark matter in the cluster Abell 2390 using the weak gravitational distortion of background galaxies. We find that the cluster light and total mass distributions are quite similar over an angular scale of $\simeq 7^\prime \;(1 \Mpc$). The cluster galaxy and mass distributions are centered on the cluster cD galaxy and exhibit elliptical isocontours in the central $\simeq 2^\prime \; (280 \kpc)$. The major axis of the ellipticity is aligned with the direction defined by the cluster cD and a ``straight arc'' located $\simeq 38^{\prime\prime}$ to the northwest. We determined the radial mass-to-light profile for this cluster and found a constant value of $(320 \pm 90) h\; M_\odot/L_{\odot V}$, which is consistent with other published determinations. We also compared our weak lensing azimuthally averaged radial mass profile with a spherical mass model proposed by the CNOC group on the basis of their detailed dynamical study of the cluster. We find good agreement between the two profiles, although there are weak indications that the CNOC density profile may be falling more steeply for $\theta\geq 3^\prime$ $(420\kpc)$.

We report the detection of weak gravitational lensing of faint, distant background galaxies by the rich, X-ray luminous cluster of galaxies MS1054-03 at z=0.83. This is the first measurement of weak lensing by a bona fide cluster at such a high redshift. We detect tangential shear at the 5% - 10% level over a range of radii 50'' < r < 250'' centered on the optical position of the cluster. Two-dimensional mass reconstruction using galaxies with 21.5 < I < 25.5 shows a strong peak which coincides with the peak of the smoothed cluster light distribution. Splitting this sample by magnitude (at I = 23.5) and color (at R-I = 0.7), we find that the brighter and redder subsamples are only very weakly distorted, indicating that the faint blue galaxies (FBG's), which dominate the shear signal, are relatively more distant. The derived cluster mass is quite sensitive to the N(z) for the FBG's. At one extreme, if all the FBG's are at z_s = 3, then the mass within a $0.5h^{-1}$Mpc aperture is $(5.9 \pm 1.24)\times 10^{14}$\h1 $M_\odot$, and the mass-to-light ratio is $M/L_V = 350 \pm 70 h$ in solar units. For $z_s = 1.5$ the derived mass is $\sim$70\% higher and $M/L \simeq 580 h$. If $N(z)$ follows the no evolution model (in shape) then $M/L \simeq 800h$, and if all the FBG's lie at $z_s\la 1$ the required $M/L$ exceeds $1600h$. These data provide clear evidence that large, dense mass concentrations existed at early epochs; that they can be weighed efficiently by weak lensing observations; and that most of the FBG's are at high redshift.

Weak lensing observations measure the shear field and hence the gradient of the dimensionless surface density $\kappa$. We present several new algorithms to recover $\kappa$ {}from shear estimates on a finite region and compare how they perform with realistically noisy data. The reconstruction methods studied here are divided into 2 classes: direct reconstruction and regularized inversion. Direct reconstruction techniques express $\kappa$ as a 2D integral of the shear field. This yields an estimator for $\kappa$ as a discrete sum over background galaxy ellipticities which is straightforward to implement, and allows a simple estimate of the noise. We study 3 types of direct reconstruction methods: 1) $\kappa$-estimators that measure the surface density at any given target point relative to the mean value in a reference region 2) a method that explicitly attempts to minimize the rotational part of $\nabla \kappa$ that is due to noise and 3) a novel, exact Fourier-space inverse gradient operator. We also develop two `regularized maximum likelihood' methods, one of which employs the conventional discrete Laplacian operator as a regularizer and the other uses regularization of all components in Fourier space. We compare the performance of all the estimators by means of simulations and noise power analysis. A general feature of these unbiased methods is an enhancement of the low frequency power which, for some of the methods, can be quite severe. We find the best performance is provided by the maximum likelihood method with Fourier space regularization although some of the other methods perform almost as well.

In this review I will describe progress that has been made in determining masses of galaxy clusters using `weak lensing' and how this technique my be applied in the future to determine the dark matter distribution both on supercluster scales and on the scale of galaxy haloes.

We apply a likelihood analysis to the data of \markcite{Lauer \& Postman 1994} With $P(k)$ parametrized by $(\sigma_8, \Gamma)$, the likelihood function peaks at $\sigma_8\simeq0.3$, $\Gamma\lesssim0.025$, indicating at face value very strong large-scale power, though at a level incompatible with COBE\@. There is, however, a ridge of likelihood such that more conventional power spectra do not seem strongly disfavored. The likelihood calculated using as data only the components of the bulk flow solution peaks at higher $\sigma_8$, in agreement with other analyses, but is rather broad. The likelihood incorporating both bulk flow and shear gives a different picture. The components of the shear are all low, and this pulls the peak to lower amplitudes as a compromise. The Lauer \& Postman velocity data alone are therefore {\em consistent}\/ with models with very strong large scale power which generates a large bulk flow, but the small shear (which also probes fairly large scales) requires that the power would have to be at {\em very}\/ large scales, which is strongly disfavored by COBE\@. The velocity data also seem compatible with more conventional $P(k)$ with $0.2\lesssim\Gamma\lesssim0.5$, and the likelihood is peaked around $\sigma_8\sim1$, in which case the bulk flow is a moderate, but not extreme, statistical fluctuation. Applying the same techniques to the data of \markcite{Riess, Press, \& Kirshner 1995}, the results are quite different. The flow is not inconsistent with the microwave dipole and we derive only an upper limit to the amplitude of the power spectrum: $\sigma_8\lesssim1.5$ at roughly 99\%.

We report on the detection of dark matter in the cluster Abell 2218 using the weak gravitational distortion of background galaxies. We find a highly significant, coherent detection of the distortion in the images of the background galaxies. The inferred 2D mass distribution has a peak that is coincident with the optical and X-ray centroid. The qualitative distributions of the cluster light, the X-ray emission and the dark matter are similar and the projected total mass, gas, and light surface densities are consistent with a $r^{-1}$ profile at distance of $r > 180^{\prime\prime}$ from the cluster cD galaxy. Using the weak lensing technique, we determine a lower bound for the total mass in A2218 of $(3.9 \pm 0.7) \times 10^{14}$~h$^{-1}$~M$_\odot$ within a fiducial aperture of radius 0.4~h$^{-1}$Mpc. The associated cluster mass-to-light ratio is $(440 \pm 80)$~h~$M_\odot/L_{\odot B}$. The mass estimated by the weak lensing method is consistent with that inferred from the X-ray data under the assumption of hydrostatic equilibrium and we derive an upper bound for the gas-to-total mass ratio at 400~h$^{-1}$kpc of $M_{gas}/M_{tot} = (0.04 \pm 0.02)$h$^{-3/2}$.

Measurement of the gravitational distortion of images of distant galaxies is rapidly becoming established as a powerful probe of the dark mass distribution in clusters of galaxies. With the advent of large mosaics of CCD's these methods should provide a composite total mass profile for galaxy haloes and also measure the power spectrum of mass fluctuations on supercluster scales. We describe how HST observations have been used to place the observational measurement of the shear on a quantitative footing. By artifically stretching and then degrading WFPC2 images to simulate ground based observing it is now possible to directly calibrate the effect of atmospheric seeing. Similar experiments show that one can remove the effect of artificial image anisotropy arising in the atmosphere or telescope. There have also been important advances in the interpretation of the shear: reconstruction techniques have been extended to encompass the strong distortion regime of giant arcs etc., progress has been made in removing a bias present in earlier reconstruction techniques, and we describe new techniques for `aperture densitometry'. We present some new results on clusters of galaxies, and discuss the intimate connections between weak lensing and deep spectroscopy.

We develop and test a method for measuring the gravitational lensing induced distortion of faint background galaxies. We first describe how we locate the galaxies and measure a 2-component `polarisation' or ellipticity statistic $e_\alpha$ whose expectation value should be proportional to the gravitational shear $\gamma_\alpha$. We then show that an anisotropic instrumental psf perturbs the polarisation by $\delta e_\alpha = P^s_{\alpha\beta} p_\beta$, where $p_\alpha$ is a measure of the psf anisotropy and $P^s_{\alpha\beta}$ is the `linearised smear polarisability tensor'. By estimating $P^s_{\alpha\beta}$ for each object we can determine $p_\alpha$ from the foreground stars and apply a correction $-P^s_{\alpha\beta}p_\beta$ to the galaxies. We test this procedure using deep high-resolution images from HST which are smeared with an anisotropic psf and then have noise added to simulate ground-based observations. We find that the procedure works very well. A similar analysis yields a linear shear polarisability tensor $P^\gamma_{\alpha\beta}$ which describes the response to a gravitational shear. This calibrates the polarisation-shear relation, but only for galaxies which are well resolved. To empirically calibrate the effect of seeing on the smaller galaxies we artificially stretch HST images to simulate lensing and then degrade them as before. These experiments provide a rigorous and exacting test of the method under realistic conditions. They show that it is possible to remove the effect of instrumental psf anisotropy, and that the method provides an efficient and quantitative measurement of the gravitational shear.

We develop a method for general non-linear cluster lens reconstruction using the observable distortion of background galaxies. The distortion measures the combination $\gamma/(1-\kappa)$ of shear $\gamma$ and surface density $\kappa$. From this we obtain an expression for the gradient of $\log (1 - \kappa)$ in terms of directly measurable quantities. This allows one to reconstruct $1 - \kappa$ up to an arbitrary constant multiplier. Recent work has emphasised an ambiguity in the relation between the distortion and $\gamma/(1-\kappa)$. Here we show that the functional relation depends only on the parity of the images, so if one has data extending to large radii, and if the critical lines can be visually identified (as lines along which the distortion diverges), this ambiguity is resolved. Moreover, we show that for a generic 2-dimensional lens it is possible to locally determine the parity from the distortion. The arbitrary multiplier, which may in fact take a different value in each region bounded by the contour $\kappa = 1$, can be determined by requiring that the mean surface excess vanish at large radii and that the gradient of $\kappa$ should be continuous across $\kappa = 1$. We show how these ideas might be implemented to reconstruct the surface density, if necessary without use of the data in regions where determination of the parity is insecure, and we show how one can measure the mass contained within an aperture, again, if necessary, without using data within the aperture.

Massive clusters of galaxies gravitationally shear the images of faint background galaxies. At large impact parameters the shear is weak, but can still be measured, to a reasonable degree of significance, as a statistical anisotropy of the faint galaxy images. We describe techniques for measuring the shear and discuss the interpretation of the shear field. We have applied this analysis to ms1224+007. We find a clear detection of the shear, but, puzzlingly, we find a mass about three times that obtained by application of the virial theorem, and obtain a very large mass-to-light ratio. Similar results have been obtained by other groups, and we discuss their implications.

We explore the dark matter distribution in ms1224.7+2007 using the gravitational distortion of the images of faint background galaxies. Projected mass image reconstruction reveals a highly significant concentration coincident with the X-ray and optical location. The concentration is seen repeatably in reconstructions from independent subsamples and the azimuthally averaged tangential shear pattern is also clearly seen in the data. The projected mass within a $2.76'$ radius aperture is $\simeq 3.5 \times 10^{14} h^{-1} M_{\sun}$. This is $\simeq 3$ times larger than that predicted if mass traces light with $M/L = 275 h$ as derived from virial analysis. It is very hard to attribute the discrepancy to a statistical fluctuation, and a further indication of a significant difference between the mass and the light comes from a second mass concentration which is again seen in independent subsamples but which is not seen at all in the cluster light. We find a mass per galaxy visible to $I = 22$ of $\simeq 8 \times 10^{12}h^{-1}M_{\sun}$ which, if representative of the universe, implies a density parameter $\Omega \sim 2$. We find a null detection of any net shear from large-scale structure with a precision of 0.9\% per component. This is much smaller than the possible detection in a recent comparable study, and the precision here is comparable to to the minimum level of rms shear fluctuations implied by observed large-scale structure.

We investigate the production of dark matter with non-standard momentum distributions arising from the decay of a relativistic massive neutrino into a lighter fermion and boson H$\rightarrow$F$+$B. We develop the Boltzmann equations for this process, paying particular attention to spin polarisation effects. Such decays can, under certain circumstances, lead to a process of runaway stimulated decays or `neutrino-lasing' resulting, for suitable boson mass, in the production of mixed-dark-matter. We discuss a number of variations on the model; we show how the yield of cold particles can be increased by a simple modification to our original proposal. We discuss the singlet-majoron model as a possible framework for these decays.

I review a number of tests of theories for structure formation. Large-scale flows and IRAS galaxies indicate a high density parameter $\Omega \simeq 1$, in accord with inflationary predictions, but it is not clear how this meshes with the uniformly low values obtained from virial analysis on scales $\sim$ 1Mpc. Gravitational distortion of faint galaxies behind clusters allows one to construct maps of the mass surface density, and this should shed some light on the large vs small-scale $\Omega$ discrepancy. Power spectrum analysis reveals too red a spectrum (compared to standard CDM) on scales $\lambda \sim 10-100$ $h^{-1}$Mpc, but the gaussian fluctuation hypothesis appears to be in good shape. These results suggest that the problem for CDM lies not in the very early universe --- the inflationary predictions of $\Omega = 1$ and gaussianity both seem to be OK; furthermore, the COBE result severely restricts modifications such as tilting the primordial spectrum --- but in the assumed matter content. The power spectrum problem can be solved by invoking a cocktail of mixed dark matter. However, if gravitational lensing fails to reveal extended dark mass around clusters then we may be forced to explore more radical possibilities for the dark matter.

We develop a general method for power spectrum analysis of three dimensional redshift surveys. We present rigorous analytical estimates for the statistical uncertainty in the power and we are able to derive a rigorous optimal weighting scheme under the reasonable (and largely empirically verified) assumption that the long wavelength Fourier components are Gaussian distributed. We apply the formalism to the updated 1-in-6 QDOT IRAS redshift survey, and compare our results to data from other probes: APM angular correlations; the CfA and the Berkeley 1.2Jy IRAS redshift surveys. Our results bear out and further quantify the impression from e.g.\ counts-in-cells analysis that there is extra power on large scales as compared to the standard CDM model with $\Omega h\simeq 0.5$. We apply likelihood analysis using the CDM spectrum with $\Omega h$ as a free parameter as a phenomenological family of models; we find the best fitting parameters in redshift space and transform the results to real space. Finally, we calculate the distribution of the estimated long wavelength power. This agrees remarkably well with the exponential distribution expected for Gaussian fluctuations, even out to powers of ten times the mean. Our results thus reveal no trace of periodicity or other non-Gaussian behavior.

Recently, Madsen has argued that relativistic decays of massive neutrinos into lighter fermions and bosons may lead, via thermalization, to the formation of a Bose condensate. If correct, this could generate mixed hot and cold dark matter, with important consequences for structure formation. From a detailed study of such decays, we arrive at substantially different conclusions; for a wide range of masses and decay times, we find that stimulated emission of bosons dominates the decay. This phenomenon can best be described as a neutrino laser, pumped by the QCD phase transition. We discuss the implications for structure formation and the dark-matter problem.