Search

The second brightest gamma-ray burst, GRB 230307A (with a duration T _90 ∼ 40 s), exhibited characteristics indicative of a magnetar engine during the prompt emission phase. Notably, a suspected kilonova was identified in its follow-up optical and infrared observations. Here we propose that the origin of GRB 230307A is a neutron star–white dwarf (NS–WD) merger as this could naturally explain the long duration and the large physical offset from the center of its host galaxy. In the framework of such an NS–WD merger event, the late-time kilonova-like emission is very likely to be powered by the spin-down of the magnetar and the radioactive decay of ^56 Ni, rather than by the decay of r -process elements as these heavy elements may not easily be synthesized in an NS-WD merger. It is demonstrated that the above scenario can be supported by our fit to the late-time observational data, where a mass of ∼10 ^−3 M _⊙ ^56 Ni is involved in the ejecta of a mass of ∼0.1 M _⊙ . Particularly, the magnetar parameters required by the fit are consistent with those derived from the early X-ray observation.

Wolf–Rayet stars in close binary systems can be tidally spun up by their companions, potentially leaving behind fast-spinning, highly magnetized neutron stars, known as “magnetars,” after core collapse. These newborn magnetars can transfer rotational energy into heating and accelerating the ejecta, producing hydrogen-poor superluminous supernovae (SLSNe). In this Letter, we propose that the magnetar wind of the newborn magnetar could significantly evaporate its companion star, typically a main-sequence or helium star, if the binary system is not disrupted by the abrupt mass loss and supernova (SN) kick. The subsequent heating and acceleration of the evaporated star material along with the SN ejecta by the magnetar wind can produce a postpeak bump in the SLSN lightcurve. Our model can reproduce the primary peaks and postpeak bumps of four example observed multiband SLSN lightcurves, revealing that the mass of the evaporated material could be ∼0.4–0.6 M _⊙ if the material is hydrogen-rich. We propose that the magnetar could induce strongly enhanced evaporation from its companion star near the pericenter if the orbit of the post-SN binary is highly eccentric, ultimately generating multiple postpeak bumps in the SLSN lightcurves. This “magnetar–star binary engine” model may offer a possible explanation for the evolution of polarization, along with the origin and velocity broadening of late-time hydrogen or helium broad spectral features observed in some bumpy SLSNe. The diversity in the lightcurves and spectra of SLSNe may be attributed to the wide variety of companion stars and post-SN binary systems.

In this work, we investigate an alternative channel for the formation of fast-spinning black hole–neutron star (BHNS) binaries, in which super-Eddington accretion is expected to occur in accreting BHs during the stable mass transfer phase within BH-stripped helium (BH–He-rich) star binary systems. We evolve intensive MESA grids of close-orbit BH–He-rich star systems to systematically explore the projected aligned spins of BHs in BHNS binaries, as well as the impact of different accretion limits on the tidal disruption probability and electromagnetic (EM) signature of BHNS mergers. Most of the BHs in BHNS mergers cannot be effectively spun up through accretion if the accretion rate is limited to $\lesssim 10\,{\dot{M}}_{\mathrm{Edd}}$ , where ${\dot{M}}_{\mathrm{Edd}}$ is the standard Eddington accretion limit. In order to reach high spins (e.g., χ _BH ≳ 0.5), the BHs are required to be born less massive (e.g., ≲3.0 M _⊙ ) in binary systems with initial periods of ≲0.2–0.3 days and accrete material at $\sim 100\,{\dot{M}}_{\mathrm{Edd}}$ . However, even under this high accretion limit, ≳6 M _⊙ BHs are typically challenging to significantly spin up and generate detectable associated EM signals. Our population simulations suggest that different accretion limits have a slight impact on the ratio of tidal disruption events. However, as the accretion limit increases, the EM counterparts from the cosmological BHNS population can become bright overall.

In this paper, we explore the formation of the mass-gap black hole-neutron star (mgBHNS) merger detected in gravitational wave (GW) event, i.e., GW230529, from the isolated binary evolution channel, and study potential signatures of its electromagnetic counterparts. By adopting the “delayed” supernova prescription and reasonable model realizations, our population synthesis simulation results can simultaneously match the rate densities of mgBHNS and total BHNS mergers inferred from the population analyses, along with the population distribution of the BH mass in BHNS mergers reported by the LIGO–Virgo–KAGRA Collaboration. Because GW230529 contributes significantly to the inferred mgBHNS rate densities, we suggest that GW230529 can be explained through the isolated binary evolution channel. Considering the AP4 (DD2) equation of state, the probability that GW230529 can make tidal disruption is 12.8% (63.2%). If GW230529 is a disrupted event, its kilonova peak apparent magnitude is predicted ∼23–24 mag, and hence, can be detected by the present survey projects and Large Synoptic Survey Telescope. Since GW230529 could be an off-axis event inferred from the GW observation, its associated gamma-ray burst (GRB) might be too dim to be observed by γ -ray detectors, interpreting the lack of GRB observations. Our study suggests the existence of mgBHNS mergers formed through the isolated binary evolution channel due to the discovery of GW230529, indicating that BHNS mergers are still likely to be multimessenger sources that emit GWs, GRBs, and kilonovae. Although mgBHNS mergers account for ∼50% of the cosmological BHNS population, we find that ≳90% of disrupted BHNS mergers are expected to originate from mgBHNS mergers.

The accretion disks of supermassive black holes (SMBHs) harboring in active galactic nuclei (AGN) are considered to be an ideal site for producing different types of gamma-ray bursts (GRBs). The detectability of these GRB phenomena hidden in AGN disks is highly dependent on the dynamical evolution of the GRB relativistic jets. By investigating the reverse- and forward-shock dynamics due to the interaction between the jets and AGN disk material, we find that the relativistic jets can successfully break out from the disks only for a sufficiently high luminosity and a long enough duration. In comparison, relatively normal GRB jets are inclined to be choked in the disks unless the GRBs occur near an SMBH with relatively low mass (e.g., ∼10 ^6 M _⊙ ). For the choked jets, unlike normal GRB prompt and afterglow emission, we can only expect to detect emission from the forward shock when the shock is very close to the edge of the disks, i.e., at the shock breakout emission and subsequent cooling of the shock.

In the second work of this series, we explore the optimal search strategy for serendipitous and gravitational-wave-triggered target-of-opportunity (ToO) observations of kilonovae and optical short-duration gamma-ray burst (sGRB) afterglows from binary neutron star (BNS) mergers, assuming that cosmological kilonovae are AT2017gfo-like (but with viewing-angle dependence) and that the properties of afterglows are consistent with those of cosmological sGRB afterglows. A one-day cadence serendipitous search strategy with an exposure time of ∼30 s can always achieve an optimal search strategy of kilonovae and afterglows for various survey projects. We show that the optimal detection rates of the kilonovae (afterglows) are ∼0.3/0.6/1/20 yr ^−1 (∼50/60/100/800 yr ^−1 ) for Zwicky the Transient Facility (ZTF)/Multi-channel Photometric Survey Telescope (Mephisto)/Wide Field Survey Telescope (WFST)/Large Synoptic Survey Telescope (LSST), respectively. A better search strategy for SiTian than the current design is to increase the exposure time. In principle, a fully built SiTian can detect ∼7(2000) yr ^−1 kilonovae (afterglows). Population properties of electromagnetic (EM) signals detected by serendipitous observations are studied in detail. For ToO observations, we predict that one can detect ∼11 yr ^−1 BNS gravitational wave (GW) events during the fourth observing run (O4) by considering an exact duty cycle of the third observing run. The median GW sky localization area is expected to be ∼10 deg ^2 for detectable BNS GW events. For O4, we predict that ZTF/Mephisto/WFST/LSST can detect ∼5/4/3/3 kilonovae (∼1/1/1/1 afterglows) per year, respectively. The GW detection rates, GW population properties, GW sky localizations, and optimistic ToO detection rates of detectable EM counterparts for BNS GW events at the Advanced Plus, LIGO Voyager, and ET&CE eras are detailedly simulated in this paper.

Despite the fact that the progenitor of fast blue optical transients (FBOTs) is still a subject of debate, FBOTs are sometimes suggested to originate from the core collapse of ultra-stripped stars and be powered by a spinning-down neutron star. Following this consideration, it is expected that the late-time evolution of the progenitor stars can leave important imprints in the circumstellar material (CSM) of FBOTs, due to the strong mass loss of the stars. The interaction of the FBOT ejecta with the CSM can drive a long-lasting shock to generate radio emission, which thus enables us to probe the CSM properties through radio observation although such observations are still rare. Within the framework of the magnetar-powered model, Liu et al. fitted the multi-band optical light curves of 40 FBOTs, and hence, the statistical distributions of the FBOT magnetar and ejecta parameters were obtained. Based on these FBOT population results, we investigate the dependence of the radio emission on the mass-loss rate of the progenitors and evaluate the detectability of radio emission from FBOTs with current and future telescopes. It is found that the distribution of the peak time and peak luminosity of the emission at 8.4 GHz are primarily in the regions of t _peak, _ν = 10 ^2.12±0.63 days and L _peak, _ν = 10 ^28.73±0.83 erg s ^−1 Hz ^−1 , respectively. A joint detection of the Zwicky Transient Facility and Very Large Array could achieve success in about 8.7% FBOTs of z ≤ 1. Furthermore, if considering a joint detection of the Chinese Space Station Telescope and the Square Kilometer Array, this rate of success could be increased to about 23.9%.

Various supernovae, compact object coalescences, and tidal disruption events are widely believed to occur embedded in active galactic nucleus (AGN) accretion disks and generate detectable electromagnetic signals. We collectively refer to them as AGN disk transients. The inelastic hadronuclear ( pp ) interactions between shock-accelerated cosmic rays and AGN disk materials shortly after the ejecta shock breaks out of the disk can produce high-energy neutrinos. However, the expected efficiency of neutrino production would decay rapidly by adopting a pure Gaussian density atmosphere profile applicable for stable gas-dominated disks. On the other hand, AGN outflows and disk winds are commonly found around AGN accretion disks. In this paper, we show that the circum-disk medium would further consume the shock kinetic energy to more efficiently produce high-energy neutrinos, especially for ∼ TeV−PeV neutrinos that IceCube detects. Thanks to the existence of the circum-disk medium, we find that the neutrino production will be enhanced significantly and make a much higher contribution to the diffuse neutrino background. Optimistically, ∼20% of the diffuse neutrino background can be contributed by AGN disk transients.

The formation of black hole–neutron star (BH–NS) or BH–BH systems may be accompanied by special supernova (SN) signals, due to the accretion feedback from the companion BH. The additional heating, which is mainly attributed to the Blandford–Payne mechanism, would disrupt the isotropic nature of the luminosity distribution on the surface of the SN ejecta, leading to the appearance of polarization. Here we develop a three-dimensional Monte Carlo polarization simulation code to conduct simulations for these special SNe. We find that the maximum polarization level of ∼2% occurs at the peak time of SN emission in the “close-binary” scenario, while in the “faraway-binary” case maximum polarization (i.e., ∼0.7%) is observed at a considerably later time than the peak of the SN. The magnitude of polarization is dependent on the degree of unevenness in the luminosity distribution and the angle between the line of sight and the equatorial direction. When considering the geometric distortion of SN ejecta at the same time, the magnitude of polarization may either increase (for a oblate ellipsoidal shape) or decrease (for a prolate ellipsoidal shape). The polarization signatures represent a promising auxiliary instrument to facilitate the identification of the companion-fed SNe. Moreover, comparing the event rate of these special SNe with the event rate density of LIGO-Virgo-detected BH–NS/BH systems could further help to distinguish the BH–NS/BH formation channel.

The "Breakthrough Starshot" program is planning to send transrelativistic probes to travel to nearby stellar systems within decades. Because the probe velocity is designed to be a good fraction of the light speed, Zhang & Li recently proposed that these transrelativistic probes can be used to study astronomical objects and to test special relativity. In this work, we propose some methods to test special relativity and constrain photon mass using the Doppler effect with the images and spectral features of astronomical objects as observed in the transrelativistic probes. We introduce more general theories to set up the framework of testing special relativity, including the parametric general Doppler effect and the Doppler effect with massive photons. We find that by comparing the spectra of a certain astronomical object, one can test Lorentz invariance and constrain photon mass. Additionally, using the imaging and spectrograph capabilities of transrelativistic probes, one can test time dilation and constrain photon mass. For a transrelativistic probe with velocity v ∼ 0.2c, aperture D ∼ 3.5 cm, and spectral resolution R ∼ 100 (or 1000), we find that the probe velocity uncertainty can be constrained to σv ∼ 0.01c (or 0.001c), and the time dilation factor uncertainty can be constrained to (or 0.001), where is the time dilation factor and γ is the Lorentz factor. Meanwhile, the photon mass limit is set to mγ ≲ 10−33 g, which is slightly lower than the energy of the optical photon. [ABSTRACT FROM AUTHOR]

By fitting the bolometric light curves of 31 superluminous supernovae (SLSNe) with the magnetar engine model, we derive the ejecta masses and magnetar parameters for these SLSNe. The lower boundary of magnetic field strengths of SLSN magnetars can be set just around the critical field strength of the electron Landau quantization. In more details, SLSN magnetars can further be divided into two subclasses of magnetic fields of and , respectively. It is revealed that these two subclasses of magnetars are just associated with the slow-evolving and fast-evolving bolometric light curves of SLSNe. In comparison, the magnetars harbored in gamma-ray bursts (GRBs) and associated hypernovae are usually inferred to have much higher magnetic fields with a lower boundary of about . This robustly suggests that it is the magnetic fields that play the crucial role in distinguishing SLSNe from GRBs/hypernovae. The rotational energy of SLSN magnetars is found to be correlated with the masses of supernova ejecta. This possible correlation provides a clue to explore the nature of their progenitors. Moreover, the distribution of ejecta masses of SLSNe is basically intermediate between those of normal core-collapse supernovae and hypernovae. This could indicate an intrinsic connection among these different stellar explosions. [ABSTRACT FROM AUTHOR]