This paper describes the organization of the database and the catalog data products from the Pan-STARRS1 $3\pi$ Steradian Survey. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI. The database is described in detail, including the construction of the database, the provenance of the data, the schema, and how the database tables are related. Examples of queries for a range of science goals are included. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI.

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

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.

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

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

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

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

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

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

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

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.

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

We present the early UV and optical light curve of Type IIP supernova (SN) 2010aq at z=0.0862, and compare it to analytical models for thermal emission following SN shock breakout in a red supergiant star. SN 2010aq was discovered in joint monitoring between the Galaxy Evolution Explorer (GALEX) Time Domain Survey (TDS) in the NUV and the Pan-STARRS1 Medium Deep Survey (PS1 MDS) in the g, r, i, and z bands. The GALEX and Pan-STARRS1 observations detect the SN less than 1 day after shock breakout, measure a diluted blackbody temperature of 31,000 +/- 6,000 K 1 day later, and follow the rise in the UV/optical light curve over the next 2 days caused by the expansion and cooling of the SN ejecta. The high signal-to-noise ratio of the simultaneous UV and optical photometry allows us to fit for a progenitor star radius of 700 +/- 200 R_sun, the size of a red supergiant star. An excess in UV emission two weeks after shock breakout compared to SNe well fitted by model atmosphere-code synthetic spectra with solar metallicity, is best explained by suppressed line blanketing due to a lower metallicity progenitor star in SN 2010aq. Continued monitoring of PS1 MDS fields by the GALEX TDS will increase the sample of early UV detections of Type II SNe by an order of magnitude, and probe the diversity of SN progenitor star properties.

We present photometric and spectroscopic observations of a luminous type IIP Supernova 2009kf discovered by the Pan-STARRS 1 (PS1) survey and detected also by GALEX. The SN shows a plateau in its optical and bolometric light curves, lasting approximately 70 days in the rest frame, with absolute magnitude of M_V = -18.4 mag. The P-Cygni profiles of hydrogen indicate expansion velocities of 9000km/s at 61 days after discovery which is extremely high for a type IIP SN. SN 2009kf is also remarkably bright in the near-ultraviolet (NUV) and shows a slow evolution 10-20 days after optical discovery. The NUV and optical luminosity at these epochs can be modelled with a black-body with a hot effective temperature (T ~16,000 K) and a large radius (R ~1x10^{15} cm). The bright bolometric and NUV luminosity, the lightcurve peak and plateau duration, the high velocities and temperatures suggest that 2009kf is a type IIP SN powered by a larger than normal explosion energy. Recently discovered high-z SNe (0.7 < z < 2.3) have been assumed to be IIn SNe, with the bright UV luminosities due to the interaction of SN ejecta with a dense circumstellar medium (CSM). UV bright SNe similar to SN 2009kf could also account for these high-z events, and its absolute magnitude M_NUV = -21.5 +/- 0.5 mag suggests such SNe could be discovered out to z ~2.5 in the PS1 survey.

We report the discovery of a massive (Mpsini = 13.02 +/- 0.64 Mjup; total mass 13.25 +/- 0.64 Mjup), large (1.95 +/- 0.16 Rjup) planet in a transiting, eccentric orbit (e = 0.260 +/- 0.017) around a 10th magnitude F5V star in the constellation Camelopardalis. We designate the planet XO-3b, and the star XO-3, also known as GSC 03727-01064. The orbital period of XO-3b is 3.1915426 +/- 0.00014 days. XO-3 lacks a trigonometric distance; we estimate its distance to be 260 +/- 23 pc. The radius of XO-3 is 2.13 +/- 0.21 Rsun, its mass is 1.41 +/- 0.08 Msun, its vsini = 18.54 +/- 0.17 km/s, and its metallicity is [Fe/H] = -0.177 +/- 0.027. This system is unusual for a number of reasons. XO-3b is one of the most massive planets discovered around any star for which the orbital period is less than 10 days. The mass is near the deuterium burning limit of 13 Mjup, which is a proposed boundary between planets and brown dwarfs. Although Burrows et al. (2001) propose that formation in a disk or formation in the interstellar medium in a manner similar to stars is a more logical way to differentiate planets and brown dwarfs, our current observations are not adequate to address this distinction. XO-3b is also unusual in that its eccentricity is large given its relatively short orbital period. Both the planetary radius and the inclination are functions of the spectroscopically determined stellar radius. Analysis of the transit light curve of XO-3b suggests that the spectroscopically derived parameters may be over estimated. Though relatively noisy, the light curves favor a smaller radius in order to better match the steepness of the ingress and egress. The light curve fits imply a planetary radius of 1.25 +/- 0.15 Rjup, which would correspond to a mass of 12.03 +/- 0.46 Mjup.

We report on a V=11.2 early K dwarf, XO-2 (GSC 03413-00005), that hosts a Rp=0.98+0.03/-0.01 Rjup, Mp=0.57+/-0.06 Mjup transiting extrasolar planet, XO-2b, with an orbital period of 2.615857+/-0.000005 days. XO-2 has high metallicity, [Fe/H]=0.45+/-0.02, high proper motion, mu_tot=157 mas/yr, and has a common proper motion stellar companion with 31" separation. The two stars are nearly identical twins, with very similar spectra and apparent magnitudes. Due to the high metallicity, these early K dwarf stars have a mass and radius close to solar, Ms=0.98+/-0.02 Msolar and Rs=0.97+0.02/-0.01 Rsolar. The high proper motion of XO-2 results from an eccentric orbit (Galactic pericenter, Rper<4 kpc) well confined to the Galactic disk (Zmax~100 pc). In addition, the phase space position of XO-2 is near the Hercules dynamical stream, which points to an origin of XO-2 in the metal-rich, inner Thin Disk and subsequent dynamical scattering into the solar neighborhood. We describe an efficient Markov Chain Monte Carlo algorithm for calculating the Bayesian posterior probability of the system parameters from a transit light curve.

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

A planet transits an 11th magnitude, G1V star in the constellation Corona Borealis. We designate the planet XO-1b, and the star, XO-1, also known as GSC 02041-01657. XO-1 lacks a trigonometric distance; we estimate it to be 200+-20 pc. Of the ten stars currently known to host extrasolar transiting planets, the star XO-1 is the most similar to the Sun in its physical characteristics: its radius is 1.0+-0.08 R_Sun, its mass is 1.0+-0.03 M_Sun, V sini < 3 km/s, and its metallicity [Fe/H] is 0.015+-0.04. The orbital period of the planet XO-1b is 3.941534+-0.000027 days, one of the longer ones known. The planetary mass is 0.90+-0.07 M_Jupiter, which is marginally larger than that of other transiting planets with periods between 3 and 4 days. Both the planetary radius and the inclination are functions of the spectroscopically determined stellar radius. If the stellar radius is 1.0+-0.08 R_Sun, then the planetary radius is 1.30+-0.11 R_Jupiter and the inclination of the orbit is 87.7+-1.2 degrees. We have demonstrated a productive international collaboration between professional and amateur astronomers that was important to distinguishing this planet from many other similar candidates.

The XO project's first objective is to find hot Jupiters transiting bright stars, i.e. V < 12, by precision differential photometry. Two XO cameras have been operating since September 2003 on the 10,000-foot Haleakala summit on Maui. Each XO camera consists of a 200-mm f/1.8 lens coupled to a 1024x1024 pixel, thinned CCD operated by drift scanning. In its first year of routine operation, XO has observed 6.6% of the sky, within six 7 deg-wide strips scanned from 0 deg to +63 deg of declination and centered at RA=0, 4, 8, 12, 16, and 20 hours. Autonomously operating, XO records 1 billion pixels per clear night, calibrates them photometrically and astrometrically, performs aperture photometry, archives the pixel data and transmits the photometric data to STScI for further analysis. From the first year of operation, the resulting database consists of photometry of 100,000 stars at more than 1000 epochs per star with differential photometric precision better than 1% per epoch. Analysis of the light curves of those stars produces transiting-planet candidates requiring detailed follow up, described elsewhere, culminating in spectroscopy to measure radial-velocity variation in order to differentiate genuine planets from the more numerous impostors, primarily eclipsing binary and multiple stars.