Search

The neutrino-driven wind cooling phase of proto-neutron stars (PNSs) follows successful supernovae. Wind models without magnetic fields or rotation fail to achieve the necessary conditions for production of the third r -process peak, but robustly produce a weak r- process in neutron-rich winds. Using 2D magnetohydrodynamic simulations with magnetar-strength magnetic fields and rotation, we show that the PNS rotation rate significantly affects the thermodynamic conditions of the wind. We show that high-entropy material is quasiperiodically ejected from the closed zone of the PNS magnetosphere with the required thermodynamic conditions to produce heavy elements. We show that maximum entropy S of the material ejected depends systematically on the magnetar spin period P ⋆ and scales as S ∝ P ⋆ − 5 / 6 for sufficiently rapid rotation. We present results from simulations at a constant neutrino luminosity representative of ∼1–2 s after the onset of cooling for P ⋆ ranging from 5–200 ms and a few simulations with evolving neutrino luminosity where we follow the evolution of the magnetar wind until 10–14 s after the onset of cooling. We estimate at magnetar polar magnetic field strength B 0 = 3 × 1015 G and 1015 G that neutron-rich magnetar winds can, respectively, produce at least ∼1–5 × 10−5 M ⊙ and ∼1–4 × 10−7 M ⊙ of material with the required parameters for synthesis of the third r -process peak, within 1–2 s and 10 s, respectively, in that order after the onset of cooling. We show that proton-rich magnetar winds can have favorable conditions for production of p- nuclei, even at a modest B 0 = 5 × 1014 G. [ABSTRACT FROM AUTHOR]

The fallback disk model is widely accepted to explain long-period neutron stars (NSs) that cannot be simulated by magnetic dipole radiation. However, there was no confirmed detection of disks from the newly discovered long-period pulsars GLEAM-X 162759.5-523504.3 and GPM J1839-10 and the slowest known isolated NS 1E 161348-5055. This might be because the disks have either been in a noninteracting/inactive state where their emission is too weak to be detected or have been disrupted. In this work, we conduct simulations to examine the lifetime of supernova fallback disks around isolated NSs. We assume that the disks' mass varies in a self-similar way, and their interaction with the NS occurs only in an interacting/active state. Our results reveal that nearly all the interacting lifetimes for the disks are shorter than 105 yr, while the existence lifetimes are considerably longer. [ABSTRACT FROM AUTHOR]

We present the semi-analytical light curve modelling of 13 supernovae associated with gamma-ray bursts (GRB-SNe) along with two relativistic broad-lined (Ic-BL) SNe without GRB association (SNe 2009bb and 2012ap), considering millisecond magnetars as central-engine-based power sources for these events. The bolometric light curves of all 15 SNe in our sample are well-regenerated utilizing a χ2-minimization code, MINIM , and numerous parameters are constrained. The median values of ejecta mass (M ej), magnetar's initial spin period (P i), and magnetic field (B) for GRB-SNe are determined to be ≈5.2 M⊙, 20.5 ms, and 20.1 × 1014 G, respectively. We leverage machine learning (ML) algorithms to comprehensively compare the three-dimensional parameter space encompassing M ej, P i, and B for GRB-SNe determined herein to those of H-deficient superluminous SNe (SLSNe-I), fast blue optical transients (FBOTs), long GRBs (LGRBs), and short GRBs (SGRBs) obtained from the literature. The application of unsupervized ML clustering algorithms on the parameters M ej, P i, and B for GRB-SNe, SLSNe-I, and FBOTs yields a classification accuracy of ∼95 per cent. Extending these methods to classify GRB-SNe, SLSNe-I, LGRBs, and SGRBs based on P i and B values results in an accuracy of ∼84 per cent. Our investigations show that GRB-SNe and relativistic Ic-BL SNe presented in this study occupy different parameter spaces for M ej, P i, and B than those of SLSNe-I, FBOTs, LGRBs, and SGRBs. This indicates that magnetars with different P i and B can give birth to distinct types of transients. [ABSTRACT FROM AUTHOR]

Core-collapse supernovae (CCSNe) are considered the primary magnetar formation channel, with 15 magnetars associated with supernova remnants (SNRs). A large fraction of these should occur in massive stellar binaries that are disrupted by the explosion, meaning that |$\sim 45~{{\ \rm per\ cent}}$| of magnetars should be nearby high-velocity stars. Here, we conduct a multiwavelength search for unbound stars, magnetar binaries, and SNR shells using public optical (uvgrizy bands), infrared (J, H, K , and Ks bands), and radio (888 MHz, 1.4 GHz, and 3 GHz) catalogues. We use Monte Carlo analyses of candidates to estimate the probability of association with a given magnetar based on their proximity, distance, proper motion, and magnitude. In addition to recovering a proposed magnetar binary, a proposed unbound binary, and 13 of 15 magnetar SNRs, we identify two new candidate unbound systems: an OB star from the Gaia catalogue we associate with SGR J1822.3−1606, and an X-ray pulsar we associate with 3XMM J185246.6 + 003317. Using a Markov Chain Monte Carlo simulation that assumes all magnetars descend from CCSNe, we constrain the fraction of magnetars with unbound companions to |$5\lesssim f_u \lesssim 24~{{\ \rm per\ cent}}$|⁠ , which disagrees with neutron star population synthesis results. Alternate formation channels are unlikely to wholly account for the lack of unbound binaries as this would require |$31\lesssim f_{nc} \lesssim 66~{{\ \rm per\ cent}}$| of magnetars to descend from such channels. Our results support a high fraction (⁠|$48\lesssim f_m \lesssim 86~{{\ \rm per\ cent}}$|⁠) of pre-CCSN mergers, which can amplify fossil magnetic fields to preferentially form magnetars. [ABSTRACT FROM AUTHOR]

SN 2020zbf is a hydrogen-poor superluminous supernova (SLSN) at z = 0.1947 that shows conspicuous C II features at early times, in contrast to the majority of H-poor SLSNe. Its peak magnitude is Mg = −21.2 mag and its rise time (≲26.4 days from first light) places SN 2020zbf among the fastest rising type I SLSNe. We used spectra taken from ultraviolet (UV) to near-infrared wavelengths to identify spectral features. We paid particular attention to the C II lines as they present distinctive characteristics when compared to other events. We also analyzed UV and optical photometric data and modeled the light curves considering three different powering mechanisms: radioactive decay of 56Ni, magnetar spin-down, and circumstellar medium (CSM) interaction. The spectra of SN 2020zbf match the model spectra of a C-rich low-mass magnetar-powered supernova model well. This is consistent with our light curve modeling, which supports a magnetar-powered event with an ejecta mass Mej = 1.5 M⊙. However, we cannot discard the CSM-interaction model as it may also reproduce the observed features. The interaction with H-poor, carbon-oxygen CSM near peak light could explain the presence of C II emission lines. A short plateau in the light curve around 35–45 days after peak, in combination with the presence of an emission line at 6580 Å, can also be interpreted as being due to a late interaction with an extended H-rich CSM. Both the magnetar and CSM-interaction models of SN 2020zbf indicate that the progenitor mass at the time of explosion is between 2 and 5 M⊙. Modeling the spectral energy distribution of the host galaxy reveals a host mass of 108.7M⊙, a star formation rate of 0.24−0.12+0.41M⊙ yr−1, and a metallicity of ∼0.4 Z⊙. [ABSTRACT FROM AUTHOR]

The development of the latest generation of Imaging Atmospheric Cherenkov Telescopes (IACTs) over recent decades has led to the discovery of new extreme astrophysical phenomena in the very-high-energy (VHE, E > 100 GeV) gamma-ray regime. Time-domain and multi-messenger astronomy are inevitably connected to the physics of transient VHE emitters, which show unexpected (and mostly unpredictable) flaring or exploding episodes at different timescales. These transients often share the physical processes responsible for the production of the gamma-ray emission, through cosmic-ray acceleration, magnetic reconnection, jet production and/or outflows, and shocks interactions. In this review, we present an up-to-date overview of the VHE transients field, spanning from novae to supernovae, neutrino counterparts or fast radio bursts, among others, and we outline the expectations for future facilities. [ABSTRACT FROM AUTHOR]

After a successful supernova, a proto-neutron star (PNS) cools by emitting neutrinos on ∼1–100 s time-scales. Provided that there are neutrino emission 'hotspots' or 'cold-spots' on the surface of the rotating PNS, we can expect a periodic modulation in the number of neutrinos observable by detectors. We show that Fourier transform techniques can be used to determine the PNS rotation rate from the neutrino arrival times. Provided there is no spin-down, a 1-parameter Discrete Fourier Transform (DFT) is sufficient to determine the spin period of the PNS. If the PNS is born as a magnetar with polar magnetic field strength B 0 ≳ 1015 G and is 'slowly' rotating with an initial spin period ≳100 ms, then it can spin-down to periods of the order of seconds during the cooling phase. We propose a modified DFT technique with three frequency parameters to detect spin-down. Due to lack of neutrino data from a nearby supernova except the ∼20 neutrinos detected from SN1987A, we use toy models and one physically motivated modulating function to generate neutrino arrival times. We use the false alarm rate (FAR) to quantify the significance of the Fourier power spectrum peaks. We show that PNS rotation and spin-down are detected with |$\rm FAR\,\lt\, 2~{{\ \rm per\ cent}}$| (2σ) for periodic signal content |$\rm M\gtrsim 13-15~{{\ \rm per\ cent}}$| if 5 × 103 neutrinos are detected in ∼3 s and with |$\rm FAR\,\lt\, 1{{\ \rm per\ cent}}$| for |$\rm M\,\ge 5{{\ \rm per\ cent}}$| if 5 × 104 neutrinos are detected in ∼3 s. Since we can expect ∼104−105 neutrino detections from a supernova at 10 kpc, detection of PNS rotation and spin-down is possible using the neutrinos from the next Galactic supernova. [ABSTRACT FROM AUTHOR]

SN 2019neq was a very fast evolving superluminous supernova. At a redshift z  = 0.1059, its peak absolute magnitude was −21.5 ± 0.2 mag in g band. In this work, we present data and analysis from an extensive spectrophotometric follow-up campaign using multiple observational facilities. Thanks to a nebular spectrum of SN 2019neq, we investigated some of the properties of the host galaxy at the location of SN 2019neq and found that its metallicity and specific star formation rate are in a good agreement with those usually measured for SLSNe-I hosts. We then discuss the plausibility of the magnetar and the circumstellar interaction scenarios to explain the observed light curves, and interpret a nebular spectrum of SN 2019neq using published sumo radiative-transfer models. The results of our analysis suggest that the spin-down radiation of a millisecond magnetar with a magnetic field |$B\simeq 6\times 10^{14}\, \mathrm{G}$| could boost the luminosity of SN 2019neq. [ABSTRACT FROM AUTHOR]

Several types of energetic supernovae, such as superluminous supernovae (SLSNe) and broad-line Ic supernovae (Ic-BL SNe), could be powered by the spin-down of a rapidly rotating magnetar. Currently, most models used to infer the parameters for potential magnetar-driven supernovae make several unsuitable assumptions that likely bias the estimated parameters. In this work, we present a new model for magnetar-driven supernovae that relaxes several of these assumptions and an inference workflow that enables accurate estimation of parameters from light curves of magnetar-driven supernovae. In particular, in this model, we include the dynamical evolution of the ejecta, coupling it to the energy injected by the magnetar itself while also allowing for non-dipole spin down. We show that the model can reproduce SLSN and Ic-BL SN light curves consistent with the parameter space from computationally expensive numerical simulations. We also show the results of parameter inference on four well-known example supernovae, demonstrating the model's effectiveness at capturing the considerable diversity in magnetar-driven supernova light curves. The model fits each light curve well and recovers parameters broadly consistent with previous works. This model will allow us to explore the full diversity of magnetar-driven supernovae under one theoretical framework, more accurately characterize these supernovae from only photometric data, and make more accurate predictions of future multiwavelength emission to test the magnetar-driven scenario better. [ABSTRACT FROM AUTHOR]

We present a first 3D magnetohydrodynamic (MHD) simulation of oxygen, neon, and carbon shell burning in a rapidly rotating |$16\hbox{-}\mathrm{M}_\odot$| core-collapse supernova progenitor. We also run a purely hydrodynamic simulation for comparison. After |$\mathord \approx 180\mathrm{s}$| (⁠|$\mathord \approx$| 15 and 7 convective turnovers, respectively), the magnetic fields in the oxygen and neon shells achieve saturation at 1011 and 5 × 1010 G. The strong Maxwell stresses become comparable to the radial Reynolds stresses and eventually suppress convection. The suppression of mixing by convection and shear instabilities results in the depletion of fuel at the base of the burning regions, so that the burning shell eventually move outward to cooler regions, thus reducing the energy generation rate. The strong magnetic fields efficiently transport angular momentum outwards, quickly spinning down the rapidly rotating convective oxygen and neon shells and forcing them into rigid rotation. The hydrodynamic model shows complicated redistribution of angular momentum and develops regions of retrograde rotation at the base of the convective shells. We discuss implications of our results for stellar evolution and for the subsequent core-collapse supernova. The rapid redistribution of angular momentum in the MHD model casts some doubt on the possibility of retaining significant core angular momentum for explosions driven by millisecond magnetars. However, findings from multidimensional models remain tentative until stellar evolution calculations can provide more consistent rotation profiles and estimates of magnetic field strengths to initialize multidimensional simulations without substantial numerical transients. We also stress the need for longer simulations, resolution studies, and an investigation of non-ideal effects. [ABSTRACT FROM AUTHOR]

Rapidly rotating magnetars have been associated with gamma-ray bursts (GRBs) and superluminous supernovae (SLSNe). Using a suite of two-dimensional magnetohydrodynamic simulations at fixed neutrino luminosity and a couple of evolutionary models with evolving neutrino luminosity and magnetar spin period, we show that magnetars are viable central engines for powering GRBs and SLSNe. We also present analytical estimates of the energy outflow rate from the proto-neutron star (PNS) as a function of polar magnetic field strength B 0, PNS angular velocity Ω⋆, PNS radius R ⋆, and mass outflow rate |$\dot{M}$|⁠. We show that rapidly rotating magnetars with spin periods P ⋆ ≲ 4 ms and polar magnetic field strength B 0 ≳ 1015 G can release 1050 to 5 × 1051 erg of energy during the first ∼2 s of the cooling phase. Based on this result, it is plausible that sustained energy injection by magnetars through the relativistic wind phase can power GRBs. We also show that magnetars with moderate field strengths of B 0 ≲ 5 × 1014 G do not release a large fraction of their rotational kinetic energy during the cooling phase and, hence, are not likely to power GRBs. Although we cannot simulate to times greater than ∼3–5 s after a supernova, we can hypothesize that moderate field strength magnetars can brighten the supernova light curves by releasing their rotational kinetic energy via magnetic dipole radiation on time-scales of days to weeks, since these do not expend most of their rotational kinetic energy during the early cooling phase. [ABSTRACT FROM AUTHOR]

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.

AT2018cow is the most extensively observed and widely studied fast blue optical transient to date; its unique observational properties challenge all existing standard models. In this paper, we model the luminosity evolution of the optical, soft X-ray, and hard X-ray emission, as well as the X-ray spectrum of AT2018cow with a magnetar-centered engine model. We consider a two-zone model with a striped magnetar wind in the interior and an expanding ejecta outside. The soft and hard X-ray emission of AT2018cow can be explained by the leakage of high-energy photons produced by internal gradual magnetic dissipation in the striped magnetar wind, while the luminous thermal UV/optical emission results from the thermalization of the ejecta by the captured photons. The two-component energy spectrum yielded by our model with a quasi-thermal component from the optically thick region of the wind superimposed on an optically thin synchrotron component well reproduces the X-ray spectral shape of AT2018cow. The Markov Chain Monte Carlo fitting results suggest that in order to explain the very short rise time to peak of the thermal optical emission, a low ejecta mass M _ej ≈ 0.1 M _⊙ and high ejecta velocity ${v}_{\mathrm{SN}}\approx 0.17c$ are required. A millisecond magnetar with P _0 ≈ 3.7 ms and B _p ≈ 2.4 × 10 ^14 G is needed to serve as the central engine of AT2018cow.

The emission from shock breakouts (SBOs) represents the earliest electromagnetic (EM) signal emitted by cataclysmic events involving the formation or the merger of neutron stars (NSs). As such, SBOs carry unique information on the structure of their progenitors and on the explosion energy. The characteristic SBO emission is expected in the UV range, and its detection is one of the key targets of the ULTRASAT satellite. Among SBO sources, we focus on a specific class involving the formation of fast-spinning magnetars in the core-collapse of massive stars. Fast-spinning magnetars are expected to produce a specific signature in the early UV supernova light curve, powered by the extra spin energy quickly released by the NS. Moreover, they are considered as optimal candidates for the emission of long-transient gravitational wave (GW) signals, the detection of which requires early EM triggers to boost the sensitivity of dedicated GW search pipelines. We calculate early supernova UV light curves in the presence of a magnetar central engine, as a function of the explosion energy, ejecta mass, and magnetar parameters. We then estimate the ULTRASAT detection horizon (z < 0.15) as a function of the same physical parameters, and the overall expected detection rate, finding that magnetar-powered SBOs may represent up to 1/5 of the total events detected by ULTRASAT. Moreover, at the expected sensitivity of the LIGO/Virgo/Kagra O5 science run, one such event occurring within 5 Mpc will provide an ideal trigger for a GW long-transient search. Future GW detectors like the Einstein Telescope will push the horizon for joint EM-GW detections to 35–40 Mpc. [ABSTRACT FROM AUTHOR]

We present observations of a peculiar hydrogen- and helium-poor stripped-envelope (SE) supernova (SN) 2020wnt, primarily in the optical and near-infrared (near-IR). Its peak absolute bolometric magnitude of −20.9 mag (L bol, peak = (6.8 ± 0.3) × 1043 erg s−1) and a rise time of 69 days are reminiscent of hydrogen-poor superluminous SNe (SLSNe I), luminous transients potentially powered by spinning-down magnetars. Before the main peak, there is a brief peak lasting <10 days post explosion, likely caused by interaction with circumstellar medium (CSM) ejected ∼years before the SN explosion. The optical spectra near peak lack a hot continuum and O ii absorptions, which are signs of heating from a central engine; they quantitatively resemble those of radioactivity-powered hydrogen/helium-poor Type Ic SESNe. At ∼1 yr after peak, nebular spectra reveal a blue pseudo-continuum and narrow O i recombination lines associated with magnetar heating. Radio observations rule out strong CSM interactions as the dominant energy source at +266 days post peak. Near-IR observations at +200–300 days reveal carbon monoxide and dust formation, which causes a dramatic optical light-curve dip. Pair-instability explosion models predict slow light curve and spectral features incompatible with observations. SN 2020wnt is best explained as a magnetar-powered core-collapse explosion of a 28 M ⊙ pre-SN star. The explosion kinetic energy is significantly larger than the magnetar energy at peak, effectively concealing the magnetar-heated inner ejecta until well after peak. SN 2020wnt falls into a continuum between normal SNe Ic and SLSNe I, and demonstrates that optical spectra at peak alone cannot rule out the presence of a central engine. [ABSTRACT FROM AUTHOR]

Context. Fast radio bursts (FRBs) are millisecond radio transients observed at cosmological distances. The nature of their progenitors is still a matter of debate, although magnetars are invoked by most models. The proposed FRB–magnetar connection was strengthened by the discovery of an FRB-like event from the Galactic magnetar SGR J1935+2154. Aims. In this work we aim to investigate how prevalent magnetars such as SGR J1935+2154 are within FRB progenitors. Methods. To this end, we carried out an FRB search in a sample of seven nearby (< 12 Mpc) galaxies with the Northern Cross Radio Telescope for a total of 692 h. Results. We detected one 1.8 ms burst in the direction of M 101 with a fluence of 58 ± 5 Jy ms. Its dispersion measure of 303 pc cm−3 places it most likely beyond M 101. Considering that no significant detection comes indisputably from the selected galaxies, we place a 38 yr−1 upper limit on the total burst rate (i.e. including the whole sample) at the 95% confidence level. This upper limit constrains the event rate per magnetar to λmag < 0.42 magnetar−1 yr−1 or, if combined with literature observations of a similar sample of nearby galaxies, it yields a joint constraint of λmag < 0.25 magnetar−1 yr−1. We also provide the first constraints on the expected rate of FRBs hypothetically originating from ultra-luminous X-ray (ULX) sources, since some of the galaxies observed during our observational campaign host confirmed ULXs. We obtain < 13 yr−1 per ULX for the total sample of galaxies observed. Conclusions. Our results indicate that bursts with energies E > 1034 erg from magnetars such as SGR J1935+2154 appear more rarely compared to previous observations and further disfavour them as unique progenitors for the cosmological FRB population. This provides support to the idea that there is a greater contribution from a population of more exotic magnetars not born via core-collapsed supernovae. [ABSTRACT FROM AUTHOR]

We present near-infrared (NIR) and optical observations of the Type Ic supernova (SN Ic) SN 2021krf obtained between days 13 and 259 at several ground-based telescopes. The NIR spectrum at day 68 exhibits a rising K -band continuum flux density longward of ∼2.0 μ m, and a late-time optical spectrum at day 259 shows strong [O i ] 6300 and 6364 Å emission-line asymmetry, both indicating the presence of dust, likely formed in the SN ejecta. We estimate a carbon-grain dust mass of ∼2 × 10−5 M ⊙ and a dust temperature of ∼900–1200 K associated with this rising continuum and suggest the dust has formed in SN ejecta. Utilizing the one-dimensional multigroup radiation-hydrodynamics code STELLA, we present two degenerate progenitor solutions for SN 2021krf, characterized by C–O star masses of 3.93 and 5.74 M ⊙, but with the same best-fit 56Ni mass of 0.11 M ⊙ for early times (0–70 days). At late times (70–300 days), optical light curves of SN 2021krf decline substantially more slowly than those expected from 56Co radioactive decay. Lack of H and He lines in the late-time SN spectrum suggests the absence of significant interaction of the ejecta with the circumstellar medium. We reproduce the entire bolometric light curve with a combination of radioactive decay and an additional powering source in the form of a central engine of a millisecond pulsar with a magnetic field smaller than that of a typical magnetar. [ABSTRACT FROM AUTHOR]

New optical photometric, spectroscopic, and imaging polarimetry data are combined with publicly available data to study some of the physical properties of the two hydrogen-poor superluminous supernovae (SLSNe) SN 2021bnw and SN 2021fpl. For each SLSN, the best-fitting parameters obtained from the magnetar model with Modular Open-Source Fitter for Transients do not depart from the range of parameter obtained on other SLSNe discussed in the literature. A spectral analysis with SYN++  shows that SN 2021bnw is a W type, fast evolver, while SN 2021fpl is a 15bn type, slow evolver. The analysis of the polarimetry data obtained on SN 2021fpl at four epochs (+1.8, +20.6, +34.1, and +43.0 d, rest frame) shows >3σ polarization detections in the range of 0.8–1 per cent. A comparison of the spectroscopy data suggests that SN 2021fpl underwent a spectral transition a bit earlier than SN 2015bn, during which, similarly, it could have underwent a polarization transition. The analysis of the polarimetry data obtained on SN 2021bnw does not show any departure from symmetry of the photosphere at an empirical diffusion time-scale of ≈2 (+81.1 d rest frame). This result is consistent with those on the sample of W-type SLSN observed at empirical diffusion time-scale ≤ 1 with that technique, even though it is not clear the effect of limited spectral windows varying from one object to the other. Measurements at higher empirical diffusion time-scale may be needed to see any departure from symmetry as it is discussed in the literature for SN 2017egm. [ABSTRACT FROM AUTHOR]

Fast X-ray Transients (FXTs) are X-ray flares with durations ranging from a few hundred seconds to a few hours. Possible origins include the tidal disruption of a white dwarf by an intermediate-mass black hole, a supernova shock breakout, or a binary neutron star merger. We present the X-ray light curve and spectrum as well as deep optical imaging of the FXT XRT 210423, which has been suggested to be powered by a magnetar produced in a binary neutron star merger. Our Very Large Telescope and Gran Telescopio Canarias (GTC) observations began on 2021 May 6, thirteen days after the onset of the flare. No transient optical counterpart is found in the 1.″ (3 σ) X-ray uncertainty region of the source to a depth g s = 27.0 AB mag. (We use the word "counterpart" for any transient light in a wave band other than the original X-ray detection wave band, whereas the word "host" refers to the host galaxy.) A candidate host lies within the 1.″ X-ray uncertainty region with a magnitude of 25.9 ± 0.1 in the GTC/HiPERCAM g s filter. Due to its faintness, it was not detected in other bands, precluding a photometric redshift determination. We detect two additional candidate host galaxies: one with z spec = 1.5082 ± 0.0001 and an offset of 4.″2 ± 1.″ (37 ± 9 kpc) from the FXT, and another one with z phot = 1.04 − 0.14 + 0.22 and an offset of 3.″6 ± 1.″ (30 ± 8 kpc). Based on the properties of all the prospective hosts, we favor a binary neutron star merger, as previously suggested in the literature, as the explanation for XRT 210423. [ABSTRACT FROM AUTHOR]

The gravitational collapse of rapidly rotating massive stars can lead to the onset of the low T /| W | instability within the central proto-neutron star (PNS), which leaves strong signatures in both the gravitational wave (GW) and neutrino emission. Strong large-scale magnetic fields are usually invoked to explain outstanding stellar explosions of rapidly rotating progenitors, but their impact on the growth of such instability has not yet been cleared. We analyse a series of three-dimensional magnetohydrodynamic models to characterize the effects of different magnetic configurations on the development of the low T /| W | and the related multimessenger features. In the absence of magnetic fields, we observe the growth on dynamical time-scales of the low T /| W |, associated with a strong burst of GW and a correlated modulation of the neutrino emission. However, models with a strong magnetic field show a quenching of the low T /| W |, due to a flattening of the rotation profile in the first ∼100 ms after shock formation caused by the magnetic transport of angular momentum. The associated GW emission is weakened by an order of magnitude, exhibits a broader spectral shape, and has no dominant feature associated with the PNS large-scale oscillation modes. Neutrino luminosities are damped along the equatorial plane due to a more oblate PNS, and the only clear modulation in the signal is due to Standing Accretion Shock Instability activity. Finally, magnetized models produce lower luminosities for ν e than for |$\bar{\nu }_e$|⁠ , which is connected to a higher concentration of neutron-rich material in the PNS surroundings. [ABSTRACT FROM AUTHOR]

The neutrino-driven wind cooling phase of proto-neutron stars (PNSs) follows successful supernovae. Wind models without magnetic fields or rotation fail to achieve the necessary conditions for production of the third r -process peak, but robustly produce a weak r- process in neutron-rich winds. Using 2D magnetohydrodynamic simulations with magnetar-strength magnetic fields and rotation, we show that the PNS rotation rate significantly affects the thermodynamic conditions of the wind. We show that high-entropy material is quasiperiodically ejected from the closed zone of the PNS magnetosphere with the required thermodynamic conditions to produce heavy elements. We show that maximum entropy S of the material ejected depends systematically on the magnetar spin period P _⋆ and scales as $S\propto {P}_{\star }^{-5/6}$ for sufficiently rapid rotation. We present results from simulations at a constant neutrino luminosity representative of ∼1–2 s after the onset of cooling for P _⋆ ranging from 5–200 ms and a few simulations with evolving neutrino luminosity where we follow the evolution of the magnetar wind until 10–14 s after the onset of cooling. We estimate at magnetar polar magnetic field strength B _0 = 3 × 10 ^15 G and 10 ^15 G that neutron-rich magnetar winds can, respectively, produce at least ∼1–5 × 10 ^−5 M _⊙ and ∼1–4 × 10 ^−7 M _⊙ of material with the required parameters for synthesis of the third r -process peak, within 1–2 s and 10 s, respectively, in that order after the onset of cooling. We show that proton-rich magnetar winds can have favorable conditions for production of p- nuclei, even at a modest B _0 = 5 × 10 ^14 G.

Rapidly spinning magnetars are potential candidates for the energy source of supernovae (SNe) and gamma-ray bursts and the most promising sources for continuous gravitational waves (GWs) detected by ground-based GW detectors. Continuous GWs can be radiated from magnetars due to magnetic-induced deformation or fluid oscillations, compatible with magnetic dipole (MD) radiation for spin-down energy. In this paper, we investigate the diverse light curves of magnetar-driven SNe in the scenario that the spin-down is dominated by GW radiation and/or MD radiation. By simulating the light curves of SNe and employing the Markov Chain Monte Carlo method, we constrain the parameters of the magnetars and SN explosions and show that the signature of GW radiation may be indicated by the bolometric luminosity curves of SNe Ic-BL 2007ru and 2009bb. We find that the ellipticity of magnetars in the order of 10 ^−3 can be induced by the magnetic field of ∼10 ^16 G. If such continuous GWs associated with SNe can be detected in the future by the Advanced LIGO and Virgo detectors, this would be a smoking gun for a magnetar engine powering SNe.

We investigate the impact of strong initial magnetic fields in core-collapse supernovae of non-rotating progenitors by simulating the collapse and explosion of a |$16.9\, \mathrm{M}_\odot$| star for a strong- and weak-field case assuming a twisted-torus field with initial central field strengths of |${\approx }10^{12}$| and |${\approx }10^{6}\, \mathrm{G}$|⁠. The strong-field model has been set up with a view to the fossil-field scenario for magnetar formation and emulates a pre-collapse field configuration that may occur in massive stars formed by a merger. This model undergoes shock revival already |$100\, \mathrm{ms}$| after bounce and reaches an explosion energy of |$9.3\times 10^{50}\, \mathrm{erg}$| at |$310\, \mathrm{ms}$|⁠ , in contrast to a more delayed and less energetic explosion in the weak-field model. The strong magnetic fields help trigger a neutrino-driven explosion early on, which results in a rapid rise and saturation of the explosion energy. Dynamically, the strong initial field leads to a fast build-up of magnetic fields in the gain region to 40 per cent of kinetic equipartition and also creates sizable pre-shock ram pressure perturbations that are known to be conducive to asymmetric shock expansion. For the strong-field model, we find an extrapolated neutron star kick of |${\approx }350\, \mathrm{km}\, \mathrm{s}^{-1}$|⁠ , a spin period of |${\approx }70\, \mathrm{ms}$|⁠ , and no spin-kick alignment. The dipole field strength of the proto-neutron star is |$2\times 10^{14}\, \mathrm{G}$| by the end of the simulation with a declining trend. Surprisingly, the surface dipole field in the weak-field model is stronger, which argues against a straightforward connection between pre-collapse fields and the birth magnetic fields of neutron stars. [ABSTRACT FROM AUTHOR]

Stripped-envelope core-collapse supernovae can be divided into two broad classes: the common Type Ib/c supernovae (SNe Ib/c), powered by the radioactive decay of 56Ni, and the rare superluminous supernovae (SLSNe), most likely powered by the spin-down of a magnetar central engine. Up to now, the intermediate regime between these two populations has remained mostly unexplored. Here, we present a comprehensive study of 40 luminous supernovae (LSNe), SNe with peak magnitudes of M r = −19 to −20 mag, bound by SLSNe on the bright end and by SNe Ib/c on the dim end. Spectroscopically, LSNe appear to form a continuum between Type Ic SNe and SLSNe. Given their intermediate nature, we model the light curves of all LSNe using a combined magnetar plus radioactive decay model and find that they are indeed intermediate, not only in terms of their peak luminosity and spectra, but also in their rise times, power sources, and physical parameters. We subclassify LSNe into distinct groups that are either as fast evolving as SNe Ib/c or as slow evolving as SLSNe, and appear to be either radioactively or magnetar powered, respectively. Our findings indicate that LSNe are powered by either an overabundant production of 56Ni or by weak magnetar engines, and may serve as the missing link between the two populations. [ABSTRACT FROM AUTHOR]

In the seconds following their formation in core-collapse supernovae, 'proto'-magnetars drive neutrino-heated magnetocentrifugal winds. Using a suite of two-dimensional axisymmetric magnetohydrodynamic simulations, we show that relatively slowly rotating magnetars with initial spin periods of P ⋆0 = 50–500 ms spin down rapidly during the neutrino Kelvin–Helmholtz cooling epoch. These initial spin periods are representative of those inferred for normal Galactic pulsars, and much slower than those invoked for gamma-ray bursts and superluminous supernovae. Since the flow is non-relativistic at early times, and because the Alfvén radius is much larger than the proto-magnetar radius, spin-down is millions of times more efficient than the typically used dipole formula. Quasi-periodic plasmoid ejections from the closed zone enhance spin-down. For polar magnetic field strengths B 0 ≳ 5 × 1014 G, the spin-down time-scale can be shorter than the Kelvin–Helmholtz time-scale. For B 0 ≳ 1015 G, it is of the order of seconds in early phases. We compute the spin evolution for cooling proto-magnetars as a function of B 0, P ⋆0, and mass (M). Proto-magnetars born with B 0 greater than |$\simeq 1.3\times 10^{15}\, {\rm \, G}\, (P_{\star 0}/{400\, \rm \, ms})^{-1.4}(M/1.4\, {\rm M}_\odot)^{2.2}$| spin down to periods >1 s in just the first few seconds of evolution, well before the end of the cooling epoch and the onset of classic dipole spin-down. Spin-down is more efficient for lower M and for larger P ⋆0. We discuss the implications for observed magnetars, including the discrepancy between their characteristic ages and supernova remnant ages. Finally, we speculate on the origin of 1E 161348−5055 in the remnant RCW 103, and the potential for other ultra-slowly rotating magnetars. [ABSTRACT FROM AUTHOR]

A non-negligible fraction of binary neutron star mergers are expected to form long-lived neutron star remnants, dramatically altering the multimessenger signatures of a merger. Here, we extend existing models for magnetar-driven kilonovae and explore the diversity of kilonovae and kilonova afterglows. Focusing on the role of the (uncertain) magnetic field strength, we study the resulting electromagnetic signatures as a function of the external dipolar and internal toroidal fields. These two parameters govern, respectively, the competition between magnetic-dipole spin-down and gravitational-wave spin-down (due to magnetic-field deformation) of the rapidly rotating remnant. We find that even in the parameter space where gravitational-wave emission is dominant, a kilonova with a magnetar central engine will be significantly brighter than one without an engine, as this parameter space is where more of the spin-down luminosity is thermalized. In contrast, a system with minimal gravitational-wave emission will produce a kilonova that may be difficult to distinguish from ordinary kilonovae unless early epoch observations are available. However, as the bulk of the energy in this parameter space goes into accelerating the ejecta, such a system will produce a brighter kilonova afterglow that will peak in shorter times. To effectively hide the presence of the magnetar from the kilonova and kilonova afterglow, the rotational energy inputted into the ejecta must be ≲10−3to 10−2 E rot. We discuss the different diagnostics available to identify magnetar-driven kilonovae in serendipitous observations and draw parallels to other potential magnetar-driven explosions, such as superluminous supernovae and broad-line supernovae Ic. [ABSTRACT FROM AUTHOR]

The neutrino dynamics in hot and dense magnetized matter, which corresponds with protoneutron star envelopes in the core collapse supernova explosions, is considered. The kinetic equation for a neutrino phase space distribution function is obtained, taking into account inelastic scattering by nuclear particles. The transfer component in a momentum space using transport properties is studied. The energy transfer coefficient is shown to change from positive to negative values when the neutrino energy exceeds four times the matter temperature. In the vicinity of a neutrino sphere, such effects are illustrated to lead to the energy strengthening in the neutrino spectra. As this paper demonstrates, such a property is favorable for the possibility of observing supernova neutrino fluxes using Large Volume Neutrino Telescopes.

We present a sample of 14 hydrogen-rich superluminous supernovae (SLSNe II) from the Zwicky Transient Facility (ZTF) between 2018 and 2020. We include all classified SLSNe with peaks Mg < −20 mag with observed broad but not narrow Balmer emission, corresponding to roughly 20 per cent of all hydrogen-rich SLSNe in ZTF phase I. We examine the light curves and spectra of SLSNe II and attempt to constrain their power source using light-curve models. The brightest events are photometrically and spectroscopically similar to the prototypical SN 2008es, while others are found spectroscopically more reminiscent of non-superluminous SNe II, especially SNe II-L. 56Ni decay as the primary power source is ruled out. Light-curve models generally cannot distinguish between circumstellar interaction (CSI) and a magnetar central engine, but an excess of ultraviolet (UV) emission signifying CSI is seen in most of the SNe with UV data, at a wide range of photometric properties. Simultaneously, the broad H α profiles of the brightest SLSNe II can be explained through electron scattering in a symmetric circumstellar medium (CSM). In other SLSNe II without narrow lines, the CSM may be confined and wholly overrun by the ejecta. CSI, possibly involving mass lost in recent eruptions, is implied to be the dominant power source in most SLSNe II, and the diversity in properties is likely the result of different mass loss histories. Based on their radiated energy, an additional power source may be required for the brightest SLSNe II, however – possibly a central engine combined with CSI. [ABSTRACT FROM AUTHOR]

Long-term periodicity in the rate of flares is observed for two repeating sources of fast radio bursts (FRBs). In this paper, we present a hydrodynamical modelling of a massive binary consisting of a magnetar and an early-type star. We model the interaction of the pulsar wind from the magnetar with an intense stellar wind. It is shown that only during a fraction of the orbital period, radio emission can escape the system. This explains the duty cycle of the two repeating FRB sources with periodic activity. The width of the transparency window depends on the eccentricity, stellar wind properties, and the viewing angle. To describe the properties of the known sources, it is necessary to assume large eccentricities ≳0.5. We apply the maser cyclotron mechanism of the radio emission generation to model spectral properties of the sources. The produced spectrum is not wide: Δν/ν ∼ 0.2 and the typical frequency depends on the radius of the shock where the emission is generated. The shock radius changes along the orbit. This, together with changing parameters of the medium, allows us to explain the frequency drift during the phase of visibility. Frequency dependence of the degree of polarization at few GHz can be a consequence of a small-scale turbulence in the shocked stellar wind. It is much more difficult to explain huge (∼105 rad/m2) and variable value of the rotation measure observed for FRB 121102. We suggest that this can be explained if the supernova explosion, which produced the magnetar happened near a dense interstellar cloud with n ∼ 100 cm−3. [ABSTRACT FROM AUTHOR]

A millisecond magnetar engine has been widely suggested to exist in gamma-ray burst (GRB) phenomena, in view of its substantial influences on the GRB afterglow emission. In this paper, we investigate the effects of the magnetar engine on the supernova (SN) emission, which is associated with long GRBs and, specifically, confront the model with the observational data of SN 2006aj/GRB 060218. SN 2006aj is featured by its remarkable double-peaked ultraviolet-optical (UV-opt) light curves. By fitting these light curves, we demonstrate that the first peak can be well accounted for by the breakout emission of the shock driven by the magnetar wind, while the primary supernova emission is also partly powered by the energy injection from the magnetar. The magnetic field strength of the magnetar is constrained to be ∼1015 G, which is in good agreement with the common results inferred from the afterglow emission of long GRBs. In more detail, it is further suggested that the UV excess in the late emission of the supernova could also be due to the leakage of the nonthermal emission of the pulsar wind nebula, if some ad hoc conditions can be satisfied. The consistency between the model and the SN 2006aj observation indicates that the magnetar engine is likely to be ubiquitous in the GRB phenomena and even further intensify their connection with the phenomena of superluminous supernovae. [ABSTRACT FROM AUTHOR]

We present optical photometry and spectroscopy of the superluminous SN 2002gh from maximum light to +204 d, obtained as part of the Carnegie Type II Supernova (CATS) project. SN 2002gh is among the most luminous discovered supernovae ever, yet it remained unnoticed for nearly two decades. Using Dark Energy Camera archival images we identify the potential supernova (SN) host galaxy as a faint dwarf galaxy, presumably having low metallicity, and in an apparent merging process with other nearby dwarf galaxies. We show that SN 2002gh is among the brightest hydrogen-poor SLSNe with MV  = −22.40 ± 0.02, with an estimated peak bolometric luminosity of 2.6 ± 0.1 × 1044 erg s−1. We discount the decay of radioactive nickel as the main SN power mechanism, and assuming that the SN is powered by the spin-down of a magnetar we obtain two alternative solutions. The first case, is characterized by significant magnetar power leakage, and M ej between 0.6 and 3.2 M⊙, P spin = 3.2 ms, and B  = 5 × 1013 G. The second case does not require power leakage, resulting in a huge ejecta mass of about 30 M⊙, a fast spin period of P spin ∼ 1 ms, and B ∼ 1.6 × 1014 G. We estimate a zero-age main-sequence mass between 14 and 25 M⊙ for the first case and of about 135 M⊙ for the second case. The latter case would place the SN progenitor among the most massive stars observed to explode as an SN. [ABSTRACT FROM AUTHOR]

Most neutron stars (NSs) and black holes (BHs) are believed to be the final remnants in the evolution of massive stars. In this study, we propose a new formation channel for the formation of BHs and peculiar NSs [specifically, magnetars and Thorne–Żytkow objects (T |$\dot{\rm Z}$| Os)], which we refer to as the core-merger-induced collapse (CMIC) model. This model involves the merger during a common-envelope phase of an oxygen/neon/magnesium composition white dwarf and the core of a hydrogen-rich or helium-rich non-degenerate star, leading to the creation of peculiar new types of objects. The results of binary population synthesis simulations show that the CMIC channel could make important contributions to the populations of (millisecond) pulsars, T |$\dot{\rm Z}$| Os, magnetars, and BHs. The possibility of superluminous supernovae powered by T |$\dot{\rm Z}$| Os, magnetars, and BHs formed through the CMIC model is also being investigated. Magnetars with immediate matter surroundings formed after the CMIC might be good sources for fast radio bursts. [ABSTRACT FROM AUTHOR]

Luminosity evolution of some stripped-envelope supernovae such as Type I superluminous supernovae is difficult to explain by the canonical 56Ni nuclear decay heating. A popular alternative heating source is rapid spin-down of strongly magnetized rapidly rotating neutron stars (magnetars). Recent observations have indicated that Type I superluminous supernovae often have bumpy light curves with multiple luminosity peaks. The cause of bumpy light curves is unknown. In this study, we investigate the possibility that the light-curve bumps are caused by variations of the thermal energy injection from magnetar spin-down. We find that a temporal increase in the thermal energy injection can lead to multiple luminosity peaks. The multiple luminosity peaks caused by the variable thermal energy injection is found to be accompanied by significant increase in photospheric temperature, and photospheric radii are not significantly changed. We show that the bumpy light curves of SN 2015bn and SN 2019stc can be reproduced by temporarily increasing magnetar spin-down energy input by a factor of 2–3 for 5–20 d. However, not all the light-curve bumps are accompanied by the clear photospheric temperature increase as predicted by our synthetic models. In particular, the secondary light-curve bump of SN 2019stc is accompanied by a temporal increase in photospheric radii rather than temperature, which is not seen in our synthetic models. We therefore conclude that not all the light-curve bumps observed in luminous supernovae are caused by the variable thermal energy injection from magnetar spin-down and some bumps are likely caused by a different mechanism. [ABSTRACT FROM AUTHOR]

Recent work has revealed that the light curves of hydrogen-poor (Type I) superluminous supernovae (SLSNe), thought to be powered by magnetar central engines, do not always follow the smooth decline predicted by a simple magnetar spin-down model. Here we present the first systematic study of the prevalence and properties of “bumps” in the post-peak light curves of 34 SLSNe. We find that the majority (44%â€"76%) of events cannot be explained by a smooth magnetar model alone. We do not find any difference in supernova properties between events with and without bumps. By fitting a simple Gaussian model to the light-curve residuals, we characterize each bump with an amplitude, temperature, phase, and duration. We find that most bumps correspond with an increase in the photospheric temperature of the ejecta, although we do not see drastic changes in spectroscopic features during the bump. We also find a moderate correlation (ρ ≈ 0.5; p ≈ 0.01) between the phase of the bumps and the rise time, implying that such bumps tend to happen at a certain “evolutionary phase,” (3.7 ± 1.4) t rise. Most bumps are consistent with having diffused from a central source of variable luminosity, although sources further out in the ejecta are not excluded. With this evidence, we explore whether the cause of these bumps is intrinsic to the supernova (e.g., a variable central engine) or extrinsic (e.g., circumstellar interaction). Both cases are plausible, requiring low-level variability in the magnetar input luminosity, small decreases in the ejecta opacity, or a thin circumstellar shell or disk. [ABSTRACT FROM AUTHOR]

The neutrino dynamics in hot and dense magnetized matter, which corresponds with protoneutron star envelopes in the core collapse supernova explosions, is considered. The kinetic equation for a neutrino phase space distribution function is obtained, taking into account inelastic scattering by nuclear particles. The transfer component in a momentum space using transport properties is studied. The energy transfer coefficient is shown to change from positive to negative values when the neutrino energy exceeds four times the matter temperature. In the vicinity of a neutrino sphere, such effects are illustrated to lead to the energy strengthening in the neutrino spectra. As this paper demonstrates, such a property is favorable for the possibility of observing supernova neutrino fluxes using Large Volume Neutrino Telescopes. [ABSTRACT FROM AUTHOR]

SN 2018hti was a very nearby (z = 0.0614) superluminous supernova with an exceedingly bright absolute magnitude of −21.7 mag in r band at maximum. The densely sampled pre-maximum light curves of SN 2018hti show a slow luminosity evolution and constrain the rise time to ∼50 rest-frame d. We fitted synthetic light curves to the photometry to infer the physical parameters of the explosion of SN 2018hti for both the magnetar and the CSM-interaction scenarios. We conclude that one of two mechanisms could be powering the luminosity of SN 2018hti; interaction with ∼10 M⊙ of circumstellar material or a magnetar with a magnetic field of B p∼ 1.3 × 1013 G, and initial period of P spin∼ 1.8 ms. From the nebular spectrum modelling we infer that SN 2018hti likely results from the explosion of a |${\sim}40\, \mathrm{M}_\odot$| progenitor star. [ABSTRACT FROM AUTHOR]

Optical imaging polarimetry was conducted on the hydrogen poor superluminous supernova 2020znr during three phases after maximum light (≈ +34 d, +288 d, and  +289 d). After instrumental and interstellar polarization correction, all measurements are consistent with null-polarization detection. Modelling the light curve with a magnetar spin-down model shows that SN2020znr has similar magnetar and ejecta parameters to other SLSNe. A comparison of the best-fitting values discussed in the literature on SN 2017egm and SN 2015bn, two hydrogen poor SLSNe showing an increase of polarization after maximum light, suggests that SN 2020znr has higher mass ejecta that may prevent access to the geometry of the inner ejecta with optical polarimetry. The combined information provided by spectroscopy and light-curve analysis of type I SLSNe may be an interesting avenue to categorize the polarization properties of this class of transients. This approach would require to expand the sample of SLSNe polarimetry data currently available with early and late time epochs new measurements. [ABSTRACT FROM AUTHOR]

In this paper, we fit the spectral energy distributions of iPTF 16asu, which has so far been classified as a luminous rapidly evolving broad-lined Ic supernova (SN Ic-BL), and reconstruct its postpeak bolometric light curve. We find that the luminosity of the postpeak bolometric light curve of iPTF 16asu is about 3 times that of the pseudobolometric light curve derived in the literature, and the extrapolated peak luminosity exceeds ∼1044 erg sâˆ'1, which is higher than the threshold of superluminous supernovae (SLSNe). We then use the 56Ni model and the magnetar plus 56Ni model to fit the multiband light curves of iPTF 16asu, and construct the theoretical bolometric light curve using the best-fitting theoretical multiband light curves. We find that the magnetar plus 56Ni model can account for the photometry of iPTF 16asu, and the peak luminosity of its theoretical bolometric light curve is ∼1.06 Ă— 1044 erg sâˆ'1. We suggest that iPTF 16asu and similar SNe (e.g., SN 2018gep) constitute the class of rapidly evolving SLSNe Ic-BL. [ABSTRACT FROM AUTHOR]

We perform a comprehensive study of the physical properties of SN 2018gk, which is a luminous Type IIb supernova (SN). We find that the early-time photospheric velocity varies from a larger value to a smaller one before the photosphere reaches a temperature floor. We generalize the photosphere modulus and fit the multiband light curves (LCs) of SN 2018gk. We find that the ^56 Ni mass model requires ∼0.90 M _⊙ of ^56 Ni, which is larger than the derived ejecta mass (∼0.10 M _⊙ ). Alternatively, we use the magnetar plus ^56 Ni and the fallback plus ^56 Ni models to fit the LCs of SN 2018gk, finding that the two models can fit the LCs. We favor the magnetar plus ^56 Ni model since the parameters are rather reasonable ( M _ej = 1.65 M _⊙ , M _Ni = 0.05 M _⊙ , which is smaller than the upper limit of the ^56 Ni mass that can by synthesized by neutrino-powered core-collapse SNe, B = 6.52 × 10 ^14 G, which is comparable to magnetic fields in luminous and superluminous SNe studied in the literature, and P _0 = 10.42 ms, which is comparable to initial periods for luminous SNe), while the validity of the fallback plus ^56 Ni model depends on the accretion efficiency ( η ). Therefore, we suggest that SN 2018gk might be an SN IIb powered mainly by a central engine. Finally, we confirm the near-IR excesses of the spectral energy distributions of SN 2018gk at some epochs and constrain the physical properties of the putative dust using the blackbody plus dust emission model.

Recent observations indicate that hydrogen-poor superluminous supernovae (SLSNe-I) often display bumpy declining light curves. However, the cause of these undulations remains unclear. In this paper, we have improved the magnetar model, which includes flare activities. We present a systematic analysis of a well-observed SLSN-I sample with bumpy light curves in their late phase. These SLSNe-I were identified from multiple transient surveys, such as the Pan-STARRS1 Medium Deep Survey and the Zwicky Transient Facility. Our study provides a set of magnetar-powered model light curve fits for five SLSNe-I, which accurately reproduce the observed light curves using reasonable physical parameters. By extracting the essential characteristics of both the explosions and central engines, these fits provide valuable insights into investigating their potential association with gamma-ray burst engines. We found that the SLSN flares tend to be the dim and long extension of GRB flares in the peak luminosity versus peak time plane. Conducting large-scale, high-cadence surveys in the near future could enhance our comprehension of both SLSN undulating light curve properties and their potential relationship with GRBs by modeling their light curve characteristics.

Fast X-ray Transients (FXTs) are X-ray flares with durations ranging from a few hundred seconds to a few hours. Possible origins include the tidal disruption of a white dwarf by an intermediate-mass black hole, a supernova shock breakout, or a binary neutron star merger. We present the X-ray light curve and spectrum as well as deep optical imaging of the FXT XRT 210423, which has been suggested to be powered by a magnetar produced in a binary neutron star merger. Our Very Large Telescope and Gran Telescopio Canarias (GTC) observations began on 2021 May 6, thirteen days after the onset of the flare. No transient optical counterpart is found in the 1.″ (3 σ ) X-ray uncertainty region of the source to a depth g _s = 27.0 AB mag. (We use the word “counterpart” for any transient light in a wave band other than the original X-ray detection wave band, whereas the word “host” refers to the host galaxy.) A candidate host lies within the 1.″ X-ray uncertainty region with a magnitude of 25.9 ± 0.1 in the GTC/HiPERCAM g _s filter. Due to its faintness, it was not detected in other bands, precluding a photometric redshift determination. We detect two additional candidate host galaxies: one with z _spec = 1.5082 ± 0.0001 and an offset of 4.″2 ± 1.″ (37 ± 9 kpc) from the FXT, and another one with ${z}_{\mathrm{phot}}={1.04}_{-0.14}^{+0.22}$ and an offset of 3.″6 ± 1.″ (30 ± 8 kpc). Based on the properties of all the prospective hosts, we favor a binary neutron star merger, as previously suggested in the literature, as the explanation for XRT 210423.

Superluminous supernovae (SLSNe) are a rare class of stellar explosions with luminosities ∼ 10–100 times greater than ordinary core-collapse supernovae. One popular model to explain the enhanced optical output of hydrogen-poor (Type I) SLSNe invokes energy injection from a rapidly spinning magnetar. A prediction in this case is that high-energy gamma-rays, generated in the wind nebula of the magnetar, could escape through the expanding supernova ejecta at late times (months or more after optical peak). This paper presents a search for gamma-ray emission in the broad energy band from 100 MeV to 30 TeV from two Type I SLSNe, SN2015bn, and SN2017egm, using observations from Fermi-LAT and VERITAS. Although no gamma-ray emission was detected from either source, the derived upper limits approach the putative magnetar’s spin-down luminosity. Prospects are explored for detecting very-high-energy (VHE; 100 GeV–100 TeV) emission from SLSNe-I with existing and planned facilities such as VERITAS and CTA.

The Nuclear Spectroscopic Telescope Array (NuSTAR) experiment detected a hard X-ray emission (10–70 keV) with a period of 8.68917 s and a pulse-phase modulation at 55 ks, or half this value, from the anomalous X-ray pulsar (AXP) 4U 0142+61. It is shown here that this evidence is naturally explained by the precession of a Keplerian supernova fallback disc surrounding this AXP. It is also found that the precession of discs formed around young neutron stars at distances larger than those considered in the past, may constitute almost neglected sources of gravitational waves with frequencies belonging to the sensitivity bands of the future space interferometers: Laser Interferometer Space Antenna (LISA), Advanced Laser Interferometer Antenna (ALIA), DECi-hertz Interferometer Gravitational wave Observatory (DECIGO), and Big Bang Observer (BBO). In this work, the gravitational wave emission from precessing fallback discs possibly formed around young pulsars such as Crab in a region extending beyond 8 × 107 m from the pulsar surface is estimated. It is also evaluated the role that infrared radiation emission from circumpulsar discs may play in contributing to inverse Compton scattering of TeV energy positrons and electrons. Extensive observational campaigns of disc formation around young and middle-aged pulsars may also contribute to solve the long-standing problem of a pulsar origin for the excess of positrons in cosmic rays observed near the Earth above 7 GeV. In the near future the James Webb Space Telescope , with unprecedented near- and mid-infrared observation capabilities, may provide direct evidence of a large sample of supernova fallback discs. [ABSTRACT FROM AUTHOR]

We perform a systematic study of the 56Ni mass (MNi) of 27 stripped-envelope supernovae (SESNe) by modeling their light-curve tails, highlighting that use of "Arnett's rule" overestimates MNi for SESNe by a factor of ∼2. Recently, Khatami & Kasen presented a new model relating the peak time (tp) and luminosity (Lp) of a radioactively powered supernova to its MNi that addresses several limitations of Arnett-like models, but depends on a dimensionless parameter, β. Using observed tp, Lp, and tail-measured MNi values for 27 SESNe, we observationally calibrate β for the first time. Despite scatter, we demonstrate that the model of Khatami & Kasen with empirically calibrated β values provides significantly improved measurements of MNi when only photospheric data are available. However, these observationally constrained β values are systematically lower than those inferred from numerical simulations, primarily because the observed sample has significantly higher (0.2–0.4 dex) Lp for a given MNi. While effects due to composition, mixing, and asymmetry can increase Lp none can explain the systematically low β values. However, the discrepancy can be alleviated if ∼7%–50% of Lp for the observed sample comes from sources other than radioactive decay. Either shock cooling or magnetar spin-down could provide the requisite luminosity. Finally, we find that even with our improved measurements, the MNi values of SESNe are still a factor of ∼3 larger than those of hydrogen-rich Type II SNe, indicating that these supernovae are inherently different in terms of the initial mass distributions of their progenitors or their explosion mechanisms. [ABSTRACT FROM AUTHOR]

Superluminous supernovae (SLSNe) are massive star explosions that are too luminous to be powered by traditional energy sources, such as the radioactive decay of 56Ni. Instead, they may be powered by a central engine, such as a millisecond pulsar or magnetar, whose relativistic wind inflates a nebula of high-energy particles and radiation behind the expanding supernova ejecta. We present three-dimensional Monte Carlo radiative transfer calculations of SLSNe that follow the production of high-energy emission in the nebula and its subsequent thermalization into optical radiation within the surrounding ejecta and, conversely, determine the gamma-ray emission that escapes the ejecta without thermalizing. We identify a novel mechanism by which γγ pair creation in the upstream pulsar wind regulates the mean energy of particles entering the nebula over the first several years after the explosion, rendering our results on this timescale insensitive to the (uncertain) intrinsic wind pair multiplicity. To explain the observed late-time steepening of SLSN optical light curves as being the result of gamma-ray leakage, we find that the nebular magnetization must be very low,. For higher , the more efficiently thermalized lower-energy synchrotron emission would overproduce the late-time (≳1 yr) optical radiation, inconsistent with observations. For magnetars to remain as viable contenders for powering SLSNe, we conclude that either magnetic dissipation in the wind/nebula is extremely efficient or the spin-down luminosity decays significantly faster than the canonical dipole rate ∝t−2 in a way that coincidentally mimics gamma-ray escape. [ABSTRACT FROM AUTHOR]

We numerically investigate the dynamics of a supernova fallback accretion confronting with a relativistic wind from a newborn neutron star (NS). The time evolution of the accretion shock in the radial direction is basically characterized by the encounter radius of the flow renc and a dimensionless parameter , where L is the NS wind luminosity and is the fallback mass accretion rate. We find that the critical condition for the fallback matter to reach near the NS surface can be simply described as or independent of the wind Lorentz factor, where M* is the NS mass. With combining the condition for the fallback matter to bury the surface magnetic field under the NS crust, we discuss the possibility that the trifurcation of NSs into rotation-powered pulsars, central compact objects, and magnetars can be induced by supernova fallback. [ABSTRACT FROM AUTHOR]

We present the photometric and spectroscopic evolution of supernova (SN) 2019cad during the first ∼100 d from explosion. Based on the light-curve morphology, we find that SN 2019cad resembles the double-peaked Type Ib/c SN 2005bf and the Type Ic PTF11mnb. Unlike those two objects, SN 2019cad also shows the initial peak in the redder bands. Inspection of the g- band light curve indicates the initial peak is reached in ∼8 d, while the r -band peak occurred ∼15 d post-explosion. A second and more prominent peak is reached in all bands at ∼45 d past explosion, followed by a fast decline from ∼60 d. During the first 30 d, the spectra of SN 2019cad show the typical features of a Type Ic SN, however, after 40 d, a blue continuum with prominent lines of Si  ii λ6355 and C  ii λ6580 is observed again. Comparing the bolometric light curve to hydrodynamical models, we find that SN 2019cad is consistent with a pre-SN mass of 11 M⊙, and an explosion energy of 3.5 × 1051 erg. The light-curve morphology can be reproduced either by a double-peaked 56Ni distribution with an external component of 0.041 M⊙, and an internal component of 0.3 M⊙ or a double-peaked 56Ni distribution plus magnetar model (P ∼ 11 ms and B ∼ 26 × 1014 G). If SN 2019cad were to suffer from significant host reddening (which cannot be ruled out), the 56Ni model would require extreme values, while the magnetar model would still be feasible. [ABSTRACT FROM AUTHOR]

We address the issue of the post-maximum bump observed in the light curve of some superluminous supernovae. We rule out the popular mechanism of a circumstellar interaction suggested for the bump explanation. Instead, we propose that the post-maximum bump is caused by the magnetar dipole field enhancement several months after the explosion. The modelling of SN 2019stc light curve based on the thin shell approximation implies that at the age of ∼90 d, the initial dipole magnetic field should be amplified by a factor of 2.8 to account for the post-maximum bump. The specific mechanism for the field amplification of the newborn magnetar on the time-scale of several months has yet to be identified. [ABSTRACT FROM AUTHOR]

The discovery of early bumps in some type-I superluminous supernovae (SLSNe-I) before the main peaks offers an important clue to their energy source mechanisms. In this paper, we updated an analytic magnetar-powered model for fitting the multiband light curves of double-peaked SLSNe-I. The early bump is powered by magnetar-driven shock-breakout thermal emission, and the main peak is powered by a radiative diffusion through the supernova (SN) ejecta as in the standard magnetar-powered model. Generally, the diffusive luminosity is greater than the shock-breakout luminosity at the early time, which usually makes the shock-breakout bumps unclear to observe. To obtain a clear double-peaked light curve, inefficient magnetar heating at early times is required. This model is applied to three well-observed double-peaked SLSNe-I (i.e., SN2006oz, LSQ14bdq, and DES14Xtaz). We find that a relatively massive SN ejecta with Mej ≃ 10.2–18.1M⊙ and relatively large kinetic energy of SN ejecta erg are required, and the thermalization efficiency of the magnetar heating is suppressed before tdelay, which is in the range of ≃15–43 days. The model can reproduce the observed light curves well, with a reasonable and similar set of physical parameters for both the early bump and the main peak, strengthening support for the magnetar-powered model. In the future, modeling of the double-peaked SLSNe-I will become more feasible as more events are discovered before the early bump. [ABSTRACT FROM AUTHOR]