Electromagnetic waves in a dynamical axion background exhibit superluminal group velocities at high frequencies and instabilities at low frequencies, altering how photons propagate through space. Local disturbances propagate causally, but unlike in ordinary Maxwell theory, propagation occurs inside as well as on the lightcone. For the unstable modes, the energy density in the electromagnetic field grows exponentially along timelike displacements. In this paper we derive retarded Green functions in axion electrodynamics in various limits and study the time-domain properties of propagating signals.

Despite considerable experimental progress large parts of the axion-like particle (ALP) parameter space remain difficult to probe in terrestrial experiments. In some cases, however, small-scale structure of the ALP dark matter (DM) distribution is strongly enhanced, offering opportunities for astrophysical tests. Such an enhancement can be produced by a period of pre-nucleosynthesis early matter domination (EMD). This cosmology arises in many ultraviolet completions and generates the correct relic abundance for weak coupling $f_a\sim 10^{16}$ GeV, ALP masses in the range $10^{-13}$ eV $<m_a < 1$ eV, and without fine-tuning of the initial misalignment angle. This range includes the QCD axion around $10^{-9}-10^{-8}$ eV. EMD enhances the growth of ALP small-scale structure, leading to the formation of dense ALP miniclusters. We study the interplay between the initial ALP oscillation, reheating temperature, and effective pressure to provide analytic estimates of the minicluster abundance and properties. ALP miniclusters in the EMD cosmology are denser and more abundant than in $\Lambda$CDM. While enhanced substructure generically reduces the prospects of direct detection experiments, we show that pulsar timing and lensing observations can discover these minihalos over a large range of ALP masses and reheating temperatures.

Cohen, Kaplan, and Nelson (CKN) conjectured that the UV and IR cutoffs of effective quantum field theories coupled to gravity are not independent, but are connected by the physics of black holes. We interpret the CKN bound as a scale-dependent depletion of the QFT density of states and discuss various aspects of the bound on small and large scales. For laboratory experiments, we argue that the bound provides small corrections to ordinary quantum field theory, which we estimate to be of order $m_e/M_p$ for $g-2$ of the electron. On large scales, we suggest a modification of the CKN bound due to the presence of cosmological horizons and discuss the connection with entropy bounds.

Many existing and proposed experiments targeting QCD axion dark matter (DM) can also search for a broad class of axion-like particles (ALPs). We analyze the experimental sensitivities to electromagnetically-coupled ALP DM in different cosmological scenarios with the relic abundance set by the misalignment mechanism. We obtain benchmark DM targets for the standard thermal cosmology, a pre-nucleosynthesis period of early matter domination, and a period of kination. These targets are theoretically simple and assume $\mathcal{O}(1)$ misalignment angles, avoiding fine-tuning of the initial conditions. We find that some experiments will have sensitivity to these ALP DM targets before they are sensitive to the QCD axion, and others can potentially reach interesting targets below the QCD band. The ALP DM abundance also depends on the origin of the ALP mass. Temperature-dependent masses that are generated by strong dynamics (as for the QCD axion) correspond to DM candidates with smaller decay constants, resulting in even better detection prospects.

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.

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.

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.

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

We present the details of the photometric and astrometric calibration of the Pan-STARRS1 $3\pi$ Survey. The photometric goals were to reduce the systematic effects introduced by the camera and detectors, and to place all of the observations onto a photometric system with consistent zero points over the entire area surveyed, the ~30,000 square degrees north of $\delta$ = -30 degrees. The astrometric calibration compensates for similar systematic effects so that positions, proper motions, and parallaxes are reliable as well. The Pan-STARRS Data Release 2 (DR2) astrometry is tied to the Gaia DR1 release.

Cosmological simulations require the use of a multiple time-stepping scheme. Without such a scheme, cosmological simulations would be impossible due to their high level of dynamic range; over eleven orders of magnitude in density. Such a large dynamic range leads to a range of over four orders of magnitude in time-step, which presents a significant load-balancing challenge. In this work, the extreme adaptivity that cosmological simulations present is tackled in three main ways through the use of the code SWIFT. First, an adaptive mesh is used to ensure that only the relevant particles are interacted in a given time-step. Second, task-based parallelism is used to ensure efficient load-balancing within a single node, using pthreads and SIMD vectorisation. Finally, a domain decomposition strategy is presented, using the graph domain decomposition library METIS, that bisects the work that must be performed by the simulation between nodes using MPI. These three strategies are shown to give SWIFT near-perfect weak-scaling characteristics, only losing 25% performance when scaling from 1 to 4096 cores on a representative problem, whilst being more than 30x faster than the de-facto standard Gadget-2 code.

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

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

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

New QCD-like `hypercolor' sectors can generate a broad class of new signatures at hadron colliders and furnish a variety of dark matter candidates. Paired diboson resonances are a particularly important collider signature, arising both from $CP$-conserving vector hypermeson decays of the form $\widetilde\rho \rightarrow\widetilde\pi\widetilde\pi\rightarrow 4V$ and from $CP$-violating pseudoscalar hypermeson decays of the form $\widetilde\eta\rightarrow \widetilde\pi\widetilde\pi\rightarrow 4V$. The latter are sensitive to the vacuum angle $\widetilde \theta$ in the hypercolor sector. We study single- and paired-diboson resonance signatures in final states involving gluons and photons at the LHC and a future 100 TeV $pp$ collider, illustrating the discovery potential at both colliders in simple benchmark models. We also describe some of the theoretical and cosmological consequences of $\widetilde \theta$. If $CP$-violating hypermeson decays are observable at hadron colliders, ordinary QCD must have an axion. Such scenarios also provide a natural setting for a dark pion component of dark matter, with its relic abundance set by $CP$-violating annihilations. If the new vacuum angle is relaxed to zero by a dark axion, the relic density can instead be a mixture of axions and dark axions. Overproduction of dark axions is most easily avoided if the universe underwent a period of early matter domination.

If the QCD axion is a significant component of dark matter, and if the universe was once hotter than a few hundred MeV, the axion relic abundance depends on the function $\chi(T)$, the temperature-dependent topological susceptibility. Uncertainties in this quantity induce uncertainties in the axion mass as a function of the relic density, or vice versa. At high temperatures, theoretical uncertainties enter through the dilute instanton gas computation, while in the intermediate and strong coupling regime, only lattice QCD can determine $\chi(T)$ precisely. We reassess the uncertainty on the instanton contribution, arguing that it amounts to less than $20\%$ in the effective action, or a factor of 20 in $\chi$ at $T=1.5$ GeV. We then combine the instanton uncertainty with a range of models for $\chi(T)$ at intermediate temperatures and determine the impact on the axion relic density. We find that for a given relic density and initial misalignment angle, the combined uncertainty amounts to a factor of 2-3 in the zero-temperature axion mass.

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.

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

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

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.

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

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

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

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

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

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

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

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

We present a new open-source cosmological code, called SWIFT, designed to solve the equations of hydrodynamics using a particle-based approach (Smooth Particle Hydrodynamics) on hybrid shared/distributed-memory architectures. SWIFT was designed from the bottom up to provide excellent strong scaling on both commodity clusters (Tier-2 systems) and Top100-supercomputers (Tier-0 systems), without relying on architecture-specific features or specialized accelerator hardware. This performance is due to three main computational approaches: (1) Task-based parallelism for shared-memory parallelism, which provides fine-grained load balancing and thus strong scaling on large numbers of cores. (2) Graph-based domain decomposition, which uses the task graph to decompose the simulation domain such that the work, as opposed to just the data, as is the case with most partitioning schemes, is equally distributed across all nodes. (3) Fully dynamic and asynchronous communication, in which communication is modelled as just another task in the task-based scheme, sending data whenever it is ready and deferring on tasks that rely on data from other nodes until it arrives. In order to use these approaches, the code had to be re-written from scratch, and the algorithms therein adapted to the task-based paradigm. As a result, we can show upwards of 60% parallel efficiency for moderate-sized problems when increasing the number of cores 512-fold, on both x86-based and Power8-based architectures.

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.

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.

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

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.

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

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

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.

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.

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

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

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.

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.

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.

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 extensible N-Dimensional Data Format (NDF) was designed and developed in the late 1980s to provide a data model suitable for use in a variety of astronomy data processing applications supported by the UK Starlink Project. Starlink applications were used extensively, primarily in the UK astronomical community, and form the basis of a number of advanced data reduction pipelines today. This paper provides an overview of the historical drivers for the development of NDF and the lessons learned from using a defined hierarchical data model for many years in data reduction software, data pipelines and in data acquisition systems.

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

SPLAT-VO is a powerful graphical tool for displaying, comparing, modifying and analyzing astronomical spectra, as well as searching and retrieving spectra from services around the world using Virtual Observatory (VO) protocols and services. The development of SPLAT-VO started in 1999, as part of the Starlink StarJava initiative, sometime before that of the VO, so initial support for the VO was necessarily added once VO standards and services became available. Further developments were supported by the Joint Astronomy Centre, Hawaii until 2009. Since end of 2011 development of SPLAT-VO has been continued by the German Astrophysical Virtual Observatory, and the Astronomical Institute of the Academy of Sciences of the Czech Republic. From this time several new features have been added, including support for the latest VO protocols, along with new visualization and spectra storing capabilities. This paper presents the history of SPLAT-VO, it's capabilities, recent additions and future plans, as well as a discussion on the motivations and lessons learned up to now.

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

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

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

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

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

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

The discovery of a large tensor-to-scalar ratio by the BICEP2 experiment points to large field excursions during inflation. One framework that predicts large $r$ is monodromy inflation. While discussed mainly in the context of string theory, the phenomenon can be illustrated and studied in the well-understood framework of SUSY QCD with a large number of colors. We discuss the requirements for viable inflation as well as various difficulties for model building, including tunneling, tuning, and the species problem.

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.

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.

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.

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

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

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

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

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

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

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

We determine the prospects for direct and indirect detection of thermal relic neutralinos in supersymmetric theories with multi-TeV squarks and sleptons. We consider the concrete example of the focus point region of minimal supergravity, but our results are generically valid for all models with decoupled scalars and mixed Bino-Higgsino or Higgsino-like dark matter. We determine the parameter space consistent with a 125 GeV Higgs boson including 3-loop corrections in the calculation of the Higgs mass. These corrections increase m_h by 1-3 GeV, lowering the preferred scalar mass scale and decreasing the fine-tuning measure in these scenarios. We then systematically examine prospects for dark matter direct and indirect detection. Direct detection constraints do not exclude these models, especially for \mu < 0. At the same time, the scenario generically predicts spin-independent signals just beyond current bounds. We also consider indirect detection with neutrinos, gamma rays, anti-protons, and anti-deuterons. Current IceCube neutrino constraints are competitive with direct detection, implying bright prospects for complementary searches with both direct and indirect detection.

We study a limit of the nearly-Peccei-Quinn-symmetric Next-to-Minimal Supersymmetric Standard Model possessing novel Higgs and dark matter (DM) properties. In this scenario, there naturally co-exist three light singlet-like particles: a scalar, a pseudoscalar, and a singlino-like DM candidate, all with masses of order 0.1-10 GeV. The decay of a Standard Model-like Higgs boson to pairs of the light scalars or pseudoscalars is generically suppressed, avoiding constraints from collider searches for these channels. For a certain parameter window annihilation into the light pseudoscalar and exchange of the light scalar with nucleons allow the singlino to achieve the correct relic density and a large direct detection cross section consistent with the CoGeNT and DAMA/LIBRA preferred region simultaneously. This parameter space is consistent with experimental constraints from LEP, the Tevatron, and Upsilon- and flavor physics.

The Planck Collaboration measured the Sunyaev-Zel'dovich (SZ) decrement of optically selected clusters from the Sloan Digital Sky Survey, finding that it falls significantly below expectations based on existing mass calibration of the maxBCG galaxy clusters. Resolving this tension requires either the data to go up, or the theoretical expectations to come down. Here, we use data from the Wilkinson Microwave Anisotropy Probe (WMAP) to perform an independent estimate of the SZ decrement of maxBCG clusters. The recovered signal is consistent with that obtained using Planck, though with larger error bars due to WMAP's larger beam size and smaller frequency range. Nevertheless, this detection serves as an independent confirmation of the magnitude of the effect, and demonstrates that the observed discrepancy must be theoretical in origin.