Source author record

K. C. Chambers

K. C. Chambers appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

97works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

97 published item(s)

preprint2026arXiv

Can tidal disruption event models reliably measure black hole masses?

Tidal disruption event (TDE) light curves are increasingly used to infer the masses of quiescent supermassive black holes ($M_{\rm{BH}}$), offering a powerful probe of low-mass black hole demographics independent of host-galaxy scaling relations. However, the reliability of most semi-analytic TDE models assume full stellar disruption, despite theoretical expectations that partial disruptions dominate the TDE population. In this work we test the robustness of current TDE models using three repeating partial TDEs (rpTDEs), in which the multiple flares produced by the same surviving stellar core must yield consistent black hole masses. We present spectroscopic observations establishing AT 2023adr as a rpTDE, making it the third such spectroscopically confirmed event. We independently model the flares of the three rpTDEs; 2020vdq, 2022dbl, and 2023adr, applying fallback-accretion fits, stream-stream collision scaling relations, luminosity-based empirical relations, and cooling-envelope fits. After accounting for statistical and model-specific systematics, we find that all TDE models generally return self-consistent $M_{\rm{BH}}$ values between flares, and are broadly consistent with host-galaxy $M_{\rm{BH}}$ proxies, recovering $M_{\rm{BH}}$ to within 0.3-0.5 dex. However, the convergence of fallback models towards unphysical stellar masses and impact parameters reveals limitations in the existing fallback model grids. We also show that light curve coverage, particularly in the near-UV, is critical for constraining model parameters. This has direct implications for interpreting the thousands of TDE light curves expected from upcoming surveys such as the Rubin Observatory's Legacy Survey of Space and Time, where from simulations, we find that $M_{\rm{BH}}$ may be underestimated on average by 0.5 dex without additional follow-up.

preprint2026arXiv

Euclid preparation. Calibrated intrinsic galaxy alignments in the Euclid Flagship simulation

Intrinsic alignments of galaxies are potentially a major contaminant of cosmological analyses of weak gravitational lensing. We construct a semi-analytic model of galaxy ellipticities and alignments in the \Euclid Flagship simulation to predict this contamination in Euclid's weak lensing observations. Galaxy shapes and orientations are determined by the corresponding properties of the host haloes in the underlying $N$-body simulation, as well as the relative positions of galaxies within their halo. Alignment strengths are moderated via stochastic misalignments, separately for central and satellite galaxies and conditional on the galaxy's redshift, luminosity, and rest-frame colour. The resulting model is calibrated against galaxy ellipticity statistics from the COSMOS Survey, selected alignment measurements based on Sloan Digital Sky Survey samples, and galaxy orientations extracted from the Horizon-AGN hydrodynamic simulation at redshift $z=1$. The best-fit model has a total of 12 alignment parameters and generally reproduces the calibration data sets well within the $1σ$ statistical uncertainties of the observations and the \flagship simulation, with notable exceptions for the most luminous sub-samples on small physical scales. The statistical power of the calibration data and the volume of the single \flagship realisation are still too small to provide informative prior ranges for intrinsic alignment amplitudes in relevant galaxy samples. As a first application, we predict that \Euclid end-of-mission tomographic weak gravitational lensing two-point statistics are modified by up to order $10\,\%$ due to intrinsic alignments.

preprint2026arXiv

Euclid preparation. Galaxy 2-point correlation function modelling in redshift space

The Euclid satellite will measure spectroscopic redshifts for tens of millions of emission-line galaxies. In the context of Stage-IV surveys, the 3-dimensional clustering of galaxies plays a key role in providing cosmological constraints. In this paper, we conduct a model comparison for the multipole moments of the galaxy 2-point correlation function (2PCF) in redshift space. We test state-of-the-art models, in particular the effective field theory of large-scale structure (EFT), one based on the velocity difference generating function (VDG$_{\infty}$), and different variants of Lagrangian perturbation theory (LPT) models, such as convolutional LPT (CLPT) and its effective-field-theory extension (CLEFT). We analyse the first three even multipoles of the 2PCF in the Flagship 1 simulation, which consists of four snapshots at $z\in\{0.9,1.2,1.5,1.8\}$. We study both template-fitting and full-shape approaches and find that with the template-fitting approach, only the VDG$_{\infty}$ model is able to reach a minimum fitting scale of $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ at $z=0.9$ without biasing the recovered parameters. Indeed, the EFT model becomes inaccurate already at $s_{\rm min}=30\,h^{-1}\,{\rm Mpc}$. Conversely, in the full-shape analysis, the CLEFT and VDG$_{\infty}$ models perform similarly well, but only the CLEFT model can reach $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ while the VDG$_{\infty}$ model is unbiased down to $s_{\rm min}=25\,h^{-1}\,{\rm Mpc}$ at the lowest redshift. Overall, in order to achieve the accuracy required by Euclid, non-perturbative modelling such as in the VDG$_{\infty}$ or CLEFT models should be considered. At $z=1.8$, the CLPT model is sufficient to describe the data with high figure of merit. This comparison selects baseline models that perform best in ideal conditions and sets the stage for an optimal analysis of Euclid data in configuration space.

preprint2026arXiv

Euclid preparation. Testing analytic models of galaxy intrinsic alignments in the Euclid Flagship simulation

We model intrinsic alignments (IA) in Euclid's Flagship simulation to investigate its impact on Euclid's weak lensing signal. Our IA implementation in the Flagship simulation takes into account photometric properties of galaxies as well as their dark matter host halos. We compare simulations against theory predictions, determining the parameters of two of the most widely used IA models: the Non Linear Alignment (NLA) and the Tidal Alignment and Tidal Torquing (TATT) models. We measure the amplitude of the simulated IA signal as a function of galaxy magnitude and colour in the redshift range $0.1<z<2.1$. We find that both NLA and TATT can accurately describe the IA signal in the simulation down to scales of $6$-$7 \,h^{-1}\,$Mpc. We measure alignment amplitudes for red galaxies comparable to those of the observations, with samples not used in the calibration procedure. For blue galaxies, our constraints are consistent with zero alignments in our first redshift bin $0.1 < z < 0.3$, but we detect a non-negligible signal at higher redshift, which is, however, consistent with the upper limits set by observational constraints. Additionally, several hydrodynamical simulations predict alignment for spiral galaxies, in agreement with our findings. Finally, the evolution of alignment with redshift is realistic and comparable to that determined in the observations. However, we find that the commonly adopted redshift power-law for IA fails to reproduce the simulation alignments above $z=1.1$. A significantly better agreement is obtained when a luminosity dependence is included, capturing the intrinsic luminosity evolution with redshift in magnitude-limited surveys. We conclude that the Flagship IA simulation is a useful tool for translating current IA constraints into predictions for IA contamination of Euclid-like samples.

preprint2026arXiv

Euclid: Galaxy SED reconstruction in the PHZ processing function: impact on the PSF and the role of medium-band filters

Weak lensing surveys require accurate correction for the point spread function (PSF) when measuring galaxy shapes. For a diffraction-limited PSF, as arises in space-based missions, this correction depends on each galaxy SED. In the Euclid mission, galaxy SED reconstruction, a tasks of the photometric-redshift processing function (PHZ PF), relies on broad- and medium-band ancillary photometry. The limited wavelength sampling of the Euclid VIS passband and signal-to-noise ratio may affect the reconstruction accuracy and translate into biases in the weak lensing measurements. In this study, we present the methodology, which is employed in the Euclid PHZ PF, for reconstructing galaxy SEDs at 55 wavelengths, sampling the VIS passband every 10 nm, and we assess whether it fulfils the accuracy requirements imposed on the Euclid PSF model. We employ both physics- and data-driven methods, focusing on a new approach of template-based flux correction and Gaussian processes, and we introduce an SED metric whose bias propagates into PSF quadrupole moment errors. Our findings demonstrate that Gaussian processes and template fitting meet the requirements only in specific, but complementary, redshift intervals. We therefore propose a hybrid approach, which leverages both methods. This solution proves to be effective in meeting the Euclid accuracy requirements for most of the redshift range of the survey. Finally, we investigate the impact on the SED reconstruction of a new set of 16 evenly-spaced medium-band filters for the Subaru telescope, providing quasi-spectroscopic coverage of the VIS passband. This study shows promising results, ensuring accurate SED reconstruction and meeting the mission PSF requirements. This work thus provides not only the methodological foundation of galaxy SED reconstruction in the Euclid PHZ PF, but also a roadmap for future improvements using a new medium-band survey.

preprint2022arXiv

A long life of excess: The interacting transient SN 2017hcc

In this study we present the results of a five-year follow-up campaign of the long-lived type IIn supernova SN 2017hcc, found in a spiral dwarf host of near-solar metallicity. The long rise time (57 $\pm$ 2 days, ATLAS $o$ band) and high luminosity (peaking at $-$20.78 $\pm$ 0.01 mag in the ATLAS $o$ band) point towards an interaction of massive ejecta with massive and dense circumstellar material (CSM). The evolution of SN 2017hcc is slow, both spectroscopically and photometrically, reminiscent of the long-lived type IIn, SN 2010jl. An infrared (IR) excess was apparent soon after the peak, and blueshifts were noticeable in the Balmer lines starting from a few hundred days, but appeared to be fading by around +1200 days. We posit that an IR light echo from pre-existing dust dominates at early times, with some possible condensation of new dust grains occurring at epochs >$\sim$+800 days.

preprint2022arXiv

Characterizing crosstalk within the Pan-STARRS GPC1 camera

Using data from a year-long dedicated campaign to observe bright stars, we study the crosstalk channels present in the GPC1 camera. By analyzing these data, we construct a dataset that checks source stars on almost every CCD of every chip within the camera against all possible crosstalk destinations. We use a clustering algorithm to find potential crosstalk occurrences, and then also check all possible combinations (driven by the hardware layout) by eye. This results in a total of 640 rules, with a flux attenuation factor ranging from 2.5x10$^{2}$ for the bright end to 2.5$\times$10$^{4}$ at the faint end. The average value of m$_{cross}$-m$_{src}\approx$-10.25 corresponds to an attenuating factor of 1.25x10$^{4}$, which produces crosstalk ghosts with an average signal-to-noise ratio of 0.64$\pm$0.1 on the bright images. We find no evidence of crosstalk signals between CCDs not connected in the hardware setup. The distribution of attenuation factors is also found to be dependent on crosstalk movement. A clear dependence on cell column offsets is found, consistent with the idea that the source star charge is progressively attenuated during the traversal of cell readout lines. While we can see the trends, the uncertainties on aperture magnitude measurements are large at this stage.

preprint2022arXiv

SN 2020kyg and the rates of faint Iax Supernovae from ATLAS

We present multi-wavelength follow-up observations of the ATLAS discovered faint Iax supernova SN 2020kyg that peaked at an absolute magnitude of $M_g \approx -14.9 \pm 0.2$, making it another member of the faint Iax supernova population. The bolometric light curve requires only $\approx 7 \times 10^{-3}$ M$_{\odot}$ of radioactive $^{56}$Ni, with an ejected mass of $M_{\rm ej} \sim 0.4$ M$_{\odot}$ and a low kinetic energy of $E \approx 0.05 \pm 0.02 \times 10^{51}$ erg. We construct a homogeneous volume-limited sample of 902 transients observed by ATLAS within 100 Mpc during a 3.5 year span. Using this sample, we constrain the rates of faint Iax ($M_r \gtrsim -16$) events within 60 Mpc at $12^{+14}_{-8}\%$ of the SN Ia rate. The overall Iax rate, at $15^{+17}_{-9}\%$ of the Ia rate, is dominated by the low-luminosity events, with luminous SNe Iax ($M_r \lesssim -17.5$) like 2002cx and 2005hk accounting for only $0.9^{+1.1}_{-0.5}\%$ of the Ia rate (a 2$σ$ upper limit of approximately 3\%). We favour the hybrid CONe WD + He star progenitor channel involving a failed deflagration of a near Chandrasekhar mass white dwarf, expected to leave a bound remnant and a surviving secondary companion, as a candidate explanation for faint Iax explosions. This scenario requires short delay times, consistent with the observed environments of SNe Iax. Furthermore, binary population synthesis calculations have suggested rates of $1-18\%$ of the SN Ia rate for this channel, consistent with our rate estimates.

preprint2021arXiv

The Foundation Supernova Survey: Photospheric Velocity Correlations in Type Ia Supernovae

The ejecta velocities of type-Ia supernovae (SNe Ia), as measured by the Si II $λ6355$ line, have been shown to correlate with other supernova properties, including color and standardized luminosity. We investigate these results using the Foundation Supernova Survey, with a spectroscopic data release presented here, and photometry analyzed with the SALT2 light-curve fitter. We find that the Foundation data do not show significant evidence for an offset in color between SNe Ia with high and normal photospheric velocities, with $Δc = 0.005 \pm 0.014$. Our SALT2 analysis does show evidence for redder high-velocity SN Ia in other samples, including objects from the Carnegie Supernova Project, with a combined sample yielding $Δc = 0.017 \pm 0.007$. When split on velocity, the Foundation SN Ia also do not show a significant difference in Hubble diagram residual, $ΔHR = 0.015 \pm 0.049$ mag. Intriguingly, we find that SN Ia ejecta velocity information may be gleaned from photometry, particularly in redder optical bands. For high-redshift SN Ia, these rest-frame red wavelengths will be observed by the Nancy Grace Roman Space Telescope. Our results also confirm previous work that SN Ia host-galaxy stellar mass is strongly correlated with ejecta velocity: high-velocity SN Ia are found nearly exclusively in high-stellar-mass hosts. However, host-galaxy properties alone do not explain velocity-dependent differences in supernova colors and luminosities across samples. Measuring and understanding the connection between intrinsic explosion properties and supernova environments, across cosmic time, will be important for precision cosmology with SNe Ia.

preprint2021arXiv

The Young Supernova Experiment: Survey Goals, Overview, and Operations

Time domain science has undergone a revolution over the past decade, with tens of thousands of new supernovae (SNe) discovered each year. However, several observational domains, including SNe within days or hours of explosion and faint, red transients, are just beginning to be explored. Here, we present the Young Supernova Experiment (YSE), a novel optical time-domain survey on the Pan-STARRS telescopes. Our survey is designed to obtain well-sampled $griz$ light curves for thousands of transient events up to $z \approx 0.2$. This large sample of transients with 4-band light curves will lay the foundation for the Vera C. Rubin Observatory and the Nancy Grace Roman Space Telescope, providing a critical training set in similar filters and a well-calibrated low-redshift anchor of cosmologically useful SNe Ia to benefit dark energy science. As the name suggests, YSE complements and extends other ongoing time-domain surveys by discovering fast-rising SNe within a few hours to days of explosion. YSE is the only current four-band time-domain survey and is able to discover transients as faint $\sim$21.5 mag in $gri$ and $\sim$20.5 mag in $z$, depths that allow us to probe the earliest epochs of stellar explosions. YSE is currently observing approximately 750 square degrees of sky every three days and we plan to increase the area to 1500 square degrees in the near future. When operating at full capacity, survey simulations show that YSE will find $\sim$5000 new SNe per year and at least two SNe within three days of explosion per month. To date, YSE has discovered or observed 8.3% of the transient candidates reported to the International Astronomical Union in 2020. We present an overview of YSE, including science goals, survey characteristics and a summary of our transient discoveries to date.

preprint2020arXiv

Luminous Red Nova AT 2019zhd, a new merger in M 31

We present the follow-up campaign of the luminous red nova (LRN) AT~2019zhd, the third event of this class observed in M 31. The object was followed by several sky surveys for about five months before the outburst, during which it showed a slow luminosity rise. In this phase, the absolute magnitude ranged from M_r=-2.8+-0.2 mag to M_r=-5.6+-0.1 mag. Then, over a four-five day period, AT 2019zhd experienced a major brightening, reaching at peak M_r=-9.61+-0.08 mag, and an optical luminosity of 1.4x10^39 erg/s. After a fast decline, the light curve settled onto a short-duration plateau in the red bands. Although less pronounced, this feature is reminiscent of the second red maximum observed in other LRNe. This phase was followed by a rapid linear decline in all bands. At maximum, the spectra show a blue continuum with prominent Balmer emission lines. The post-maximum spectra show a much redder continuum, resembling that of an intermediate-type star. In this phase, Halpha becomes very weak, Hbeta is no longer detectable and a forest of narrow absorption metal lines now dominate the spectrum. The latest spectra, obtained during the post-plateau decline, show a very red continuum (T_eff ~ 3000 K) with broad molecular bands of TiO, similar to those of M-type stars. The long-lasting, slow photometric rise observed before the peak resembles that of LRN V1309 Sco, which was interpreted as the signature of the common-envelope ejection. The subsequent outburst is likely due to the gas outflow following a stellar merging event. The inspection of archival HST images taken 22 years before the LRN discovery reveals a faint red source (M_F555W=0.21+-0.14 mag, with F555W-F814W = 2.96+-0.12 mag) at the position of AT 2019zhd, which is the most likely quiescent precursor. The source is consistent with expectations for a binary system including a predominant M5-type star.

preprint2019arXiv

Asteroid Discovery and Light Curve Extraction Using the Hough Transform -- A Rotation Period Study for Sub-Kilometer Main-Belt Asteroids

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

preprint2016arXiv

A Global Astrometric Solution for Pan-STARRS referenced to ICRF2

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

preprint2016arXiv

A Multi-Wavelength Photometric Census of AGN and Star Formation Activity in the Brightest Cluster Galaxies of X-ray Selected Clusters

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.

preprint2016arXiv

A search for an optical counterpart to the gravitational wave event GW151226

We present a search for an electromagnetic counterpart of the gravitational wave source GW151226. Using the Pan-STARRS1 telescope we mapped out 290 square degrees in the optical i_ps filter starting 11.5hr after the LIGO information release and lasting for a further 28 days. The first observations started 49.5hr after the time of the GW151226 detection. We typically reached sensitivity limits of i_ps = 20.3-20.8 and covered 26.5% of the LIGO probability skymap. We supplemented this with ATLAS survey data, reaching 31% of the probability region to shallower depths of m~19. We found 49 extragalactic transients (that are not obviously AGN), including a faint transient in a galaxy at 7Mpc (a luminous blue variable outburst) plus a rapidly decaying M-dwarf flare. Spectral classification of 20 other transient events showed them all to be supernovae. We found an unusual transient, PS15dpn, with an explosion date temporally coincident with GW151226 which evolved into a type Ibn supernova. The redshift of the transient is secure at z=0.1747 +/- 0.0001 and we find it unlikely to be linked, since the luminosity distance has a negligible probability of being consistent with that of GW151226. In the 290 square degrees surveyed we therefore do not find a likely counterpart. However we show that our survey strategy would be sensitive to NS-NS mergers producing kilonovae at D < 100 Mpc which is promising for future LIGO/Virgo searches.

preprint2016arXiv

Brightness variation distributions among main belt asteroids from sparse light curve sampling with Pan-STARRS 1

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

preprint2016arXiv

Dust in three dimensions in the Galactic Plane

We present three dimensional maps in monochromatic extinction $A_{\rm 0}$ and the extinction parameter $R_0$ within a few degrees of the Galactic plane. These are inferred using photometry from the Pan-STARRS1 and Spitzer Glimpse surveys of nearly $20$ million stars located in the region $l = 0-250^{\circ}$ and from $b = -4.5^{\circ}$ to $b=4.5^{\circ}$. Given the available stellar number density, we use an angular resolution of $7' \times 7'$ and steps of $1{\rm mag}$ in distance modulus. We simultaneously estimate distance modulus and effective temperature $T_{\rm eff}$ alongside the other parameters for stars individually using the method of \citet{Hanson2014} before combining these estimates to a complete map. The full maps are available via the MNRAS website.

preprint2016arXiv

Localization and broadband follow-up of the gravitational-wave transient GW150914

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.

preprint2016arXiv

Multiwavelength Observations of the Redback Millisecond Pulsar J1048+2339

We report on radio timing and multiwavelength observations of the 4.66 ms redback pulsar J1048+2339, which was discovered in an Arecibo search targeting the Fermi-LAT source 3FGLJ1048.6+2338. Two years of timing allowed us to derive precise astrometric and orbital parameters for the pulsar. PSR J1048+2339 is in a 6-hour binary, and exhibits radio eclipses over half the orbital period and rapid orbital period variations. The companion has a minimum mass of 0.3 solar masses, and we have identified a $V \sim 20$ variable optical counterpart in data from several surveys. The phasing of its $\sim 1$~mag modulation at the orbital period suggests highly efficient and asymmetric heating by the pulsar wind, which may be due to an intrabinary shock that is distorted near the companion, or to the companion's magnetic field channeling the pulsar wind to specific locations on its surface. We also present gamma-ray spectral analysis of the source and preliminary results from searches for gamma-ray pulsations using the radio ephemeris.

preprint2016arXiv

Of Genes and Machines: application of a combination of machine learning tools to astronomy datasets

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.

preprint2016arXiv

Pan-STARRS and PESSTO search for an optical counterpart to the LIGO gravitational wave source GW150914

We searched for an optical counterpart to the first gravitational wave source discovered by LIGO (GW150914), using a combination of the Pan-STARRS1 wide-field telescope and the PESSTO spectroscopic follow-up programme. As the final LIGO sky maps changed during analysis, the total probability of the source being spatially coincident with our fields was finally only 4.2 per cent. Therefore we discuss our results primarily as a demonstration of the survey capability of Pan-STARRS and spectroscopic capability of PESSTO. We mapped out 442 square degrees of the northern sky region of the initial map. We discovered 56 astrophysical transients over a period of 41 days from the discovery of the source. Of these, 19 were spectroscopically classified and a further 13 have host galaxy redshifts. All transients appear to be fairly normal supernovae and AGN variability and none is obviously linked with GW150914. We illustrate the sensitivity of our survey by defining parameterised lightcurves with timescales of 4, 20 and 40 days and use the sensitivity of the Pan-STARRS1 images to set limits on the luminosities of possible sources. The Pan-STARRS1 images reach limiting magnitudes of i = 19.2, 20.0 and 20.8 respectively for the three timescales. For long timescale parameterised lightcurves (with FWHM=~40d) we set upper limits of M_i <= -17.2 -0.9/+1.4 if the distance to GW150914 is D = 400 +/- 200Mpc. The number of type Ia SN we find in the survey is similar to that expected from the cosmic SN rate, indicating a reasonably complete efficiency in recovering supernova like transients out to D = 400 +/- 200 Mpc.

preprint2016arXiv

PS1-14bj: A Hydrogen-Poor Superluminous Supernova With a Long Rise and Slow Decay

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] $λ$5007 and [O III] $λ$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.

preprint2016arXiv

SN 2015bn: a detailed multi-wavelength view of a nearby superluminous supernova

We present observations of SN 2015bn (= PS15ae = CSS141223-113342+004332 = MLS150211-113342+004333), a Type I superluminous supernova (SLSN) at redshift $z=0.1136$. As well as being one of the closest SLSNe I yet discovered, it is intrinsically brighter ($M_U\approx-23.1$) and in a fainter galaxy ($M_B\approx-16.0$) than other SLSNe at $z\sim0.1$. We used this opportunity to collect the most extensive dataset for any SLSN I to date, including densely-sampled spectroscopy and photometry, from the UV to the NIR, spanning $-$50 to +250 days from optical maximum. SN 2015bn fades slowly, but exhibits surprising undulations in the light curve on a timescale of 30-50 days, especially in the UV. The spectrum shows extraordinarily slow evolution except for a rapid transformation between +7 and +20-30 days. No narrow emission lines from slow-moving material are observed at any phase. We derive physical properties including the bolometric luminosity, and find slow velocity evolution and non-monotonic temperature and radial evolution. A deep radio limit rules out a healthy off-axis gamma-ray burst, and places constraints on the pre-explosion mass loss. The data can be consistently explained by a $\gtrsim10\,{\rm M}_\odot$ stripped progenitor exploding with $\sim 10^{51}\,$erg kinetic energy, forming a magnetar with a spin-down timescale of $\sim20$ days (thus avoiding a gamma-ray burst) that reheats the ejecta and drives ionization fronts. The most likely alternative scenario -- interaction with $\sim20\,{\rm M}_\odot$ of dense, inhomogeneous circumstellar material -- can be tested with continuing radio follow-up.

preprint2016arXiv

Superluminous supernova 2015bn in the nebular phase: evidence for the engine-powered explosion of a stripped massive star

We present nebular-phase imaging and spectroscopy for the hydrogen-poor superluminous supernova SN 2015bn, at redshift z=0.1136, spanning +250-400 d after maximum light. The light curve exhibits a steepening in the decline rate from 1.4 mag/(100 d) to 1.7 mag/(100 d), suggestive of a significant decrease in the opacity. This change is accompanied by a transition from a blue continuum superposed with photospheric absorption lines to a nebular spectrum dominated by emission lines of oxygen, calcium and magnesium. There are no obvious signatures of circumstellar interaction or large nickel mass. We show that the spectrum at +400 d is virtually identical to a number of energetic Type Ic supernovae such as SN 1997dq, SN 2012au, and SN 1998bw, indicating similar core conditions and strengthening the link between `hypernovae'/long gamma-ray bursts and superluminous supernovae. A single explosion mechanism may unify these events that span absolute magnitudes of -22 < M_B < -17. Both the light curve and spectrum of SN 2015bn are consistent with an engine-driven explosion ejecting 7-30 M$_\odot$ of oxygen-dominated ejecta (for reasonable choices in temperature and opacity). A strong and relatively narrow O II $λ$7774 line, seen in a number of these energetic events but not in normal supernovae, may point to an inner shell that is the signature of a central engine.

preprint2016arXiv

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.

preprint2016arXiv

The Optical-Infrared Extinction Curve and its Variation in the Milky Way

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.

preprint2016arXiv

The Pan-STARRS1 distant z>5.6 quasar survey: more than 100 quasars within the first Gyr of the universe

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.

preprint2015arXiv

A supernova distance to the anchor galaxy NGC 4258

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

preprint2015arXiv

Constraining the radio-loud fraction of quasars at z>5.5

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_{ν,5{\rm GHz}}/f_{ν,4400Å}>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.

preprint2015arXiv

Evidence of Fanning in the Ophiuchus Stream

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

preprint2015arXiv

Galaxy Interactions in Compact Groups II: abundance and kinematic anomalies in HCG 91c

Galaxies in Hickson Compact Group 91 (HCG 91) were observed with the WiFeS integral field spectrograph as part of our ongoing campaign targeting the ionized gas physics and kinematics inside star forming members of compact groups. Here, we report the discovery of HII regions with abundance and kinematic offsets in the otherwise unremarkable star forming spiral HCG 91c. The optical emission line analysis of this galaxy reveals that at least three HII regions harbor an oxygen abundance ~0.15 dex lower than expected from their immediate surroundings and from the abundance gradient present in the inner regions of HCG 91c. The same star forming regions are also associated with a small kinematic offset in the form of a lag of 5-10 km/s with respect to the local circular rotation of the gas. HI observations of HCG 91 from the Very Large Array and broadband optical images from Pan-STARRS suggest that HCG 91c is caught early in its interaction with the other members of HCG 91. We discuss different scenarios to explain the origin of the peculiar star forming regions detected with WiFeS, and show that evidence point towards infalling and collapsing extra-planar gas clouds at the disk-halo interface, possibly as a consequence of long-range gravitational perturbations of HCG 91c from the other group members. As such, HCG 91c provides evidence that some of the perturbations possibly associated with the early phase of galaxy evolution in compact groups impact the star forming disk locally, and on sub-kpc scales.

preprint2015arXiv

GALEX Detection of Shock Breakout in Type II-P Supernova PS1-13arp: Implications for the Progenitor Star Wind

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

preprint2015arXiv

Machine learning for transient discovery in Pan-STARRS1 difference imaging

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

preprint2015arXiv

On the triple peaks of SNHunt248 in NGC 5806

We present our findings on a supernova (SN) impostor, SNHunt248, based on optical and near-IR data spanning $\sim$15 yrs before discovery, to $\sim$1 yr post-discovery. The light curve displays three distinct peaks, the brightest of which is at $M_{R} \sim -15.0$ mag. The post-discovery evolution is consistent with the ejecta from the outburst interacting with two distinct regions of circumstellar material. The 0.5 - 2.2 $μ$m spectral energy distribution at -740 d is well-matched by a single 6700 K blackbody with $\log(L/L_\odot) \sim 6.1$. This temperature and luminosity support previous suggestions of a yellow hypergiant progenitor; however, we find it to be brighter than the brightest and most massive Galactic late-F to early-G spectral type hypergiants. Overall the historical light curve displays variability of up to $\sim \pm1$ mag. At current epochs ($\sim$1 yr post-outburst), the absolute magnitude ($M_{R} \sim -9$ mag) is just below the faintest observed historical absolute magnitude $\sim$10 yrs before discovery.

preprint2015arXiv

Optical Confirmation and Redshift Estimation of the Planck Cluster Candidates overlapping the Pan-STARRS Survey

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 $σ_{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.

preprint2015arXiv

Pan-Planets: Searching for hot Jupiters around cool dwarfs

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.

preprint2015arXiv

Pan-STARRS1 variability of XMM-COSMOS AGN. I. Impact on photometric redshifts

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

preprint2015arXiv

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.

preprint2015arXiv

Selection of Burst-like Transients and Stochastic Variables Using Multi-Band Image Differencing in the Pan-STARRS1 Medium-Deep Survey

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

preprint2015arXiv

The Globular Cluster System of NGC 6822

We present a comprehensive analysis of the globular cluster (GC) system of the Local Group dwarf irregular galaxy NGC 6822. Our study is based on homogeneous optical and near-IR photometry, as well as long-slit spectroscopic observations which are used to determine new radial velocities for 6 GCs, two of which had no previous spectroscopic information. We construct optical-near IR colour-colour diagrams and through comparison to simple stellar population models infer that the GCs have old ages consistent with being 9 Gyr or older, while their metallicities are in the range between -1.6 < [Fe/H] < -0.4. We conduct a kinematic analysis of the GC population and find tentative evidence for weak net rotation of the GC system, in the same sense as that exhibited by the underlying spheroid. The most likely amplitude of rotation is ~10 km/s, approximately half the magnitude of the observed velocity dispersion. Finally, we use the GCs to estimate the dynamical mass of NGC 6822 within 11 kpc and we formally find it to be in the range between (3-4)10^9 Msun. This implies an overall mass-to-light ratio in the range of ~ 30-40 and indicates that NGC 6822 is highly dark matter dominated. The mass and the corresponding mass-to-light ratio estimates are affected by various additional systematic effects due to limitations of the data and the model that are not necessary reflected in the formal uncertainties.

preprint2015arXiv

The identification of z-dropouts in Pan-STARRS1: three quasars at 6.5<z<6.7

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.

preprint2015arXiv

The Nature and Orbit of the Ophiuchus Stream

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

preprint2015arXiv

Total eclipse of the heart: The AM CVn Gaia14aae / ASSASN-14cn

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

preprint2015arXiv

Zooming In on the Progenitors of Superluminous Supernovae With the HST

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

preprint2014arXiv

3D Dust Mapping Reveals that Orion Forms Part of a Large Ring of Dust

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.

preprint2014arXiv

A Large Catalog of Accurate Distances to Molecular Clouds from PS1 Photometry

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.

preprint2014arXiv

A Map of Dust Reddening to 4.5 kpc from Pan-STARRS1

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

preprint2014arXiv

A search for planetary eclipses of white dwarfs in the Pan-STARRS1 medium-deep fields

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

preprint2014arXiv

Characterization of the Praesepe Star Cluster by Photometry and Proper Motions with 2MASS, PPMXL, and Pan-STARRS

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.

preprint2014arXiv

Cosmological Constraints from Measurements of Type Ia Supernovae discovered during the first 1.5 years of the Pan-STARRS1 Survey

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 $Ω_{\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$σ$ 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σ$. We cannot conclude whether the tension with flat $Λ$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.

preprint2014arXiv

Discovery of eight z ~ 6 quasars from Pan-STARRS1

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

preprint2014arXiv

Hydrogen-Poor Superluminous Supernovae and Long-Duration Gamma-Ray Bursts Have Similar Host Galaxies

We present optical spectroscopy and optical/near-IR photometry of 31 host galaxies of hydrogen-poor superluminous supernovae (SLSNe), including 15 events from the Pan-STARRS1 Medium Deep Survey. Our sample spans the redshift range 0.1 < z < 1.6 and is the first comprehensive host galaxy study of this specific subclass of cosmic explosions. Combining the multi-band photometry and emission-line measurements, we determine the luminosities, stellar masses, star formation rates and metallicities. We find that as a whole, the hosts of SLSNe are a low-luminosity (<M_B> ~ -17.3 mag), low stellar mass (<M_*> ~ 2 x 10^8 M_sun) population, with a high median specific star formation rate (<sSFR> ~ 2 Gyr^-1). The median metallicity of our spectroscopic sample is low, 12 + log(O/H}) ~ 8.35 ~ 0.45 Z_sun, although at least one host galaxy has solar metallicity. The host galaxies of H-poor SLSNe are statistically distinct from the hosts of GOODS core-collapse SNe (which cover a similar redshift range), but resemble the host galaxies of long-duration gamma-ray bursts (LGRBs) in terms of stellar mass, SFR, sSFR and metallicity. This result indicates that the environmental causes leading to massive stars forming either SLSNe or LGRBs are similar, and in particular that SLSNe are more effectively formed in low metallicity environments. We speculate that the key ingredient is large core angular momentum, leading to a rapidly-spinning magnetar in SLSNe and an accreting black hole in LGRBs.

preprint2014arXiv

Measuring Quasar Variability with Pan-STARRS1 and SDSS

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

preprint2014arXiv

Observational Constraints on the Catastrophic Disruption Rate of Small Main Belt Asteroids

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.

preprint2014arXiv

Pan-STARRS 1 observations of the unusual active Centaur P/2011 S1(Gibbs)

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

preprint2014arXiv

Properties of M31. IV: Candidate Luminous Blue Variables from PAndromeda

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.

preprint2014arXiv

Rapidly-Evolving and Luminous Transients from Pan-STARRS1

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

preprint2014arXiv

Systematic Uncertainties Associated with the Cosmological Analysis of the First Pan-STARRS1 Type Ia Supernova Sample

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 $Ω_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}$.

preprint2013arXiv

A Search for Fast Optical Transients in the Pan-STARRS1 Medium-Deep Survey: M Dwarf Flares, Asteroids, Limits on Extragalactic Rates, and Implications for LSST

[Abridged] We present a search for fast optical transients (~0.5 hr-1 day) using repeated observations of the Pan-STARRS1 Medium-Deep Survey (PS1/MDS) fields. Our search takes advantage of the consecutive g/r-band observations (16.5 min in each filter), by requiring detections in both bands, with non-detections on preceding and subsequent nights. We identify 19 transients brighter than 22.5 AB mag (S/N>10). Of these, 11 events exhibit quiescent counterparts in the deep PS1/MDS templates that we identify as M4-M9 dwarfs. The remaining 8 transients exhibit a range of properties indicative of main-belt asteroids near the stationary point of their orbits. With identifications for all 19 transients, we place an upper limit of R_FOT(0.5hr)<0.12 deg^-2 d^-1 (95% confidence level) on the sky-projected rate of extragalactic fast transients at <22.5 mag, a factor of 30-50 times lower than previous limits; the limit for a timescale of ~day is R_FOT<2.4e-3 deg^-2 d^-1. To convert these sky-projected rates to volumetric rates, we explore the expected peak luminosities of fast optical transients powered by various mechanisms, and find that non-relativistic events are limited to M~-10 mag (M~-14 mag) for a timescale of ~0.5 hr (~day), while relativistic sources (e.g., GRBs, magnetar-powered transients) can reach much larger luminosities. The resulting volumetric rates are <13 (M~-10 mag), <0.05 (M~-14 mag) and <1e-6 Mpc^-3 yr^-1 (M~-24 mag), significantly above the nova, supernova, and GRB rates, respectively, indicating that much larger surveys are required to provide meaningful constraints. Motivated by the results of our search we discuss strategies for identifying fast optical transients in the LSST main survey, and reach the optimistic conclusion that the veil of foreground contaminants can be lifted with the survey data, without the need for expensive follow-up observations.

preprint2013arXiv

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.

preprint2013arXiv

Main-Belt Comet P/2012 T1 (PANSTARRS)

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

preprint2013arXiv

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π$ 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π$ survey for studies of large scale structure.

preprint2013arXiv

Properties of M31. II: A Cepheid disk sample derived from the first year of PS1 PAndromeda data

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

preprint2013arXiv

Properties of M31. III: Candidate Beat Cepheids from PS1 PAndromeda Data and Their Implication on Metallicity Gradient

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

preprint2013arXiv

PS1-10afx at z=1.388: Pan-STARRS1 Discovery of a New Type of Superluminous Supernova

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

preprint2013arXiv

PS1-10bzj: A Fast, Hydrogen-Poor Superluminous Supernova in a Metal Poor Host Galaxy

We present observations and analysis of PS1-10bzj, a superluminous supernova (SLSN) discovered in the Pan-STARRS Medium Deep Survey at a redshift z = 0.650. Spectroscopically, PS1-10bzj is similar to the hydrogen-poor SLSNe 2005ap and SCP 06F6, though with a steeper rise and lower peak luminosity (M_bol = -21.4 mag) than previous events. We construct a bolometric light curve, and show that while PS1-10bzj's energetics were less extreme than previous events, its luminosity still cannot be explained by radioactive nickel decay alone. We explore both a magnetar spin-down and circumstellar interaction scenario and find that either can fit the data. PS1-10bzj is located in the Extended Chandra Deep Field South and the host galaxy is imaged in a number of surveys, including with the Hubble Space Telescope. The host is a compact dwarf galaxy (M_B ~ -18 mag, diameter < 800 pc), with a low stellar mass (M_* ~ 2.4 * 10^7 M_sun), young stellar population (τ_* ~ 5 Myr), and a star formation rate of ~ 2-3 M_sun/yr. The specific star formation rate is the highest seen in a SLSN host so far (~ 100 Gyr^{-1}). We detect the [O III]λ4363 line, and find a low metallicity: 12+(O/H) = 7.8 +/- 0.2 (~ 0.1 Z_sun). Together, this indicates that at least some of the progenitors of SLSNe come from young, low-metallicity populations.

preprint2013arXiv

PS1-12sk is a Peculiar Supernova From a He-rich Progenitor System in a Brightest Cluster Galaxy Environment

We report on our discovery and observations of the Pan-STARRS1 supernova (SN) PS1-12sk, a transient with properties that indicate atypical star formation in its host galaxy cluster or pose a challenge to popular progenitor system models for this class of explosion. The optical spectra of PS1-12sk classify it as a Type Ibn SN (c.f. SN 2006jc), dominated by intermediate-width (3x10^3 km/s) and time variable He I emission. Our multi-wavelength monitoring establishes the rise time dt = 9-23 days and shows an NUV-NIR SED with temperature > 17x10^3 K and a peak rise magnitude of Mz = -18.9 mag. SN Ibn spectroscopic properties are commonly interpreted as the signature of a massive star (17 - 100 M_sun) explosion within a He-enriched circumstellar medium. However, unlike previous Type Ibn supernovae, PS1-12sk is associated with an elliptical brightest cluster galaxy, CGCG 208-042 (z = 0.054) in cluster RXC J0844.9+4258. The expected probability of an event like PS1-12sk in such environments is low given the measured infrequency of core-collapse SNe in red sequence galaxies compounded by the low volumetric rate of SN Ibn. Furthermore, we find no evidence of star formation at the explosion site to sensitive limits (Sigma Halpha < 2x10^-3 M_sun/yr/kpc^2). We therefore discuss white dwarf binary systems as a possible progenitor channel for SNe Ibn. We conclude that PS1-12sk represents either a fortuitous and statistically unlikely discovery, evidence for a top-heavy IMF in galaxy cluster cooling flow filaments, or the first clue suggesting an alternate progenitor channel for Type Ibn SNe.

preprint2013arXiv

Stacked reverberation lags at high redshift

Over the past 20years reverberation mapping has proved one of the most successful techniques for studying the local (<1pc) environment of super-massive black holes that drive active galactic nuclei. Key successes of reverberation mapping have been direct black hole mass estimates, the radius-luminosity relation for the Hbeta line and the calibration of single-epoch mass estimators commonly employed up to z~7. However, observing constraints mean that few studies have been successful at z>0.1, or for the more-luminous quasars that make up the majority of current spectroscopic samples, or for the rest-frame ultra-violet emission lines available in optical spectra of z>0.5 objects. Previously we described a technique for stacking cross correlations to obtain reverberation mapping results at high z. Here we present the first results from a campaign designed for this purpose. We construct stacked cross-correlation functions for the CIV and MgII lines and find a clear peak in both. We find the peak in the MgII correlation is at longer lags than CIV consistent with previous results at low redshift. For the CIV sample we are able to bin by luminosity and find evidence for increasing lags for more-luminous objects. This CIV radius-luminosity relation is consistent with previous studies but with a fraction of the observational cost.

preprint2013arXiv

Super Luminous Ic Supernovae: catching a magnetar by the tail

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

preprint2013arXiv

The Extremely Red, Young L Dwarf PSO J318-22: A Free-Floating Planetary-Mass Analog to Directly Imaged Young Gas-Giant Planets

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.

preprint2013arXiv

The GALEX Time Domain Survey I. Selection and Classification of Over a Thousand UV Variable Sources

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

preprint2013arXiv

The Pan-STARRS 1 Photometric Reference Ladder, Release 12.0

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

preprint2013arXiv

The Pan-STARRS Moving Object Processing System

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.

preprint2013arXiv

The Pan-STARRS1 Medium-Deep Survey: The role of galaxy group environment in the star formation rate versus stellar mass relation and quiescent fraction out to $z \sim 0.8$

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$σ$ 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)

preprint2013arXiv

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

preprint2013arXiv

The superluminous supernova PS1-11ap: bridging the gap between low and high redshift

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.

preprint2013arXiv

The UV-bright, Slowly Declining Transient PS1-11af as a Partial Tidal Disruption Event

We present the Pan-STARRS1 discovery of the long-lived and blue transient PS1-11af, which was also detected by GALEX with coordinated observations in the near-ultraviolet (NUV) band. PS1-11af is associated with the nucleus of an early-type galaxy at redshift z=0.4046 that exhibits no evidence for star formation or AGN activity. Four epochs of spectroscopy reveal a pair of transient broad absorption features in the UV on otherwise featureless spectra. Despite the superficial similarity of these features to P-Cygni absorptions of supernovae (SNe), we conclude that PS1-11af is not consistent with the properties of known types of SNe. Blackbody fits to the spectral energy distribution are inconsistent with the cooling, expanding ejecta of a SN, and the velocities of the absorption features are too high to represent material in homologous expansion near a SN photosphere. However, the constant blue colors and slow evolution of the luminosity are similar to previous optically-selected tidal disruption events (TDEs). The shape of the optical light curve is consistent with models for TDEs, but the minimum accreted mass necessary to power the observed luminosity is only ~0.002M_sun, which points to a partial disruption model. A full disruption model predicts higher bolometric luminosities, which would require most of the radiation to be emitted in a separate component at high energies where we lack observations. In addition, the observed temperature is lower than that predicted by pure accretion disk models for TDEs and requires reprocessing to a constant, lower temperature. Three deep non-detections in the radio with the VLA over the first two years after the event set strict limits on the production of any relativistic outflow comparable to Swift J1644+57, even if off-axis.

preprint2013arXiv

X-ray selected galaxy clusters in the Pan-STARRS Medium-Deep Survey

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

preprint2012arXiv

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.

preprint2012arXiv

An ultraviolet-optical flare from the tidal disruption of a helium-rich stellar core

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.

preprint2012arXiv

Discovery of Main-Belt Comet P/2006 VW139 by 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.

preprint2012arXiv

HIP 38939B: A New Benchmark T Dwarf in the Galactic Plane Discovered with Pan-STARRS1

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

preprint2012arXiv

Observational and Dynamical Characterization of Main-Belt Comet P/2010 R2 (La Sagra)

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.

preprint2012arXiv

PAndromeda - first results from the high-cadence monitoring of M31 with Pan-STARRS 1

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.

preprint2012arXiv

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.

preprint2012arXiv

SN 2010ay is a Luminous and Broad-lined Type Ic Supernova within a Low-metallicity Host Galaxy

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.

preprint2012arXiv

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.

preprint2011arXiv

A Search for High Proper Motion T Dwarfs with Pan-STARRS1 + 2MASS + WISE

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

preprint2011arXiv

First Results from Pan-STARRS1: Faint, High Proper Motion White Dwarfs in the Medium-Deep Fields

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

preprint2011arXiv

Four new T dwarfs identified in PanSTARRS 1 commissioning data

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

preprint2011arXiv

Pan-STARRS1 Discovery of Two Ultra-Luminous Supernovae at z ~ 0.9

We present the discovery of two ultra-luminous supernovae (SNe) at z ~ 0.9 with the Pan-STARRS1 Medium-Deep Survey. These SNe, PS1-10ky and PS1-10awh, are amongst the most luminous SNe ever discovered, comparable to the unusual transients SN 2005ap and SCP 06F6. Like SN 2005ap and SCP 06F6, they show characteristic high luminosities (M_bol ~ -22.5 mag), blue spectra with a few broad absorption lines, and no evidence for H or He. We have constructed a full multi-color light curve sensitive to the peak of the spectral energy distribution in the rest-frame ultraviolet, and we have obtained time-series spectroscopy for these SNe. Given the similarities between the SNe, we combine their light curves to estimate a total radiated energy over the course of explosion of (0.9-1.4) x 10^51 erg. We find photospheric velocities of 12,000-19,000 km/s with no evidence for deceleration measured across ~3 rest-frame weeks around light-curve peak, consistent with the expansion of an optically-thick massive shell of material. We show that, consistent with findings for other ultra-luminous SNe in this class, radioactive decay is not sufficient to power PS1-10ky, and we discuss two plausible origins for these events: the initial spin-down of a newborn magnetar in a core-collapse SN, or SN shock breakout from the dense circumstellar wind surrounding a Wolf-Rayet star.

preprint2011arXiv

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.

preprint2010arXiv

Displaying the Heterogeneity of the SN 2002cx-like Subclass of Type Ia Supernovae with Observations of the Pan-STARRS-1 Discovered SN2009ku

SN2009ku, discovered by Pan-STARRS-1, is a Type Ia supernova (SNIa), and a member of the distinct SN2002cx-like class of SNeIa. Its light curves are similar to the prototypical SN2002cx, but are slightly broader and have a later rise to maximum in g. SN2009ku is brighter (~0.6 mag) than other SN2002cx-like objects, peaking at M_V = -18.4 mag - which is still significantly fainter than typical SNeIa. SN2009ku, which had an ejecta velocity of ~2000 kms^-1 at 18 days after maximum brightness is spectroscopically most similar to SN2008ha, which also had extremely low-velocity ejecta. However, SN2008ha had an exceedingly low luminosity, peaking at M_V = -14.2 mag, ~4 mag fainter than SN2009ku. The contrast of high luminosity and low ejecta velocity for SN2009ku is contrary to an emerging trend seen for the SN2002cx class. SN2009ku is a counter-example of a previously held belief that the class was more homogeneous than typical SNeIa, indicating that the class has a diverse progenitor population and/or complicated explosion physics. As the first example of a member of this class of objects from the new generation of transient surveys, SN2009ku is an indication of the potential for these surveys to find rare and interesting objects.

preprint2010arXiv

GALEX and Pan-STARRS1 Discovery of SN IIP 2010aq: The First Few Days After Shock Breakout in a Red Supergiant Star

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.

preprint2010arXiv

SN 2009kf : a UV bright type IIP supernova discovered with Pan-STARRS 1 and GALEX

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.

preprint2010arXiv

Ultra-bright optical transients are linked with type Ic supernovae

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.

preprint2009arXiv

58 Radio Sources Near Bright Natural Guide Stars

We present a preliminary survey of 58 radio sources within the isoplanatic patches (r < 25") of bright (11<R<12) stars suitable for use as natural guide stars with high-order adaptive optics (AO). An optical and near-infrared imaging survey was conducted utilizing tip-tilt corrections in the optical and AO in the near-infrared. Spectral Energy Distributions (SEDs) were fit to the multi-band data for the purpose of obtaining photometric redshifts using the Hyperz code. Several of these photometric redshifts were confirmed with spectroscopy, a result that gives more confidence to the redshift distribution for the whole sample. Additional long-wavelength data from Spitzer, SCUBA, SHARC2, and VLA supplement the optical and near-infrared data. We find the sample generally follows and extends the magnitude-redshift relation found for more powerful local radio galaxies. The survey has identified several reasonably bright (H=19-20) objects at significant redshifts (z>1) that are now within the capabilities of the current generation of AO-fed integral-field spectrographs. These objects constitute a unique sample that can be used for detailed ground-based AO studies of galactic structure, evolution, and AGN formation at high redshift.