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 similar to 32 000 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 20 x 20 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 per cent 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 per cent, the classifier initially suggests a missed detection rate of around 10 per cent. 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 per cent.

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 five bands (grizyP1) 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 superluminous supernovae (SLSNe). A total of 249 hostless transients were discovered down to a limiting magnitude of MAB ∼ 23.5, of which 76 were classified as Type Ia supernovae (SNe Ia). There were 57 SNe with complete light curves that are likely core-collapse SNe (CCSNe) or type Ic 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 and 1.4. This illustrates that the discovery rate of type Ic SLSNe can be maximized by concentrating on hostless transients and removing normal SNe Ia. We present data for two possible SLSNe; PS1-10pm (z = 1.206) and PS1-10ahf (z = 1.1), and estimate the rate of type Ic SLSNe to be between 3+3−2×10−53+3−2×10−53+3−2×10−5 and 8+2−1×10−58+2−1×10−58+2−1×10−5 that of the CCSN rate within 0.3 ≤ z ≤ 1.4 by applying a Monte Carlo technique. The rate of slowly evolving, type Ic SLSNe (such as SN2007bi) is estimated as a factor of 10 lower than this range.

Selection of burst-like transients and stochastic variables using multi-brand 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 (SVs) 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 gP1, rP1, iP1, and zP1. We use three deterministic light-curve models to fit BL transients; a Gaussian, a Gamma distribution, and an analytic supernova (SN) model, and one stochastic light-curve model, the Ornstein-Uhlenbeck process, in order to fit variability that is characteristic of active galactic nuclei (AGNs). We assess the quality of fit of the models band-wise and 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, resulting in two measures of classification quality, from the averages across the photometric filters of (1) the classifications determined from the closest K-means cluster centers, and (2) the square distances from the clustering centers in the K-means clustering spaces. For a verification set of AGNs and SNe, we show that SV and BL occupy distinct regions in the plane constituted by these measures. We use our clustering method to characterize 4361 extragalactic image difference detected sources, in the first 2.5 yr of the PS1 MDS, into 1529 BL, and 2262 SV, with a purity of 95.00% for AGNs, 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 AGNs and 812 SNe. With these two samples, we characterize their variability and host galaxy properties, and identify simple photometric priors that would enable their real-time identification in future wide-field synoptic surveys.

Potentially harmful excipients in neonatal medicines: A pan-European observational study

Objectives We aimed to describe administration of eight potentially harmful excipients of interest (EOI)-parabens, polysorbate 80, propylene glycol, benzoates, saccharin sodium, sorbitol, ethanol and benzalkonium chloride-to hospitalised neonates in Europe and to identify risk factors for exposure. Methods All medicines administered to neonates during 1 day with individual prescription and demographic data were registered in a web-based point prevalence study. Excipients were identified from the Summaries of Product Characteristics. Determinants of EOI administration (geographical region, gestational age (GA), active pharmaceutical ingredient, unit level and hospital teaching status) were identified using multivariable logistical regression analysis. Results Overall 89 neonatal units from 21 countries participated. Altogether 2095 prescriptions for 530 products administered to 726 neonates were recorded. EOI were found in 638 (31%) prescriptions and were administered to 456 (63%) neonates through a relatively small number of products (n=142; 27%). Parabens, found in 71 (13%) products administered to 313 (43%) neonates, were used most frequently. EOI administration varied by geographical region, GA and route of administration. Geographical region remained a significant determinant of the use of parabens, polysorbate 80, propylene glycol and saccharin sodium after adjustment for the potential covariates including anatomical therapeutic chemical class of the active ingredient. Conclusions European neonates receive a number of potentially harmful pharmaceutical excipients. Regional differences in EOI administration suggest that EOI-free products are available and provide the potential for substitution to avoid side effects of some excipients.

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

We probe the systematic uncertainties from the 113 Type Ia supernovae (SN Ia) 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. 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 ± 0.031 mag for host galaxies with high and low masses. Assuming flatness and including systematic uncertainties in our analysis of only SNe measurements, we find w = -1.120^+0.360_-0.206(Stat)^+0.269_-0.291(Sys). With additional constraints from Baryon acoustic oscillation, cosmic microwave background (CMB) (Planck) and H₀ 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). The significance of the inconsistency with w = -1 depends on whether we use Planck or Wilkinson Microwave Anisotropy Probe measurements of the CMB: w_{BAO+H0+SN+WMAP} = -1.124^+0.083_-0.065.

Cosmological Constraints from Measurements of Type Ia Supernovae Discovered During the First 1.5 Yr of the Pan-STARRS1 Survey

We present griz_P1 light curves of 146 spectroscopically confirmed Type Ia supernovae (SNe Ia; 0.03 < z < 0.65) discovered during the first 1.5 yr 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 Hubble Space Telescope 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 supernovae (SNe) and assuming a constant dark energy equation of state and flatness, yields w = -1.120^+0.360_-0.206(Stat)^+0.269_0.291(Sys). When combined with BAO+CMB(Planck)+H₀, the analysis yields Ω_M = 0.280^+0.013_0.012 and w = -1.166^+0.072_-0.069 including all identified systematics. The value of w is inconsistent with the cosmological constant value of -1 at the 2.3σ level. Tension endures after removing either the baryon acoustic oscillation (BAO) or the H₀ constraint, though it is strongest when including the H₀ constraint. If we include WMAP9 cosmic microwave background (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 SN sample with ∼three times as many SNe should provide more conclusive results.

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 (SNe) models. The sample size of these objects has remained small due, at least in part, to the challenges 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 (t_1/2) 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_P1-rP1 ≲ -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⁴³erg s^-1), and lack of UV line blanketing observed in many of these transients with an explosion powered mainly by the radioactive decay of ⁵⁶Ni. 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 that ejected very little (<0.03 M) 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^-1Gpc^-3(4%-7% of the core-collapse SN rate at z = 0.2).

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 redshift z = 1.388. The light curve peaked at z P1 = 21.7 mag, making PS1-10afx comparable to the most luminous known SNe, with Mu = -22.3 mag. Our extensive optical and near-infrared observations indicate that the bolometric light curve of PS1-10afx rose on the unusually fast timescale of ~12 days to the extraordinary peak luminosity of 4.1 × 1044 erg s-1 (M bol = -22.8 mag) and subsequently faded rapidly. Equally important, the spectral energy distribution is unusually red for an SLSN, with a color temperature of ~6800 K near maximum light, in contrast to previous hydrogen-poor SLSNe, which are bright in the ultraviolet (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-1 and evidence for line blanketing in the rest-frame UV. Despite the fast rise, these parameters imply a very large emitting radius (gsim 5 × 1015 cm). We demonstrate that no existing theoretical model can satisfactorily explain this combination of properties: (1) a nickel-powered light curve cannot match the combination of high peak luminosity with the fast timescale; (2) models powered by the spindown energy of a rapidly rotating magnetar predict significantly hotter and faster ejecta; and (3) 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 ⊙ yr-1, and is fairly massive (~2 × 1010 M ⊙), with a stellar population age of ~108 yr, also in contrast to the young dwarf hosts of known hydrogen-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.

Alpha particles induce pan-nuclear phosphorylation of H2AX in primary human lymphocytes mediated through ATM

The use of high linear energy transfer radiations in the form of carbon ions in heavy ion beam lines or alpha particles in new radionuclide treatments has increased substantially over the past decade and will continue to do so due to the favourable dose distributions they can offer versus conventional therapies. Previously it has been shown that exposure to heavy ions induces pan-nuclear phosphorylation of several DNA repair proteins such as H2AX and ATM in vitro. Here we describe similar effects of alpha particles on ex vivo irradiated primary human peripheral blood lymphocytes. Following alpha particle irradiation pan-nuclear phosphorylation of H2AX and ATM, but not DNA-PK and 53BP1, was observed throughout the nucleus. Inhibition of ATM, but not DNA-PK, resulted in the loss of pan-nuclear phosphorylation of H2AX in alpha particle irradiated lymphocytes. Pan-nuclear gamma-H2AX signal was rapidly lost over 24h at a much greater rate than foci loss. Surprisingly, pan-nuclear gamma-H2AX intensity was not dependent on the number of alpha particle induced double strand breaks, rather the number of alpha particles which had traversed the cell nucleus. This distinct fluence dependent damage signature of particle radiation is important in both the fields of radioprotection and clinical oncology in determining radionuclide biological dosimetry and may be indicative of patient response to new radionuclide cancer therapies.

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 ∼ 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 ± 0.13: 0.71 ± 0.13, with larger objects more likely to have spin axes perpendicular to the orbital plane.

SN 2011dh: DISCOVERY OF A TYPE IIb SUPERNOVA FROM A COMPACT PROGENITOR IN THE NEARBY GALAXY M51

On 2011 May 31 UT a supernova (SN) exploded in the nearby galaxy M51 (the Whirlpool Galaxy). We discovered this event using small telescopes equipped with CCD cameras and also detected it with the Palomar Transient Factory survey, rapidly confirming it to be a Type II SN. Here, we present multi-color ultraviolet through infrared photometry which is used to calculate the bolometric luminosity and a series of spectra. Our early-time observations indicate that SN 2011dh resulted from the explosion of a relatively compact progenitor star. Rapid shock-breakout cooling leads to relatively low temperatures in early-time spectra, compared to explosions of red supergiant stars, as well as a rapid early light curve decline. Optical spectra of SN 2011dh are dominated by H lines out to day 10 after explosion, after which He I lines develop. This SN is likely a member of the cIIb (compact IIb) class, with progenitor radius larger than that of SN 2008ax and smaller than the eIIb (extended IIb) SN 1993J progenitor. Our data imply that the object identified in pre-explosion Hubble Space Telescope images at the SN location is possibly a companion to the progenitor or a blended source, and not the progenitor star itself, as its radius (~1013 cm) would be highly inconsistent with constraints from our post-explosion spectra.