A novel method for transient detection in high-cadence optical surveys: Its application for a systematic search for novae in M31

[abridged] In large-scale time-domain surveys, the processing of data, from procurement up to the detection of sources, is generally automated. One of the main challenges is contamination by artifacts, especially in regions of strong unresolved emission. We present a novel method for identifying candidates for variables and transients from the outputs of such surveys' data pipelines. We use the method to systematically search for novae in iPTF observations of the bulge of M31. We demonstrate that most artifacts produced by the iPTF pipeline form a locally uniform background of false detections approximately obeying Poissonian statistics, whereas genuine variables and transients as well as artifacts associated with bright stars result in clusters of detections, whose spread is determined by the source localization accuracy. This makes the problem analogous to source detection on images produced by X-ray telescopes, enabling one to utilize tools developed in X-ray astronomy. In particular, we use a wavelet-based source detection algorithm from the Chandra data analysis package CIAO. Starting from ~2.5x10^5 raw detections made by the iPTF data pipeline, we obtain ~4000 unique source candidates. Cross-matching these candidates with the source-catalog of a deep reference image, we find counterparts for ~90% of them. These are either artifacts due to imperfect PSF matching or genuine variable sources. The remaining ~400 detections are transient sources. We identify novae among these candidates by applying selection cuts based on the expected properties of nova lightcurves. Thus, we recovered all 12 known novae registered during the time span of the survey and discovered three nova candidates. Our method is generic and can be applied for mining any target out of the artifacts in optical time-domain data. As it is fully automated, its incompleteness can be accurately computed and corrected for.Comment: 16 pages, 8 figures, accepted to A&

Variability of Red Supergiants in M31 from the Palomar Transient Factory

Most massive stars end their lives as Red Supergiants (RSGs), a short-lived evolution phase when they are known to pulsate with varying amplitudes. The RSG period-luminosity (PL) relation has been measured in the Milky Way, the Magellanic Clouds and M33 for about 120 stars in total. Using over 1500 epochs of R-band monitoring from the Palomar Transient Factory (PTF) survey over a five-year period, we study the variability of 255 spectroscopically cataloged RSGs in M31. We find that all RGSs brighter than M_K~ -10 mag (log(L/L_sun)>4.8) are variable at dm_R>0.05 mag. Our period analysis finds 63 with significant pulsation periods. Using the periods found and the known values of M_K for these stars, we derive the RSG PL relation in M31 and show that it is consistent with those derived earlier in other galaxies of different metallicities. We also detect, for the first time, a sequence of likely first-overtone pulsations. Comparison to stellar evolution models from MESA confirms the first overtone hypothesis and indicates that the variable stars in this sample have 12 M_sun<M<24 M_sun. As these RSGs are the immediate progenitors to Type II-P core-collapse supernovae (SNe), we also explore the implication of their variability in the initial-mass estimates for SN progenitors based on archival images of the progenitors. We find that this effect is small compared to the present measurement errors.Comment: 17 pages, 10 figure

Multiwavelength approach to classifying transient events in the direction of M31

Context. In the hunt for rare time-domain events, it is important to consider confusing exotic extragalactic phenomena with more common Galactic foreground events. Aims. We show how observations from multiple wavebands, in this case optical and X-ray observations, can be used to facilitate the distinction between the two. Methods. We discovered an extremely bright and rapid transient event during optical observations of the M 31 galaxy taken by the intermediate Palomar Transient Factory (iPTF). The persistent optical counterpart of this transient was previously thought to be a variable star in M 31 without any dramatic flux excursions. The iPTF event initially appeared to be an extraordinarily rapid and energetic extragalactic transient, which had a ≈3 mag positive flux excursion in less than a kilosecond; one of the exciting possibilities was this event could be a very fast nova in M 31. The nature of the source was resolved with the help of Chandraarchival data, where we found an X-ray counterpart and obtained its X-ray spectrum. Results. We find the X-ray spectrum of the quiescent emission can be described by a model of optically thin plasma emission with a temperature of ≈7 MK, typical for coronal emission from an active star. The combination of the X-ray luminosity, which is calculated assuming the source is located in M 31 (~3 × 10^(36) erg s^(−1)), and the color temperature exclude any type of known accreting compact object or active star in M 31. We argue instead that the optical transient source is an M-type main-sequence, active star located in the disk of the Milky Way at a distance of ~0.5–1 kpc. Its persistent X-ray luminosity is in the ≈1.3–5 × 10^(30)erg s^(−1) range and it has the absolute optical magnitude of 9.5–11.0 mag in the R band. The observed optical flare has the equivalent duration of ≈95 min and total energy of ≈(0.3–1) × 10^(35) erg in the R band, which places it among the brightest flares ever observed from an M-type star. This case can serve as an example for the classification of Galactic and extragalactic events in upcoming high-cadence time-domain projects, such as the Zwicky Transient Facility and the Large Synoptic Survey Telescope

Variability of massive stars in M31 from the Palomar Transient Factory

Using data from the (intermediate) Palomar Transient Factory (iPTF), we characterize the time variability of ≈500 massive stars in M31. Our sample is those stars that are spectrally typed by Massey and collaborators, including Luminous Blue Variables, Wolf–Rayets, and warm and cool supergiants. We use the high-cadence, long-baseline (≈5 yr) data from the iPTF survey, coupled with data-processing tools that model complex features in the light curves. We find widespread photometric (R-band) variability in the upper Hertzsprung Russell diagram (or CMD) with an increasing prevalence of variability with later spectral types. Red stars (V − I > 1.5) exhibit larger amplitude fluctuations than their bluer counterparts. We extract a characteristic variability timescale, t_(ch), via wavelet transformations that are sensitive to both continuous and localized fluctuations. Cool supergiants are characterized by longer timescales (>100 days) than the hotter stars. The latter have typical timescales of tens of days but cover a wider range, from our resolution limit of a few days to longer than 100 days. Using a 60 night block of data straddling two nights with a cadence of around 2 minutes, we extracted t_(ch) in the range 0.1–10 days with amplitudes of a few percent for 13 stars. Though there is broad agreement between the observed variability characteristics in the different parts of the upper CMD with theoretical predictions, detailed comparison requires models with a more comprehensive treatment of the various physical processes operating in these stars, such as pulsation, subsurface convection, and the effect of binary companions

The SN 2023ixf Progenitor in M101: II. Properties

We follow our first paper with an analysis of the ensemble of the extensive pre-explosion ground- and space-based infrared observations of the red supergiant (RSG) progenitor candidate for the nearby core-collapse supernova SN 2023ixf in Messier 101, together with optical data prior to explosion obtained with the Hubble Space Telescope (HST). We have confirmed the association of the progenitor candidate with the SN, as well as constrained the metallicity at the SN site, based on SN observations with instruments at Gemini-North. The internal host extinction to the SN has also been confirmed from a high-resolution Keck spectrum. We fit the observed spectral energy distribution (SED) for the star, accounting for its intrinsic variability, with dust radiative-transfer modeling, which assume a silicate-rich dust shell ahead of the underlying stellar photosphere. The star is heavily dust-obscured, likely the dustiest progenitor candidate yet encountered. We found maximum-likelihood estimates of the star's effective temperature and luminosity of 3450 K and 9.3e4 L_Sun, with 68% credible intervals of 2370--3700 K and (7.6--10.8)e4 L_sun. The candidate may have a Galactic RSG analog, IRC -10414, with a strikingly similar SED and luminosity. Via comparison with single-star evolutionary models we have constrained the initial mass of the progenitor candidate from 12 M_sun to as high as 15 M_sun. We have had available to us an extraordinary view of the SN 2023ixf progenitor candidate, which should be further followed up in future years with HST and the James Webb Space Telescope.Comment: 29 pages, submitted to AAS Journal

The ANTARES Astronomical Time-Domain Event Broker

We describe the Arizona-NOIRLab Temporal Analysis and Response to Events System (ANTARES), a software instrument designed to process large-scale streams of astronomical time-domain alerts. With the advent of large-format CCDs on wide-field imaging telescopes, time-domain surveys now routinely discover tens of thousands of new events each night, more than can be evaluated by astronomers alone. The ANTARES event broker will process alerts, annotating them with catalog associations and filtering them to distinguish customizable subsets of events. We describe the data model of the system, the overall architecture, annotation, implementation of filters, system outputs, provenance tracking, system performance, and the user interface.Comment: 24 Pages, 8 figures, Accepted by A

Variability of Red Supergiants in M31 from the Palomar Transient Factory

Most massive stars end their lives as red supergiants (RSGs), a short-lived evolutionary phase when they are known to pulsate with varying amplitudes. The RSG period–luminosity (PL) relation has been measured in the Milky Way, the Magellanic Clouds and M33 for about 120 stars in total. Using over 1500 epochs of R-band monitoring from the Palomar Transient Factory survey over a five-year period, we study the variability of 255 spectroscopically cataloged RSGs in M31. We find that all RGSs brighter than M_K ≈ −10 mag (log(L/L⊙) > 4.8) are variable at Δm_R > 0.05 mag. Our period analysis finds 63 with significant pulsation periods. Using the periods found and the known values of M K for these stars, we derive the RSG PL relation in M31 and show that it is consistent with those derived earlier in other galaxies of different metallicities. We also detect, for the first time, a sequence of likely first-overtone pulsations. Comparison to stellar evolution models from MESA confirms the first-overtone hypothesis and indicates that the variable stars in this sample have 12 M⊙ < M < 24 M⊙. As these RSGs are the immediate progenitors to Type II-P core-collapse supernovae (SNe), we also explore the implication of their variability in the initial-mass estimates for SN progenitors based on archival images of the progenitors. We find that this effect is small compared to the present measurement errors