Search

Abstract The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic emission, and dense matter reliably. In addition, one needs efficient and accurate computational tools to ensure a correct cross-correlation between the models and the observational data. For this purpose, we have developed the Nuclear-physics and Multi-Messenger Astrophysics framework NMMA. The code allows incorporation of nuclear-physics constraints at low densities as well as X-ray and radio observations of isolated neutron stars. In previous works, the NMMA code has allowed us to constrain the equation of state of supranuclear dense matter, to measure the Hubble constant, and to compare dense-matter physics probed in neutron-star mergers and in heavy-ion collisions, and to classify electromagnetic observations and perform model selection. Here, we show an extension of the NMMA code as a first attempt of analyzing the gravitational-wave signal, the kilonova, and the gamma-ray burst afterglow simultaneously. Incorporating all available information, we estimate the radius of a 1.4M ⊙ neutron star to be $$R=11.9{8}_{-0.40}^{+0.35}$$ R = 11.9 8 − 0.40 + 0.35 km.

We present initial results from a James Webb Space Telescope (JWST) survey of the youngest Galactic core-collapse supernova remnant, Cassiopeia A (Cas A), made up of NIRCam and MIRI imaging mosaics that map emission from the main shell, interior, and surrounding circumstellar/interstellar material (CSM/ISM). We also present four exploratory positions of MIRI Medium Resolution Spectrograph integral field unit spectroscopy that sample ejecta, CSM, and associated dust from representative shocked and unshocked regions. Surprising discoveries include (1) a weblike network of unshocked ejecta filaments resolved to ∼0.01 pc scales exhibiting an overall morphology consistent with turbulent mixing of cool, low-entropy matter from the progenitor’s oxygen layer with hot, high-entropy matter heated by neutrino interactions and radioactivity; (2) a thick sheet of dust-dominated emission from shocked CSM seen in projection toward the remnant’s interior pockmarked with small (∼1″) round holes formed by ≲0.″1 knots of high-velocity ejecta that have pierced through the CSM and driven expanding tangential shocks; and (3) dozens of light echoes with angular sizes between ∼0.″1 and 1′ reflecting previously unseen fine-scale structure in the ISM. NIRCam observations place new upper limits on infrared emission (≲20 nJy at 3 μ m) from the neutron star in Cas A’s center and tightly constrain scenarios involving a possible fallback disk. These JWST survey data and initial findings help address unresolved questions about massive star explosions that have broad implications for the formation and evolution of stellar populations, the metal and dust enrichment of galaxies, and the origin of compact remnant objects.

JWST observations of the young Galactic supernova remnant Cassiopeia A revealed an unexpected structure seen as a green emission feature in colored composite MIRI F1130W and F1280W images—hence dubbed the Green Monster—that stretches across the central parts of the remnant in projection. Combining the kinematic information from NIRSpec and the MIRI Medium Resolution Spectrograph with the multiwavelength imaging from NIRCam and MIRI, we associate the Green Monster with circumstellar material (CSM) that was lost during an asymmetric mass-loss phase. MIRI images are dominated by dust emission, but their spectra show emission lines from Ne, H, and Fe with low radial velocities indicative of a CSM nature. An X-ray analysis of this feature in a companion paper supports its CSM nature and detects significant blueshifting, thereby placing the Green Monster on the nearside, in front of the Cas A supernova remnant. The most striking features of the Green Monster are dozens of almost perfectly circular 1″–3″ sized holes, most likely created by interaction between high-velocity supernova ejecta material and the CSM. Further investigation is needed to understand whether these holes were formed by small 8000–10,500 km s ^−1 N-rich ejecta knots that penetrated and advanced out ahead of the remnant’s 5000–6000 km s ^−1 outer blast wave or by narrow ejecta fingers that protrude into the forward-shocked CSM. The detection of the Green Monster provides further evidence of the highly asymmetric mass loss that Cas A’s progenitor star underwent prior to its explosion.

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.5 M _⊙ . Based on the nebular spectra and lightcurve properties, we estimate that the progenitor masses for these are less than ∼12 M _⊙ . The rest have an ejecta mass >2.4 M _⊙ 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.

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 ( Q that 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 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 the r band of ${\overline{{M}}}_{{\rm{r}},\max }=-18.51\pm 0.15$ mag. Using spectra obtained around peak light, we compute expansion velocities from the Fe ii 5169 Å line for each event with high enough signal-to-noise ratio spectra, and find an average value of $\bar{{v}_{\mathrm{ph}}}=16,100\pm 1100$ km s ^−1 . 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{{M}}}_{\mathrm{Ni}}={0.37}_{-0.06}^{+0.08}\,{{M}}_{\odot }$ , ${\overline{{M}}}_{\mathrm{ej}}={2.45}_{-0.41}^{+0.47}\,{{M}}_{\odot }$ , and ${\overline{{E}}}_{{\rm{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.

We analyze a sample of 45 Type II supernovae from the Zwicky Transient Facility public survey using a grid of hydrodynamical models in order to assess whether theoretically driven forecasts can intelligently guide follow-up observations supporting all-sky survey alert streams. We estimate several progenitor properties and explosion physics parameters, including zero-age main-sequence (ZAMS) mass, mass-loss rate, kinetic energy, ^56 Ni mass synthesized, host extinction, and the time of the explosion. Using complete light curves we obtain confident characterizations for 34 events in our sample, with the inferences of the remaining 11 events limited either by poorly constraining data or the boundaries of our model grid. We also simulate real-time characterization of alert stream data by comparing our model grid to various stages of incomplete light curves (Δ t < 25 days, Δ t < 50 days, all data), and find that some parameters are more reliable indicators of true values at early epochs than others. Specifically, ZAMS mass, time of the explosion, steepness parameter β , and host extinction are reasonably constrained with incomplete light-curve data, whereas mass-loss rate, kinetic energy, and ^56 Ni mass estimates generally require complete light curves spanning >100 days. We conclude that real-time modeling of transients, supported by multi-band synthetic light curves tailored to survey passbands, can be used as a powerful tool to identify critical epochs of follow-up observations. Our findings are relevant to identifying, prioritizing, and coordinating efficient follow-up of transients discovered by the Vera C. Rubin Observatory.

We present proper motion measurements of the oxygen-rich ejecta of the LMC supernova remnant N132D using two epochs of Hubble Space Telescope Advanced Camera for Surveys data spanning 16 years. The proper motions of 120 individual knots of oxygen-rich gas were measured and used to calculate a center of expansion (CoE) of α = 5 ^h 25 ^m 01.ˢ71 and δ = −69°38′41.″64 (J2000) with a 1 σ uncertainty of 2.″90. This new CoE measurement is 9.″2 and 10.″8 from two previous CoE estimates based on the geometry of the optically emitting ejecta. We also derive an explosion age of 2770 ± 500 yr, which is consistent with recent age estimates of ≈2500 yr made from 3D ejecta reconstructions. We verified our estimates of the CoE and age using a new automated procedure that detected and tracked the proper motions of 137 knots, with 73 knots that overlap with the visually identified knots. We find that the proper motions of the ejecta are still ballistic, despite the remnant’s age, and are consistent with the notion that the ejecta are expanding into an interstellar medium cavity. Evidence for explosion asymmetry from the parent supernova is also observed. Using the visually measured proper motion measurements and corresponding CoE and age, we compare N132D to other supernova remnants with proper motion ejecta studies.

Narrow transient emission lines (flash-ionization features) in early supernova (SN) spectra trace the presence of circumstellar material (CSM) around the massive progenitor stars of core-collapse SNe. The lines disappear within days after the SN explosion, suggesting that this material is spatially confined, and originates from enhanced mass loss shortly (months to a few years) prior to the explosion. We performed a systematic survey of H-rich (Type II) SNe discovered within less than 2 days from the explosion during the first phase of the Zwicky Transient Facility survey (2018–2020), finding 30 events for which a first spectrum was obtained within 36% at the 95% confidence level) confirms that elevated mass loss in massive stars prior to SN explosion is common. We find that SNe II showing flash-ionization features are not significantly brighter, nor bluer, nor more slowly rising than those without. This implies that CSM interaction does not contribute significantly to their early continuum emission, and that the CSM is likely optically thin. We measured the persistence duration of flash-ionization emission and find that most SNe show flash features for ≈5 days. Rarer events, with persistence timescales >10 days, are brighter and rise longer, suggesting these may be intermediate between regular SNe II and strongly interacting SNe IIn.

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 _0 ), the Galactic enrichment from r -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).