Recent observations presented in Cendes et al. (2024a) show that optically selected tidal disruption events (TDEs) commonly produce delayed radio emission that can peak years post-disruption. Here, we explore the multi-wavelength properties of a sample of radio-observed optically selected TDEs to shed light on the physical process(es) responsible for the late-rising radio emission. We combine new late-time X-ray observations with archival optical, UV, X-ray, and radio data to conclude that a diversity of accretion-driven outflows may power delayed radio emission in TDEs. Our analysis suggests that some late radio outflows may be launched by a delayed phase of super-Eddington accretion onto the central supermassive black hole (SMBH), while others may result from a state transition to a “low-hard” radiatively inefficient accretion flow or the deceleration of an off-axis relativistic jet. We find that TDEs with delayed radio emission are less likely to exhibit helium emission lines at early times (p=0.002) and may have larger optical/UV photospheric radii (p=0.026) than other TDEs, possibly also indicating that the onset of SMBH accretion is delayed in these systems. Our results have implications for our understanding of state changes in SMBH accretion flows, the circularization timescale for TDE debris, and the prevalence of off-axis jets in TDEs, and motivate systematic, long-term monitoring of these unique transients. The objects in our sample with the brightest radio emission are also detected in the VLA Sky Survey (VLASS), demonstrating that all-sky radio surveys can play an important role in discovering unexpected properties of the TDE population.

We present multiband ATLAS and ZTF photometry for SN 2019tsf, a Type Ib stripped-envelope supernova (SESN). The slow spectral evolution could be associated with an uncommon explosion mechanism specific to this SN. Possible explanations include fallback accretion onto a compact remnant or a long-lived central engine, both of which could provide extended energy injection responsible for the late-time rebrightening and unusual spectral features. The rebrightening observations represent the latest photometric measurements of a multipeaked Type Ib SN. As late-time photometry and spectroscopy suggest no hydrogen, the potential circumstellar material (CSM) must be H-poor. The absence of a nebular phase and the lack of narrow emission lines in the late-time spectra (>142 days) of the SNe suggest that any CSM interaction is likely asymmetric and enveloped by the SN ejecta. However, an extended CSM structure is evident through a follow-up radio campaign with the Karl G. Jansky Very Large Array (VLA), indicating a source of bright optically thick radio emission at late times, which is highly unusual among H-poor SESNe. We attribute this phenomenology to an interaction of the supernova ejecta with asymmetric CSM, potentially disk-like, and we present several models that may explain the origin of this rare Type Ib supernova. We propose a warped disk model in which a tertiary companion—commonly present around massive stars—perturbs the progenitor’s CSM, producing density enhancements that may explain the observed multipeaked SN 2019tsf light curve. This SN 2019tsf is a unique SN Type Ib among the recently discovered class of SNe that undergo mass transfer at the moment of explosion.

The Legacy Survey of Space and Time (LSST) at Vera C. Rubin Observatory is planned to begin in the Fall of 2025. The LSST survey cadence has been designed via a community-driven process regulated by the Survey Cadence Optimization Committee (SCOC), which recommended up to 3 Target of Opportunity (ToO) observations. Experts from the scientific community, Rubin Observatory personnel, and members of the SCOC were brought together to deliver a recommendation for the implementation of the ToO program during a workshop held in March 2024. Four main science cases were identified: gravitational wave multi-messenger astronomy, high energy neutrinos, Galactic supernovae, and small potentially hazardous asteroids possible impactors. Additional science cases were identified and briefly addressed in the documents, including lensed or poorly localized gamma-ray bursts and twilight discoveries. Trigger prioritization, automated response, and detailed strategies were discussed for each science case. This document represents the outcome of the Rubin ToO 2024 workshop, with additional contributions from members of the Rubin Science Collaborations. The implementation of the selection criteria and strategies presented in this document has been endorsed in the SCOC Phase 3 Recommendations document (PSTN-056). Although the ToO program is still to be finalized, this document serves as a baseline plan for ToO observations with the Rubin Observatory.

We introduce a novel method that utilizes the longitude-velocity (l - v) envelope to constrain the Milky Way (MW) bar potential. Previous work used the l - v diagram to explain the distribution of the observed high-velocity stars. We successfully reproduce their results but find that their method is limited to only one type of periodic orbits. In contrast, we propose that the l - v envelope provides much more comprehensive constraints. We compare the properties of test particles in the M. Portail et al. MW potential model (P17) with the observed SiO maser stars from the Bulge Asymmetries and Dynamical Evolution (BAaDE) survey. We find that the l - v envelope generated by the bar potential demonstrates reasonable agreement with the observational data, albeit with slight discrepancies near the Galactic center. The inconsistencies suggest that the P17 potential yields a lower central rotation curve, a slightly larger quadrupole strength, or a possibly underestimated pattern speed. We also adopt an updated version of the P17 potential with a modified central mass component proposed by G. H. Hunter et al. (H24). The fitting of the l - v envelope suggests that the H24 potential does not completely address the existing challenges and may hint at a possible underestimation of the central bar mass. Our study demonstrates that the l - v envelope can be used as a valuable tool for constraining the Galactic potential and provides insights into the MW bar potential.

We present a method to estimate distances to asymptotic giant branch (AGB) stars in the Galaxy, using spectral energy distributions (SEDs) in the near- and mid-infrared. By assuming that a given set of source properties (initial mass, stellar temperature, composition, and evolutionary stage) will provide a typical SED shape and brightness, sources are color matched to a distance-calibrated template and thereafter scaled to extract the distance. The method is tested by comparing the distances obtained to those estimated from very long baseline interferometry or Gaia parallax measurements, yielding a strong correlation in both cases. Additional templates are formed by constructing a source sample likely to be close to the Galactic center, and thus with a common, typical distance for calibration of the templates. These first results provide statistical distance estimates to a set of almost 15,000 Milky Way AGB stars belonging to the Bulge Asymmetries and Dynamical Evolution (BAaDE) survey, with typical distance errors of ±35%. With these statistical distances, a map of the intermediate-age population of stars traced by AGBs is formed, and a clear bar structure can be discerned, consistent with the previously reported inclination angle of 30° to the GC–Sun direction vector. These results motivate deeper studies of the AGB population to tease out the intermediate-age stellar distribution throughout the Galaxy, as well as determining statistical properties of the AGB population luminosity and mass-loss-rate distributions.

The Vera C. Rubin Legacy Survey of Space and Time (LSST) holds the potential to revolutionize time domain astrophysics, reaching completely unexplored areas of the Universe and mapping variability time scales from minutes to a decade. To prepare to maximize the potential of the Rubin LSST data for the exploration of the transient and variable Universe, one of the four pillars of Rubin LSST science, the Transient and Variable Stars Science Collaboration, one of the eight Rubin LSST Science Collaborations, has identified research areas of interest and requirements, and paths to enable them. While our roadmap is ever-evolving, this document represents a snapshot of our plans and preparatory work in the final years and months leading up to the survey’s first light.

We present multi-band ATLAS photometry for SN 2019tsf, a stripped-envelope Type Ib supernova (SESN). The SN shows a triple-peaked light curve and a late (re-)brightening, making it unique among stripped-envelope systems. The re-brightening observations represent the latest photometric measurements of a multi-peaked Type Ib SN to date. As late-time photometry and spectroscopy suggest no hydrogen, the potential circumstellar material (CSM) must be H-poor. Moreover, late (>150 days) spectra show no signs of narrow emission lines, further disfavouring CSM interaction. On the contrary, an extended CSM structure is seen through a follow-up radio campaign with Karl G. Jansky Very Large Array (VLA), indicating a source of bright optically thick radio emission at late times, which is highly unusual among H-poor SESNe. We attribute this phenomenology to an interaction of the supernova ejecta with spherically-asymmetric CSM, potentially disk-like, and we present several models that can potentially explain the origin of this rare Type Ib supernova. The warped disc model paints a novel picture, where the tertiary companion perturbs the progenitors CSM, that can explain the multi-peaked light curves of SNe, and here we apply it to SN 2019tsf. This SN 2019tsf is likely a member of a new sub-class of Type Ib SNe and among the recently discovered class of SNe that undergo mass transfer at the moment of explosion

We present the results from our 7 yr long broadband X-ray observing campaign of SN 2014C with Chandra and NuSTAR. These coordinated observations represent the first look at the evolution of a young extragalactic SN in the 0.3–80 keV energy range in the years after core collapse. We find that the spectroscopic metamorphosis of SN 2014C from an ordinary type Ib SN into an interacting SN with copious hydrogen emission is accompanied by luminous X-rays reaching L x ≈ 5.6 × 10 40 erg s −1 (0.3–100 keV) at ∼1000 days post-explosion and declining as L x ∝ t −1 afterwards. The broadband X-ray spectrum is of thermal origin and shows clear evidence for cooling after peak, with T ( t ) ≈ 20 keV ( t / t pk ) − 0.5 . Soft X-rays of sub-keV energy suffer from large photoelectric absorption originating from the local SN environment with NH int ( t ) ≈ 3 × 10 22 ( t / 400 days ) − 1.4 cm − 2 . We interpret these findings as the result of the interaction of the SN shock with a dense ( n ≈ 10 5 − 10 6 cm −3 ), H-rich disk-like circumstellar medium (CSM) with inner radius ∼2 × 10 16 cm and extending to ∼10 17 cm. Based on the declining NH int ( t ) and X-ray luminosity evolution, we infer a CSM mass of ∼(1.2 f –2.0 f ) M ⊙ , where f is the volume filling factor. We place SN 2014C in the context of 121 core-collapse SNe with evidence for strong shock interaction with a thick circumstellar medium. Finally, we highlight the challenges that the current mass-loss theories (including wave-driven mass loss, binary interaction, and line-driven winds) face when interpreting the wide dynamic ranges of CSM parameters inferred from observations.

We present extensive multifrequency Karl G. Jansky Very Large Array (VLA) and Very Long Baseline Array (VLBA) observations of the radio-bright supernova (SN) IIb SN 2004C that span ∼40–2793 days post-explosion. We interpret the temporal evolution of the radio spectral energy distribution in the context of synchrotron self-absorbed emission from the explosion’s forward shock as it expands in the circumstellar medium (CSM) previously sculpted by the mass-loss history of the stellar progenitor. VLBA observations and modeling of the VLA data point to a blastwave with average velocity ∼0.06 c that carries an energy of ≈10 49 erg. Our modeling further reveals a flat CSM density profile ρ CSM ∝ R −0.03±0.22 up to a break radius R br ≈ (1.96 ± 0.10) × 10 16 cm, with a steep density gradient following ρ CSM ∝ R −2.3±0.5 at larger radii. We infer that the flat part of the density profile corresponds to a CSM shell with mass ∼0.021 M ☉ , and that the progenitor’s effective mass-loss rate varied with time over the range (50–500) × 10 −5 M ☉ yr −1 for an adopted wind velocity v w = 1000 km s −1 and shock microphysical parameters ϵ e = 0.1, ϵ B = 0.01. These results add to the mounting observational evidence for departures from the traditional single-wind mass-loss scenarios in evolved, massive stars in the centuries leading up to core collapse. Potentially viable scenarios include mass loss powered by gravity waves and/or interaction with a binary companion.

We present cm-band and mm-band afterglow observations of five long-duration $\gamma$-ray bursts (GRBs; GRB 130131A, 130420B, 130609A, 131229A, 140713A) with dust-obscured optical afterglow emission, known as "dark" GRBs. We detect the radio afterglow of two of the dark GRBs (GRB 130131A and 140713A), along with a tentative detection of a third (GRB 131229A) with the Karl G. Jansky Very Large Array (VLA). Supplemented by three additional VLA-detected dark GRBs from the literature, we present uniform modeling of their broadband afterglows. We derive high line-of-sight dust extinctions of $A_{V, \rm GRB} \gtrsim 2.2 - 10.6~{\rm mag}$. Additionally, we model the host galaxies of the six bursts in our sample, and derive host galaxy dust extinctions of $A_{V, \rm Host} \approx 0.3-4.7~{\rm mag}$. Across all tested $\gamma$-ray (fluence and duration) and afterglow properties (energy scales, geometries and circumburst densities), we find dark GRBs to be representative of more typical unobscured long GRBs, except in fluence, for which observational biases and inconsistent classification may influence the dark GRB distribution. Additionally, we find that $A_{V, \rm GRB}$ is not related to a uniform distribution of dust throughout the host, nor to the extremely local environment of the burst, indicating that a larger scale patchy dust distribution is the cause of the high line-of-sight extinction. Since radio observations are invaluable to revealing heavily dust-obscured GRBs, we make predictions for the detection of radio emission from host star formation with the next generation VLA.

We conduct X-ray spectral fits on 184 likely counterparts to Fermi-LAT 3FGL unassociated sources. Characterization and classification of these sources allows for more complete population studies of the high-energy sky. Most of these X-ray spectra are well fit by an absorbed power-law model, as expected for a population dominated by blazars and pulsars. A small subset of seven X-ray sources have spectra unlike the power law expected from a blazar or pulsar and may be linked to coincident stars or background emission. We develop a multiwavelength machine learning classifier to categorize unassociated sources into pulsars and blazars using gamma-ray and X-ray observations. Training a random forest (RF) procedure with known pulsars and blazars, we achieve a cross-validated classification accuracy of 98.6%. Applying the RF routine to the unassociated sources returned 126 likely blazar candidates (defined as P-bzr >= 90%) and five likely pulsar candidates (P-bzr <= 10%). Our new X-ray spectral analysis does not drastically alter the RF classifications of these sources compared to previous works, but it builds a more robust classification scheme and highlights the importance of X-ray spectral fitting. Our procedure can be further expanded with UV, visual, or radio spectral parameters or by measuring flux variability.

We utilize machine-learning methods to distinguish BL Lacertae (BL Lac) objects from flat spectrum radio quasars (FSRQs) within a sample of likely X-ray blazar counterparts to Fermi 3FGL unassociated gamma-ray sources. From our previous work, we have extracted 84 sources that were classified as ≥99% likely to be blazars. We then utilize Swift-XRT, Fermi, and The Wide-field Infrared Survey Explorer (WISE) data together to distinguish the specific type of blazar, FSRQs, or BL Lac objects. Various X-ray and gamma-ray parameters can be used to differentiate between these subclasses. These are also known to occupy different parameter space on the WISE color–color diagram. Using all these data together would provide more robust results for the classified sources. We utilized a random forest classifier to calculate the probability for each blazar to be associated with a BL Lac object or an FSRQ. Based on P bll, which is the probability for each source to be a BL Lac object, we placed our sources into five different categories based on this value as follows: P bll ≥ 99%: highly likely BL Lac object, P bll ≥ 90%: likely BL Lac object, P bll ≤ 1%: highly likely FSRQ, P bll ≤ 10%: likely FSRQ, and 90% < P bll < 10%: ambiguous. Our results categorize the 84 blazar candidates as 50 likely BL Lac objects and the other 34 as being ambiguous. A small subset of these sources have been listed as associated sources in the most recent Fermi catalog, 4FGL, and in these cases our results are in agreement on the classification.

The Fermi-LAT unassociated sources represent some of the most enigmatic gamma-ray sources in the sky. Observations with the Swift-XRT and -UVOT telescopes have identified hundreds of likely X-ray and UV/optical counterparts in the uncertainty ellipses of the unassociated sources. In this work we present spectral fitting results for 205 possible X-ray/UV/optical counterparts to 4FGL unassociated targets. Assuming that the unassociated sources contain mostly pulsars and blazars, we develop a neural network classifier approach that applies gamma-ray, X-ray, and UV/optical spectral parameters to yield a descriptive classification of unassociated spectra into pulsars and blazars. From our primary sample of 174 Fermi sources with a single X-ray/UV/optical counterpart, we present 132 P bzr > 0.99 likely blazars and 14 P bzr < 0.01 likely pulsars, with 28 remaining ambiguous. These subsets of the unassociated sources suggest a systematic expansion to catalogs of gamma-ray pulsars and blazars. Compared to previous classification approaches our neural network classifier achieves significantly higher validation accuracy and returns more bifurcated P bzr values, suggesting that multiwavelength analysis is a valuable tool for confident classification of Fermi unassociated sources.

We present a population of 19 radio-luminous supernovae (SNe) with emission reaching L-v similar to 10(26)-10(29) erg s(-1) Hz(-1) in the first epoch of the Very Large Array Sky Survey (VLASS) at 2-4 GHz. Our sample includes one long gamma-ray burst, SN 2017iuk/GRB 171205A, and 18 core-collapse SNe detected at approximate to 1-60 yr after explosion. No thermonuclear explosion shows evidence for bright radio emission, and hydrogen-poor progenitors dominate the subsample of corecollapse events with spectroscopic classification at the time of explosion (79%). We interpret these findings in the context of the expected radio emission from the forward shock interaction with the circumstellar medium (CSM). We conclude that these observations require a departure from the single wind-like density profile (i.e., rho(csm) proportional to r(-2)) that is expected around massive stars and/or from a spherical Newtonian shock. Viable alternatives include the shock interaction with a detached, dense shell of CSM formed by a large effective progenitor mass-loss rate, (M) over dot similar to 10(-4)-10(-1) M-circle dot yr(-1) (for an assumed wind velocity of 1000 km s(-1)); emission from an off-axis relativistic jet entering our line of sight; or the emergence of emission from a newly born pulsar-wind nebula. The relativistic SN 2012ap that is detected 5.7 and 8.5 yr after explosion with L-v similar to 10(28) erg s(-1) Hz(-1) might constitute the first detections of an off-axis jet+cocoon system in a massive star. However, none of the VLASS SNe with archival data points are consistent with our model off-axis jet light curves. Future multiwavelength observations will distinguish among these scenarios. Our VLASS source catalogs, which were used to perform the VLASS cross-matching, are publicly available at https://doi.org/10.5281/zenodo.4895112.

Detections of SiO masers from the Bulge Asymmetries and Dynamical Evolution (BAaDE) survey more tightly define the region where Oxygen-rich (O) Asymptotic Giant Branch (AGB) stars reside in multiple infrared (IR) color-color diagrams. Using MSX and 2MASS data along with radio spectra from the BAaDE survey we find that three main populations were observed in the BAaDE survey: O-rich AGB stars of which about $73\%$ host SiO masers, Carbon-rich (C) AGB stars which do not host these masers, and a small contaminating set of possible Young Stellar Objects (YSOs). The distinction between YSOs and AGB stars can be drawn using only MSX data, specifically the $[D]-[E]$ color, while the shorter wavelengths provided by 2MASS are necessary to divide potential C- and O-rich AGB stars. Divisions similar to these have been seen in multiple earlier IR-studies, but BAaDE currently provides a sample of $\sim$ 15000 sources which far exceeds previous studies in sample size, and therefore provides much more distinct divisions. With these IR distinctions in place, we discuss the sources which are exceptions in either their molecular detections or IR colors, as well as the distribution of the three populations in Galactic coordinates.

We utilize machine learning methods to distinguish BL Lacertae objects (BL Lac) from Flat Spectrum Radio Quasars (FSRQ) within a sample of likely X-ray blazar counterparts to Fermi 3FGL unassociated gamma-ray sources. From our previous work, we have extracted 84 sources that were classified as ≥ 99 likley to be blazars. We then utilize Swift-XRT, Fermi, and WISE (The Wide-field Infrared Survey Explorer) data together to distinguish the specific type of blazar, FSRQs or BL Lacs. Various X-ray and Gamma-ray parameters can be used to differentiate between these subclasses. These are also known to occupy different parameter space on the WISE color-color diagram. Using all these data together would provide more robust results for the classified sources. We utilized a Random Forest Classifier to calculate the probability for each blazar to be associated with a BL Lac or an FSRQ. Based on P_bll, which is the probability for each source to be a BL Lac, we placed our sources into five different categories based on this value as follows; P_bll ≥ 99 highly likely BL Lac, P_bll ≥ 90 highly likely FSRQ, P_bll ≤ 10 < 10 likely BL Lacs and the rest 34 being ambiguous. A small subset of these sources have been listed as associated sources in the most recent Fermi catalog, 4FGL, and in these cases our results are in agreement on the classification.

Radio campaigns using maser stellar beacons have provided crucial information to characterize Galactic stellar populations. Currently, the Bulge Asymmetries and Dynamical Evolution (BAaDE) project is surveying infrared (IR) color-selected targets for SiO masers. This provides a sample of evolved stars that can be used to study the inner, optically obscured Galaxy using line-of-sight velocities and possibly very long baseline interferometry proper motions. In order to use the BAaDE sample for kinematic studies, the stellar population should be characterized. In this study, the BAaDE targets have been cross-matched with IR (2MASS) and optical (Gaia) samples. By exploring the synergies of this cross-match together with Gaia parallaxes and extinction maps, the local (d < 2 kpc) asymptotic giant branch (AGB) stars can be characterized. We have defined a BAaDE-Gaia sample of 20,111 sources resulting from cross-matching BAaDE targets with IR and optical surveys. From this sample, a local sample of 1812 evolved stars with accurate parallax measurements, confirmed evolved stellar evolution stage, and within 2 kpc distance around the Sun was selected, for which absolute (bolometric) magnitudes are estimated. The evolved stellar population with Gaia counterparts that are variable seems to be predominantly associated with AGB stars with moderate luminosity (1500(-500)(+3000) L-circle dot) and periods between 250 and 1250 days.

We report on the first 1432 sources observed using the Atacama Large Millimeter/submillimeter Array, from the Bulge Asymmetries and Dynamical Evolution survey, which aims to obtain tens of thousands of line-of-sight velocities from SiO masers in Asymptotic Giant Branch (AGB) stars in the Milky Way. A 71% detection rate of 86 GHz SiO masers is obtained from the infrared color-selected sample and increases to 80% when considering the likely oxygen-rich stars using Midcourse Space Experiment colors isolated in a region where [D]-[E] <= 1.38. Based on Galactic distributions, the presence of extended CS emission, and likely kinematic associations, the population of sources with [D]-[E] > 1.38 probably consists of young stellar objects, or alternatively, of planetary nebulae. For the SiO detections, we examined whether individual SiO transitions provide comparable stellar line-of-sight velocities and found that any SiO transition is suitable for determining a stellar AGB line-of-sight velocity. Finally, we discuss the relative SiO detection rates and line strengths in the context of current pumping models.

Abstract The BAaDE (Bulge Asymmetries and Dynamical Evolution) project is an SiO maser survey of the Galactic Plane. About 19,000 sources have been observed at 43 GHz with the VLA, and the production of spectra for each of these sources is well underway. The primary goal of the project is to collect line-of-sight velocities for all the detected masers in the sample to probe Galactic dynamics. With an expected detection rate of over 60% we should collect over 11,000 velocities to probe the Galactic potential. The survey is also a large sample of infrared sources to explore the different evolved stellar populations within the Milky Way. So far we discern three distinct groups in the BAaDE sample: the main group containing oxygen-rich, evolved stars with a high SiO maser detection rate, a much smaller population of carbon-rich evolved stars, and finally a group of likely young stellar objects with no maser emission. These populations are separated using 2MASS and MSX color-color diagrams, and we find a particularly useful cut between the young and evolved objects using the MSX [D] –[E] color. Identification of these populations will isolate BAaDE’s evolved star sample, and will more tightly define the region in IR color-color diagrams where SiO masers occur yielding a better understanding of these kinematical probes. Using our color-divisions we can also study the distribution of each of the populations within the Galactic Plane.