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.

Massive stars exploding in a He-rich circumstellar medium - V. Observations of the slow-evolving SN Ibn OGLE-2012-SN-006

We present optical observations of the peculiar Type Ibn supernova (SN Ibn) OGLE-2012-SN-006, discovered and monitored by the Optical Gravitational Lensing Experiment-IV survey, and spectroscopically followed by Public ESO Spectroscopic Survey of Transient Objects (PESSTO) at late phases. Stringent pre-discovery limits constrain the explosion epoch with fair precision to JD = 245 6203.8 +/- 4.0. The rise time to the I-band light-curve maximum is about two weeks. The object reaches the peak absolute magnitude M-I = -19.65 +/- 0.19 on JD = 245 6218.1 +/- 1.8. After maximum, the light curve declines for about 25 d with a rate of 4 mag (100 d)(-1). The symmetric I-band peak resembles that of canonical Type Ib/c supernovae (SNe), whereas SNe Ibn usually exhibit asymmetric and narrower early-time light curves. Since 25 d past maximum, the light curve flattens with a decline rate slower than that of the Co-56-Fe-56 decay, although at very late phases it steepens to approach that rate. However, other observables suggest that the match with the Co-56 decay rate is a mere coincidence, and the radioactive decay is not the main mechanism powering the light curve of OGLE-2012-SN-006. An early-time spectrum is dominated by a blue continuum, with only a marginal evidence for the presence of He I lines marking this SN type. This spectrum shows broad absorptions bluewards than 5000 angstrom, likely O II lines, which are similar to spectral features observed in superluminous SNe at early epochs. The object has been spectroscopically monitored by PESSTO from 90 to 180 d after peak, and these spectra show the typical features observed in a number of SN 2006jc-like events, including a blue spectral energy distribution and prominent and narrow (v(FWHM) approximate to 1900 km s(-1)) He I emission lines. This suggests that the ejecta are interacting with He-rich circumstellar material. The detection of broad (10(4) km s(-1)) O I and Ca II features likely produced in the SN ejecta (including the [OI] lambda lambda 6300,6364 doublet in the latest spectra) lends support to the interpretation of OGLE-2012-SN-006 as a core-collapse event.

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.

OGLE-2013-SN-079: A Lonely Supernova Consistent with a Helium Shell Detonation

We present observational data for a peculiar supernova discovered by the OGLE-IV survey and followed by the Public ESO Spectroscopic Survey for Transient Objects. The inferred redshift of z = 0.07 implies an absolute magnitude in the rest-frame I-band of M-1 similar to -17.6 mag. This places it in the luminosity range between normal Type Ia SNe and novae. Optical and near infrared spectroscopy reveal mostly Ti and Ca lines, and an unusually red color arising from strong depression of flux at rest wavelengths

PESSTO: Survey description and products from the first data release by the Public ESO Spectroscopic Survey of Transient Objects

Context. The Public European Southern Observatory Spectroscopic Survey of Transient Objects (PESSTO) began as a public spectroscopic survey in April 2012. PESSTO classifies transients from publicly available sources and wide-field surveys, and selects science targets for detailed spectroscopic and photometric follow-up. PESSTO runs for nine months of the year, January - April and August - December inclusive, and typically has allocations of 10 nights per month. 

Aims. We describe the data reduction strategy and data products that are publicly available through the ESO archive as the Spectroscopic Survey data release 1 (SSDR1). 

Methods. PESSTO uses the New Technology Telescope with the instruments EFOSC2 and SOFI to provide optical and NIR spectroscopy and imaging. We target supernovae and optical transients brighter than 20.5^m for classification. Science targets are selected for follow-up based on the PESSTO science goal of extending knowledge of the extremes of the supernova population. We use standard EFOSC2 set-ups providing spectra with resolutions of 13-18 Å between 3345-9995 Å. A subset of the brighter science targets are selected for SOFI spectroscopy with the blue and red grisms (0.935-2.53 μm and resolutions 23-33 Å) and imaging with broadband JHK<inf>s</inf> filters. 

Results. This first data release (SSDR1) contains flux calibrated spectra from the first year (April 2012-2013). A total of 221 confirmed supernovae were classified, and we released calibrated optical spectra and classifications publicly within 24 h of the data being taken (via WISeREP). The data in SSDR1 replace those released spectra. They have more reliable and quantifiable flux calibrations, correction for telluric absorption, and are made available in standard ESO Phase 3 formats. We estimate the absolute accuracy of the flux calibrations for EFOSC2 across the whole survey in SSDR1 to be typically ∼15%, although a number of spectra will have less reliable absolute flux calibration because of weather and slit losses. Acquisition images for each spectrum are available which, in principle, can allow the user to refine the absolute flux calibration. The standard NIR reduction process does not produce high accuracy absolute spectrophotometry but synthetic photometry with accompanying JHK<inf>s</inf> imaging can improve this. Whenever possible, reduced SOFI images are provided to allow this. 

Conclusions. Future data releases will focus on improving the automated flux calibration of the data products. The rapid turnaround between discovery and classification and access to reliable pipeline processed data products has allowed early science papers in the first few months of the survey.

PESSTO monitoring of SN 2012hn: further heterogeneity among faint Type I supernovae

We present optical and infrared monitoring data of SN 2012hn collectedby the Public European Southern Observatory Spectroscopic Survey forTransient Objects. We show that SN 2012hn has a faint peak magnitude(MR ˜ -15.65) and shows no hydrogen and no clearevidence for helium in its spectral evolution. Instead, we detectprominent Ca II lines at all epochs, which relates this transient topreviously described `Ca-rich' or `gap' transients. However, thephotospheric spectra (from -3 to +32 d with respect to peak) of SN2012hn show a series of absorption lines which are unique and a redcontinuum that is likely intrinsic rather than due to extinction. Linesof Ti II and Cr II are visible. This may be a temperature effect, whichcould also explain the red photospheric colour. A nebular spectrum at+150 d shows prominent Ca II, O I, C I and possibly Mg I lines whichappear similar in strength to those displayed by core-collapsesupernovae (SNe). To add to the puzzle, SN 2012hn is located at aprojected distance of 6 kpc from an E/S0 host and is not close to anyobvious star-forming region. Overall SN 2012hn resembles a group offaint H-poor SNe that have been discovered recently and for which aconvincing and consistent physical explanation is still missing. Theyall appear to explode preferentially in remote locations offset from amassive host galaxy with deep limits on any dwarf host galaxies,favouring old progenitor systems. SN 2012hn adds heterogeneity to thissample of objects. We discuss potential explosion channels includingHe-shell detonations and double detonations of white dwarfs as well aspeculiar core-collapse SNe.

Interpreting signals from astrophysical transient experiments

Time domain astronomy has come of age with astronomers now able to monitor the sky at high cadence both across the electromagnetic spectrum and using neutrinos and gravitational waves. The advent of new observing facilities permits new science, but the ever increasing throughput of facilities demands efficient communication of coincident detections and better subsequent coordination among the scientific community so as to turn detections into scientific discoveries. To discuss the revolution occurring in our ability to monitor the Universe and the challenges it brings, on 2012 April 25-26 a group of scientists from observational and theoretical teams studying transients met with representatives of the major international transient observing facilities at the Kavli Royal Society International Centre, UK. This immediately followed the Royal Society Discussion meeting "New windows on transients across the Universe" held in London. Here we present a summary of the Kavli meeting at which the participants discussed the science goals common to the transient astronomy community and analysed how to better meet the challenges ahead as ever more powerful observational facilities come on stream.

Synchrophasor-Based Islanding Detection for Distributed Generation Systems Using Systematic Principal Component Analysis Approaches

Systematic principal component analysis (PCA) methods are presented in this paper for reliable islanding detection for power systems with significant penetration of distributed generations (DGs), where synchrophasors recorded by Phasor Measurement Units (PMUs) are used for system monitoring. Existing islanding detection methods such as Rate-of-change-of frequency (ROCOF) and Vector Shift are fast for processing local information, however with the growth in installed capacity of DGs, they suffer from several drawbacks. Incumbent genset islanding detection cannot distinguish a system wide disturbance from an islanding event, leading to mal-operation. The problem is even more significant when the grid does not have sufficient inertia to limit frequency divergences in the system fault/stress due to the high penetration of DGs. To tackle such problems, this paper introduces PCA methods for islanding detection. Simple control chart is established for intuitive visualization of the transients. A Recursive PCA (RPCA) scheme is proposed as a reliable extension of the PCA method to reduce the false alarms for time-varying process. To further reduce the computational burden, the approximate linear dependence condition (ALDC) errors are calculated to update the associated PCA model. The proposed PCA and RPCA methods are verified by detecting abnormal transients occurring in the UK utility network.

Locations of peculiar supernovae as a diagnostic of their origins

We put constraints on the properties of the progenitors of peculiar calcium-rich transients using the distribution of locations within their host galaxies. We confirm that this class of transients do not follow the galaxy stellar mass profile and are more likely to be found in remote locations of their apparent hosts. We test the hypothesis that these transients are from low-metallicity progenitors by comparing their spatial distributions with the predictions of self-consistent cosmological simulations that include star formation and chemical enrichment. We find that while metal-poor stars and our transient sample show a consistent preference for large offsets, metallicity alone cannot explain the extreme cases. Invoking a lower age limit on the progenitor helps to improve the match, indicating these events may result from a very old metal-poor population. We also investigate the radial distribution of globular cluster systems, and show that they too are consistent with the class of calcium-rich transients. Because photometric upper limits exist for globular clusters for some members of the class, a production mechanism related to the dense environment of globular clusters is not favoured for the calcium-rich events. However, the methods developed in this paper may be used in the future to constrain the effects of low metallicity on radially distant core-collapse events or help establish a correlation with globular clusters for other classes of peculiar explosions.

Human odontoblasts express functional TRPV2 channels

Background: The transient receptor potential (TRP) ion channels play a critical role in sensory physiology, where they act as transducers of thermal, mechanical and chemical stimuli. We have previously shown the functional expression of several TRP channels by human odontoblast-like cells and proposed their significance in odontoblast sensory perception. Functional expression of the mechano-sensitiveTRPV2 channel by human odontoblasts would further support a role for TRP channels in odontoblast physiology. Objective: The objective of the current study was to determine the functional expression of TRPV2 by human odontoblasts. Methods: Human dental pulp cells were cultured in the presence of 2 mM β-glycerophoshate to induce an odontoblast phenotype. TRPV2 gene expression was determined by qPCR employing custom designed FAM TRPV2 specific primers and probes (Roche, UK) and the Light Cycler 480 Probes Master (Roche). TRPV2 protein expression was determined following SDS-PAGE and Western blotting of cell lysate preparations. Functional expression of TRPV2 was investigated by Ca2+ microfluorimetry. Results: qPCR data indicated robust expression of TRPV2 in odontoblast-like cells. Western blotting revealed a discrete immunoreactive protein band indicating expression of TRPV2 in cell lysates. In functional assays, the chemical agonist of TRPV2, cannabidiol, was shown to elicit [Ca2+]i transients, that were reduced to baseline in the presence of the TRPV2 antagonist Tranilast, suggesting channel functionality in odontoblast-like cells. Conclusion: These results provide the first evidence for the functional expression of TRPV2 in human odontoblast-like cells, providing further support for the role of TRP channels in odontoblast physiology.

SN 2011A: A low-luminosity interacting transient with a double plateau and strong sodium absorpton

We present optical photometry and spectroscopy of the optical transient SN 2011A. Our data span 140 days after discovery including BVRI u′g′r′i′z′ photometry and 11 epochs of optical spectroscopy. Originally classified as a type IIn supernova (SN IIn) due to the presence of narrow Hα emission, this object shows exceptional characteristics. First, the light curve shows a double plateau, a property only observed before in the impostor SN 1997bs. Second, SN 2011A has a very low luminosity (M_V=-15.72), placing it between normal luminous SNe IIn and SN impostors. Third, SN 2011A shows low velocity and high equivalent width absorption close to the sodium doublet, which increases with time and is most likely of circumstellar origin. This evolution is also accompanied by a change in line profile; when the absorption becomes stronger, a P Cygni profile appears. We discuss SN 2011A in the context of interacting SNe IIn and SN impostors, which appears to confirm the uniqueness of this transient. While we favor an impostor origin for SN 2011A, we highlight the difficulty in differentiating between terminal and non-terminal interacting transients.

Spectral analysis of four 'hypervariable' AGN: a micro-needle in the haystack?

We analyze four extreme AGN transients to explore the possibility that they are caused by rare, high-amplitude microlensing events. These previously unknown type-I AGN are located in the redshift range 0.6-1.1 and show changes of > 1.5 magnitudes in the g-band on a timescale of ~years. Multi-epoch optical spectroscopy, from the William Herschel Telescope, shows clear differential variability in the broad line fluxes with respect to the continuum changes and also evolution in the line profiles. In two cases a simple point-source, point-lens microlensing model provides an excellent match to the long-term variability seen in these objects. For both models the parameter constraints are consistent with the microlensing being due to an intervening stellar mass object but as yet there is no confirmation of the presence of an intervening galaxy. The models predict a peak amplification of 10.3/13.5 and an Einstein timescale of 7.5/10.8 years respectively. In one case the data also allow constraints on the size of the CIII] emitting region, with some simplifying assumptions, to to be ~1.0-6.5 light-days and a lower limit on the size of the MgII emitting region to be > 9 light-days (half-light radii). This CIII] radius is perhaps surprisingly small. In the remaining two objects there is spectroscopic evidence for an intervening absorber but the extra structure seen in the lightcurves requires a more complex lensing scenario to adequately explain.