Search for: All records

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

 

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.

 

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.

 

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. 

We report observations of the optical counterpart of the long gamma-ray burst (GRB) GRB 230812B and its associated supernova (SN) SN 2023pel. The proximity (z= 0.36) and high energy (E_γ,iso∼ 10⁵³erg) make it an important event to study as a probe of the connection between massive star core collapse and relativistic jet formation. With a phenomenological power-law model for the optical afterglow, we find a late-time flattening consistent with the presence of an associated SN. SN 2023pel has an absolute peakr-band magnitude ofM_r= −19.46 ± 0.18 mag (about as bright as SN 1998bw) and evolves on quicker timescales. Using a radioactive heating model, we derive a nickel mass powering the SN ofM_Ni= 0.38 ± 0.01M_⊙and a peak bolometric luminosity ofL_bol∼ 1.3 × 10⁴³erg s⁻¹. We confirm SN 2023pel’s classification as a broad-line Type Ic SN with a spectrum taken 15.5 days after its peak in therband and derive a photospheric expansion velocity ofv_ph= 11,300 ± 1600 km s⁻¹at that phase. Extrapolating this velocity to the time of maximum light, we derive the ejecta massM_ej= 1.0 ± 0.6M_⊙and kinetic energy$E_{\mathrm{KE}}={1.3}_{-1.2}^{+3.3}\times {10}^{51}\mathrm{erg}$. We find that GRB 230812B/SN 2023pel has SN properties that are mostly consistent with the overall GRB-SN population. The lack of correlations found in the GRB-SN population between SN brightness andE_γ,isofor their associated GRBs across a broad range of 7 orders of magnitude provides further evidence that the central engine powering the relativistic ejecta is not coupled to the SN powering mechanism in GRB-SN systems.

 

We report observations of the optical counterpart of the long gamma-ray burst GRB 221009A. Due to the extreme rarity of being both nearby (z= 0.151) and highly energetic (E_γ,iso≥ 10⁵⁴erg), GRB 221009A offers a unique opportunity to probe the connection between massive star core collapse and relativistic jet formation across a very broad range ofγ-ray properties. Adopting a phenomenological power-law model for the afterglow and host galaxy estimates from high-resolution Hubble Space Telescope imaging, we use Bayesian model comparison techniques to determine the likelihood of an associated supernova (SN) contributing excess flux to the optical light curve. Though not conclusive, we find moderate evidence (K_Bayes= 10^1.2) for the presence of an additional component arising from an associated SN, SN 2022xiw, and find that it must be substantially fainter (<67% as bright at the 99% confidence interval) than SN 1998bw. Given the large and uncertain line-of-sight extinction, we attempt to constrain the SN parameters (M_Ni,M_ej, andE_KE) under several different assumptions with respect to the host galaxy’s extinction. We find properties that are broadly consistent with previous GRB-associated SNe:M_Ni= 0.05–0.25M_⊙,M_ej= 3.5–11.1M_⊙, andE_KE= (1.6–5.2) × 10⁵²erg. We note that these properties are weakly constrained due to the faintness of the SN with respect to the afterglow and host emission, but we do find a robust upper limit onM_NiofM_Ni< 0.36M_⊙. Given the tremendous range in isotropic gamma-ray energy release exhibited by GRBs (seven orders of magnitude), the SN emission appears to be decoupled from the central engine in these systems.

 

Luminous red novae (LRNe) are transients characterized by low luminosities and expansion velocities, and they are associated with mergers or common-envelope ejections in stellar binaries. Intermediate-luminosity red transients (ILRTs) are an observationally similar class with unknown origins, but they are generally believed to be either electron-capture supernovae in super-asymptotic giant branch stars or outbursts in dusty luminous blue variables (LBVs). In this paper, we present a systematic sample of eight LRNe and eight ILRTs detected as part of the Census of the Local Universe (CLU) experiment on the Zwicky Transient Facility (ZTF). The CLU experiment spectroscopically classifies ZTF transients associated with nearby (<150 Mpc) galaxies, achieving 80% completeness form_r< 20 mag. Using the ZTF-CLU sample, we derive the first systematic LRNe volumetric rate of${7.8}_{-3.7}^{+6.5}\times {10}^{-5}$Mpc⁻³yr⁻¹in the luminosity range −16 ≤M_r≤ −11 mag. We find that, in this luminosity range, the LRN rate scales as$\mathit{dN}/\mathit{dL}\propto L^{-2.5\pm 0.3}$—significantly steeper than the previously derived scaling ofL^−1.4±0.3for lower-luminosity LRNe (M_V≥ −10 mag). The steeper power law for LRNe at high luminosities is consistent with the massive merger rates predicted by binary population synthesis models. We find that the rates of the brightest LRNe (M_r≤ −13 mag) are consistent with a significant fraction of them being progenitors of double compact objects that merge within a Hubble time. For ILRTs, we derive a volumetric rate of${2.6}_{-1.4}^{+1.8}\times {10}^{-6}$Mpc⁻³yr⁻¹forM_r≤ −13.5 mag, which scales as$\mathit{dN}/\mathit{dL}\propto L^{-2.5\pm 0.5}$. This rate is ∼1%–5% of the local core-collapse supernova rate and is consistent with theoretical ECSN rate estimates.

 

Searches for electromagnetic counterparts of gravitational-wave signals have redoubled since the first detection in 2017 of a binary neutron star merger with a gamma-ray burst, optical/infrared kilonova, and panchromatic afterglow. Yet, one LIGO/Virgo observing run later, there has not yet been a second, secure identification of an electromagnetic counterpart. This is not surprising given that the localization uncertainties of events in LIGO and Virgo’s third observing run, O3, were much larger than predicted. We explain this by showing that improvements in data analysis that now allow LIGO/Virgo to detect weaker and hence more poorly localized events have increased the overall number of detections, of which well-localized,gold-platedevents make up a smaller proportion overall. We present simulations of the next two LIGO/Virgo/KAGRA observing runs, O4 and O5, that are grounded in the statistics of O3 public alerts. To illustrate the significant impact that the updated predictions can have, we study the follow-up strategy for the Zwicky Transient Facility. Realistic and timely forecasting of gravitational-wave localization accuracy is paramount given the large commitments of telescope time and the need to prioritize which events are followed up. We include a data release of our simulated localizations as a public proposal planning resource for astronomers.

 

The Zwicky Transient Facility (ZTF) performs a systematic neutrino follow-up programme, searching for optical counterparts to high-energy neutrinos with dedicated Target-of-Opportunity (ToO) observations. Since first light in March 2018, ZTF has taken prompt observations for 24 high-quality neutrino alerts from the IceCube Neutrino Observatory, with a median latency of 12.2 h from initial neutrino detection. From two of these campaigns, we have already reported tidal disruption event (TDE) AT 2019dsg and likely TDE AT 2019fdr as probable counterparts, suggesting that TDEs contribute >7.8 per cent of the astrophysical neutrino flux. We here present the full results of our programme through to December 2021. No additional candidate neutrino sources were identified by our programme, allowing us to place the first constraints on the underlying optical luminosity function of astrophysical neutrino sources. Transients with optical absolutes magnitudes brighter that −21 can contribute no more than 87 per cent of the total, while transients brighter than −22 can contribute no more than 58 per cent of the total, neglecting the effect of extinction and assuming they follow the star formation rate. These are the first observational constraints on the neutrino emission of bright populations such as superluminous supernovae. None of the neutrinos were coincident with bright optical AGN flares comparable to that observed for TXS 0506+056/IC170922A, with such optical blazar flares producing no more than 26 per cent of the total neutrino flux. We highlight the outlook for electromagnetic neutrino follow-up programmes, including the expected potential for the Rubin Observatory.