Search

Tidal disruption events (TDEs), where stars are captured or tidally disrupted by supermassive black holes, are potential sources of high-energy neutrinos. We report the discovery of a potential neutrino flare that is spatially and temporally associated with X-ray emission from TDE ATLAS17jrp. The best-fit spectrum of the neutrino flare follows a power-law with an index of $\rm{\gamma=2.7\pm0.4}$ and a flux normalization of $\rm{\Phi_0 =1.7^{+6.3}_{-1.5}\times 10^{-18}\;GeV^{-1} cm^{-2} s^{-1}}$ at 100 TeV based on an analysis of 10-year track data from IceCube, and the flare duration is 61 days. We calculate that the probability of this association occurring by chance is $0.17\%$. Therefore, ATLAS17jrp is the second TDE (not including candidates) associated with high-energy neutrinos, following TDE AT2019dsg associated with an IceCube neutrino alert. This association can be attributed to the interaction of X-ray photons produced by the hot corona with high-energy particles accelerated by disk winds or outflows, resulting in the production of neutrinos., Comment: 8 pages, 2 figures, submitted to ApJL

Magnetars are young, highly magnetized neutron stars that are associated with magnetar short bursts (MSBs), magnetar giant flares (MGFs), and at least some fast radio bursts (FRBs). In this work, we consider a magnetar and a main sequence star in a binary system and analyze the properties of the electromagnetic signals generated by the interaction between the magnetar bursts and the companion star. During the preburst period, persistent radiation could be generated by the interaction between the $e^+e^-$-pair wind from the magnetar and the companion or its stellar wind. We find that for a newborn magnetar, the persistent preburst radiation from the strong magnetar wind can be dominant, and it is mainly at the optical and ultraviolet (UV) bands. For relatively old magnetars, the re-emission from a burst interacting with the companion is larger than the persistent preburst radiation and the luminosity of the companion itself. The transient re-emission produced by the heating process has a duration of $0.1 - 10^5 {\rm~s}$ at the optical, UV, and X-ray bands. Additionally, we find that if these phenomena occur in nearby galaxies within a few hundred kiloparsecs, they could be detected by current or future optical telescopes., Comment: 15 pages, 8 figures; Published in A&A

Gamma-ray bursts (GRBs) have long been considered potential sources of ultra-high-energy cosmic rays (UHECRs; with energy $\gtrsim 10^{18} {\rm~eV}$). In this work, we propose a novel model generating MeV emission lines in GRB, which can constrain the properties of heavy nuclei that potentially exist in GRB jets. Specifically, we find that relativistic hydrogen-like high-atomic-number ions originating from the $\beta$ decay of unstable nuclei and/or the recombination entrained in the GRB jet can generate narrow MeV emission lines through the de-excitation of excited-electrons. This model can successfully explain the MeV emission line observed in the most luminous GRB ever recorded, GRB~221009A, with suitable parameters including a Lorentz factor $\gamma \sim 820-1700$ and a total mass of heavy nuclei $M_{\rm tot} \sim 10^{23} - 10^{26}$~g. Especially, the emission line broadening can be reasonably attributed to both the expansion of the jet shell and the thermal motion of nuclei, naturally resulting in a narrow width ($\sigma_{\rm line} / E_{\rm line} \lesssim 0.2$) consistent with the observation. Furthermore, we predict that different GRBs can exhibit lines in different bands with various evolving behaviors, which might be confirmed with further observations. Finally, our model provides indirect evidence that GRBs may be one of the sources of UHECRs., Comment: 13 pages, 4 figures; Published in ApJL, https://doi.org/10.3847/2041-8213/ad4ce1

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 $\sim 10^{16} \mbox{ 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., Comment: 13pages, 9 figures, Accepted for publication in ApJ

In 2019, Neutron star Interior Composition ExploreR (NICER) mission released its findings on the mass and radius of the isolated neutron star (INS) PSR J0030+0451, revealing a mass of approximately 1.4 solar masses ($M_{\odot}$) and a radius near 13 kilometers. However, the recent re-analysis by the NICER collaboration \citep{vinciguerra2024updated} suggests that the available data primarily yields a precise inference of the compactness for this source while the resulting mass and radius are strongly model-dependent and diverse (the 68.3\% credible regions just overlap slightly for the ST+PDT and PDT-U models). By integrating this compactness data with the equation of state (EoS) refined by our latest investigations, we have deduced the mass and radius for PSR J0030+0451, delivering estimates of $M=1.48^{+0.09}_{-0.10}~M_\odot$ and $R=12.38_{-0.70}^{+0.51}~{\rm km}$ for the compactness found in ST+PDT model, alongside $M=1.47^{+0.14}_{-0.20}~M_\odot$ and $R=12.37_{-0.69}^{+0.50}~{\rm km}$ for the compactness in PDT-U model. These two groups of results are well consistent with each other and the direct X-ray data inference within the ST+PDT model seems to be favored. Additionally, we have calculated the tidal deformability, moment of inertia, and gravitational binding energy for this NS. Furthermore, employing these refined EoS models, we have updated mass-radius estimates for three INSs with established gravitational redshifts., Comment: 10 pages,6 figures, Accepted for publication in ApJ

The measurements of the bulk properties of most isolated neutron stars (INSs) are challenging tasks. Tang et al. (2020) have developed a new method, based on the equation of state (EoS) of neutron star (NS) material constrained by the observational data, to infer the gravitational masses of a few INSs whose gravitational redshifts are available. However, in that work, the authors only considered the constraints on the EoS from nuclear experiments/theories and the gravitational wave data of GW170817; the possible phase transition has not been taken into account. In this work, we adopt three EoS models (including the one incorporates a first-order strong phase transition) that are constrained by the latest multimessenger NS data, including in particular the recent mass\textendash radius measurements of two NSs by Neutron Star Interior Composition Explorer, to update the estimation of the gravitational masses of RBS 1223, RX J0720.4-3125, and RX J1856.5-3754. In comparison to our previous approach, the new constraints are tighter, and the gravitational masses are larger by about $0.1M_\odot$. All the inferred gravitational masses are within the range of the NS masses measured in other ways. We have also calculated the radius, tidal-deformability, and moment of inertia of these sources. The inclusion of the first-order strong phase transition has little influence on modifying the results., Comment: 8pages,4figures

A stellar-mass black hole (BH) or a millisecond magnetar is believed to be born as the central engine of Gamma-ray bursts (GRBs). The presence of plateaus in the X-ray extended emission or afterglow of GRBs is widely accepted as an indicator of magnetar central engine, particularly those with a sharp decay (faster than $t^{-3}$), so-called internal plateau. However, an alternative model, by taking the evolution of the magnetic flux at the BH horizon into account, suggests that an internal plateau can also arise from a Blandford-Znajek (BZ) mechanism powered jet (hereafter referred to as the BZ jet). In this study, we propose that a precessional BZ jet would manifest a Quasi-Periodic Oscillation (QPO) signature on the internal plateau and the subsequent sharp decay. Such lightcurves cannot be readily explained by the activity of a short-lived, supermassive magnetar, thus favoring a Kerr BH as the central engine. The X-ray afterglow of GRB 050904, comprising nine flares, is characterized by a QPO-modulated plateau and sharp decay, which can be well reproduced by a precessional BZ jet model. Therefore, one potential clue for distinguishing between these two engines lies in whether QPO signature throughout the entire plateau and the subsequent sharp decay, as the magnetar scenario suggests a collapse at the end of the plateau., Comment: Accepted for publication in ApJ

Recently, an association of GW190425 and FRB 20190425A had been claimed and a highly magnetized neutron star (NS) remnant was speculated. Given the $\sim 2.5$-h delay of the occurrence of FRB 20190425A, a uniformly rotating supramassive magnetar is favored since the differential rotation would have been promptly terminated by the magnetic braking. The required maximum gravitational mass ($M_{\rm TOV}$) of the nonrotating NS is $\approx 2.77M_\odot$, which is strongly in tension with the relatively low $M_{\rm TOV}\approx 2.25M_\odot$ obtained in current equation of state (EOS) constraints incorporating perturbative quantum chromodynamics (pQCD) information. However, the current mass-radius and mass-tidal deformability measurements of NSs alone do not convincingly exclude the high $M_{\rm TOV}$ possibility. By performing EOS constraints with mock measurements, we find that with a $2\%$ determination for the radius of PSR J0740+6620-like NS it is possible to distinguish between the low and high $M_{\rm TOV}$ scenarios. We further explore the prospect to resolve the issue of the appropriate density to impose the pQCD constraints with future massive NS observations or determinations of $M_{\rm TOV}$ and/or $R_{\rm TOV}$. It turns out that measuring the radius of a PSR J0740+6620-like NS is insufficient to probe the EOSs around 5 nuclear saturation density, where the information from pQCD becomes relevant. The additional precise $M_{\rm TOV}$ measurements anyhow could provide insights into the EOS at such a density. Indeed, supposing the central engine of GRB 170817A is a black hole formed via the collapse of a supramassive NS, the resulting $M_{\rm TOV}\approx 2.2M_\odot$ considerably softens the EOS at the center of the most massive NS, which is in favor of imposing the pQCD constraint at density beyond the one achievable in the NSs., Comment: 12 pages, 5 figures, PRD published

The long gamma-ray burst GRB 191221B has abundant observations in X-ray, optical and radio bands. In the literature, the observed optical light curve of GRB 191221B displays a plateau around 0.1-day, which is rather peculiar in gamma-ray bursts. Here we performed detailed analysis of the observational data from Swift/UVOT, VLT and LCO, obtained the light curve of the multi-band afterglow of GRB 191221B. By examining optical, ultraviolet, X-ray, and radio data for this event, we demonstrate that an on-axis two-component jet model can explain the observations. Our analysis suggests that the narrow component has an initial Lorentz factor of 400 and a jet opening half-angle of $1.4^{\circ}$, while the wide component has an initial Lorentz factor of 25 and a jet opening half-angle of $2.8^{\circ}$. The narrow jet dominates the early decay, whereas the wider jet causes the optical plateau and dominates late decay. According to this model, the reason for the absence of the X-ray plateau is due to the steeper spectral index of the wide component, resulting in a less significant flux contribution from the wide jet in the X-ray bands than in the optical bands. Moreover, we have explained the inconsistency in the decay indices of the UVOT and Rc-band data around 2000 seconds using reverse shock emission., Comment: 10 pages, 7 figures, Accepted by the Monthly Notices of the Royal Astronomical Society

The mass ($M_{\rm TOV}$) and radius ($R_{\rm TOV}$) of the maximum-mass nonrotating neutron star (NS) play a crucial role in constraining the elusive equation of state of cold dense matter and in predicting the fate of remnants from binary neutron star (BNS) mergers. In this study, we introduce a novel method to deduce these parameters by examining the mergers of second-generation (2G) black holes (BHs) with NSs. These 2G BHs are assumed to originate from supramassive neutron stars (SMNSs) formed in BNS mergers. Since the properties of the remnant BHs arising from the collapse of SMNSs follow a universal relation governed by $M_{\rm TOV}$ and $R_{\rm TOV}$, we anticipate that by analyzing a series ($\sim 100$ detections) of mass and spin measurements of the 2G BHs using the third-generation ground-based gravitational-wave detectors, $M_{\rm TOV}$ and $R_{\rm TOV}$ can be determined with a precision of $\sim 0.01M_\odot$ and $\sim 0.6$ km, respectively., Comment: 11 pages, 5 figures, ApJ published

For bright transients such as Gamma-Ray Bursts (GRBs), the Ultra-Violet/Optical Telescope (UVOT) operates under event mode at early phases, which records incident positions and arrival time for each photon. The event file is able to be screened into many exposures to study the early light curve of GRBs with a high time resolution, including in particular the rapid brightening of the UV/Optical emission. Such a goal, however, is hampered for some extremely bright GRBs by the saturation in UVOT event images. For moderately saturated UVOT sources, in this work we develop the method proposed in Jin et al. (2023) to recover their photometries. The basic idea is to assume a stable point spread function (PSF) of UVOT images, for which the counts in the core region (i.e., an aperture of a radius of 5 arcsec) and the wing region (i.e., an annulus ranging from 15 arcsec to 25 arcsec) should be a constant and the intrinsic flux can be reliably inferred with data in the ring. We demonstrate that in a given band, a tight correlation does hold among the background-removed count rates in the core and the wing. With the new method, the bright limit of measuring range for UVOT V and B bands increases ~ 1.7 mag, while only ~ 0.7 mag for U band due to the lack of bright calibration sources. Systematic uncertainties are ~ 0.2 mag for V, B and U bands., Comment: 27 pages, 19 figures, 3 tables, accepted to ApJS

Studies have proposed that there is evidence for cosmological coupling of black holes (BHs) with an index of $k\approx 3$; hence, BHs serve as the astrophysical source of dark energy. However, the data sample is limited for the redshifts of $\leq 2.5$. In recent years, the James Webb Space Telescope (JWST) has detected many high-redshift active galactic nuclei (AGNs) and quasars. Among the JWST NIRSpec-/NIRCam-resolved AGNs, three are determined to be in early-type host galaxies with a redshift of $z\sim 4.5--7$. However, their $M_{\star}$ and $M_{\rm BH}$ are in tension with the predicted cosmological coupling of black holes with $k = 3$ at a confidence level of $\sim 2\sigma$, which challenges the hypothesis that BHs serve as the origin of dark energy. Future work on high-redshift AGNs using the JWST will further assess such a hypothesis by identifying more early-type host galaxies in the higher mass range., Comment: 9 pages, 3 figures, 2 tables; Comments are welcome!

Hyperluminous optical-ultraviolet fares have been detected in gamma-ray bursts and the luminosity record was held by naked-eye event GRB 080319B. Such fares are widely attributed to internal shock or external reverse shock radiation. Here, with a new method developed to derive reliable photometry from saturated sources of Swift/UVOT, we carry out time-resolved analysis of the initial white-band 150s exposure of GRB 220101A, a burst at a redshift of 4.618, and report a rapidly evolving optical-ultraviolet fare with a high absolute AB magnitude of -39.4$\pm$0.2. In contrast to GRB 080319B, the temporal behaviour of this new fare does not trace the gamma-ray activity. Instead of either internal shocks or reverse shock, this extremely energetic optical-ultraviolet fare is most likely to originate from the refreshed shocks induced by the late-ejected extremely energetic material catching up with the earlier-launched decelerating outflow. This finding reveals the diverse origins of the extremely energetic optical-ultraviolet fares and demonstrates the necessity of high-time-resolution observations at early times., Comment: Published in Nature Astronomy, the definitive version is available at https://www.nature.com/articles/s41550-023-02005-w

The temporal behavior of the very dim optical afterglow of GRB 080503 is at odds with the regular forward shock afterglow model and a sole kilonova component responsible for optical emission has been speculated in some literature. Here we analyze the optical afterglow data available in archive and construct time-resolved spectra. The significant detection by Keck-I in {\it G/R} bands at $t\sim 3$ day, which has not been reported before, as well as the simultaneous Gemini-North {\it r} band measurement, are in favor of a power-law spectrum that is well consistent with the optical to X-ray spectrum measured at $t\sim 4.5$ day. However, for $t\leq 2$ day, the spectra are thermal-like and a straightforward interpretation is a kilonova emission from a neutron star merger, making it, possibly, the first detection of a very early kilonova signal at $t\sim 0.05$ day. A non-thermal nature of optical emission at late times ($t\sim 2$ day), anyhow, can not be ruled out because of the large uncertainty of the {\it g}-band data. We also propose to classify the neutron star merger induced optical transients, according to the temporal behaviors of the kilonova and the non-thermal afterglow emission, into four types. GRB 080503 would then represent the first observation of a sub-group of neutron star merger driven optical transients (i.e., Type IV) consisting of an early blue kilonova and an adjacent non-thermal afterglow radiation., Comment: Accepted to ApJ. 15 pages, 4 figures, 5 tables

The LIGO, Virgo, and KAGRA (LVK) collaboration has announced 90 coalescing binary black holes (BBHs) with $p_{\rm astro} > 50\%$ to date, however, the origin of their formation channels is still an open scientific question. Given various properties of BBHs (BH component masses and individual spins) inferred using the default priors by the LVK, independent groups have been trying to explain the formation of the BBHs with different formation channels. Of all formation scenarios, the chemically homogeneous evolution (CHE) channel has stood out with distinguishing features, namely, nearly-equal component masses and preferentially high individual spins aligned with the orbital angular momentum. We perform Bayesian inference on the BBH events officially reported in GWTC-3 with astrophysically-predicted priors representing different formation channels of the isolated binary evolution (CEE: common-envelope evolution channel; CHE; SMT: stable mass transfer). Given assumed models, we report strong evidence for GW190517\_055101 being most likely to have formed through the CHE channel. Assuming the BBH events in the subsample are all formed through one of the isolated binary evolution channels, we obtain the lower limits on the local merger rate density of these channels at $11.45 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (CEE), $0.18 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (CHE), and $0.63 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (SMT) at $90\%$ credible level., Comment: 13 pages, 4 figures, 1 table

Hawking's area theorem is one of the fundamental laws of black holes (BHs), which has been tested at a confidence level of $\sim 95\%$ with gravitational wave (GW) observations by analysing the inspiral and ringdown portions of GW signals, independently. In this work, we propose to carry out the test in a new way with the hierarchical triple merger (i.e., two successive BH mergers occurred sequentially within the observation window of GW detectors), for which the properties of the progenitor BHs and the remnant BH of the first coalescence can be inferred from the inspiral portions of the two mergers. As revealed in our simulations, the BH area law can be well confirmed for some plausible hierarchical triple merger events detectable in LIGO/Virgo/KAGRA/LIGO-India's O4/O5 runs. Our proposed method provides significant facilitation for testing the area law and complements previous methods., Comment: 8 pages, 8 figures, 2 tables, published in MNRAS

The origins of coalescing binary black holes (BBHs) detected by the advanced LIGO/Virgo are still under debate, and clues may be present in the joint mass-spin distribution of these merger events. Here we construct phenomenological models containing two sub-populations to investigate the BBH population detected in gravitational wave observations. We find that our models can explain the GWTC-3 data rather well, and several constraints to our model are required by the data: first, the maximum mass for the component with a stellar-origin, $m_{\rm max}$, is $39.1^{+2.4}_{-2.7}M_{\odot}$ at 90\% credibility; second, about $15\%$ of the mergers happen in dynamical environments, in which $7-16\%$ of events are hierarchical mergers, and these BHs have an average spin magnitude significantly larger than the first-generation mergers, with ${\rm d}\mu_{\rm a} > 0.4 $ at $99\%$ credibility; third, the dynamical component BHs tend to pair with each other with larger total mass and higher mass ratio. An independent analysis focusing on spins is also carried out, and we find that the spin amplitude of component BHs can be divided into two groups according to a division mass $m_{\rm d} = 46.1^{+5.6}_{-5.1}M_{\odot}$. These constraints can be naturally explained by current formation channels, and our results suggest that some of the observed events were likely from AGN disks., Comment: 21 pages, 8 figures. Published in The Astrophysical Journal Letters

The X-ray plateau emission observed in many Long gamma-ray bursts (LGRBs) has been usually interpreted as the spin-down luminosity of a rapidly spinning, highly magnetized neutron star (millisecond magnetar). If this is true, then the magnetar may emit extended gravitational wave (GW) emission associated with the X-ray plateau due to non-axisymmetric deformation or various stellar oscillations. The advanced LIGO and Virgo detectors have searched for long-duration GW transients for several years, no evidence of GWs from any magnetar has been found until now. In this work, we attempt to search for signature of GW radiation in the electromagnetic observation of 30 LGRBs under the assumption of the magnetar model. We utilize the observations of the LGRB plateau to constrain the properties of the new-born magnetar, including the initial spin period $P_0$, diploe magnetic field strength $B_p$ and the ellipticity $\epsilon$. We find that there are some tight relations between magnetar parameters, e.g., $\epsilon \propto B_p^{1.29}$ and $B_p\propto P_0^{1.14}$. In addition, we derive the GW strain for magnetar sample via their spin-down processes, and find that the GWs from these objects may not be detectable by the aLIGO and ET detectors. For a rapidly spinning magnetar ($P\sim 1\mbox{ ms}$, $ B \sim10^{15}\mbox{ G}$), the detection horizon for advanced LIGO O5 detector is $\sim 180\mbox{ Mpc}$. The detection of such GW signal associated with the X-ray plateau would be a smoking gun that the central engine of GRB is a magnetar., Comment: 15 pages, 7 figures,accepted for publication in ApJ

Coalescing binary black hole (BBH) systems are likely formed via several channels, and it is challenging to understand their formation / evolutionary processes. Some features in the mass function of the primary components ($m_1$), such as the distinct Gaussian-like peak located at $\sim 34M_\odot$, have been previously found. In this work, we investigate the possible dependence of the mass ratio ($q=m_2/m_1$) distribution on the primary mass. We find a Bayesian odds ratio of 18.1 in favor of divergence in the mass ratio distributions between the low- and high-mass ranges over an invariable mass ratio distribution. The BBHs with $m_1\gtrsim29M_{\odot}$ have a stronger preference to be symmetric compared to those with $m_1\lesssim29M_{\odot}$ at a 97.6\% credible level. Additionally, we find mild evidence that the BBHs with $m_1$ located in the Gaussian-like peak have a mass ratio distribution different from that of other BBHs. Our findings may be in favor of some formation channels, such as the chemically homogeneous evolution and the dynamical assembly in globular clusters/nuclear star clusters, which are more likely to provide symmetric BBHs in the high-mass range., Comment: 11 pages, 6 figures, 2 tables, published In ApJL

The DArk Matter Particle Explorer (DAMPE) is well suitable for searching for monochromatic and sharp $\gamma$-ray structures in the GeV$-$TeV range thanks to its unprecedented high energy resolution. In this work, we search for $\gamma$-ray line structures using five years of DAMPE data. To improve the sensitivity, we develop two types of dedicated data sets (including the BgoOnly data which is the first time to be used in the data analysis for the calorimeter-based gamma-ray observatories) and adopt the signal-to-noise ratio optimized regions of interest (ROIs) for different DM density profiles. No line signals or candidates are found between 10 and 300 GeV in the Galaxy. The constraints on the velocity-averaged cross section for $\chi\chi \to \gamma\gamma$ and the decay lifetime for $\chi \to \gamma\nu$, both at 95% confidence level, have been calculated and the systematic uncertainties have been taken into account. Comparing to the previous Fermi-LAT results, though DAMPE has an acceptance smaller by a factor of $\sim 10$, similar constraints on the DM parameters are achieved and below 100 GeV the lower limits on the decay lifetime are even stronger by a factor of a few. Our results demonstrate the potential of high-energy-resolution observations on dark matter detection., Comment: 14 pages, 8 figures. Update the content to keep up with the published version