Search for: All records

Joint observations in electromagnetic and gravitational waves shed light on the physics of objects and surrounding environments with extreme gravity that are otherwise unreachable via siloed observations in each messenger. However, such detections remain challenging due to the rapid and faint nature of counterparts. Protocols for discovery and inference still rely on human experts manually inspecting survey alert streams and intuiting optimal usage of limited follow-up resources. Strategizing an optimal follow-up program requires adaptive sequential decision-making given evolving light curve data that (i) maximizes a global objective despite incomplete information and (ii) is robust to stochasticity introduced by detectors/observing conditions. Reinforcement learning (RL) approaches allow agents to implicitly learn the physics or detector dynamics and the behavior policy that maximizes a designated objective through experience. To demonstrate the utility of such an approach for the kilonova follow-up problem, we train a toy RL agent with the goal of maximizing follow-up photometry for the true kilonova among several contaminant transient light curves. In a simulated environment where the agent learns online, it achieves 3× higher accuracy compared to a random strategy. However, it is surpassed by human agents by up to a factor of 2. This is likely because our hypothesis function (Qthat is linear in state-action features) is an insufficient representation of the optimal behavior policy. More complex agents could perform at par or surpass human experts. Agents like these could pave the way for machine-directed software infrastructure to efficiently respond to next generation detectors, for conducting science inference and optimally planning expensive follow-up observations, scalably and with demonstrable performance guarantees.

 

We present results from a systematic infrared (IR) census of R Coronae Borealis (RCB) stars in the Milky Way, using data from the Palomar Gattini IR (PGIR) survey. RCB stars are dusty, erratic variable stars presumably formed from the merger of a He-core and a CO-core white dwarf (WD). PGIR is a 30 cmJ-band telescope with a 25 deg²camera that surveys 18,000 deg²of the northern sky (δ> −28°) at a cadence of 2 days. Using PGIRJ-band lightcurves for ∼60 million stars together with mid-IR colors from WISE, we selected a sample of 530 candidate RCB stars. We obtained near-IR spectra for these candidates and identified 53 RCB stars in our sample. Accounting for our selection criteria, we find that there are a total of$\approx {350}_{-100}^{+150}$RCB stars in the Milky Way. Assuming typical RCB lifetimes, this corresponds to an RCB formation rate of 0.8–5 × 10⁻³yr⁻¹, consistent with observational and theoretical estimates of the He-CO WD merger rate. We searched for quasi-periodic pulsations in the PGIR lightcurves of RCB stars and present pulsation periods for 16 RCB stars. We also examined high-cadenced TESS lightcurves for RCB and the chemically similar, but dustless hydrogen-deficient carbon (dLHdC) stars. We find that dLHdC stars show variations on timescales shorter than RCB stars, suggesting that they may have lower masses than RCB stars. Finally, we identified 3 new spectroscopically confirmed and 12 candidate Galactic DY Per type stars—believed to be colder cousins of RCB star—doubling the sample of Galactic DY Per type stars.

 

We study a magnitude-limited sample of 36 broad-lined type Ic supernovae (SNe Ic-BL) from the Zwicky Transient Facility Bright Transient Survey (detected between 2018 March and 2021 August), which is the largest systematic study of SNe Ic-BL done in literature thus far. We present the light curves (LCs) for each of the SNe and analyze the shape of the LCs to derive empirical parameters, along with the explosion epochs for every event. The sample has an average absolute peak magnitude in therband of${\overline{\mathit{M}}}_{\mathrm{r},\max }=-18.51\pm 0.15$mag. Using spectra obtained around peak light, we compute expansion velocities from the Feii5169 Å line for each event with high enough signal-to-noise ratio spectra, and find an average value of$\overline{v_{\mathrm{ph}}}=16,100\pm 1100$km s⁻¹. We also compute bolometric LCs, study the blackbody temperature and radii evolution over time, and derive the explosion properties of the SNe. The explosion properties of the sample have average values of${\overline{\mathit{M}}}_{\mathrm{Ni}}={0.37}_{-0.06}^{+0.08}{\mathit{M}}_{\odot }$,${\overline{\mathit{M}}}_{\mathrm{ej}}={2.45}_{-0.41}^{+0.47}{\mathit{M}}_{\odot }$, and${\overline{\mathit{E}}}_{\mathrm{K}}=({4.02}_{-1.00}^{+1.37})\times {10}^{51}$erg. Thirteen events have radio observations from the Very Large Array, with eight detections and five non-detections. We find that the populations that have radio detections and radio non-detections are indistinct from one another with respect to their optically inferred explosion properties, and there are no statistically significant correlations present between the events’ radio luminosities and optically inferred explosion properties. This provides evidence that the explosion properties derived from optical data alone cannot give inferences about the radio properties of SNe Ic-BL and likely their relativistic jet formation mechanisms.

 

An advanced LIGO and Virgo’s third observing run brought another binary neutron star merger (BNS) and the first neutron-star black hole mergers. While no confirmed kilonovae were identified in conjunction with any of these events, continued improvements of analyses surrounding GW170817 allow us to project constraints on the Hubble Constant (H₀), the Galactic enrichment fromr-process nucleosynthesis, and ultra-dense matter possible from forthcoming events. Here, we describe the expected constraints based on the latest expected event rates from the international gravitational-wave network and analyses of GW170817. We show the expected detection rate of gravitational waves and their counterparts, as well as how sensitive potential constraints are to the observed numbers of counterparts. We intend this analysis as support for the community when creating scientifically driven electromagnetic follow-up proposals. During the next observing run O4, we predict an annual detection rate of electromagnetic counterparts from BNS of${0.43}_{-0.26}^{+0.58}$(${1.97}_{-1.2}^{+2.68}$) for the Zwicky Transient Facility (Rubin Observatory).

 

Eruptive mass loss of massive stars prior to supernova (SN) explosion is key to understanding their evolution and end fate. An observational signature of pre-SN mass loss is the detection of an early, short-lived peak prior to the radioactive-powered peak in the lightcurve of the SN. This is usually attributed to the SN shock passing through an extended envelope or circumstellar medium. Such an early peak is common for double-peaked Type IIb SNe with an extended hydrogen envelope but uncommon for normal Type Ibc SNe with very compact progenitors. In this paper, we systematically study a sample of 14 double-peaked Type Ibc SNe out of 475 Type Ibc SNe detected by the Zwicky Transient Facility. The rate of these events is ∼3%–9% of Type Ibc SNe. A strong correlation is seen between the peak brightness of the first and the second peak. We perform a holistic analysis of this sample’s photometric and spectroscopic properties. We find that six SNe have ejecta mass less than 1.5M_⊙. Based on the nebular spectra and lightcurve properties, we estimate that the progenitor masses for these are less than ∼12M_⊙. The rest have an ejecta mass >2.4M_⊙and a higher progenitor mass. This sample suggests that the SNe with low progenitor masses undergo late-time binary mass transfer. Meanwhile, the SNe with higher progenitor masses are consistent with wave-driven mass loss or pulsation-pair instability-driven mass-loss simulations.

 

During the first half of the fourth observing run (O4a) of the International Gravitational Wave Network, the Zwicky Transient Facility (ZTF) conducted a systematic search for kilonova (KN) counterparts to binary neutron star (BNS) and neutron star–black hole (NSBH) merger candidates. Here, we present a comprehensive study of the five high-significance (False Alarm Rate less than 1 yr⁻¹) BNS and NSBH candidates in O4a. Our follow-up campaigns relied on both target-of-opportunity observations and re-weighting of the nominal survey schedule to maximize coverage. We describe the toolkit we have been developing,Fritz, an instance ofSkyPortal, instrumental in coordinating and managing our telescope scheduling, candidate vetting, and follow-up observations through a user-friendly interface. ZTF covered a total of 2841 deg²within the skymaps of the high-significance GW events, reaching a median depth ofg≈ 20.2 mag. We circulated 15 candidates, but found no viable KN counterpart to any of the GW events. Based on the ZTF non-detections of the high-significance events in O4a, we used a Bayesian approach,nimbus, to quantify the posterior probability of KN model parameters that are consistent with our non-detections. Our analysis favors KNe with initial absolute magnitude fainter than −16 mag. The joint posterior probability of a GW170817-like KN associated with all our O4a follow-ups was 64%. Additionally, we use a survey simulation software,simsurvey, to determine that our combined filtered efficiency to detect a GW170817-like KN is 36%, when considering the 5 confirmed astrophysical events in O3 (1 BNS and 4 NSBH events), along with our O4a follow-ups. Following Kasliwal et al., we derived joint constraints on the underlying KN luminosity function based on our O3 and O4a follow-ups, determining that no more than 76% of KNe fading at 1 mag day⁻¹can peak at a magnitude brighter than −17.5 mag.

 

Electromagnetic follow-up of gravitational-wave detections is very resource intensive, taking up hours of limited observation time on dozens of telescopes. Creating more efficient schedules for follow-up will lead to a commensurate increase in counterpart location efficiency without using more telescope time. Widely used in operations research and telescope scheduling, mixed-integer linear programming is a strong candidate to produce these higher-efficiency schedules, as it can make use of powerful commercial solvers that find globally optimal solutions to provided problems. We detail a new target-of-opportunity scheduling algorithm designed with Zwicky Transient Facility in mind that uses mixed-integer linear programming. We compare its performance togwemopt, the tuned heuristic scheduler used by the Zwicky Transient Facility and other facilities during the third LIGO–Virgo gravitational-wave observing run. This new algorithm uses variable-length observing blocks to enforce cadence requirements and to ensure field observability, along with having a secondary optimization step to minimize slew time. We show that by employing a hybrid method utilizing both this scheduler andgwemopt, the previous scheduler used, in concert, we can achieve an average improvement in detection efficiency of 3%–11% overgwemoptalone for a simulated binary neutron star merger data set consistent with LIGO–Virgo’s third observing run, highlighting the potential of mixed-integer target of opportunity schedulers for future multimessenger follow-up surveys.

 

Abstract The detonation of a thin (≲0.03 M ⊙ ) helium shell (He-shell) atop a ∼1 M ⊙ white dwarf (WD) is a promising mechanism to explain normal Type Ia supernovae (SNe Ia), while thicker He-shells and less massive WDs may explain some recently observed peculiar SNe Ia. We present observations of SN 2020jgb, a peculiar SN Ia discovered by the Zwicky Transient Facility (ZTF). Near maximum brightness, SN 2020jgb is slightly subluminous (ZTF g -band absolute magnitude −18.7 mag ≲ M g ≲ −18.2 mag depending on the amount of host-galaxy extinction) and shows an unusually red color (0.2 mag ≲ g ZTF − r ZTF ≲ 0.4 mag) due to strong line-blanketing blueward of ∼5000 Å. These properties resemble those of SN 2018byg, a peculiar SN Ia consistent with an He-shell double detonation (DDet) SN. Using detailed radiative transfer models, we show that the optical spectroscopic and photometric evolution of SN 2020jgb is broadly consistent with a ∼0.95–1.00 M ⊙ (C/O core + He-shell) progenitor ignited by a ≳0.1 M ⊙ He-shell. However, one-dimensional radiative transfer models without non-local-thermodynamic-equilibrium treatment cannot accurately characterize the line-blanketing features, making the actual shell mass uncertain. We detect a prominent absorption feature at ∼1 μ m in the near-infrared (NIR) spectrum of SN 2020jgb, which might originate from unburnt helium in the outermost ejecta. While the sample size is limited, we find similar 1 μ m features in all the peculiar He-shell DDet candidates with NIR spectra obtained to date. SN 2020jgb is also the first peculiar He-shell DDet SN discovered in a star-forming dwarf galaxy, indisputably showing that He-shell DDet SNe occur in both star-forming and passive galaxies, consistent with the normal SN Ia population. 

The fate of stars in the zero-age main-sequence (ZAMS) range ≈8–12M_⊙is unclear. They could evolve to form white dwarfs or explode as electron-capture supernovae (SNe) or iron core-collapse SNe (CCSNe). Even though the initial mass function indicates that this mass range should account for over 40% of all CCSN progenitors, few have been observationally confirmed, likely due to the faintness and rapid evolution of some of these transients. In this paper, we present a sample of nine Ca-rich/O-poor Type IIb SNe detected by the Zwicky Transient Facility with progenitors likely in this mass range. These sources have a [Caii]λλ7291, 7324/[Oi]λλ6300, 6364 flux ratio of ≳2 in their nebular spectra. Comparing the measured [Oi] luminosity (≲10³⁹erg s⁻¹) and derived oxygen mass (≈0.01M_⊙) with theoretical models, we infer that the progenitor ZAMS mass for these explosions is less than 12M_⊙. The ejecta properties (M_ej≲ 1M_⊙andE_kin∼ 10⁵⁰erg) are also consistent. The low ejecta mass of these sources indicates a class of strongly-stripped SNe that is a transition between the regular stripped-envelope SNe and ultra-stripped SNe. The progenitor could be stripped by a main-sequence companion and result in the formation of a neutron star−main sequence binary. Such binaries have been suggested to be progenitors of neutron star−white dwarf systems that could merge within a Hubble time and be detectable with LISA.

 

One of the open questions following the discovery of GW170817 is whether neutron star (NS) mergers are the only astrophysical sites capable of producingr-process elements. Simulations have shown that 0.01–0.1M_⊙ofr-process material could be generated in the outflows originating from the accretion disk surrounding the rapidly rotating black hole that forms as a remnant to both NS mergers and collapsing massive stars associated with long-duration gamma-ray bursts (collapsars). The hallmark signature ofr-process nucleosynthesis in the binary NS merger GW170817 was its long-lasting near-infrared (NIR) emission, thus motivating a systematic photometric study of the light curves of broad-lined stripped-envelope (Ic-BL) supernovae (SNe) associated with collapsars. We present the first systematic study of 25 SNe Ic-BL—including 18 observed with the Zwicky Transient Facility and 7 from the literature—in the optical/NIR bands to determine what quantity ofr-process material, if any, is synthesized in these explosions. Using semi-analytic models designed to account forr-process production in SNe Ic-BL, we perform light curve fitting to derive constraints on ther-process mass for these SNe. We also perform independent light curve fits to models without ther-process. We find that ther-process-free models are a better fit to the light curves of the objects in our sample. Thus, we find no compelling evidence ofr-process enrichment in any of our objects. Further high-cadence infrared photometric studies and nebular spectroscopic analysis would be sensitive to smaller quantities ofr-process ejecta mass or indicate whether all collapsars are completely devoid ofr-process nucleosynthesis.