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 which 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 years) 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, tch, 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 timescales. Using a 60-night block of data straddling two nights with a cadence of around 2 minutes, we extracted tch 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.

With the advent of the Large Synoptic Survey Telescope (LSST), time-domain astronomy will be faced with an unprecedented volume and rate of data. Real-time processing of variables and transients detected by such large-scale surveys is critical to identifying the more unusual events and allocating scarce follow-up resources efficiently. We develop an algorithm to identify these novel events within a given population of variable sources. We determine the distributions of magnitude changes (dm) over time intervals (dt) for a given passband f, pf(dm|dt), and use these distributions to compute the likelihood of a test source being consistent with the population, or an outlier. We demonstrate our algorithm by applying it to the DECam multi-band time-series data of more than 2000 variable stars identified by Saha et al. (2019) in the Galactic Bulge that are largely dominated by long-period variables and pulsating stars. Our algorithm discovers 18 outlier sources in the sample, including a microlensing event, a dwarf nova, and two chromospherically active RS CVn stars, as well as sources in the Blue Horizontal Branch region of the color-magnitude diagram without any known counterparts. We compare the performance of our algorithm for novelty detection with multivariate KDE and Isolation Forest on the simulated PLAsTiCC dataset. We find that our algorithm yields comparable results despite its simplicity. Our method provides an efficient way for flagging the most unusual events in a real-time alert-broker system.

We have obtained repeated images of 6 fields towards the Galactic bulge in 5 passbands (u, g, r, i, z) with the DECam imager on the Blanco 4m telescope at CTIO. From over 1.6 billion individual photometric measurements in the field centered on Baade's window, we have detected 4877 putative variable stars. 474 of these have been confirmed as fundamental mode RR Lyrae stars, whose colors at minimum light yield line-of-sight reddening determinations as well as a reddening law towards the Galactic Bulge which differs significantly from the standard R_V = 3.1 formulation. Assuming that the stellar mix is invariant over the 3 square-degree field, we are able to derive a line-of-sight reddening map with sub-arcminute resolution, enabling us to obtain de-reddened and extinction corrected color-magnitude diagrams (CMD's) of this bulge field using up to 2.5 million well-measured stars. The corrected CMD's show unprecedented detail and expose sparsely populated sequences: e.g., delineation of the very wide red giant branch, structure within the red giant clump, the full extent of the horizontal branch, and a surprising bright feature which is likely due to stars with ages younger than 1 Gyr. We use the RR Lyrae stars to trace the spatial structure of the ancient stars, and find an exponential decline in density with Galactocentric distance. We discuss ways in which our data products can be used to explore the age and metallicity properties of the bulge, and how our larger list of all variables is useful for learning to interpret future LSST alerts.

We propose a new strategy of finding strongly-lensed supernovae (SNe) by monitoring known galaxy-scale strong-lens systems. Strongly lensed SNe are potentially powerful tools for the study of cosmology, galaxy evolution, and stellar populations, but they are extremely rare. By targeting known strongly lensed starforming galaxies, our strategy significantly boosts the detection efficiency for lensed SNe compared to a blind search. As a reference sample, we compile the 128 galaxy-galaxy strong-lens systems from the Sloan Lens ACS Survey (SLACS), the SLACS for the Masses Survey, and the Baryon Oscillation Spectroscopic Survey Emission-Line Lens Survey. Within this sample, we estimate the rates of strongly-lensed Type Ia SN (SNIa) and core-collapse SN (CCSN) to be $1.23 \pm 0.12$ and $10.4 \pm 1.1$ events per year, respectively. The lensed SN images are expected to be widely separated with a median separation of 2 arcsec. Assuming a conservative fiducial lensing magnification factor of 5 for the most highly magnified SN image, we forecast that a monitoring program with a single-visit depth of 24.7 mag (5$\sigma$ point source, $r$ band) and a cadence of 5 days can detect 0.49 strongly-lensed SNIa event and 2.1 strongly-lensed CCSN events per year within this sample. Our proposed targeted-search strategy is particularly useful for prompt and efficient identifications and follow-up observations of strongly-lensed SN candidates. It also allows telescopes with small field of views and limited time to efficiently discover strongly-lensed SNe with a pencil-beam scanning strategy.

(abridged) In the hunt for rare time-domain events, confusing exotic extragalactic phenomena with more common Galactic foreground events is an important consideration. 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. We discovered an extremely bright and rapid transient event in optical observations of M31 with the intermediate Palomar Transient Factory (iPTF). The persistent optical counterpart of this transient was previously thought to be a variable star in M31 without any dramatic flux excursions. The iPTF event initially appeared to be an extraordinarily rapid and energetic extragalactic transient, a very fast nova in M31 being one of the exciting possibilities, with a ~3mag positive flux excursion in less than a kilosecond. The nature of the source was resolved with the help of Chandra archival data, where we found an X-ray counterpart and obtained its X-ray spectrum. 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 ~7MK, typical for coronal emission from an active star. The combination of the X-ray luminosity, calculated assuming the source is located in M31 (~3.10^36 erg/s), and the color temperature exclude any type of known accreting compact object or active star in M31. 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 range and it has the absolute optical magnitude of 9.5-11.0 mag in R-band. The observed optical flare has the equivalent duration of ~95min and total energy of ~(0.3-1).10^35 erg in R-band, which places it among the brightest flares ever observed from an M-type star.

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.

The unprecedented volume and rate of transient events that will be discovered by the Large Synoptic Survey Telescope (LSST) demands that the astronomical community update its followup paradigm. Alert-brokers -- automated software system to sift through, characterize, annotate and prioritize events for followup -- will be critical tools for managing alert streams in the LSST era. The Arizona-NOAO Temporal Analysis and Response to Events System (ANTARES) is one such broker. In this work, we develop a machine learning pipeline to characterize and classify variable and transient sources only using the available multiband optical photometry. We describe three illustrative stages of the pipeline, serving the three goals of early, intermediate and retrospective classification of alerts. The first takes the form of variable vs transient categorization, the second, a multi-class typing of the combined variable and transient dataset, and the third, a purity-driven subtyping of a transient class. While several similar algorithms have proven themselves in simulations, we validate their performance on real observations for the first time. We quantitatively evaluate our pipeline on sparse, unevenly sampled, heteroskedastic data from various existing observational campaigns, and demonstrate very competitive classification performance. We describe our progress towards adapting the pipeline developed in this work into a real-time broker working on live alert streams from time-domain surveys.

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

Novae undergo a supersoft X-ray phase of varying duration after the optical outburst. Such transient post-nova supersoft X-ray sources (SSSs) are the majority of the observed SSSs in M31. In this paper, we use the post-nova evolutionary models of Wolf et al. to compute the expected population of post-nova SSSs in M31. We predict that depending on the assumptions about the WD mass distribution in novae, at any instant there are about 250-600 post-nova SSSs in M31 with (unabsorbed) 0.2-1.0 keV luminosity L_x>10^36 erg/s. Their combined unabsorbed luminosity is of the order of ~10^39 erg/s. Their luminosity distribution shows significant steepening around log(L_x)~37.7-38 and becomes zero at L_x~2x10^38 erg/s, the maximum L_x achieved in the post-nova evolutionary tracks. Their effective temperature distribution has a roughly power law shape with differential slope of ~4-6 up to the maximum temperature of T_eff~1.5x10^6 K. We compare our predictions with the results of the XMM-Newton monitoring of the central field of M31 between 2006 and 2009. The predicted number of post-nova SSSs exceed the observed number by a factor of ~2-5, depending on the assumed WD mass distribution in novae. This is good agreement, considering the number and magnitude of uncertainties involved in calculations of the post-nova evolutionary models and their X-ray output. Furthermore, only a moderate circumstellar absorption, with hydrogen column density of the order of ~10^21 cm^-2, will remove the discrepancy.

(Abridged) White dwarfs (WDs) undergoing unstable nuclear burning on their surfaces, resulting in the nova phenomenon, have been considered as one of the prospective candidates for the still elusive progenitors of SNeIa. We propose that statistics of novae in nearby galaxies may be a powerful tool to gauge the role of such systems in producing SNeIa. Using multicycle nova evolutionary models, we compute the number and temporal distribution of novae that would be produced by a typical (unstably burning) SNIa progenitor before reaching the Chandrasekhar mass limit (M_ch) and exploding. Combined with the observed nova rates, this can be used to constrain the maximal contribution of the nova channel to the SNIa rate in nearby galaxies. The M31 nova rate measured by the POINT-AGAPE survey is 65/yr. Assuming that all these novae will reach M_ch, we obtain the maximal SNIa rate novae may produce to be (1-5)x10^-4/yr. This constrains the overall contribution of the nova channel to the SNIa rate in M31 at <2-7%. Furthermore, ~80% of novae generated by a successful SNIa progenitor have short decay times (t2<10 days), and are produced when the WD mass is close to M_ch. We point out that statistics of such fast novae can provide a sensitive diagnostic of the contribution of the nova channel to the final stages of mass accumulation by the single degenerate (SD) SNIa progenitors. To explore the prospects of their use, we investigate the efficiency of detecting fast novae in an M31 nova survey of the PTF class. We find that a survey with limiting magnitude of m_R~22 observing at least every 2nd night will catch ~90% of fast novae expected in the SD scenario. Such surveys should be detecting fast novae in M31 at a rate of the order of >10^3xf per yr, where f is the fraction of SNeIa which accreted in the unstable nuclear burning regime while accumulating the final ~0.1M_sun before the SNIa explosion.