Search

Luminous interacting supernovae are a class of stellar explosions whose progenitors underwent vigorous mass loss in the years prior to core-collapse. While the mechanism by which this material is ejected is still debated, obtaining the full density profile of the circumstellar medium (CSM) could reveal more about this process. Here, we present an extensive multi-wavelength study of PS1-11aop, a luminous and slowly declining Type IIn SN discovered by the PanSTARRS Medium Deep Survey. PS1-11aop had a peak r-band magnitude of $-$20.5\,mag, a total radiated energy $>$ 8$\times$10$^{50}$\,erg, and it exploded near the center of a star-forming galaxy with super-solar metallicity. We obtained multiple detections at the location of PS1-11aop in the radio and X-ray bands between 4 and 10\,years post-explosion, and if due to the SN, it is one of the most luminous radio supernovae identified to date. Taken together, the multiwavelength properties of PS1-11aop are consistent with a CSM density profile with multiple zones. The early optical emission is consistent with the supernova blastwave interacting with a dense and confined CSM shell which contains multiple solar masses of material that was likely ejected in the final $<$10-100 years prior to the explosion,($\sim$0.05$-$1.0 M$_{\odot}$yr$^{-1}$ at radii of $\lesssim$10$^{16}$\,cm). The radio observations, on the other hand, are consistent with a sparser environment ($\lesssim$2$\times 10^{-3}$ M$_{\odot}$yr$^{-1}$ at radii of $\sim$0.5-1$\times$10$^{17}$\,cm) -- thus probing the history of the progenitor star prior to its final mass loss episode., Comment: 37 pages, 17 figures

The identification of persistent radio sources (PRSs) coincident with two repeating fast radio bursts (FRBs) supports FRB theories requiring a compact central engine. However, deep non-detections in other cases highlight the diversity of repeating FRBs and their local environments. Here, we perform a systematic search for radio sources towards 37 CHIME/FRB repeaters using their arcminute localizations and a combination of archival surveys and targeted observations. Through multi-wavelength analysis of individual radio sources, we identify two (20181030A-S1 and 20190417A-S1) for which we disfavor an origin of either star formation or an active galactic nucleus in their host galaxies and thus consider them candidate PRSs. We do not find any associated PRSs for the majority of the repeating FRBs in our sample. For 8 FRB fields with Very Large Array imaging, we provide deep limits on the presence of PRSs that are 2--4 orders of magnitude fainter than the PRS associated with FRB\,20121102A. Using Very Large Array Sky Survey imaging of all 37 fields, we constrain the rate of luminous ($\gtrsim$10$^{40}$ erg s$^{-1}$) PRSs associated with repeating FRBs to be low. Within the context of FRB-PRS models, we find that 20181030A-S1 and 20190417A-S1 can be reasonably explained within the context of magnetar, hypernebulae, gamma-ray burst afterglow, or supernova ejecta models -- although we note that both sources follow the radio luminosity versus rotation measure relationship predicted in the nebula model framework. Future observations will be required to both further characterize and confirm the association of these PRS candidates with the FRBs., Comment: 37 pages, 12 figures

We conduct a systematic analysis of the early multi-band light curves and colors of 19 Type Ia Supernovae (SNe) from the Korea Microlensing Telescope Network SN Program, including 16 previously unpublished events. Seven are detected $\lesssim$ 1 day since the estimated epoch of first light and the rest within $\lesssim$ 3 days. Some show excess emission within $<$ 0.5 days to $\sim$ 2 days, but most show pure power-law rises. The colors are initially diverse before $\sim$ 5 days, but converge to a similar color at $\sim$ 10 days. We identify at least three populations based on 2--5-day color evolution: (1) "early-blues" exhibit slowly-evolving colors consistent with a $\sim$ 17,000 K blackbody; (2) "early-reds" have initially blue $B-V$ and red $V-i$ colors that cannot simultaneously be fit with a blackbody -- likely due to suppression of $B$- and $i$-band flux by Fe II/III and Ca II -- and evolve more rapidly; and (3) "early-yellows" evolve blueward, consistent with thermal heating from $\sim$ 8,000 to 13,000 K. The distributions of early-blue and early-red colors are compatible with them being either distinct populations -- with early-reds comprising (60 $\pm$ 15)% of them -- or extreme ends of one continuous population; whereas the early-yellow population identified here is clearly distinct. Compared to the other populations, early-blues in our sample differ by exhibiting excess emission within 1--2 days, nearly constant peak brightness regardless of $\Delta M_{15}(B)$ after standardization, and shallower Si II features. Early-blues also prefer star-forming host environments, while early-yellows and, to a lesser extent, early-reds prefer quiescent ones. These preferences appear to indicate at least two Type Ia SN production channels based on stellar population age, while early-reds and early-blues may still share a common origin., Comment: Submitted for publication in ApJ. 48 pages, 29 figures, 7 tables

Close binary systems are the progenitors to both Type Ia supernovae and the compact object mergers that can be detected via gravitational waves. To achieve a binary with a small radial separation, it is believed that the system likely undergoes common envelope (CE) evolution. Despite its importance, CE evolution may be one of the largest uncertainties in binary evolution due to a combination of computational challenges and a lack of observed benchmarks where both the post-CE and pre-CE conditions are known. Identifying post-CE systems in star clusters can partially circumvent this second issue by providing an independent age constraint on the system. For the first time, we conduct a systematic search for white dwarf (WD) and main-sequence (MS) binary systems in 299 Milky Way open star clusters. Coupling Gaia DR3 photometry and kinematics with multi-band photometry from Pan-STARRS1 and 2MASS, we apply a machine learning based approach and find 52 high-probability candidates in 38 open clusters. For a subset of our systems, we present follow-up spectroscopy from the Gemini and Lick Observatories and archival light curves from TESS, Kepler/K2 and the Zwicky Transient Facility. Examples of M-dwarfs with hot companions are spectroscopically observed, along with regular system variability. While the kinematics of our candidates are consistent with their host clusters, some systems have spatial positions offset relative to their hosts, potentially indicative of natal kicks. Ultimately, this catalogue is a first step to obtaining a set of observational benchmarks to better link post-CE systems to their pre-CE progenitors., Comment: 35 pages, 18 figures. Submitted to ApJ

Recent works have constrained the binary fraction of evolved populations of massive stars in local galaxies such as red supergiants and Wolf-Rayet stars, but the binary fraction of yellow supergiants (YSGs) in the Hertzsprung Gap remains unconstrained. Binary evolution theory predicts that the Hertzsprung Gap is home to multiple populations of binary systems with varied evolutionary histories. In this paper, we develop a method to distinguish single YSGs from YSG plus O- or B-type main sequence binaries using optical and ultraviolet photometry, and then apply this method to identify candidate YSG binaries in the Magellanic Clouds. After constructing a set of combined stellar atmosphere models, we find that optical photometry is, given typical measurement and reddening uncertainties, sufficient to discern single YSGs from YSG+OB binaries if the OB-star is at least $\sim5$M$_{\odot}$ for T$_{\mathrm{eff,YSG}}\sim$ 4000 K, but requires a $\sim$20M$_{\odot}$ OB star for YSGs up to T$_{\mathrm{eff,YSG}}\sim$ 9000 K. For these hotter YSG temperatures, ultraviolet photometry allows binaries with OB companions as small as $\sim$7M$_{\odot}$ to be identified. We use color-color spaces developed from these models to search for evidence of excess blue or ultraviolet light in a set of $\sim$1000 YSG candidates in the Magellanic Clouds. We identify hundreds of candidate YSG binary systems and report a preliminary fraction of YSGs that show a blue/UV color excess of 20-60\%. Spectroscopic follow-up is now required to confirm the true nature of this population., Comment: 21 pages, 10 figures, 3 tables. Submitted to The Astrophysical Journal on June 24th 2024, Accepted on September 3rd

We present a systematic search for tidal disruption events (TDEs) using radio data from the Variables and Slow Transients (VAST) Pilot Survey conducted using the Australian Square Kilometre Array Pathfinder (ASKAP). Historically, TDEs have been identified using observations at X-ray, optical, and ultraviolet wavelengths. After discovery, a few dozen TDEs have been shown to have radio counterparts through follow-up observations. With systematic time-domain radio surveys becoming available, we can now identify new TDEs in the radio regime. A population of radio-discovered TDEs has the potential to provide several key insights including an independent constraint on their volumetric rate. We conducted a search to select variable radio sources with a single prominent radio flare and a position consistent within 2$\sigma$ of the nucleus of a known galaxy. While TDEs were the primary target of our search, sources identified in this search may also be consistent with active galactic nuclei exhibiting unusual flux density changes at the timescales probed, uncharacteristically bright supernovae, or a population of gamma-ray bursts. We identify a sample of 12 radio-bright candidate TDEs. The timescales and luminosities range from ~6 to 230 days and ~10$^{38}$ to 10$^{41}$ erg s$^{-1}$, consistent with models of radio emission from TDEs that launch relativistic jets. After calculating the detection efficiency of our search using a Monte Carlo simulation of TDEs, and assuming all 12 sources are jetted TDEs, we derive a volumetric rate for jetted TDEs of 0.80$^{+0.31}_{-0.23}$ Gpc$^{-3}$ yr$^{-1}$, consistent with previous empirically estimated rates., Comment: 37 pages, 37 figures, accepted to ApJ

The nearby type II supernova, SN2023ixf in M101 exhibits signatures of early-time interaction with circumstellar material in the first week post-explosion. This material may be the consequence of prior mass loss suffered by the progenitor which possibly manifested in the form of a detectable pre-supernova outburst. We present an analysis of the long-baseline pre-explosion photometric data in $g$, $w$, $r$, $i$, $z$ and $y$ filters from Pan-STARRS as part of the Young Supernova Experiment, spanning $\sim$5,000 days. We find no significant detections in the Pan-STARRS pre-explosion light curve. We train a multilayer perceptron neural network to classify pre-supernova outbursts. We find no evidence of eruptive pre-supernova activity to a limiting absolute magnitude of $-7$. The limiting magnitudes from the full set of $gwrizy$ (average absolute magnitude $\approx$-8) data are consistent with previous pre-explosion studies. We use deep photometry from the literature to constrain the progenitor of SN2023ixf, finding that these data are consistent with a dusty red supergiant (RSG) progenitor with luminosity $\log\left(L/L_\odot\right)$$\approx$5.12 and temperature $\approx$3950K, corresponding to a mass of 14-20 M$_\odot$, Comment: 19 pages, 8 figures, 1 table

In kilonovae, freshly-synthesized $r$-process elements imprint features on optical spectra, as observed in AT2017gfo, the counterpart to the GW170817 binary neutron star merger. However, measuring the $r$-process compositions of the merger ejecta is computationally challenging. Vieira et al. (2023) introduced Spectroscopic $r$-Process Abundance Retrieval for Kilonovae (SPARK), a software tool to infer elemental abundance patterns of the ejecta, and associate spectral features with particular species. Previously, we applied SPARK to the 1.4 day spectrum of AT2017gfo and inferred its abundance pattern for the first time, characterized by electron fraction $Y_e=0.31$, a substantial abundance of strontium, and a dearth of lanthanides and heavier elements. This ejecta is consistent with wind from a remnant hypermassive neutron star and/or accretion disk. We now extend our inference to spectra at 2.4 and 3.4 days, and test the need for multicomponent ejecta, where we stratify the ejecta in composition. The ejecta at 1.4 and 2.4 days is described by the same single blue component. At 3.4 days, a new redder component with lower $Y_e=0.16$ and a significant abundance of lanthanides emerges. This new redder component is consistent with dynamical ejecta and/or neutron-rich ejecta from a magnetized accretion disk. As expected from photometric modelling, this component emerges as the ejecta expands, the photosphere recedes, and the earlier bluer component dims. At 3.4 days, we find an ensemble of lanthanides, with the presence of cerium most concrete. This presence of lanthanides has important implications for the contribution of kilonovae to the $r$-process abundances observed in the Universe., Comment: Accepted to ApJ; 19 pages, 8 figures in body; 10 pages, 7 figures in appendix

We present optical, ultraviolet, and infrared data of the type II supernova (SN II) 2020jfo at 14.5 Mpc. This wealth of multiwavelength data allows to compare different metrics commonly used to estimate progenitor masses of SN II for the same object. Using its early light curve, we infer SN 2020jfo had a progenitor radius of $\approx$700 $R_{\odot}$, consistent with red supergiants of initial mass $M_{\rm ZAMS}=$11-13 $M_{\odot}$. The decline in its late-time light curve is best fit by a ${}^{56}$Ni mass of 0.018$\pm$0.007 $M_{\odot}$ consistent with that ejected from SN II-P with $\approx$13 $M_{\odot}$ initial mass stars. Early spectra and photometry do not exhibit signs of interaction with circumstellar matter, implying that SN 2020jfo experienced weak mass loss within the final years prior to explosion. Our spectra at $>$250 days are best fit by models from 12 $M_{\odot}$ initial mass stars. We analyzed integral field unit spectroscopy of the stellar population near SN 2020jfo, finding its massive star population had a zero age main sequence mass of 9.7$\substack{+2.5\\-1.3} M_{\odot}$. We identify a single counterpart in pre-explosion imaging and find it has an initial mass of at most $7.2\substack{+1.2\\-0.6} M_{\odot}$. We conclude that the inconsistency between this mass and indirect mass indicators from SN 2020jfo itself is most likely caused by extinction with $A_{V}=2$-3 mag due to matter around the progenitor star, which lowered its observed optical luminosity. As SN 2020jfo did not exhibit extinction at this level or evidence for interaction with circumstellar matter between 1.6-450 days from explosion, we conclude that this material was likely confined within $\approx$3000 $R_{\odot}$ from the progenitor star., Comment: 31 pages, 12 figures, accepted to MNRAS

We present pre-explosion optical and infrared (IR) imaging at the site of the type II supernova (SN II) 2023ixf in Messier 101 at 6.9 Mpc. We astrometrically registered a ground-based image of SN 2023ixf to archival Hubble Space Telescope (HST), Spitzer Space Telescope (Spitzer), and ground-based near-IR images. A single point source is detected at a position consistent with the SN at wavelengths ranging from HST $R$-band to Spitzer 4.5 $\mu$m. Fitting to blackbody and red supergiant (RSG) spectral-energy distributions (SEDs), we find that the source is anomalously cool with a significant mid-IR excess. We interpret this SED as reprocessed emission in a 8600 $R_{\odot}$ circumstellar shell of dusty material with a mass $\sim$5$\times10^{-5} M_{\odot}$ surrounding a $\log(L/L_{\odot})=4.74\pm0.07$ and $T_{\rm eff}=3920\substack{+200\\-160}$ K RSG. This luminosity is consistent with RSG models of initial mass 11 $M_{\odot}$, depending on assumptions of rotation and overshooting. In addition, the counterpart was significantly variable in pre-explosion Spitzer 3.6 $\mu$m and 4.5 $\mu$m imaging, exhibiting $\sim$70% variability in both bands correlated across 9 yr and 29 epochs of imaging. The variations appear to have a timescale of 2.8 yr, which is consistent with $\kappa$-mechanism pulsations observed in RSGs, albeit with a much larger amplitude than RSGs such as $\alpha$ Orionis (Betelgeuse)., Comment: 14 pages, 5 figures, submitted to ApJL, comments welcome

We present a search for host galaxy associations for the third set of repeating fast radio burst (FRB) sources discovered by the CHIME/FRB Collaboration. Using the $\sim$ 1 arcmin CHIME/FRB baseband localizations and probabilistic methods, we identify potential host galaxies of two FRBs, 20200223B and 20190110C at redshifts of 0.06024(2) and 0.12244(6), respectively. We also discuss the properties of a third marginal candidate host galaxy association for FRB 20191106C with a host redshift of 0.10775(1). The three putative host galaxies are all relatively massive, fall on the standard mass-metallicity relationship for nearby galaxies, and show evidence of ongoing star formation. They also all show signatures of being in a transitional regime, falling in the ``green valley'' which is between the bulk of star-forming and quiescent galaxies. The plausible host galaxies identified by our analysis are consistent with the overall population of repeating and non-repeating FRB hosts while increasing the fraction of massive and bright galaxies. Coupled with these previous host associations, we identify a possible excess of FRB repeaters whose host galaxies have $M_{\mathrm{u}}-M_{\mathrm{r}}$ colors redder than the bulk of star-forming galaxies. Additional precise localizations are required to confirm this trend., Comment: 23 pages, 9 figures, submitted to AAS journals

\object{SN 2021aefx} is a normal Type Ia Supernova (SN) with red excess emission over the first $\sim$ 2 days. We present detailed analysis of this SN using our high-cadence KMTNet multi-band photometry, spectroscopy, and publicly available data. We provide the first measurements of its epochs of explosion (MJD 59529.32 $\pm$ 0.16) as well as ``first light'' (MJD 59529.85 $\pm$ 0.55) associated with the main ejecta ${\rm{^{56}Ni}}$ distribution. This places our first detection of SN 2021aefx at $\sim -$0.5 hours since ``first light'', indicating the presence of additional power sources. Our peak-spectrum confirms its Type Ia sub-classification as intermediate between Core-Normal and Broad-Line, and we estimate the ejecta mass to be $\sim$ 1.34 $M_{\odot}$. The pre-peak spectral evolution identifies fast-expanding material reaching $>$ 40,000 km s$^{-1}$ (the fastest ever observed in Type Ia SNe) and at least two distinct homologously-expanding ejecta components: (1) a normal-velocity (12,400 km s$^{-1}$) component consistent with the typical photospheric evolution of Chandrasekhar-mass ejecta; and (2) a high-velocity (23,500 km s$^{-1}$) component visible during the first $\sim$ 3.6 days post-explosion, which locates the component within the outer $<$ 16\% of the ejecta mass. Asymmetric, subsonic explosion processes producing a non-spherical photosphere provide an explanation for the simultaneous presence of the two components, as well as the red excess emission via a slight ${\rm{^{56}Ni}}$ enrichment in the outer $\sim$ 0.5\% of the ejecta mass. Our spectrum from 300 days post-peak advances the constraint against non-degenerate companions and further supports a near-Chandrasekhar-mass explosion origin. Off-center ignited delayed-detonations of Chandrasekhar-mass white dwarfs may be responsible for the observed features of SN 2021aefx in some normal Type Ia SNe., Comment: Submitted for publication in ApJ. 29 pages, 14 figures, 4 tables

The exact nature of the luminous fast blue optical transient AT 2018cow is still debated. In this first of a two-paper series, we present a detailed analysis of three Hubble Space Telescope (HST) observations of AT 2018cow covering $\sim$50--60 days post-discovery in combination with other observations throughout the first two months, and derive significantly improved constraints of the late thermal properties. By modeling the spectral energy distributions (SEDs), we confirm that the UV-optical emission over 50--60 days was still a smooth blackbody (i.e., optically thick) with a high temperature ($T_{\mathrm{BB}}\sim15000\,\mathrm{K}$) and small radius ($R_{\mathrm{BB}}\lesssim1000\,R_\odot$). Additionally, we report for the first time a break in the bolometric light curve: the thermal luminosity initially declined at a rate of $L_{\mathrm{BB}}\propto t^{-2.40}$, but faded much faster at $t^{-3.06}$ after day 13. Reexamining possible late-time power sources, we disfavor significant contributions from radioactive decay based on the required $^{56}$Ni mass and lack of UV line blanketing in the HST SEDs. We argue that the commonly proposed interaction with circumstellar material may face significant challenges in explaining the late thermal properties, particularly the effects of the optical depth. Alternatively, we find that continuous outflow/wind driven by a central engine can still reasonably explain the combination of a receding photosphere, optically thick and rapidly fading emission, and intermediate-width lines. However, the rapid fading may have further implications on the power output and structure of the system. Our findings may support the hypothesis that AT 2018cow and other "Cow-like transients" are powered mainly by accretion onto a central engine., Comment: 22 pages, 10 figures. This is paper I of a two-paper series. v2: updated to the accepted version

In this second of a two-paper series, we present a detailed analysis of three HST observations taken $\sim$2--4 years post-discovery, examining the evolution of a UV-bright underlying source at the precise position of AT 2018cow. While observations at $\sim$2--3 years post-discovery revealed an exceptionally blue ($L_\nu\propto \nu^{1.99}$) underlying source with relatively stable optical brightness, fading in the NUV was observed at year 4, indicating flattening in the spectrum (to $L_\nu\propto \nu^{1.64}$). The resulting spectral energy distributions can be described by an extremely hot but small blackbody, and the fading may be intrinsic (cooling) or extrinsic (increased absorption). Considering possible scenarios and explanations, we disfavor significant contributions from stellar sources and dust formation based on the observed color and brightness. By comparing the expected power and the observed luminosity, we rule out interaction with the known radio-producing circumstellar material as well as magnetar spin down with $B\sim10^{15}\,\mathrm{G}$ as possible power sources, though we cannot rule out the possible existence of a denser CSM component (e.g., previously ejected hydrogen envelope) or a magnetar with $B\lesssim10^{14}\,\mathrm{G}$. Finally, we find that a highly-inclined precessing accretion disk can reasonably explain the color, brightness, and evolution of the underlying source. However, a major uncertainty in this scenario is the mass of the central black hole (BH), as both stellar-mass and intermediate-mass BHs face notable challenges that cannot be explained by our simple disk model, and further observations and theoretical works are needed to fully constrain the nature of this underlying source., Comment: 24 pages, 12 figures. This is paper II of a two-paper series. v2: updated to the accepted version. Note a mistake in v1 in the caption of Figure 6 is now fixed: removed the word "dereddened" because the plotted observed magnitudes are directly from Table 1 which are not extinction-corrected (the analyses are not affected by this)

In previous work we identified a population of 38 cool and luminous variable stars in the Magellanic Clouds and examined 11 in detail in order to classify them as either Thorne-\.Zytkow Objects (T\.ZOs, red supergiants with a neutron star cores) or super-AGB stars (the most massive stars that will not undergo core collapse). This population includes HV\,2112, a peculiar star previously considered in other works to be either a T\.ZO or high-mass AGB star. Here we continue this investigation, using the kinematic and radio environments and local star formation history of these stars to place constraints on the age of the progenitor systems and the presence of past supernovae. These stars are not associated with regions of recent star formation, and we find no evidence of past supernovae at their locations. Finally, we also assess the presence of heavy elements and lithium in their spectra compared to red supergiants. We find strong absorption in Li and s-process elements compared to RSGs in most of the sample, consistent with super-AGB nucleosynthesis, while HV\,2112 shows additional strong lines associated with T\.ZO nucleosynthesis. Coupled with our previous mass estimates, the results are consistent with the stars being massive (~4-6.5M$_{\odot}$) or super-AGB (~6.5-12M$_{\odot}$) stars in the thermally pulsing phase, providing crucial observations of the transition between low- and high-mass stellar populations. HV\,2112 is more ambiguous; it could either be a maximally massive sAGB star, or a T\.ZO if the minimum mass for stability extends down to <13 M$_\odot$., Comment: 32 pages, 10 figures, 7 tables, spectroscopic data available at https://zenodo.org/record/7058608, accepted to The Astrophysical Journal

Freshly-synthesized r-process elements in kilonovae ejecta imprint absorption features on optical spectra, as observed in the GW170817 binary neutron star merger. These spectral features encode insights into the physical conditions of the r-process and the origins of the ejected material, but associating features with particular elements and inferring the resultant abundance pattern is computationally challenging. We introduce Spectroscopic r-Process Abundance Retrieval for Kilonovae (SPARK), a modular framework to perform Bayesian inference on kilonova spectra with the goals of inferring elemental abundance patterns and identifying absorption features at early times. SPARK inputs an atomic line list and abundance patterns from reaction network calculations into the TARDIS radiative transfer code. It then performs fast Bayesian inference on observed kilonova spectra by training a Gaussian process surrogate for the approximate posteriors of kilonova ejecta parameters, via active learning. We use the spectrum of GW170817 at 1.4 days to perform the first inference on a kilonova spectrum, and recover a complete abundance pattern. Our inference shows that this ejecta was generated by an r-process with either (1) high electron fraction Y_e ~ 0.35 and high entropy s/k_B ~ 25, or, (2) a more moderate Y_e ~ 0.30 and s/k_B ~ 14. These parameters are consistent with a shocked, polar dynamical component, and a viscously-driven outflow from a remnant accretion disk, respectively. We also recover previous identifications of strontium absorption at ~8000 AA, and tentatively identify yttrium and/or zirconium at < 4500 AA. Our approach will enable computationally-tractable inference on the spectra of future kilonovae discovered through multi-messenger observations., Comment: Accepted to ApJ; 34 pages, 13 figures

We present a low-cost ultraviolet to infrared absolute quantum efficiency detector characterization system developed using commercial off-the-shelf components. The key components of the experiment include a light source,a regulated power supply, a monochromator, an integrating sphere, and a calibrated photodiode. We provide a step-by-step procedure to construct the photon and quantum efficiency transfer curves of imaging sensors. We present results for the GSENSE 2020 BSI CMOS sensor and the Sony IMX 455 BSI CMOS sensor. As a reference for similar characterizations, we provide a list of parts and associated costs along with images of our setup., Comment: 15 pages, 14 figures, SPIE Astronomical Telescopes and Instrumentation, Conference 12191: X-ray, Optical, and Infrared Detectors for Astronomy X, Paper Number: 12191-39 (Montreal, July 2022)

SN~2018aoz is a Type Ia SN with a $B$-band plateau and excess emission in the infant-phase light curves $\lesssim$ 1 day after first light, evidencing an over-density of surface iron-peak elements as shown in our previous study. Here, we advance the constraints on the nature and origin of SN~2018aoz based on its evolution until the nebular phase. Near-peak spectroscopic features show the SN is intermediate between two subtypes of normal Type Ia: Core-Normal and Broad-Line. The excess emission could have contributions from the radioactive decay of surface iron-peak elements as well as ejecta interaction with either the binary companion or a small torus of circumstellar material. Nebular-phase limits on H$\alpha$ and He~I favour a white dwarf companion, consistent with the small companion size constrained by the low early SN luminosity, while the absence of [O~I] and He~I disfavours a violent merger of the progenitor. Of the two main explosion mechanisms proposed to explain the distribution of surface iron-peak elements in SN~2018aoz, the asymmetric Chandrasekhar-mass explosion is less consistent with the progenitor constraints and the observed blueshifts of nebular-phase [Fe~II] and [Ni~II]. The helium-shell double-detonation explosion is compatible with the observed lack of C spectral features, but current 1-D models are incompatible with the infant-phase excess emission, $B_{\rm max}-V_{\rm max}$ color, and absence of nebular-phase [Ca~II]. Although the explosion processes of SN~2018aoz still need to be more precisely understood, the same processes could produce a significant fraction of Type Ia SNe that appear normal after $\sim$ 1 day., Comment: Submitted for publication in ApJ. 35 pages, 16 figures, 7 tables

Type Ia Supernovae are thermonuclear explosions of white dwarf stars. They play a central role in the chemical evolution of the Universe and are an important measure of cosmological distances. However, outstanding questions remain about their origins. Despite extensive efforts to obtain natal information from their earliest signals, observations have thus far failed to identify how the majority of them explode. Here, we present infant-phase detections of SN 2018aoz from a brightness of -10.5 absolute AB magnitudes -- the lowest luminosity early Type Ia signals ever detected -- revealing a hitherto unseen plateau in the $B$-band that results in a rapid redward color evolution between 1.0 and 12.4 hours after the estimated epoch of first light. The missing $B$-band flux is best-explained by line-blanket absorption from Fe-peak elements in the outer 1% of the ejected mass. The observed $B-V$ color evolution of the SN also matches the prediction from an over-density of Fe-peak elements in the same outer 1% of the ejected mass, whereas bluer colors are expected from a purely monotonic distribution of Fe-peak elements. The presence of excess nucleosynthetic material in the extreme outer layers of the ejecta points to enhanced surface nuclear burning or extended sub-sonic mixing processes in some normal Type Ia Supernova explosions., Comment: Accepted for publication in Nature Astronomy. Main text = 8 pages, 4 figures, 1 table; Full document = 46 pages, with Methods, Supplementary Information, 7 Supplementary figures, 2 Supplementary tables and references. Nat Astron (2022)

We present the first ever discovery of a short-period and unusually helium-deficient dwarf nova KSP-OT-201701a by the Korea Microlensing Telescope Network Supernova Program. The source shows three superoutbursts, each led by a precursor outburst, and several normal outbursts in BVI during the span of ~2.6 years with supercycle and normal cycle lengths of about 360 and 76 days, respectively. Spectroscopic observations near the end of a superoutburst reveal the presence of strong double-peaked HI emission lines together with weak HeI emission lines. The helium-to-hydrogen intensity ratios measured by HeI{\lambda}5876 and H{\alpha} lines are 0.10 {\pm} 0.01 at a quiescent phase and 0.26 {\pm} 0.04 at an outburst phase, similar to the ratios found in long-period dwarf novae while significantly lower than those in helium cataclysmic variables (He CVs). Its orbital period of 51.91 {\pm} 2.50 minutes, which is estimated based on time series spectroscopy, is a bit shorter than the superhump period of 56.52 {\pm} 0.19 minutes, as expected from the gravitational interaction between the eccentric disk and the secondary star. We measure its mass ratio to be 0.37^{+0.32}_{-0.21} using the superhump period excess of 0.089 {\pm} 0.053. The short orbital period, which is under the period minimum, the unusual helium deficiency, and the large mass ratio suggest that KSP-OT-201701a is a transition object evolving to a He CV from a long-period dwarf nova with an evolved secondary star., Comment: 11 pages, 5 figures, published in ApJL

Periodic variables illuminate the physical processes of stars throughout their lifetime. Wide-field surveys continue to increase our discovery rates of periodic variable stars. Automated approaches are essential to identify interesting periodic variable stars for multi-wavelength and spectroscopic follow-up. Here, we present a novel unsupervised machine learning approach to hunt for anomalous periodic variables using phase-folded light curves presented in the Zwicky Transient Facility Catalogue of Periodic Variable Stars by \citet{Chen_2020}. We use a convolutional variational autoencoder to learn a low dimensional latent representation, and we search for anomalies within this latent dimension via an isolation forest. We identify anomalies with irregular variability. Most of the top anomalies are likely highly variable Red Giants or Asymptotic Giant Branch stars concentrated in the Milky Way galactic disk; a fraction of the identified anomalies are more consistent with Young Stellar Objects. Detailed spectroscopic follow-up observations are encouraged to reveal the nature of these anomalies., Comment: 26 pages, 17 figures. The full version of Table 4 and Table 5 are available upon request

V471 Tau is a post-common-envelope binary consisting of an eclipsing DA white dwarf and a K-type main-sequence star in the Hyades star cluster. We analyzed publicly available photometry and spectroscopy of V471 Tau to revise the stellar and orbital parameters of the system. We used archival K2 photometry, archival Hubble Space Telescope spectroscopy, and published radial-velocity measurements of the K-type star. Employing Gaussian processes to fit for rotational modulation of the system flux by the main-sequence star, we recovered the transits of the white dwarf in front of the main-sequence star for the first time. The transits are shallower than would be expected from purely geometric occultations owing to gravitational microlensing during transit, which places an additional constraint on the white-dwarf mass. Our revised mass and radius for the main-sequence star is consistent with single-star evolutionary models given the age and metallicity of the Hyades. However, as noted previously in the literature, the white dwarf is too massive and too hot to be the result of single-star evolution given the age of the Hyades, and may be the product of a merger scenario. We independently estimate the conditions of the system at the time of common envelope that would result in the measured orbital parameters today., Comment: 11 pages, 7 figures, Accepted to The Astronomical Journal

We present panchromatic observations and modeling of supernova (SN) 2020tlf, the first normal type II-P/L SN with confirmed precursor emission, as detected by the Young Supernova Experiment transient survey with the Pan-STARRS1 telescope. Pre-explosion emission was detected in $riz-$bands at 130 days prior to SN 2020tlf and persisted at relatively constant flux until first light. Soon after discovery, "flash" spectroscopy of SN 2020tlf revealed prominent narrow symmetric emission lines ($v_w < 300$ km s$^{-1}$) that resulted from the photo-ionization of unshocked circumstellar material (CSM) shedded in progenitor mass loss episodes in the final weeks to months before explosion. Surprisingly, this novel display of pre-SN emission and associated mass loss occurred in a RSG progenitor with ZAMS mass of only 10-12 M$_{\odot}$, as inferred from nebular spectra. Modeling of the light curve and multi-epoch spectra with the non-LTE radiative transfer code CMFGEN and radiation-hydrodynamical (RHD) code HERACLES suggests a dense CSM limited to $r \approx 10^{15}$ cm, and mass loss rate of $10^{-2}$ M$_{\odot}$ yr$^{-1}$. The subsequent luminous light-curve plateau and persistent blue excess indicates an extended progenitor, compatible with a RSG model with $R_{\star} = 1100$ R$_{\odot}$. Inferences from the limits on the shock-powered X-ray and radio luminosity are consistent with these conclusions and suggest a CSM density of $\rho < 2 \times 10^{-16}$ g cm$^{-3}$ for distances of $r \approx 5 \times 10^{15}$ cm, as well as a mass loss rate of $\dot M<1.3 \times 10^{-5}\,\rm{M_{\odot}\,yr^{-1}}$ at larger distances. A promising power source for the observed precursor emission is the ejection of stellar material following energy disposition into the stellar envelope as a result of gravity waves emitted during either neon/oxygen burning or a nuclear flash from silicon combustion., Comment: 36 pages, 20 figures. Accepted for publication in ApJ

APOGEE is a high-resolution (R sim 22,000), near-infrared, multi-epoch, spectroscopic survey of the Milky Way. The second generation of the APOGEE project, APOGEE-2, includes an expansion of the survey to the Southern Hemisphere called APOGEE-2S. This expansion enabled APOGEE to perform a fully panoramic mapping of all the main regions of the Milky Way; in particular, by operating in the H-band, APOGEE is uniquely able to probe the dust-hidden inner regions of the Milky Way that are best accessed from the Southern Hemisphere. In this paper we present the targeting strategy of APOGEE-2S, with special attention to documenting modifications to the original, previously published plan. The motivation for these changes is explained as well as an assessment of their effectiveness in achieving their intended scientific objective. In anticipation of this being the last paper detailing APOGEE targeting, we present an accounting of all such information complete through the end of the APOGEE-2S project; this includes several main survey programs dedicated to exploration of major stellar populations and regions of the Milky Way, as well as a full list of programs contributing to the APOGEE database through allocations of observing time by the Chilean National Time Allocation Committee (CNTAC) and the Carnegie Institution for Science (CIS). This work was presented along with a companion article, R. Beaton et al. (submitted; AAS29028), presenting the final target selection strategy adopted for APOGEE-2 in the Northern Hemisphere., Comment: 43 pages; 5 figures; 6 Tables; 1 Appendix; Submitted to Journal and Under Review; Posting to accompany papers using the SDSS-IV/APOGEE-2 Data Release 17 scheduled for December 2021

We present photometric and spectroscopic observations of Supernova 2020oi (SN 2020oi), a nearby ($\sim$17 Mpc) type-Ic supernova (SN Ic) within the grand-design spiral M100. We undertake a comprehensive analysis to characterize the evolution of SN 2020oi and constrain its progenitor system. We detect flux in excess of the fireball rise model $\delta t \approx 2.5$ days from the date of explosion in multi-band optical and UV photometry from the Las Cumbres Observatory and the Neil Gehrels Swift Observatory, respectively. The derived SN bolometric luminosity is consistent with an explosion with $M_{\rm ej} = 0.81 \pm 0.03 M_{\odot}$, $E_{k}= 0.79 \pm 0.09 \times 10^{51} \rm{erg} \rm{s}^{-1}$, and $M_{\rm Ni56} = 0.08 \pm 0.02 M_{\odot}$. Inspection of the event's decline reveals the highest $\Delta m_{15,\rm{bol}}$ reported for a stripped-envelope event to date. Modeling of optical spectra near event peak indicates a partially mixed ejecta comparable in composition to the ejecta observed in SN 1994I, while the earliest spectrum shows signatures of a possible interaction with material of a distinct composition surrounding the SN progenitor. Further, Hubble Space Telescope (HST) pre-explosion imaging reveals a stellar cluster coincident with the event. From the cluster photometry, we derive the mass and age of the SN progenitor using stellar evolution models implemented in the BPASS library. Our results indicate that SN 2020oi occurred in a binary system from a progenitor of mass $M_{\rm ZAMS} \approx 9.5 \pm 1.0 M_{\odot}$, corresponding to an age of $27 \pm 7$ Myr. SN 2020oi is the dimmest SN Ic event to date for which an early-time flux excess has been observed, and the first in which an early excess is unlikely to be associated with shock-cooling., Comment: 49 pages, 21 figures; accepted to ApJ

We report the early discovery and multi-color ($BVI$) high-cadence light curve analyses of a rapidly-declining sub-Chandrasekhar Type Ia supernova KSP-OT-201509b (= AT2015cx) from the KMTNet Supernova Program. The Phillips parameter and color stretch parameter of KSP-OT-201509b (= AT2015cx) are $\Delta$$M_{B,15}$ $\simeq$ 1.62 mag and $s_{BV}$ $\simeq$ 0.54, respectively, at an inferred redshift of 0.072. These, together with other measured parameters (such as the strength of the secondary $I$-band peak, colors and luminosity), identify the source to be a rapidly-declining Type Ia of transitional nature that is closer to Branch Normal than 91bg-like. Its early light curve evolution and bolometric luminosity are consistent with those of homologously expanding ejecta powered by radioactive decay and a Type Ia SN explosion with 0.32 $\pm$ 0.01 $M_{\odot}$ of synthesized $^{56}$Ni mass, 0.84 $\pm$ 0.12 $M_{\odot}$ of ejecta mass and (0.61 $\pm$ 0.14) $\times$ 10$^{51}$ erg of ejecta kinetic energy. While its $B-V$ and $V-I$ colors evolve largely synchronously with the changes in the $I$-band light curve as found in other supernovae, we also find the presence of an early redward evolution in $V-I$ prior to --10 days since peak. The bolometric light curve of the source is compatible with a stratified $^{56}$Ni distribution extended to shallow layers of the exploding progenitor. Comparisons between the observed light curves and those predicted from ejecta-companion interactions clearly disfavor Roche Lobe-filling companion stars at large separation distances, thus supporting a double-degenerate scenario for its origin. The lack of any apparent host galaxy in our deep stack images reaching a sensitivity limit of $\sim$ 28 $\rm mag\;arcsec^{-2}$ makes KSP-OT-201509b a hostless Type Ia supernova and offers new insights into supernova host galaxy environments., Comment: 39 pages, 20 figures, To appear in ApJ

We present the discovery of the fading radio transient FIRST J153350.8+272729. The source had a maximum observed 5-GHz radio luminosity of $8\times10^{39}$ erg s$^{-1}$ in 1986, but by 2019 had faded by a factor of nearly 400. It is located 0.15 arcsec from the center of a galaxy (SDSS J153350.89+272729) at 147 Mpc, which shows weak Type II Seyfert activity. We show that a tidal disruption event (TDE) is the preferred scenario for FIRST J153350.8+272729, although it could plausibly be interpreted as the afterglow of a long-duration gamma-ray burst. This is only the second TDE candidate to be first discovered at radio wavelengths. Its luminosity fills a gap between the radio afterglows of sub-relativistic TDEs in the local universe, and relativistic TDEs at high redshifts. The unusual properties of FIRST J153350.8+272729 (ongoing nuclear activity in the host galaxy, high radio luminosity) motivate more extensive TDE searches in untargeted radio surveys., Comment: 10 pages, 3 figures, 1 table, submitted to AAS Journals (ApJ)

Using {\it Gaia} Early Data Release 3 (EDR3) parallaxes and Bayesian inference, we infer a parallax of the Westerlund 1 (Wd1) cluster. We find a parallax of $0.34\pm{0.05}$ mas corresponding to a distance of $2.8^{+0.7}_{-0.6}$ kpc. The new {\it Gaia} EDR3 distance is consistent with our previous result using {\it Gaia} DR2 parallaxes. This confirms that Wd1 is less massive and older than previously assumed. Compared to DR2, the EDR3 individual parallax uncertainties for each star decreased by 30\%. However, the aggregate parallax uncertainty for the cluster remained the same. This suggests that the uncertainty is dominated by systematics, which is possibly due to crowding, motions within the cluster, or motions due to binary orbits., Comment: 3 pages, 1 figure

We present new proper motion measurements of optically emitting oxygen-rich knots of supernova remnant 1E 0102.2-7219 (E0102), which are used to estimate the remnant's center of expansion and age. Four epochs of high resolution Hubble Space Telescope images spanning 19 yr were retrieved and analyzed. We found a robust center of expansion of alpha=1:04:02.48 and delta=-72:01:53.92 (J2000) with 1-sigma uncertainty of 1.77 arcseconds using 45 knots from images obtained with the Advanced Camera for Surveys using the F475W filter in 2003 and 2013 having the highest signal-to-noise ratio. We also estimate an upper limit explosion age of 1738 +/- 175 yr by selecting knots with the highest proper motions, that are assumed to be the least decelerated. We find evidence of an asymmetry in the proper motions of the knots as a function of position angle. We conclude that these asymmetries were most likely caused by interaction between E0102's original supernova blast wave and an inhomogeneous surrounding environment, as opposed to intrinsic explosion asymmetry. The observed non-homologous expansion suggests that the use of a free expansion model inaccurately offsets the center of expansion and leads to an overestimated explosion age. We discuss our findings as they compare to previous age and center of expansion estimates of E0102 and their relevance to a recently identified candidate central compact object., Comment: 14 pages, 11 figures, revised according to referee comments and resubmitted to ApJ

We present Hubble Space Telescope imaging of a pre-explosion counterpart to SN 2019yvr obtained 2.6 years before its explosion as a type Ib supernova (SN Ib). Aligning to a post-explosion Gemini-S/GSAOI image, we demonstrate that there is a single source consistent with being the SN 2019yvr progenitor system, the second SN Ib progenitor candidate after iPTF13bvn. We also analyzed pre-explosion Spitzer/IRAC imaging, but we do not detect any counterparts at the SN location. SN 2019yvr was highly reddened, and comparing its spectra and photometry to those of other, less extinguished SNe Ib we derive $E(B-V)=0.51\substack{+0.27\\-0.16}$ mag for SN 2019yvr. Correcting photometry of the pre-explosion source for dust reddening, we determine that this source is consistent with a $\log(L/L_{\odot}) = 5.3 \pm 0.2$ and $T_{\mathrm{eff}} = 6800\substack{+400\\-200}$ K star. This relatively cool photospheric temperature implies a radius of 320$\substack{+30\\-50} R_{\odot}$, much larger than expectations for SN Ib progenitor stars with trace amounts of hydrogen but in agreement with previously identified SN IIb progenitor systems. The photometry of the system is also consistent with binary star models that undergo common envelope evolution, leading to a primary star hydrogen envelope mass that is mostly depleted but seemingly in conflict with the SN Ib classification of SN 2019yvr. SN 2019yvr had signatures of strong circumstellar interaction in late-time ($>$150 day) spectra and imaging, and so we consider eruptive mass loss and common envelope evolution scenarios that explain the SN Ib spectroscopic class, pre-explosion counterpart, and dense circumstellar material. We also hypothesize that the apparent inflation could be caused by a quasi-photosphere formed in an extended, low-density envelope or circumstellar matter around the primary star., Comment: 22 pages, 9 figures, submitted to MNRAS

We identify red supergiants (RSGs) in our spiral neighbors M31 and M33 using near-IR (NIR) photometry complete to a luminosity limit of log L/Lo=4.0. Our archival survey data cover 5 deg^2 of M31, and 3 deg^2 for M33, and are likely spatially complete for these massive stars. Gaia is used to remove foreground stars, after which the RSGs can be separated from asymptotic giant branch (AGB) stars in the color-magnitude diagram. The photometry is used to derive effective temperatures and bolometric luminosities via MARCS stellar atmosphere models. The resulting H-R diagrams show superb agreement with the evolutionary tracks of the Geneva evolutionary group. Our census includes 6400 RSGs in M31 and 2850 RSGs in M33 within their Holmberg radii; by contrast, only a few hundred RSGs are known so far in the Milky Way. Our catalog serves as the basis for a study of the RSG binary frequency being published separately, as well as future studies relating to the evolution of massive stars. Here we use the matches between the NIR-selected RSGs and their optical counterparts to show that the apparent similarity in the reddening of OB stars in M31 and M33 is the result of Malmquist bias; the average extinction in M31 is likely higher than that of M33. As expected, the distribution of RSGs follows that of the spiral arms, while the much older AGB population is more uniformly spread across each galaxy's disk., Comment: Accepted for publication in the Astronomical Journal. Revision is to fix minor typo in Table 2, where M31 and M33 temperature equations had been interchanged

It is unclear if neutron star mergers can explain the observed r-process abundances of metal-poor stars. Collapsars, defined here as rotating massive stars whose collapse results in a rapidly accreting disk around a black hole that can launch jets, are a promising alternative. We find that we can produce a self-consistent model in which a population of collapsars with stochastic europium yields synthesizes all of the r-process material in metal-poor ([Fe/H] < -2.5) stars. Our model reproduces the observed scatter and evolution of scatter of [Eu/Fe] abundances. We find that if collapsars are the dominant r-process site for metal-poor stars, r-process synthesis may be linked to supernovae that produce long gamma-ray bursts. Our results also allow for the possibility that core-collapse supernovae beyond those that launch gamma-ray bursts also produce r-process material (e.g., potentially a subset of Type Ic-BL supernovae). Furthermore, we identify collapsar jet properties (isotropic energy, engine luminosity, or engine time) which may trace r-process yield and verify that the amount of r-process yield produced per collapsar in our model (~0.07 Msun) is consistent with other independent estimates. In the future, achieving 0.05 dex precision on distribution scatter or a reliable selection function would further constrain our probe of r-process production. Our model would also hold for another prompt r-process site with a power-law yield, and work is needed to determine if, for example, fast-merging neutron stars can also explain abundance scatter., Comment: 17 pages, 8 figures. Accepted by ApJ

We perform a systematic study of the $^{56}$Ni mass ($M_{\rm Ni}$) of 27 stripped envelope supernovae (SESNe) by modeling their light-curve tails, highlighting that use of ``Arnett's rule'' overestimates $M_{\rm Ni}$ for SESN by a factor of $\sim$2. Recently, \citet{Khatami2019} presented a new model relating the peak time ($t_{\rm p}$) and luminosity ($L_{\rm p}$) of a radioactive-powered SN to its $M_{\rm Ni}$ that addresses several limitations of Arnett-like models, but depends on a dimensionless parameter, $\beta$. Using observed $t_{\rm p}$, $L_{\rm p}$, and tail-measured $M_{\rm Ni}$ values for 27 SESN, we observationally calibrate $\beta$ for the first time. Despite scatter, we demonstrate that the model of \citet{Khatami2019} with empirically-calibrated $\beta$ values provides significantly improved measurements of $M_{\rm Ni}$ when only photospheric data is available. However, these observationally-constrained $\beta$ values are systematically lower than those inferred from numerical simulations, primarily because the observed sample has significantly higher (0.2-0.4 dex) $L_{\rm p}$ for a given $M_{\rm Ni}$. While effects due to composition, mixing, and asymmetry can increase $L_{\rm p}$ current models cannot explain the systematically low $\beta$ values. However, the discrepancy can be alleviated if $\sim$7--50\% of $L_{\rm p}$ for the observed sample originates from sources other than $^{56}$Ni. Either shock cooling or magnetar spin-down could provide the requisite luminosity. Finally, we find that even with our improved measurements, the $M_{\rm Ni}$ values of SESN are still a factor of $\sim$3 larger than those of hydrogen-rich Type II SN, indicating that these supernovae are inherently different in terms of their progenitor initial mass distributions or explosion mechanisms., Comment: 27 pages, 12 figures, accepted for publication in ApJ

Stars with unusual properties can provide a wealth of information about rare stages of stellar evolution and exotic physics. However, determining the true nature of peculiar stars is often difficult. In this work, we conduct a systematic search for cool and luminous stars in the Magellanic Clouds with extreme variability, motivated by the properties of the unusual SMC star and Thorne-\.Zytkow Object (T\.ZO) candidate HV2112. Using light curves from ASAS-SN we identify 38 stars with surface temperatures T $<$ 4800K, luminosities $\log$(L/L$_\odot$) $>$ 4.3, variability periods $>$ 400 days, and variability amplitudes $\Delta$V $>$ 2.5 mag. Eleven of these stars possess the distinctive double-peaked light curve morphology of HV2112. We use the pulsation properties and derived occurrence rates for these 12 objects to constrain their nature. From comparisons to stellar populations and models, we find that one star may be a red supergiant with large amplitude pulsations. For the other 11 stars we derive current masses of $\sim$5-10 M$_{\odot}$, below the theoretical minimum mass of $\sim$15 M$_{\odot}$ for T\.ZOs to be stable, casting doubt on this interpretation. Instead, we find that the temperatures, luminosities, mass-loss rates, and periods of these stars are consistent with predictions for super-Asymptotic Giant Branch (s-AGB) stars that have begun carbon burning but have not reached the superwind phase. We infer lifetimes in this phase of $\sim($1$-$7) $\times$ 10$^{4}$ years, also consistent with an s-AGB interpretation. If confirmed, these objects would represent the first identified population of s-AGB stars, illuminating the transition between low- and high-mass stellar evolution., Comment: 29 pages, 14 figures, 4 tables, plus Appendix. Accepted for publication in the Astrophysical Journal

The classification of supernovae (SNe) and its impact on our understanding of the explosion physics and progenitors have traditionally been based on the presence or absence of certain spectral features. However, current and upcoming wide-field time-domain surveys have increased the transient discovery rate far beyond our capacity to obtain even a single spectrum of each new event. We must therefore rely heavily on photometric classification, connecting SN light curves back to their spectroscopically defined classes. Here we present Superphot, an open-source Python implementation of the machine-learning classification algorithm of Villar et al., and apply it to 2315 previously unclassified transients from the Pan-STARRS1 Medium Deep Survey for which we obtained spectroscopic host-galaxy redshifts. Our classifier achieves an overall accuracy of 82%, with completenesses and purities of >80% for the best classes (SNe Ia and superluminous SNe). For the worst performing SN class (SNe Ibc), the completeness and purity fall to 37% and 21%, respectively. Our classifier provides 1257 newly classified SNe Ia, 521 SNe II, 298 SNe Ibc, 181 SNe IIn, and 58 SLSNe. These are among the largest uniformly observed samples of SNe available in the literature and will enable a wide range of statistical studies of each class., Comment: Updated to match published version in ApJ. Companion paper to Villar et al. (2020). Full machine-readable tables available as ancillary files at right

Automated classification of supernovae (SNe) based on optical photometric light curve information is essential in the upcoming era of wide-field time domain surveys, such as the Legacy Survey of Space and Time (LSST) conducted by the Rubin Observatory. Photometric classification can enable real-time identification of interesting events for extended multi-wavelength follow-up, as well as archival population studies. Here we present the complete sample of 5,243 "SN-like" light curves (in griz) from the Pan-STARRS1 Medium-Deep Survey (PS1-MDS). The PS1-MDS is similar to the planned LSST Wide-Fast-Deep survey in terms of cadence, filters and depth, making this a useful training set for the community. Using this dataset, we train a novel semi-supervised machine learning algorithm to photometrically classify 2,315 new SN-like light curves with host galaxy spectroscopic redshifts. Our algorithm consists of a random forest supervised classification step and a novel unsupervised step in which we introduce a recurrent autoencoder neural network (RAENN). Our final pipeline, dubbed SuperRAENN, has an accuracy of 87% across five SN classes (Type Ia, Ibc, II, IIn, SLSN-I). We find the highest accuracy rates for Type Ia SNe and SLSNe and the lowest for Type Ibc SNe. Our complete spectroscopically- and photometrically-classified samples break down into: 62.0% Type Ia (1839 objects), 19.8% Type II (553 objects), 4.8% Type IIn (136 objects), 11.7% Type Ibc (291 objects), and 1.6% Type I SLSNe (54 objects). Finally, we discuss how this algorithm can be modified for online LSST data streams., Comment: Submitted to ApJ; Companion paper to Hosseinzadeh et al.; Tables 1 and 2 available on ashleyvillar.com prior to publication

The binary fraction of unevolved massive stars is thought to be 70-100% but there are few observational constraints on the binary fraction of the evolved version of a subset of these stars, the red supergiants (RSGs). Here we identify a complete sample of RSGs in the Large Magellanic Cloud (LMC) using new spectroscopic observations and archival UV, IR and broadband optical photometry. We find 4090 RSGs with log L/Lo > 3.5 with 1820 of them having log L/Lo > 4, which we believe is our completeness limit. We additionally spectroscopically confirmed 38 new RSG+B star binaries in the LMC, bringing the total known up to 55. We then estimated the binary fraction using a k-nearest neighbors algorithm that classifies stars as single or binary based on photometry with a spectroscopic sample as a training set. We take into account observational biases such as line-of-sight stars and binaries in eclipse while also calculating model-dependent corrections for RSGs with companions that our observations were not designed to detect. Based on our data, we find an initial result of 13.5 +7.56/-6.67% for RSGs with O or B-type companions. Using the Binary Population and Spectral Synthesis (BPASS) models to correct for unobserved systems, this corresponds to a total RSG binary fraction of 19.5 +7.6/-6.7%. This number is in broad agreement with what we would expect given an initial OB binary distribution of 70%, a predicted merger fraction of 20-30% and a binary interaction fraction of 40-50%., Comment: ApJ, accepted

This compilation is the fourth data release from the $R$-Process Alliance (RPA) search for $r$-process-enhanced stars, and the second release based on "snapshot" high-resolution ($R \sim 30,000$) spectra collected with the du Pont 2.5m Telescope. In this data release, we propose a new delineation between the $r$-I and $r$-II stellar classes at $\mathrm{[Eu/Fe]} = +0.7$, instead of the empirically chosen $\mathrm{[Eu/Fe]} = +1.0$ level previously in use, based on statistical tests of the complete set of RPA data released to date. We also statistically justify the minimum level of [Eu/Fe] for definition of the $r$-I stars, [Eu/Fe] $> +0.3$. Redefining the separation between $r$-I and $r$-II stars will aid in analysis of the possible progenitors of these two classes of stars and whether these signatures arise from separate astrophysical sources at all. Applying this redefinition to previous RPA data, the number of identified $r$-II and $r$-I stars changes to 51 and 121, respectively, from the initial set of data releases published thus far. In this data release, we identify 21 new $r$-II, 111 new $r$-I (plus three re-identified), and 7 new (plus one re-identified) limited-$r$ stars out of a total of 232 target stars, resulting in a total sample of 72 new $r$-II stars, 232 new $r$-I stars, and 42 new limited-$r$ stars identified by the RPA to date., Comment: 17 pages, 7 figures, 6 tables; accepted to ApJS

We report results of a novel high-energy follow-up observation of a potential Fast Radio Burst. The radio burst was detected by VLA/realfast and followed-up by the Neil Gehrels Swift Observatory in very low latency utilizing new operational capabilities of Swift (arXiv:2005.01751), with pointed soft X-ray and UV observations beginning at T0+32 minutes, and hard X-ray/gamma-ray event data saved around T0. These observations are $>10$x faster than previous X-ray/UV follow-up of any radio transient to date. No emission is seen coincident with the FRB candidate at T0, with a 0.2s fluence $5\sigma$ upper limit of $1.35\times10^{-8}$ erg cm$^{-2}$ (14-195 keV) for a SGR 1935+2154-like flare, nor at T0+32 minutes down to $3\sigma$ upper limits of 22.18 AB mag in UVOT u band, and $3.33\times10^{-13}$ erg cm$^{-2}$ s$^{-1}$ from 0.3-10 keV for the 2 ks observation. The candidate FRB alone is not significant enough to be considered astrophysical, so this note serves as a technical demonstration. These new Swift operational capabilities will allow future FRB detections to be followed up with Swift at even lower latencies than demonstrated here: 15-20 minutes should be regularly achievable, and 5-10 minutes occasionally achievable. We encourage FRB detecting facilities to release alerts in low latency to enable this science., Comment: Technical note and capability update for the community. We encourage low latency FRB alerts from relevant facilities to enable this science

We present a wide-field optical imaging search for electromagnetic counterparts to the likely neutron star - black hole (NS-BH) merger GW190814/S190814bv. This compact binary merger was detected through gravitational waves by the LIGO/Virgo interferometers, with masses suggestive of a NS-BH merger. We imaged the LIGO/Virgo localization region using the MegaCam instrument on the Canada-France-Hawaii Telescope. We describe our hybrid observing strategy of both tiling and galaxy-targeted observations, as well as our image differencing and transient detection pipeline. Our observing campaign produced some of the deepest multi-band images of the region between 1.7 and 8.7 days post-merger, reaching a 5sigma depth of g > 22.8 (AB mag) at 1.7 days and i > 23.1 and i > 23.9 at 3.7 and 8.7 days, respectively. These observations cover a mean total integrated probability of 67.0% of the localization region. We find no compelling candidate transient counterparts to this merger in our images, which suggests that either the lighter object was tidally disrupted inside of the BH's innermost stable circular orbit, the transient lies outside of the observed sky footprint, or the lighter object is a low-mass BH. We use 5sigma source detection upper limits from our images in the NS-BH interpretation of this merger to constrain the mass of the kilonova ejecta to be Mej < 0.015Msun for a 'blue' (kappa = 0.5 cm^2 g^-1) kilonova, and Mej < 0.04Msun for a 'red' (kappa = 5-10 cm^2 g^-1) kilonova. Our observations emphasize the key role of large-aperture telescopes and wide-field imagers such as CFHT MegaCam in enabling deep searches for electromagnetic counterparts to gravitational wave events., Comment: 23 pages, 13 figures including appendices; accepted to ApJ with minor revisions

We have imaged over 35 years of archival Very Large Array (VLA) observations of the nearby (d$_{\rm{L}}$ $=$ 3.15 Mpc) Type Ia supernovae SN\,1972E and SN\,1895B between 9 and 121 years post-explosion. No radio emission is detected, constraining the 8.5 GHz luminosities of SN\,1972E and SN\,1895B to be L$_{\nu,8.5\rm{GHz}}$ $<$ 6.0 $\times$ 10$^{23}$ erg s$^{-1}$ Hz$^{-1}$ 45 years post-explosion and L$_{\nu,8.5\rm{GHz}}$ $<$ 8.9 $\times$ 10$^{23}$ erg s$^{-1}$ Hz$^{-1}$ 121 years post-explosion, respectively. These limits imply a clean circumstellar medium (CSM), with $n$ $<$ 0.9 cm$^{-3}$ out to radii of a few $\times$ 10$^{18}$ cm, if the SN blastwave is expanding into uniform density material. Due to the extensive time coverage of our observations, we also constrain the presence of CSM shells surrounding the progenitor of SN\,1972E. We rule out essentially all medium and thick shells with masses of 0.05$-$0.3 M$_\odot$ at radii between $\sim$10$^{17}$ and 10$^{18}$ cm, and thin shells at specific radii with masses down to $\lesssim$0.01 M$_\odot$. These constraints rule out swaths of parameter space for a range of single and double degenerate progenitor scenarios, including recurrent nova, core-degenerate objects, ultra-prompt explosions and white dwarf (WD) mergers with delays of a few hundred years between the onset of merger and explosion. Allowed progenitors include WD-WD systems with a significant ($>$ 10$^{4}$ years) delay from the last episode of common envelope evolution and single degenerate systems undergoing recurrent nova, provided that the recurrence timescale i short and the system has been in the nova phase for $\gtrsim$10$^{4}$ yr, such that a large ($>$ 10$^{18}$ cm) cavity has been evacuated. Future multi-epoch observations of additional intermediate-aged Type Ia SNe will provide a comprehensive view of the large-scale CSM environments around these explosions., Comment: Accepted for publication in the Astrophysical Journal

The mass-loss rates of red supergiant stars (RSGs) are poorly constrained by direct measurements, and yet the subsequent evolution of these stars depends critically on how much mass is lost during the RSG phase. In 2012 the Geneva evolutionary group updated their mass-loss prescription for RSGs with the result that a 20 solar mass star now loses 10x more mass during the RSG phase than in the older models. Thus, higher mass RSGs evolve back through a second yellow supergiant phase rather than exploding as Type II-P supernovae, in accord with recent observations (the so-called "RSG Problem"). Still, even much larger mass-loss rates during the RSG phase cannot be ruled out by direct measurements of their current dust-production rates, as such mass-loss is episodic. Here we test the models by deriving a luminosity function for RSGs in the nearby spiral galaxy M31 which is sensitive to the total mass loss during the RSG phase. We carefully separate RSGs from asymptotic giant branch stars in the color-magnitude diagram following the recent method exploited by Yang and collaborators in their Small Magellanic Cloud studies. Comparing our resulting luminosity function to that predicted by the evolutionary models shows that the new prescription for RSG mass-loss does an excellent job of matching the observations, and we can readily rule out significantly larger values., Comment: ApJ, accepted

The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear astrophysics, dating back to the work of Alastair Cameron in the 1950s. Recent faculty hires have further boosted activity in the field, including transient observation and theory, survey science on galactic nucleosynthesis, and nuclear experiments. This white paper contains a brief overview of recent activity in the community, highlighting strengths in each sub-field, and provides recommendations to improve interdisciplinary collaboration. Sustaining Canadian leadership in the next decade will require, on the observational side, access to transient and non-transient surveys like LSST, SKA, or MSE, support for target-of-opportunity observing in current and future Canadian telescopes, and participation in next-generation X-ray telescopes such as ATHENA. State-of-the-art theoretical predictions will require an ambitious succession plan for the Niagara supercomputer to support large parallel jobs. We propose a funding instrument for postdoctoral training that reflects the interdisciplinary nature of nuclear astrophysics research, and the creation of a national collaborative funding program that allows for joint projects and workshop organization., Comment: White paper submitted to the Canadian Long Range Plan 2020. Minor formatting changes relative to submitted version

The relationship between galaxies and the state/chemical enrichment of the warm-hot intergalactic medium (WHIM) expected to dominate the baryon budget at low-z provides sensitive constraints on structure formation and galaxy evolution models. We present a deep redshift survey in the field of 1ES1553+113, a blazar with a unique combination of UV+X-ray spectra for surveys of the circum-/intergalactic medium (CGM/IGM). Nicastro et al. 2018 reported the detection of two O VII WHIM absorbers at $z=0.4339$ and $0.3551$ in its spectrum, suggesting that the WHIM is metal-rich and sufficient to close the missing baryons problem. Our survey indicates that the blazar is a member of a $z=0.433$ group and that the higher-$z$ O VII candidate arises from its intragroup medium. The resulting bias precludes its use in baryon censuses. The $z=0.3551$ candidate occurs in an isolated environment 630 kpc from the nearest galaxy (with stellar mass $\log M_*/M_\odot \approx 9.7$) which we show is unexpected for the WHIM. Finally, we characterize the galactic environments of broad H I Ly$\alpha$ absorbers (Doppler widths of $b=40-80$ \kms; $T\lesssim4\times10^5$ K) which provide metallicity independent WHIM probes. On average, broad Ly$\alpha$, absorbers are ${\approx}2\times$ closer to the nearest luminous ($L>0.25 L_*$) galaxy (700 kpc) than narrow ($b<30$ \kms; $T\lesssim4\times10^5$ K) ones (1300 kpc) but ${\approx}2\times$ further than O\,VI absorbers (350 kpc). These observations suggest that gravitational collapse heats portions of the IGM to form the WHIM but with feedback that does not enrich the IGM far beyond galaxy/group halos to levels currently observable in UV/X-ray metal lines., Comment: Accepted to ApJ Letters; 10 pages, 4 figures, 1 table with the full version available on the journal webpage (or by request); v2 corrects one reference

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., Comment: 20 pages + Appendix, accepted to ApJ

We present the discovery and the photometric and spectroscopic study of H-rich Type II supernova (SN) KSP-SN-2016kf (SN2017it) observed in the KMTNet Supernova Program in the outskirts of a small irregular galaxy at $z\simeq0.043$ within a day from the explosion. Our high-cadence, multi-color ($BVI$) light curves of the SN show that it has a very long rise time ($t_\text{rise}\simeq 20$ days in $V$ band), a moderately luminous peak ($M_V\simeq -$17.6 mag), a notably luminous and flat plateau ($M_V\simeq -$17.4 mag and decay slope $s\simeq0.53$ mag per 100 days), and an exceptionally bright radioactive tail. Using the color-dependent bolometric correction to the light curves, we estimate the $^{56}$Ni mass powering the observed radioactive tail to be $0.10\pm0.01$ M$_\odot$, making it a H-rich Type II SN with one of the largest $^{56}$Ni masses observed to date. The results of our hydrodynamic simulations of the light curves constrain the mass and radius of the progenitor at the explosion to be $\sim$15 M$_\odot$ (evolved from a star with an initial mass of $\sim$ 18.8 M$_\odot$) and $\sim1040$ R$_\odot$, respectively, with the SN explosion energy of $\sim 1.3\times 10^{51}$ erg s$^{-1}$. The above-average mass of the KSP-SN-2016kf progenitor, together with its low metallicity $ Z/Z_\odot \simeq0.1-0.4$ obtained from spectroscopic analysis, is indicative of a link between the explosion of high-mass red supergiants and their low-metallicity environment. The early part of the observed light curves shows the presence of excess emission above what is predicted in model calculations, suggesting there is interaction between the ejecta and circumstellar material. We further discuss the implications of the high progenitor initial mass and low-metallicity environment of KSP-SN-2016kf on our understanding of the origin of Type II SNe., Comment: 18 pages, 12 figures, accepted for publication in ApJ

Multimessenger observations of the neutron star merger GW170817 and its kilonova proved that neutron star mergers can synthesize large quantities of r-process elements. If neutron star mergers in fact dominate all r-process element production, then the distribution of kilonova ejecta compositions should match the distribution of r-process abundance patterns observed in stars. The lanthanide fraction (XLa) is a measurable quantity in both kilonovae and metal-poor stars, but it has not previously been explicitly calculated for stars. Here, we compute the lanthanide fraction distribution of metal-poor stars ([Fe/H] < -2.5) to enable comparison to current and future kilonovae. The full distribution peaks at log XLa ~ -1.8, but r-process-enhanced stars ([Eu/Fe] > 0.7) have distinctly higher lanthanide fractions; log XLa > -1.5. We review observations of GW170817 and find general consensus that the total log XLa = -2.2 +/- 0.5, somewhat lower than the typical metal-poor star and inconsistent with the most highly r-enhanced stars. For neutron star mergers to remain viable as the dominant r-process site, future kilonova observations should be preferentially lanthanide-rich (including a population of ~10% with log XLa > -1.5). These high-XLa kilonovae may be fainter and more rapidly evolving than GW170817, posing a challenge for discovery and follow-up observations. Both optical and (mid-)infrared observations will be required to robustly constrain kilonova lanthanide fractions. If such high-XLa kilonovae are not found in the next few years, that likely implies that the stars with the highest r-process enhancements have a different origin for their r-process elements., Comment: ApJ accepted after minor revisions. 18 pages + appendices, 8 figures, 5 tables. See Figure 5 for summary

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network., Comment: White paper submitted to the Astro2020 Decadal Survey

Westerlund 1 (Wd1) is potentially the largest star cluster in the Galaxy. That designation critically depends upon the distance to the cluster, yet the cluster is highly obscured, making luminosity-based distance estimates difficult. Using {\it Gaia} Data Release 2 (DR2) parallaxes and Bayesian inference, we infer a parallax of $0.35^{+0.07}_{-0.06}$ mas corresponding to a distance of $2.6^{+0.6}_{-0.4}$ kpc. To leverage the combined statistics of all stars in the direction of Wd1, we derive the Bayesian model for a cluster of stars hidden among Galactic field stars; this model includes the parallax zero-point. Previous estimates for the distance to Wd1 ranged from 1.0 to 5.5 kpc, although values around 5 kpc have usually been adopted. The {\it Gaia} DR2 parallaxes reduce the uncertainty from a factor of 3 to 18\% and rules out the most often quoted value of 5 kpc with 99\% confidence. This new distance allows for more accurate mass and age determinations for the stars in Wd1. For example, the previously inferred initial mass at the main-sequence turn-off was around 40 M$_{\odot}$; the new {\it Gaia} DR2 distance shifts this down to about 22 M$_{\odot}$. This has important implications for our understanding of the late stages of stellar evolution, including the initial mass of the magnetar and the LBV in Wd1. Similarly, the new distance suggests that the total cluster mass is about four times lower than previously calculated., Comment: 14 pages, 10 figures