High Energy Astrophysical Phenomena (astro-ph.HE)

Abstract: In this study, we analyze the post-maximum spectra of a sample of 27 Type I superluminous supernovae (SLSNe-I) in order to search for physical differences between the so-called Type W and Type 15bn sub-types. This paper is a continuation of \citet{ktr21} and \citet{ktr22}. In the former, it was revealed that not all SLSNe-I show the W-shaped absorption feature between 4000 and 5000 \AA\ in the pre-maximum spectra, and two new SLSN-subgroups were disclosed: Type W, where the W-shaped feature is present, and Type 15bn, where it is missing. In the latter, it was shown that the pre-maximum photosphere of Type W SLSNe-I tend to be hotter compared to Type 15bn objects, and they are different regarding their ion composition, their early light curves and their geometry as well. For completeness, post-maximum data are analyzed in this paper. It is concluded that in terms of photospheric temperature and velocity, Type W and Type 15bn SLSNe decrease to a similar value by the post-maximum phases, and their pseudo-nebular spectra are nearly uniform. Pseudo-equivalent width calculations show that the pEW of the wavelength range between 4166 and 5266 \AA\ evolve differently in case of the two sub-types, while the other parts of the spectra seem to evolve similarly. It was found that the host galaxies of the studied objects do not differ significantly in their star formation rate, morphology, stellar mass and absolute brightness. The main difference behind the bimodality of Type W and Type 15bn SLSNe-I therefore is in their pre-maximum evolution.

Abstract: Recent observations indicate that hydrogen-poor superluminous supernovae 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 SLSNe-I sample with bumpy light curves in the late-phase. These SLSNe-I were identified from multiple transient surveys, such as the Pan-STARRS1 Medium Deep Survey (PS1 MDS) and the Zwicky Transient Facility (ZTF). Our study provides a set of magnetar-powered model light curve fits for five SLSNe-I, which accurately reproduce observed light curves using reasonable physical parameters. By extracting essential characteristics of both 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 the 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 undulation properties and their potential relationship with GRBs by modeling their light curve characteristics.

Abstract: The origin of fast radio bursts (FRBs), the brightest cosmic explosion in radio bands, remains unknown. Magnetar-related mechanisms are currently favored. The searches for short-term periodicity that is naturally expected for such fast-spinning compact objects, however, have failed. We introduce here a novel method for a comprehensive analysis of active FRBs' behaviors in the time-energy domain. Using ``Pincus Index'' and ``Maximum Lyapunov Exponent'', we were able to quantify the stochasticity and chaos, respectively, of the bursting events and put FRBs in the context of common transient physical phenomena, such as pulsars, earthquakes, and solar flares. In the bivariate time-energy domain, repeated FRB bursts' behaviors deviate significantly (more random, less chaotic) from pulsars, earthquakes, and solar flares. FRB bursts wander in time-energy space stochastically, akin to Brownian motions. The high degree of stochasticity suggests complex and even multi-origins for FRBs.

Abstract: We conduct a spectral and timing analysis of GX 339-4 and EXO 1846-031 with the aim of studying the evolution of Type-C QPOs with spectral parameters. The high cadence data from Insight-HXMT and NICER allow us to track them. Type-C QPOs appear at the end of low-hard state and/or hard-intermediate state. The results reveal that the QPO frequency is closely related to the inner disk radius and mass accretion rate in the two sources. Such a correlation is nicely consistent with the dynamic frequency model.

Abstract: Swift has now completed 18 years of mission, during which it discovered thousands of gamma-ray bursts (GRBs) as well as new classes of high-energy transient phenomena. Its first breakthrough result was the localization of short duration GRBs, which enabled for redshift measurements and kilonova searches. Swift, in synergy with the Hubble Space Telescope and a wide array of ground-based telescopes, provided the first tantalizing evidence of a kilonova in the aftermath of a short GRB. In 2017, Swift observations of the gravitational wave event GW170817 captured the early UV photons from the kilonova AT2017gfo, opening a new window into the physics of kilonovae. Since then, Swift has continued to expand the sample of known kilonovae, leading to the surprising discovery of a kilonova in a long duration GRB. This article will discuss recent advances in the study of kilonovae driven by the fundamental contribution of Swift.

Abstract: While supermassive binary black holes (SMBBHs) inspiral toward merger they may also accrete significant amounts of matter. To study the dynamics of such a system requires simultaneously describing the evolving spacetime and the dynamics of magnetized plasma. Here we present the first relativistic calculation simulating two equal-mass, non-spinning black holes as they inspiral from an initial separation of $20M$ ($G=c=1$) almost to merger, $\simeq 9M$, while accreting gas from a surrounding disk, where $M$ is the total binary mass. We find that the accretion rate $\dot M$ onto the black holes first decreases during this period and then reaches a plateau, dropping by only a factor of $\sim 3$ despite its rapid inspiral. An estimated bolometric light curve follows the same profile. The minidisks through which the accretion reaches the black holes are very non-standard. Reynolds, not Maxwell, stresses dominate, and they oscillate between two distinct structural states. In one part of the cycle, ``sloshing" streams transfer mass from one minidisk to the other through the L1 point at a rate $\sim 0.1\times$ the accretion rate, carrying kinetic energy at a rate that can be as large as the peak minidisk bolometric luminosity. We also discover that the minidisks have time-varying tilts with respect to the orbital plane similar in magnitude to the circumbinary disk's aspect ratio. The unsigned poloidal flux on the black hole event horizon is roughly constant at a dimensionless level $\phi\sim 2-3$, but doubles just before merger; if the black holes had significant spin, this flux could support jets whose power could approach the radiated luminosity. This simulation is the first to employ our multipatch infrastructure \pwmhd, decreasing computational expense per physical time to $\sim 3\%$ of similar runs using conventional single-grid methods.

Abstract: Gravitational-wave observations have revealed sources whose unusual properties challenge our understanding of compact-binary formation. Inferring the formation processes that are best able to reproduce such events may therefore yield key astrophysical insights. A common approach is to simulate a population of mergers and count the fraction of these synthetic events that lie within a chosen region in the measured parameters of a real event. Though appealing owing to its simplicity, this approach is flawed because it neglects the full posterior information, depends on a heuristically constructed confidence region, and fails for high signal-to-noise detections. We point out that the statistically consistent solution is to compute the posterior odds between two simulated populations, which crucially is a relative measure, and show how to include the effect of observational biases by conditioning on source detectability. Applying the approach to several gravitational-wave events and astrophysical populations, we assess the degree to which we can conclude model preference not just between distinct formation pathways but also between subpopulations within a given pathway.

Abstract: The gamma-ray binary LS I+61 303 has been widely monitored at different wavelengths since its discovery more than sixty years ago. However, the nature of the compact object and the peculiar behavior of the system are still largely debated. Aimed at investigating the rapid X-ray variability of LS I+61 303, we have analysed all the archival RXTE/PCA data of the source, taken between 1996 and 2011. The timing analysis yields a periodicity of $P\sim 26.6\pm 0.3$ days, which is statistically compatible with several periodicities reported in the literature for LS I+61 303. Using this period, we performed a data phase-resolved analysis to produce a set of phase-bin-averaged energy spectra and power density spectra. These power density spectra are dominated by weak red noise below 0.1 Hz, and show no signal above this frequency. The amplitude of the red noise varies mildly with the phase, and shows a maximum that coincides with a dip of the X-ray flux and a softer photon index. Aside from low-frequency noise, this analysis does not provide any statistically significant periodic or quasi-periodic timing feature in the RXTE/PCA data of LS I+61 303.

Abstract: Aims. Fast Radio Bursts are bright radio transients whose origin has not yet explained. The search for a multi-wavelength counterpart of those events can put a tight constrain on the emission mechanism and the progenitor source. Methods. We conducted a multi-wavelength observational campaign on FRB 20180916B between October 2020 and August 2021 during eight activity cycles of the source. Observations were led in the radio band by the SRT both at 336 MHz and 1547 MHz and the uGMRT at 400 MHz. Simultaneous observations have been conducted by the optical telescopes Asiago (Galileo and Copernico), CMO SAI MSU, CAHA 2.2m, RTT-150 and TNG, and X/Gamma-ray detectors on board the AGILE, Insight-HXMT, INTEGRAL and Swift satellites. Results. We present the detection of 14 new bursts detected with the SRT at 336 MHz and seven new bursts with the uGMRT from this source. We provide the deepest prompt upper limits in the optical band fro FRB 20180916B to date. In fact, the TNG/SiFAP2 observation simultaneous to a burst detection by uGMRT gives an upper limit E_optical / E_radio < 1.3 x 10^2. Another burst detected by the SRT at 336 MHz was also co-observed by Insight-HMXT. The non-detection in the X-rays yields an upper limit (1-30 keV band) of E_X-ray / E_radio in the range of (0.9-1.3) x 10^7, depending on which model is considered for the X-ray emission.