Pan-STARRS1 Survey Research Articles

Selecting superluminous supernovae in faint galaxies from the first year of the Pan-STARRS1 Medium Deep Survey

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.

The main-belt comets: The Pan-STARRS1 perspective

We analyze a set of 760475 observations of 333026 unique main-belt objects obtained by the Pan-STARRS1 (PS1) survey telescope between 2012 May 20 and 2013 November 9, a period during which PS1 discovered two main-belt comets, P/2012 T1 (PANSTARRS) and P/2013 R3 (Catalina-PANSTARRS). PS1 comet detection procedures currently consist of the comparison of the point spread functions (PSFs) of moving objects to those of reference stars, and the flagging of objects that show anomalously large radial PSF widths for human evaluation and possible observational follow-up. Based on the number of missed discovery opportunities among comets discovered by other observers, we estimate an upper limit comet discovery efficiency rate of ∼70% for PS1. Additional analyses that could improve comet discovery yields in future surveys include linear PSF analysis, modeling of trailed stellar PSFs for comparison to trailed moving object PSFs, searches for azimuthally localized activity, comparison of point-source-optimized photometry to extended-source-optimized photometry, searches for photometric excesses in objects with known absolute magnitudes, and crowd-sourcing. Analysis of the discovery statistics of the PS1 survey indicates an expected fraction of 59 MBCs per 106 outer main-belt asteroids (corresponding to a total expected population of ∼140 MBCs among the outer main-belt asteroid population with absolute magnitudes of 12<HV<19.5), and a 95% confidence upper limit of 96 MBCs per 106 outer main-belt asteroids (corresponding to a total of ∼230 MBCs), assuming a detection efficiency of 50%. We note however that significantly more sensitive future surveys (particularly those utilizing larger aperture telescopes) could detect many more MBCs than estimated here. Examination of the orbital element distribution of all known MBCs reveals an excess of high eccentricities (0.1<e<0.3) relative to the background asteroid population. Theoretical calculations show that, given these eccentricities, the sublimation rate for a typical MBC is orders of magnitude larger at perihelion than at aphelion, providing a plausible physical explanation for the observed behavior of MBCs peaking in observed activity strength near perihelion. These results indicate that the overall rate of mantle growth should be slow, consistent with observational evidence that MBC activity can be sustained over multiple orbit passages.

STACKING THE INVISIBLES: A GUIDED SEARCH FOR LOW-LUMINOSITY MILKY WAY SATELLITES

Almost every known low-luminosity Milky Way dwarf spheroidal (dSph) satellite galaxy contains at least one RR Lyrae star. Assuming that a fraction of distant (60 < d_{helio} < 100 kpc) Galactic halo RR Lyrae stars are members of yet to be discovered low-luminosity dSph galaxies, we perform a {\em guided} search for these low-luminosity dSph galaxies. In order to detect the presence of dSph galaxies, we combine stars selected from more than 123 sightlines centered on RR Lyrae stars identified by the Palomar Transient Factory. We find that this method is sensitive enough to detect the presence of Segue 1-like galaxies (M_V= -1.5^{+0.6}_{-0.8}, r_h=30 pc) even if only ~20 sightlines were occupied by such dSph galaxies. Yet, when our method is applied to the SDSS DR10 imaging catalog, no signal is detected. An application of our method to sightlines occupied by pairs of close (<200 pc) horizontal branch stars, also did not yield a detection. Thus, we place upper limits on the number of low-luminosity dSph galaxies with half-light radii from 30 pc to 120 pc, and in the probed volume of the halo. Stronger constraints on the luminosity function may be obtained by applying our method to sightlines centered on RR Lyrae stars selected from the Pan-STARRS1 survey, and eventually, from LSST. In the Appendix, we present spectroscopic observations of an RRab star in the Bo\"{o}tes 3 dSph and a light curve of an RRab star near the Bo\"{o}tes 2 dSph.

Chasing the phantom: A closer look at type Ia supernovae and the dark energy equation of state

Some recent observations provide $> 2\sigma$ evidence for phantom dark energy -- a value of the dark energy equation of state less than the cosmological-constant value of $-1$. We focus on constraining the equation of state by combining current data from the most mature geometrical probes of dark energy: Type Ia supernovae (SNe Ia) from the Supernova Legacy Survey (SNLS3), the Supernova Cosmology Project (Union2.1), and the Pan-STARRS1 survey (PS1); cosmic microwave background measurements from Planck and WMAP9; and a combination of measurements of baryon acoustic oscillations. The combined data are consistent with $w = -1$ for the Union2.1 sample, though they present moderate ($\sim 1.9\sigma$) evidence for a phantom value when either the SNLS3 or PS1 sample is used instead. We study the dependence of the constraints on the redshift, stretch, color, and host galaxy stellar mass of SNe, but we find no unusual trends. In contrast, the constraints strongly depend on any external $H_0$ prior: a higher adopted value for the direct measurement of the Hubble constant ($H_0 \gtrsim 71~\text{km/s/Mpc}$) leads to $\gtrsim 2\sigma$ evidence for phantom dark energy. Given Planck data, we can therefore make the following statement at $2\sigma$ confidence: either the SNLS3 and PS1 data have systematics that remain unaccounted for or the Hubble constant is below 71 km/s/Mpc; else the dark energy equation of state is indeed phantom.

MEASURING DISTANCES AND REDDENINGS FOR A BILLION STARS: TOWARD A 3D DUST MAP FROM PAN-STARRS 1

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.

DETECTION OF AN OUTBURST ONE YEAR PRIOR TO THE EXPLOSION OF SN 2011ht

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.

Pan-STARRS1: Galaxy clustering in the Small Area Survey 2

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.

OUTGASSING BEHAVIOR OF C/2012 S1 (ISON) FROM 2011 SEPTEMBER TO 2013 JUNE

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

The Pan-STARRS1 Small Area Survey 2

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.

DETERMINATION OF AN UPPER LIMIT FOR THE WATER OUTGASSING RATE OF MAIN-BELT COMET P/2012 T1 (PANSTARRS)

A new Main-Belt Comet (MBC) P/2012 T1 (PANSTARRS) was discovered on 2012 October 6, approximately one month after its perihelion, by the Pan-STARRS1 survey based in Hawaii. It displayed cometary activity upon its discovery with one hypothesis being that the activity was driven by sublimation of ices; as a result, we searched for emission assumed to be driven by the sublimation of subsurface ices. Our search was of the H2O 110-101 ground state rotational line at 557 GHz from P/2012 T1 (PANSTARRS) with the Heterodyne Instrument for the Far Infrared on board the Herschel Space Observatory on 2013 January 16, when the object was at a heliocentric distance of 2.504 AU and a geocentric distance of 2.064 AU. Perihelion was in early 2012 September at a distance of 2.411 AU. While no H2O line emission was detected in our observations, we were able to derive sensitive 3σ upper limits for the water production rate and column density of <7.63 × 1025 molecules s–1 and of <1.61 × 1011 cm–2, respectively. An observation taken on 2013 January 15 using the Very Large Telescope found the MBC to be active during the Herschel observation, suggesting that any ongoing sublimation due to subsurface ice was lower than our upper limit.

LACERTA I AND CASSIOPEIA III. TWO LUMINOUS AND DISTANT ANDROMEDA SATELLITE DWARF GALAXIES FOUND IN THE 3π PAN-STARRS1 SURVEY

We report the discovery of two new dwarf galaxies, Lacerta I/Andromeda XXXI (Lac I/And XXXI) and Cassiopeia III/Andromeda XXXII (Cas III/And XXXII), in stacked Pan-STARRS1 rP1- and iP1-band imaging data. Both are luminous systems (MV ∼ −12) located at projected distances of 203 and 105 from M31. Lac I and Cas III are likely satellites of the Andromeda galaxy with heliocentric distances of and , respectively, and corresponding M31-centric distances of 275 ± 7 kpc and . The brightest of recent Local Group member discoveries, these two new dwarf galaxies owe their late discovery to their large sizes ( arcmin or for Lac I; arcmin or 1456 ± 267 pc for Cas III) and consequently low surface brightness (μ0 ∼ 26.0 mag arcsec−2), as well as to the lack of a systematic survey of regions at large radii from M31, close to the Galactic plane. This latter limitation is now alleviated by the 3π Pan-STARRS1 survey, which could lead to the discovery of other distant Andromeda satellite dwarf galaxies.

Pan-STARRS, ATLAS and optical transient searches

The Pan-STARRS1 survey is collecting multi-epoch, multi-colour observations of the sky north of declination -30°, and has designated 70 deg(2) for nightly observations that are particularly useful for transient detection. A duplicate, Pan-STARRS2, is nearing completion that offers opportunities to improve the quality of transient search and observation, as well as simply increasing the number of detections. A new system, the Asteroid Terrestrial-impact Last Alert System (ATLAS), increases the search area to all-sky in return for diminished sensitivity, and highlights tension among optimization for static sky images, optimization for faint transients and optimization for an unbiased number of transients. ATLAS gives up sub-arcsecond images and full colour information to specialize in the third category, but should detect many more transients than the Pan-STARRS1 Medium Deep fields or the Palomar Transient Factory, with examples of transient classes that are considerably closer and brighter.

PS1-10jh – a tidal disruption event with an extremely low disc temperature

Abstract The cooler than expected optical–UV transient PS1-10jh detected by the Pan-STARRS1 survey is probably related to a tidal disruption event in which a He-rich stellar core remnant is implied. The evolution of bound debris during the disc-phase is studied by solving the hydrodynamic equations. The model provides a good fit either of the raising part of the light curve in the bands gP1, rP1 and iP1 or in the early decay. The parameters characterizing this optimized model are the mass of the central black hole, i.e. 6.3 × 106 M⊙, and the critical Reynolds number ${\cal R} = 10^4$ that fixes the viscosity and the accretion time-scale. Such a high value of ${\cal R}$ explains the low disc temperature and the consequent absence of X-ray emission. The predicted bolometric peak luminosity is about 1045 erg s− 1 and the predicted total radiated energy is about Erad = 2.67 × 1051 erg.

A Pan-STARRS1 VIEW OF THE BIFURCATED SAGITTARIUS STREAM

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

Improved Asteroid Astrometry and Photometry with Trail Fitting

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

PHOTOMETRIC CALIBRATION OF THE FIRST 1.5 YEARS OF THE PAN-STARRS1 SURVEY

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

P/2006 VW139: a main-belt comet born in an asteroid collision?

In this paper, we apply different methods to examine the possibility that a small group of 24 asteroids dynamically linked to a main-belt comet P/2006 VW139, recently discovered by the Pan-STARRS1 survey telescope, shares a common physical origin. By applying the hierarchical clustering and backward integration methods, we find strong evidence that 11 of these asteroids form a sub-group which likely originated in a recent collision event, and that this group includes P/2006 VW139. The objects not found to be part of the 11-member sub-group, which we designate as the P/2006 VW139 family, were either found to be dynamically unstable or are likely interlopers which should be expected due to the close proximity of the Themis family. As we demonstrated, statistical significance of the P/2006 VW139 family is >99 per cent. We determine the age of the family to be 7.5 ± 0.3 Myr, and estimate the diameter of the parent body to be ∼11 km. Results show that the family is produced by an impact which can be best characterized as a transition from the catastrophic to the cratering regime. The dynamical environment of this family is studied as well, including the identification of the most influential mean motion and secular resonances in the region. Our findings now make P/2006 VW139 the second main-belt comet to be dynamically associated with a young asteroid family, a fact with important implications for the origin and activation mechanism of such objects.

THE Pan-STARRS1 PHOTOMETRIC SYSTEM

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

DISCOVERY OF MAIN-BELT COMET P/2006 VW139BY Pan-STARRS1

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

THE PHOTOMETRIC CLASSIFICATION SERVER FOR Pan-STARRS1

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