Source author record

H. Flewelling

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

52works
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

52 published item(s)

preprint2021arXiv

Comparison of the physical properties of the L4 and L5 Trojan asteroids from ATLAS data

Jupiter has nearly 8000~known co-orbital asteroids orbiting in the L4 and L5 Lagrange points called Jupiter Trojan asteroids. Aside from the greater number density of the L4 cloud the two clouds are in many ways considered to be identical. Using sparse photometric data taken by the Asteroid Terrestrial-impact Last Alert System (ATLAS) for 863 L4 Trojans and 380 L5 Trojans we derive the shape distribution for each of the clouds and find that, on average, the L4 asteroids are more elongated than the L5 asteroids. This shape difference is most likely due to the greater collision rate in the L4 cloud that results from its larger population. We additionally present the phase functions and $c-o$ colours of 266~objects.

preprint2020arXiv

Design and operation of the ATLAS Transient Science Server

The Asteroid Terrestrial impact Last Alert System (ATLAS) system consists of two 0.5m Schmidt telescopes with cameras covering 29 square degrees at plate scale of 1.86 arcsec per pixel. Working in tandem, the telescopes routinely survey the whole sky visible from Hawaii (above $δ> -50^{\circ}$) every two nights, exposing four times per night, typically reaching $o < 19$ magnitude per exposure when the moon is illuminated and $c < 19.5$ per exposure in dark skies. Construction is underway of two further units to be sited in Chile and South Africa which will result in an all-sky daily cadence from 2021. Initially designed for detecting potentially hazardous near earth objects, the ATLAS data enable a range of astrophysical time domain science. To extract transients from the data stream requires a computing system to process the data, assimilate detections in time and space and associate them with known astrophysical sources. Here we describe the hardware and software infrastructure to produce a stream of clean, real, astrophysical transients in real time. This involves machine learning and boosted decision tree algorithms to identify extragalactic and Galactic transients. Typically we detect 10-15 supernova candidates per night which we immediately announce publicly. The ATLAS discoveries not only enable rapid follow-up of interesting sources but will provide complete statistical samples within the local volume of 100 Mpc. A simple comparison of the detected supernova rate within 100 Mpc, with no corrections for completeness, is already significantly higher (factor 1.5 to 2) than the current accepted rates.

preprint2020arXiv

Investigating Taxonomic Diversity within Asteroid Families through ATLAS Dual-Band Photometry

We present here the c-o colors for identified Flora, Vesta, Nysa-Polana, Themis, and Koronis family members within the historic data set (2015-2018) of the Asteroid Terrestrial-impact Last Alert System (ATLAS). The Themis and Koronis families are known to be relatively pure C- and S-type Bus-DeMeo taxonomic families, respectively, and the extracted color data from the ATLAS broadband c- and o-filters of these two families is used to demonstrate that the ATLAS c-o color is a sufficient parameter to distinguish between the C- and S-type taxonomies. The Vesta and Nysa-Polana families are known to display a mixture of taxonomies possibly due to Vesta's differentiated parent body origin and Nysa-Polana actually consisting of two nested families with differing taxonomies. Our data show that the Flora family also displays a large degree of taxonomic mixing and the data reveal a substantial H-magnitude dependence on color. We propose and exclude several interpretations for the observed taxonomic mix. Additionally, we extract rotation periods of all of the targets reported here and find good agreement with targets that have previously reported periods.

preprint2019arXiv

Search for transient optical counterparts to high-energy IceCube neutrinos with Pan-STARRS1

In order to identify the sources of the observed diffuse high-energy neutrino flux, it is crucial to discover their electromagnetic counterparts. IceCube began releasing alerts for single high-energy ($E > 60$ TeV) neutrino detections with sky localisation regions of order 1 deg radius in 2016. We used Pan-STARRS1 to follow-up five of these alerts during 2016-2017 to search for any optical transients that may be related to the neutrinos. Typically 10-20 faint ($m < 22.5$ mag) extragalactic transients are found within the Pan-STARRS1 footprints and are generally consistent with being unrelated field supernovae (SNe) and AGN. We looked for unusual properties of the detected transients, such as temporal coincidence of explosion epoch with the IceCube timestamp. We found only one transient that had properties worthy of a specific follow-up. In the Pan-STARRS1 imaging for IceCube-160427A (probability to be of astrophysical origin of $\sim$50 %), we found a SN PS16cgx, located at 10.0' from the nominal IceCube direction. Spectroscopic observations of PS16cgx showed that it was an H-poor SN at z = 0.2895. The spectra and light curve resemble some high-energy Type Ic SNe, raising the possibility of a jet driven SN with an explosion epoch temporally coincident with the neutrino detection. However, distinguishing Type Ia and Type Ic SNe at this redshift is notoriously difficult. Based on all available data we conclude that the transient is more likely to be a Type Ia with relatively weak SiII absorption and a fairly normal rest-frame r-band light curve. If, as predicted, there is no high-energy neutrino emission from Type Ia SNe, then PS16cgx must be a random coincidence, and unrelated to the IceCube-160427A. We find no other plausible optical transient for any of the five IceCube events observed down to a 5$σ$ limiting magnitude of $m \sim 22$ mag, between 1 day and 25 days after detection.

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

Distant Solar System Objects identified in the Pan-STARRS1 survey

We present a method to identify distant solar system objects in long-term wide-field asteroid survey data, and conduct a search for them in the Pan-STARRS1 (PS1) image data acquired from 2010 to mid-2015. We demonstrate that our method is able to find multi-opposition orbital links, and we present the resulting orbital distributions which consist of 154 Centaurs, 255 classical Trans-Neptunian Objects (TNOs), 121 resonant TNOs, 89 Scattered Disc Objects (SDOs) and 10 comets. Our results show more than half of these are new discoveries, including a newly discovered 19th magnitude TNO. Our identified objects do not show clustering in their argument of perihelia, which if present, might support the existence of a large unknown planetary-sized object in the outer solar system.

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

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

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 host galaxy properties of variability selected AGN in the Pan-STARRS1 Medium-Deep Survey

We study the properties of 975 active galactic nuclei (AGN) selected by variability in the Pan-STARRS1 Medium-Deep Survey. Using complementary multi wavelength data from the ultraviolet to the far-infrared, we use SED fitting to determine the AGN and host properties at $z<1$, and compare to a well-matched control sample. We confirm the trend previously observed that the variability amplitude decreases with AGN luminosity, but on the other hand, we observe that the slope of this relation steepens with wavelength resulting in a "redder when brighter" trend at low luminosities. Our results show that AGN are hosted by more massive hosts than control sample galaxies, while the restframe, dust-corrected $NUV-r$ color distribution of AGN hosts is similar to control galaxies. We find a positive correlation between the AGN luminosity and star formation rate (SFR), independent of redshift. AGN hosts populate the whole range of SFRs within and outside the Main Sequence of star forming galaxies. Comparing the distribution of AGN hosts and control galaxies, we show that AGN hosts are less likely to be hosted by quiescent galaxies, but more likely to be hosted by Main Sequence or starburst galaxies.

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

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

Finding, characterizing and classifying variable sources in multi-epoch sky surveys: QSOs and RR Lyrae in PS1 3$π$ data

In area and depth, the Pan-STARRS1 (PS1) 3$π$ survey is unique among many-epoch, multi-band surveys and has enormous potential for all-sky identification of variable sources. PS1 has observed the sky typically seven times in each of its five bands ($grizy$) over 3.5 years, but unlike SDSS not simultaneously across the bands. Here we develop a new approach for quantifying statistical properties of non-simultaneous, sparse, multi-color lightcurves through light-curve structure functions, effectively turning PS1 into a $\sim 35$-epoch survey. We use this approach to estimate variability amplitudes and timescales $(ω_r, τ)$ for all point-sources brighter than $r_{\mathrm{P1}}=21.5$ mag in the survey. With PS1 data on SDSS Stripe 82 as ``ground truth", we use a Random Forest Classifier to identify QSOs and RR Lyrae based on their variability and their mean PS1 and WISE colors. We find that, aside from the Galactic plane, QSO and RR Lyrae samples of purity $\sim$75\% and completeness $\sim$92\% can be selected. On this basis we have identified a sample of $\sim 1,000,000$ QSO candidates, as well as an unprecedentedly large and deep sample of $\sim$150,000 RR Lyrae candidates with distances from $\sim$10 kpc to $\sim$120 kpc. Within the Draco dwarf spheroidal, we demonstrate a distance precision of 6\% for RR Lyrae candidates. We provide a catalog of all likely variable point sources and likely QSOs in PS1, a total of $25.8\times 10^6$ sources.

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

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 fastest unbound star in our Galaxy ejected by a thermonuclear supernova

Hypervelocity stars (HVS) travel with velocities so high, that they exceed the escape velocity of the Galaxy. Several acceleration mechanisms have been discussed. Only one HVS (US 708, HVS 2) is a compact helium star. Here we present a spectroscopic and kinematic analysis of US\,708. Travelling with a velocity of $\sim1200\,{\rm km\,s^{-1}}$, it is the fastest unbound star in our Galaxy. In reconstructing its trajectory, the Galactic center becomes very unlikely as an origin, which is hardly consistent with the most favored ejection mechanism for the other HVS. Furthermore, we discovered US\,708 to be a fast rotator. According to our binary evolution model it was spun-up by tidal interaction in a close binary and is likely to be the ejected donor remnant of a thermonuclear supernova.

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

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

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

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

Newly Discovered RR Lyrae Stars in the SDSSXPanXSTARRS1XCatalina Footprint

We present the detection of 6,371 RR Lyrae (RRL) stars distributed across ~14,000 deg^2 of the sky from the combined data of the Sloan Digital Sky Survey (SDSS), the Panoramic Survey Telescope and Rapid Response System 1 (PS1), and the second photometric catalogue from the Catalina Survey (CSDR2), out of these, ~2,021 RRL stars (~572 RRab and 1,449 RRc) are new discoveries. The RRL stars have heliocentric distances in the 4--28 kpc distance range. RRL-like color cuts from the SDSS and variability cuts from the PS1 are used to cull our candidate list. We then use the CSDR2 multi-epoch data to refine our sample. Periods were measured using the Analysis of Variance technique while the classification process is performed with the Template Fitting Method in addition to the visual inspection of the light curves. A cross-match of our RRL star discoveries with previous published catalogs of RRL stars yield completeness levels of ~50% for both RRab and RRc stars, and an efficiency of ~99% and ~87% for RRab and RRc stars, respectively. We show that our method for selecting RRL stars allows us to recover halo structures. The full lists of all the RRL stars are made publicly available.

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

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

preprint2014arXiv

Wide, Cool and Ultracool Companions to Nearby Stars from Pan-STARRS1

We present the discovery of 61 wide (>5 arcsecond) separation, low-mass (stellar and substellar) companions to stars in the solar neighborhood identified from Pan-STARRS\,1 (PS1) data and the spectral classification of 27 previously known companions. Our companions represent a selective subsample of promising candidates and span a range in spectral type of K7-L9 with the addition of one DA white dwarf. These were identified primarily from a dedicated common proper motion search around nearby stars, along with a few as serendipitous discoveries from our Pan-STARRS1 brown dwarf search. Our discoveries include 24 new L dwarf companions and one known L dwarf not previously identified as a companion. The primary stars around which we searched for companions come from a list of bright stars with well-measured parallaxes and large proper motions from the Hipparcos catalog (8583 stars, mostly A-K~dwarfs) and fainter stars from other proper motion catalogues (79170 stars, mostly M~dwarfs). We examine the likelihood that our companions are chance alignments between unrelated stars and conclude that this is unlikely for the majority of the objects that we have followed-up spectroscopically. We also examine the entire population of ultracool (>M7) dwarf companions and conclude that while some are loosely bound, most are unlikely to be disrupted over the course of $\sim$10 Gyr. Our search increases the number of ultracool M dwarf companions wider than 300 AU by 88% and increases the number of L dwarf companions in the same separation range by 96%. Finally, we resolve our new L dwarf companion to HIP 6407 into a tight (0.13 arcsecond, 7.4 AU) L1+T3 binary, making the system a hierarchical triple. Our search for these key benchmarks against which brown dwarf and exoplanet atmosphere models are tested has yielded the largest number of discoveries to date.

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

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

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

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

Panchromatic observations of the textbook GRB 110205A: constraining physical mechanisms of prompt emission and afterglow

We present a comprehensive analysis of a bright, long duration (T90 ~ 257 s) GRB 110205A at redshift z= 2.22. The optical prompt emission was detected by Swift/UVOT, ROTSE-IIIb and BOOTES telescopes when the GRB was still radiating in the gamma-ray band. Nearly 200 s of observations were obtained simultaneously from optical, X-ray to gamma-ray, which makes it one of the exceptional cases to study the broadband spectral energy distribution across 6 orders of magnitude in energy during the prompt emission phase. By fitting the time resolved prompt spectra, we clearly identify, for the first time, an interesting two-break energy spectrum, roughly consistent with the standard GRB synchrotron emission model in the fast cooling regime. Although the prompt optical emission is brighter than the extrapolation of the best fit X/gamma-ray spectra, it traces the gamma-ray light curve shape, suggesting a relation to the prompt high energy emission. The synchrotron + SSC scenario is disfavored by the data, but the models invoking a pair of internal shocks or having two emission regions can interpret the data well. Shortly after prompt emission (~ 1100 s), a bright (R = 14.0) optical emission hump with very steep rise (alpha ~ 5.5) was observed which we interpret as the emission from the reverse shock. It is the first time that the rising phase of a reverse shock component has been closely observed. The full optical and X-ray afterglow lightcurves can be interpreted within the standard reverse shock (RS) + forward shock (FS) model. In general, the high quality prompt emission and afterglow data allow us to apply the standard fireball shock model to extract valuable information about the GRB including the radiation mechanism, radius of prompt emission R, initial Lorentz factor of the outflow, the composition of the ejecta, as well as the collimation angle and the total energy budget.

preprint2012arXiv

The origin of the early time optical emission of Swift GRB 080310

We present broadband multi-wavelength observations of GRB 080310 at redshift z = 2.43. This burst was bright and long-lived, and unusual in having extensive optical and near IR follow-up during the prompt phase. Using these data we attempt to simultaneously model the gamma-ray, X-ray, optical and IR emission using a series of prompt pulses and an afterglow component. Initial attempts to extrapolate the high energy model directly to lower energies for each pulse reveal that a spectral break is required between the optical regime and 0.3 keV to avoid over predicting the optical flux. We demonstrate that afterglow emission alone is insufficient to describe all morphology seen in the optical and IR data. Allowing the prompt component to dominate the early-time optical and IR and permitting each pulse to have an independent low energy spectral indices we produce an alternative scenario which better describes the optical light curve. This, however, does not describe the spectral shape of GRB 080310 at early times. The fit statistics for the prompt and afterglow dominated models are nearly identical making it difficult to favour either. However one enduring result is that both models require a low energy spectral index consistent with self absorption for at least some of the pulses identified in the high energy emission model.

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.

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.