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.

Energy levels, radiative rates and electron impact excitation rates for transitions in AlX

Energy levels, radiative rates and lifetimes are calculated among the lowest 98 levels of the n ≤4 configurations of Be-like Al X. The GRASP (General-purpose Relativistic Atomic Structure Package) is adopted and data are provided for all E1, E2, M1 and M2 transitions. Similar data are also obtained with the FAC (Flexible Atomic Code) to assess the accuracy of the calculations. Based on comparisons between calculations with the two codes as well as with available measurements, our listed energy levels are assessed to be accurate to better than 0.3 per cent. However, the accuracy for radiative rates and lifetimes is estimated to be about 20 per cent. Collision strengths are also calculated for which the DARC (Dirac Atomic R-matrix Code) is used. A wide energy range (up to 380 Ryd) is considered and resonances resolved in a fine energy mesh in the thresholds region. The collision strengths are subsequently averaged over a Maxwellian velocity distribution to determine effective collision strengths up to a temperature of 1.6 × 10⁷ K. Our results are compared with the previous (limited) atomic data and significant differences (up to a factor of 4) are noted for several transitions, particularly those which are not allowed in jj coupling. 

Early-type stars observed in the ESO UVES Paranal observatory project - IV. studies of CN, CH+ and CH in the interstellar medium

High spectral resolution (~80 000) and signal-to-noise observations from the Ultraviolet and Visual Echelle Spectrograph Paranal Observatory Project (UVES-POP) are used to study the interstellarmolecular lines CN (3874 Å), CH⁺ (3957, 4232 Å) and CH (3886, 4300 Å) towards 74 O- and B-type stellar sightlines. Additionally, archive data are presented for 140 ELODIE early-type stellar sightlines at R = 42 000, plus 25 FEROS at R = 48 000 and 3 UVES at R > 50 000, mainly in the CH⁺ (4232 Å) and CH (3886, 4300 Å) transitions. Detection rates are ~45 per cent for CN and ~67 per cent for the other lines in the POP sample, and ~10-15 per cent for CH⁺ and CH lines in the additional sample. CH and CH⁺ are well correlated between log[N(CH) cm^-2]~12-14, implying that these clouds are CH⁺-like CH and not CN-like CH. CH is also very well correlated with Na I D in the range log[N(Na I cm^-2]) ~12.2-14.2. A few sightlines show tentative velocity shifts of ~2 km s^-1 between CH and CH⁺, which appear to be caused by differences in component strength in blends, and hence do not provide firm evidence for shocks. Finally, we describe a search for ¹³CH⁺ in a sightline towards HD 76341. No ¹³CH⁺ is detected, placing a limit on the ¹³CH⁺ to ¹²CH⁺ ratio of ~0.01. If a formal fit is attempted, the equivalent width ratio in the two isotopes is a factor ~90 but with large errors. 

Energy levels, radiative rates and electron impact excitation rates for transitions in Si II

Energies for the lowest 56 levels, belonging to the 3s² 3p, 3s 3p², 3p³, 3s² 3d, 3s 3p 3d, 3s² 4ℓ and 3s² 5ℓ configurations of Si II, are calculated using the General-purpose Relativistic Atomic Structure Package (GRASP) code. Analogous calculations have also been performed (for up to 175 levels) using the FlexibleAtomicCode (FAC). Furthermore, radiative rates are calculated for all E1, E2, M1 and M2 transitions. Extensive comparisons are made with available theoretical and experimental energy levels, and the accuracy of the present results is assessed to be better than 0.1Ryd. Similarly, the accuracy for radiative rates (and subsequently lifetimes) is estimated to be better than 20 per cent for most of the (strong) transitions. Electron impact excitation collision strengths are also calculated, with the Dirac Atomic R-matrix Code (DARC), over a wide energy range up to 13 Ryd. Finally, to determine effective collision strengths, resonances are resolved in a fine energy mesh in the thresholds region. These collision strengths are averaged over a Maxwellian velocity distribution and results listed over a wide range of temperatures, up to 10^5.5 K. Our data are compared with earlier R-matrix calculations and differences noted, up to a factor of 2, for several transitions. Although scope remains for improvement, the accuracy for our results of collision strengths and effective collision strengths is assessed to be about 20 per cent for a majority of transitions. 

Transit timing variations in WASP-10b induced by stellar activity

The hot-JupiterWASP-10bwas reported by Maciejewski et al. to showtransit timing variations (TTVs) with an amplitude of ~3.5 min. These authors proposed that the observed TTVs were caused by a 0.1M_Jup perturbing companion with an orbital period of ~5.23 d, and hence, close to the outer 5:3 mean-motion resonance with WASP-10b. To test this scenario, we present eight new transit light curves of WASP-10b obtained with the Faulkes Telescope North and the Liverpool Telescope. The new light curves, together with 22 previously published ones, were modelled with a Markov Chain Monte Carlo transit fitting code. Transit depth differences reported forWASP-10b are thought to be due to starspot-induced brightness modulation of the host star. Assuming the star is brighter at the activity minimum, we favour a small planetary radius. We find Rp = 1.039^+0.043 _-0.049R_Jup in agreement with Johnson et al. and Maciejewski et al. Recent studies find no evidence for a significant eccentricity in this system. We present consistent system parameters for a circular orbit and refine the orbital ephemeris ofWASP-10b. Our homogeneously derived transit times do not support the previous claimed TTV signal, which was strongly dependent on two previously published transits that have been incorrectly normalized. Nevertheless, a linear ephemeris is not a statistically good fit to the transit times of WASP-10b. We show that the observed transit time variations are due to spot occultation features or systematics. We discuss and exemplify the effects of occultation spot features in the measured transit times and show that despite spot occultation during egress and ingress being difficult to distinguish in the transit light curves, they have a significant effect in the measured transit times. We conclude that if we account for spot features, the transit times of WASP-10b are consistent with a linear ephemeris with the exception of one transit (epoch 143) which is a partial transit. Therefore, there is currently no evidence for the existence of a companion to WASP-10b. Our results support the lack of TTVs of hot-Jupiters reported for the Kepler sample. 

Expanded iron uta spectra - Probing the thermal stability limits in AGN clouds

The Fe unresolved transition arrays (UTAs) produce prominent features in the 15-17 Å wavelength range in the spectra of active galactic nuclei (AGNs). Here, we present new calculations of the energies and oscillator strengths of inner-shell lines from Fe XIV, Fe XV, and Fe XVI. These are crucial ions since they are dominant at inflection points in the gas thermal stability curve, and UTA excitation followed by autoionization is an important ionization mechanism for these species. We incorporate these, and data reported in previous papers, into the plasma simulation code Cloudy. This updated physics is subsequently employed to reconsider the thermally stable phases in absorbing media in AGNs. We show how the absorption profile of the Fe XIV UTA depends on density, due to the changing populations of levels within the ground configuration. © 2013. The American Astronomical Society. All rights reserved.

Early-type stars observed in the ESO UVES Paranal Observatory Project - V. Time-variable interstellar absorption

The structure and properties of the diffuse interstellar medium (ISM) on small scales, sub-au to 1 pc, are poorly understood. We compare interstellar absorption-lines, observed towards a selection of O- and B-type stars at two or more epochs, to search for variations over time caused by the transverse motion of each star combined with changes in the structure in the foreground ISM. Two sets of data were used: 83 VLT- UVES spectra with approximately 6 yr between epochs and 21 McDonald observatory 2.7m telescope echelle spectra with 6 - 20 yr between epochs, over a range of scales from 0 - 360 au. The interstellar absorption-lines observed at the two epochs were subtracted and searched for any residuals due to changes in the foreground ISM. Of the 104 sightlines investigated with typically five or more components in Na I D, possible temporal variation was identified in five UVES spectra (six components), in Ca II, Ca I and/or Na I absorption-lines. The variations detected range from 7\% to a factor of 3.6 in column density. No variation was found in any other interstellar species. Most sightlines show no variation, with 3{\sigma} upper limits to changes of the order 0.1 - 0.3 dex in Ca II and Na I. These variations observed imply that fine-scale structure is present in the ISM, but at the resolution available in this study, is not very common at visible wavelengths. A determination of the electron densities and lower limits to the total number density of a sample of the sightlines implies that there is no striking difference between these parameters in sightlines with, and sightlines without, varying components.

Energy levels, radiative rates and electron impact excitation rates for transitions in C III

We report energy levels, radiative rates (A-values) and lifetimes for the astrophysically important Be-like ion C III. For the calculations, 166 levels belonging to the n ≤ 5 configurations are considered and the GRASP (General-purpose Relativistic Atomic Structure Package) is adopted. Einstein A-coefficients are provided for all E1, E2, M1 and M2 transitions, while lifetimes are compared with available measurements as well as theoretical results, and no large discrepancies noted. Our energy levels are assessed to be accurate to better than 1 per cent for a majority of levels, and A-values to better than 20 per cent for most transitions. Collision strengths are also calculated, for which the Dirac Atomic R-matrix Code (DARC) is used. A wide energy range, up to 21 Ryd, is considered and resonances resolved in a fine energy mesh in the thresholds region. The collision strengths are subsequently averaged over a Maxwellian velocity distribution to determine effective collision strengths up to a temperature of 8.0 ×10[5]K, sufficient for most astrophysical applications. Our data are compared with the recent R-matrix calculations of Fernández-Menchero et al., and significant differences (up to over an order of magnitude) are noted for several transitions over the complete temperature range of the results.

Electron impact excitation of Be-like ions: a comparison of DARC and ICFT results

Emission lines of Be-like ions are frequently observed in astrophysical plasmas, and many are useful for density and temperature diagnostics. However, accurate atomic data for energy levels, radiative rates (A-values) and effective electron excitation collision strengths ($\Upsilon$) are required for reliable plasma modelling. In general it is reasonably straightforward to calculate energy levels and A- values to a high level of accuracy. By contrast, considerable effort is required to calculate $\Upsilon$, and hence it is not always possible to assess the accuracy of available data. Recently, two independent calculations (adopting the $R$-matrix method) but with different approaches (DARC and ICFT) have appeared for a range of Be-like ions. Therefore, in this work we compare the two sets of $\Upsilon$, highlight the large discrepancies for a significant number of transitions and suggest possible reasons for these.

Energy levels, radiative rates and electron impact excitation rates for transitions in Fe XIV

Energies and lifetimes are reported for the lowest 136 levels of Fe XIV, belonging to the (1s(2)2s(2)2p(6)) 3s(2)3p, 3s(3)p(2), 3s(2)3d, 3p(3), 3s(3)p(3)d, 3p(2)3d, 3s(3)d(2), 3p(3)d(2) and 3s(2)4l configurations. Additionally, radiative rates, oscillator strengths and line strengths are calculated for all electric dipole (E1), magnetic dipole (M1), electric quadrupole (E2) and magnetic quadrupole (M2) transitions. Theoretical lifetimes determined from these radiative rates for most levels show satisfactory agreement with earlier calculations, a swell as with measurements. Electron impact excitation collision strengths are also calculated with the Dirac atomic R-matrix code (DARC) over a wide energy range up to 260 Ryd. Furthermore, resonances have been resolved in a fine energy mesh to determine effective collision strengths, obtained after integrating the collision strengths over a Maxwellian distribution of electron velocities. Results are listed for all 9180 transitions among the 136 levels over a wide range of electron temperatures, up to 10(7.1) K. Comparisons are made with available results in the literature, and the accuracy of the present data is assessed.

Large- and small-scale structure of the intermediate- and high-velocity clouds towards the LMC and SMC.

We employ Ca II K and Na I D interstellar absorption-line spectroscopy of early-type stars in the Large and Small Magellanic Clouds (LMC, SMC) to investigate the large- and small-scale structure in foreground intermediate- and high-velocity clouds (I/HVCs). Data include FLAMES-GIRAFFE Ca II K observations of 403 stars in four open clusters, plus FEROS or UVES spectra of 156 stars in the LMC and SMC. The FLAMES observations are amongst the most extensive probes to date of Ca II structures on ∼20 arcsec scales in Magellanic I/HVCs. From the FLAMES data within a 0 ∘.∘.∘.5 field of view, the Ca II K equivalent width in the I/HVC components towards three clusters varies by factors of ≥10. There are no detections of molecular gas in absorption at intermediate or high velocities, although molecular absorption is present at LMC and Galactic velocities towards some sightlines. The FEROS/UVES data show Ca II K I/HVC absorption in ∼60 per cent of sightlines. The range in the Ca II/Na I ratio in I/HVCs is from –0.45 to +1.5 dex, similar to previous measurements for I/HVCs. In 10 sightlines we find Ca II/O I ratios in I/HVC gas ranging from 0.2 to 1.5 dex below the solar value, indicating either dust or ionization effects. In nine sightlines I/HVC gas is detected in both H I and Ca II at similar velocities, implying that the two elements form part of the same structure.

Deflagrations in hybrid CONe white dwarfs: a route to explain the faint Type Iax supernova 2008ha

Stellar evolution models predict the existence of hybrid white dwarfs (WDs) with a carbon-oxygen core surrounded by an oxygen-neon mantle. Being born with masses similar to 1.1 M-aS (TM), hybrid WDs in a binary system may easily approach the Chandrasekhar mass (M-Ch) by accretion and give rise to a thermonuclear explosion. Here, we investigate an off-centre deflagration in a near-M-Ch hybrid WD under the assumption that nuclear burning only occurs in carbon-rich material. Performing hydrodynamics simulations of the explosion and detailed nucleosynthesis post-processing calculations, we find that only 0.014 M-aS (TM) of material is ejected while the remainder of the mass stays bound. The ejecta consist predominantly of iron-group elements, O, C, Si and S. We also calculate synthetic observables for our model and find reasonable agreement with the faint Type Iax SN 2008ha. This shows for the first time that deflagrations in near-M-Ch WDs can in principle explain the observed diversity of Type Iax supernovae. Leaving behind a near-M-Ch bound remnant opens the possibility for recurrent explosions or a subsequent accretion-induced collapse in faint Type Iax SNe, if further accretion episodes occur. From binary population synthesis calculations, we find the rate of hybrid WDs approaching M-Ch to be of the order of 1 per cent of the Galactic SN Ia rate.

The impact of accretion disc winds on the optical spectra of cataclysmic variables

Many high-state non-magnetic cataclysmic variables (CVs) exhibit blueshifted absorption or P-Cygni profiles associated with ultraviolet (UV) resonance lines. These features imply the existence of powerful accretion disc winds in CVs. Here, we use our Monte Carlo ionization and radiative transfer code to investigate whether disc wind models that produce realistic UV line profiles are also likely to generate observationally significant recombination line and continuum emission in the optical waveband. We also test whether outflows may be responsible for the single-peaked emission line profiles often seen in high-state CVs and for the weakness of the Balmer absorption edge (relative to simple models of optically thick accretion discs). We find that a standard disc wind model that is successful in reproducing the UV spectra of CVs also leaves a noticeable imprint on the optical spectrum, particularly for systems viewed at high inclination. The strongest optical wind-formed recombination lines are H alpha and He ii lambda 4686. We demonstrate that a higher density outflow model produces all the expected H and He lines and produces a recombination continuum that can fill in the Balmer jump at high inclinations. This model displays reasonable verisimilitude with the optical spectrum of RW Trianguli. No single-peaked emission is seen, although we observe a narrowing of the double-peaked emission lines from the base of the wind. Finally, we show that even denser models can produce a single-peaked H alpha line. On the basis of our results, we suggest that winds can modify, and perhaps even dominate, the line and continuum emission from CVs.

Polarization spectral synthesis for Type Ia supernova explosion models

We present a Monte Carlo radiative transfer technique for calculating synthetic spectropolarimetry for multidimensional supernova explosion models. The approach utilizes 'virtual-packets' that are generated during the propagation of the Monte Carlo quanta and used to compute synthetic observables for specific observer orientations. Compared to extracting synthetic observables by direct binning of emergent Monte Carlo quanta, this virtual-packet approach leads to a substantial reduction in the Monte Carlo noise. This is not only vital for calculating synthetic spectropolarimetry (since the degree of polarization is typically very small) but also useful for calculations of light curves and spectra. We first validate our approach via application of an idealized test code to simple geometries. We then describe its implementation in the Monte Carlo radiative transfer code ARTIS and present test calculations for simple models for Type Ia supernovae. Specifically, we use the well-known one-dimensional W7 model to verify that our scheme can accurately recover zero polarization from a spherical model, and to demonstrate the reduction in Monte Carlo noise compared to a simple packet-binning approach. To investigate the impact of aspherical ejecta on the polarization spectra, we then use ARTIS to calculate synthetic observables for prolate and oblate ellipsoidal models with Type Ia supernova compositions.

Optical and near-infrared observations of SN 2011dh - The first 100 days

We present optical and near-infrared (NIR) photometry and spectroscopy of the Type IIb supernova (SN) 2011dh for the first 100 days. We complement our extensive dataset with Swift ultra-violet (UV) and Spitzer mid-infrared (MIR) data to build a UV to MIR bolometric lightcurve using both photometric and spectroscopic data. Hydrodynamical modelling of the SN based on this bolometric lightcurve have been presented in Bersten et al. (2012, ApJ, 757, 31). We find that the absorption minimum for the hydrogen lines is never seen below ~11 000 km s-1 but approaches this value as the lines get weaker. This suggests that the interface between the helium core and hydrogen rich envelope is located near this velocity in agreement with the Bersten et al. (2012) He4R270 ejecta model. Spectral modelling of the hydrogen lines using this ejecta model supports the conclusion and we find a hydrogen mass of 0.01-0.04 M⊙ to be consistent with the observed spectral evolution. We estimate that the photosphere reaches the helium core at 5-7 days whereas the helium lines appear between ~10 and ~15 days, close to the photosphere and then move outward in velocity until ~40 days. This suggests that increasing non-thermal excitation due to decreasing optical depth for the γ-rays is driving the early evolution of these lines. The Spitzer 4.5 μm band shows a significant flux excess, which we attribute to CO fundamental band emission or a thermal dust echo although further work using late time data is needed. Thedistance and in particular the extinction, where we use spectral modelling to put further constraints, is discussed in some detail as well as the sensitivity of the hydrodynamical modelling to errors in these quantities. We also provide and discuss pre- and post-explosion observations of the SN site which shows a reduction by ~75 percent in flux at the position of the yellow supergiant coincident with SN 2011dh. The B, V and r band decline rates of 0.0073, 0.0090 and 0.0053 mag day-1 respectively are consistent with the remaining flux being emitted by the SN. Hence we find that the star was indeed the progenitor of SN 2011dh as previously suggested by Maund et al. (2011, ApJ, 739, L37) and which is also consistent with the results from the hydrodynamical modelling. Figures 2, 3, Tables 3-10, and Appendices are available in electronic form at http://www.aanda.orgThe photometric tables are only available at the CDS via anonymous ftp to http://cdsarc.u-strasbg.fr (ftp://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/562/A17