A Machine Learning Method to Infer Fundamental Stellar Parameters from Photometric Light Curves

A fundamental challenge for wide-field imaging surveys is obtaining follow-up spectroscopic observations: there are > $10^9$ photometrically cataloged sources, yet modern spectroscopic surveys are limited to ~few x $10^6$ targets. As we approach the Large Synoptic Survey Telescope (LSST) era, new algorithmic solutions are required to cope with the data deluge. Here we report the development of a machine-learning framework capable of inferring fundamental stellar parameters (Teff, log g, and [Fe/H]) using photometric-brightness variations and color alone. A training set is constructed from a systematic spectroscopic survey of variables with Hectospec/MMT. In sum, the training set includes ~9000 spectra, for which stellar parameters are measured using the SEGUE Stellar Parameters Pipeline (SSPP). We employed the random forest algorithm to perform a non-parametric regression that predicts Teff, log g, and [Fe/H] from photometric time-domain observations. Our final, optimized model produces a cross-validated root-mean-square error (RMSE) of 165 K, 0.39 dex, and 0.33 dex for Teff, log g, and [Fe/H], respectively. Examining the subset of sources for which the SSPP measurements are most reliable, the RMSE reduces to 125 K, 0.37 dex, and 0.27 dex, respectively, comparable to what is achievable via low-resolution spectroscopy. For variable stars this represents a ~12-20% improvement in RMSE relative to models trained with single-epoch photometric colors. As an application of our method, we estimate stellar parameters for ~54,000 known variables. We argue that this method may convert photometric time-domain surveys into pseudo-spectrographic engines, enabling the construction of extremely detailed maps of the Milky Way, its structure, and history.

Discovery of Bright Galactic R Coronae Borealis and DY Persei Variables: Rare Gems Mined from ACVS

We present the results of a machine-learning (ML) based search for new R Coronae Borealis (RCB) stars and DY Persei-like stars (DYPers) in the Galaxy using cataloged light curves from the All-Sky Automated Survey (ASAS) Catalog of Variable Stars (ACVS). RCB stars - a rare class of hydrogen-deficient carbon-rich supergiants - are of great interest owing to the insights they can provide on the late stages of stellar evolution. DYPers are possibly the low-temperature, low-luminosity analogs to the RCB phenomenon, though additional examples are needed to fully establish this connection. While RCB stars and DYPers are traditionally identified by epochs of extreme dimming that occur without regularity, the ML search framework more fully captures the richness and diversity of their photometric behavior. We demonstrate that our ML method can use newly discovered RCB stars to identify additional candidates within the same data set. Our search yields 15 candidates that we consider likely RCB stars/DYPers: new spectroscopic observations confirm that four of these candidates are RCB stars and four are DYPers. Our discovery of four new DYPers increases the number of known Galactic DYPers from two to six; noteworthy is that one of the new DYPers has a measured parallax and is m ~ 7 mag, making it the brightest known DYPer to date. Future observations of these new DYPers should prove instrumental in establishing the RCB connection. We consider these results, derived from a machine-learned probabilistic classification catalog, as an important proof-of-concept for the efficient discovery of rare sources with time-domain surveys.

Evidence for an FU Orionis-like Outburst from a Classical T Tauri Star

We present pre- and post-outburst observations of the new FU Orionis-like young stellar object PTF 10qpf (also known as LkHa 188-G4 and HBC 722). Prior to this outburst, LkHa 188-G4 was classified as a classical T Tauri star on the basis of its optical emission-line spectrum superposed on a K8-type photosphere, and its photometric variability. The mid-infrared spectral index of LkHa 188-G4 indicates a Class II-type object. LkHa 188-G4 exhibited a steady rise by ~1 mag over ~11 months starting in Aug. 2009, before a subsequent more abrupt rise of > 3 mag on a time scale of ~2 months. Observations taken during the eruption exhibit the defining characteristics of FU Orionis variables: (i) an increase in brightness by > 4 mag, (ii) a bright optical/near-infrared reflection nebula appeared, (iii) optical spectra are consistent with a G supergiant and dominated by absorption lines, the only exception being Halpha which is characterized by a P Cygni profile, (iv) near-infrared spectra resemble those of late K--M giants/supergiants with enhanced absorption seen in the molecular bands of CO and H_2O, and (v) outflow signatures in H and He are seen in the form of blueshifted absorption profiles. LkHa 188-G4 is the first member of the FU Orionis-like class with a well-sampled optical to mid-infrared spectral energy distribution in the pre-outburst phase. The association of the PTF 10qpf outburst with the previously identified classical T Tauri star LkHa 188-G4 (HBC 722) provides strong evidence that FU Orionis-like eruptions represent periods of enhanced disk accretion and outflow, likely triggered by instabilities in the disk. The early identification of PTF 10qpf as an FU Orionis-like variable will enable detailed photometric and spectroscopic observations during its post-outburst evolution for comparison with other known outbursting objects.

SN 2008iy: An Unusual Type IIn Supernova with an Enduring 400 Day Rise Time

We present spectroscopic and photometric observations of the Type IIn supernova (SN) 2008iy. SN 2008iy showed an unprecedentedly long rise time of ~400 days, making it the first SN to take significantly longer than 100 days to reach peak optical luminosity. The peak absolute magnitude of SN 2008iy was M_r ~ -19.1 mag, and the total radiated energy over the first ~700 days was ~2 x 10^50 erg. Spectroscopically, SN 2008iy is very similar to the Type IIn SN 1988Z at late times, and, like SN 1988Z, it is a luminous X-ray source (both supernovae had an X-ray luminosity L_ X > 10^41 erg/s). The Halpha emission profile of SN 2008iy shows a narrow P Cygni absorption component, implying a pre-SN wind speed of ~100 km/s. We argue that the luminosity of SN 2008iy is powered via the interaction of the SN ejecta with a dense, clumpy circumstellar medium. The ~400 day rise time can be understood if the number density of clumps increases with distance over a radius ~1.7 x 10^16 cm from the progenitor. This scenario is possible if the progenitor experienced an episodic phase of enhanced mass-loss < 1 century prior to explosion or the progenitor wind speed increased during the decades before core collapse. We favour the former scenario, which is reminiscent of the eruptive mass-loss episodes observed for luminous blue variable (LBV) stars. The progenitor wind speed and increased mass-loss rates serve as further evidence that at least some, and perhaps all, Type IIn supernovae experience LBV-like eruptions shortly before core collapse. We also discuss the host galaxy of SN 2008iy, a subluminous dwarf galaxy, and offer a few reasons why the recent suggestion that unusual, luminous supernovae preferentially occur in dwarf galaxies may be the result of observational biases.

New Observations of the Very Luminous Supernova 2006gy: Evidence for Echoes

Supernova (SN) 2006gy was a hydrogen-rich core-collapse SN that remains one of the most luminous optical supernovae ever observed. The total energy budget (> 2 x 10^51 erg radiated in the optical alone) poses many challenges for standard SN theory. We present new ground-based near-infrared (NIR) observations of SN 2006gy, as well as a single epoch of Hubble Space Telescope (HST) imaging obtained more than two years after the explosion. Our NIR data taken around peak optical emission show an evolution that is largely consistent with a cooling blackbody, with tentative evidence for a growing NIR excess starting at day ~100. Our late-time Keck adaptive optics (AO) NIR image, taken on day 723, shows little change from previous NIR observations taken around day 400. Furthermore, the optical HST observations show a reduced decline rate after day 400, and the SN is bluer on day 810 than it was at peak. This late-time decline is inconsistent with Co56 decay, and thus is problematic for the various pair-instability SN models used to explain the nature of SN 2006gy. The slow decline of the NIR emission can be explained with a light echo, and we confirm that the late-time NIR excess is the result of a massive (>10 Msun) dusty shell heated by the SN peak luminosity. The late-time optical observations require the existence of a scattered light echo, which may be generated by the same dust that contributes to the NIR echo. Both the NIR and optical echoes originate in the proximity of the progenitor, ~10^18 cm for the NIR echo and <~10-40 pc for the optical echo, which provides further evidence that the progenitor of SN 2006gy was a very massive star.