Search

The ultraviolet/optical telescope (UVOT) onboard the Neil Gehrels Swift Observatory is capable of imaging with 7 lenticular filters and of taking slitless spectra with 2 grisms. Both image and grism data have been widely used to study gamma-ray bursts, supernovae and other ultraviolet/optical transients, and proved UVOT is a powerful instrument in time-domain astronomy. However, the second order contamination, for blue sources, strongly limits the red end of ultraviolet (UV) grism spectra. This, in turn, reduces the valid wavelength range to only about 33% of the total. However, to explore the broadband spectral energy distribution of GRBs at the early stage, a larger valid wavelength range is required. Hence based on the uvotpy package, we propose a method to remove the second order contamination from UV grism spectra (nominal mode) up to about 4000\AA, i.e., about 70% of the full wavelength range. The 1-sigma systematic uncertainty of this method is about 11.2%. In addition, if a source is red enough, the red end of the valid range could reach about 5000\AA. The source code is available on GitHub., Comment: 14 pages, 8 figures, 3 tables. Accepted by the ApJS

We present a detailed analysis of the long-duration GRB 241030A detected by {\it Swift}. Thanks to the rapid response of XRT and UVOT, the strongest part of the prompt emission of GRB 241030A has been well measured simultaneously from optical to hard X-ray band. The time-resolved WHITE band emission shows strong variability, largely tracing the activity of the prompt gamma-ray emission, may be produced by internal shocks too. The joint analysis of the XRT and BAT data reveals the presence of a thermal component with a temperature of a few keV, which can be interpreted as the photosphere radiation, and the upper limit of the Lorentz factor of this region is found to range between approximately 20 and 80. The time-resolved analysis of the initial U-band exposure data yields a very rapid rise ($ \sim t^{5.3}$) with a bright peak reaching 13.6 AB magnitude around 410 seconds, which is most likely attributed to the onset of the external shock emission. The richness and fineness of early observational data have made this burst a unique sample for studying the various radiation mechanisms of gamma-ray bursts.

Ultrahigh-energy cosmic rays (UHECRs) are the highest energy messenger from space, with energies exceeding 1 EeV. Although UHECRs were discovered over 60 years ago, their origin still remains a mystery. Pinpointing sources of UHECRs is crucial for understanding the extreme astrophysical processes that accelerate particles to such extraordinary energies. We searched for UHECR multiplets via analyzing 17 years of data with energies greater than 40 EeV from the Pierre Auger Observatory. A spatial association is found between a multiplet of $25.7^{+6.2}_{-7.0}$ cosmic rays and the Sombrero galaxy with a local (global) significance of $4.5~\sigma~(3.3~\sigma)$. The Sombrero galaxy hosts a supermassive central black hole with a mass of $\sim1\times 10^9 M_{\odot}$ and exhibits large-scale radio lobes and jets. Our finding provides critical evidence on active supermassive black holes as the source of the highest-energy cosmic rays.

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

The origins of the coalescing binary black holes (BBHs) detected by the advanced LIGO/Virgo are still in debate and clues may present in the mass and effective spin ($\chi_{\rm eff}$) distributions of these merger events. Here we analyze the GWTC-1 and GWTC-2.1 BBH events, and find evidence for a primary mass dependent $\chi_{\rm eff}$ distribution. Below (above) the mass of $\sim 50M_\odot$, a nearly zero (relatively high) effective spin is typical. In the joint analysis of the primary mass and $\chi_{\rm eff}$ distributions, we identify a sub-population of BBHs characterized by both a bump-like primary mass distribution extending to $\sim 100M_\odot$ and a relatively high $\chi_{\rm eff}\sim 0.3$, consistent with the expectation for the $\geq 2$g merger events in the disks of active galactic nuclei (AGNs). Our results suggest that the merger rate in AGN disks is around $1.7~{\rm Gpc}^{-3}{\rm yr}^{-1}$., Comment: 7 pages, 3 figures

GRB 060505 was the first well-known nearby (at redshift 0.089) "hybrid" gamma-ray burst that has a duration longer than 2 seconds but without the association of a supernova down to very stringent limits. The prompt $\gamma-$ray flash lasting $\sim 4$ sec could consist of an intrinsic short burst and its tail emission, but the sizable temporal lag ($\sim 0.35$ sec) as well as the environment properties led to the widely-accepted classification of a long duration gamma-ray burst originated from the collapse of a massive star. Here for the $ first$ time we report the convincing evidence for a thermal-like optical radiation component in the spectral energy distribution of the early afterglow emission. In comparison to AT2017gfo, the thermal radiation is $\sim 2$ times brighter and the temperature is comparable at similar epochs. The optical decline is much quicker than that in X-rays, which is also at odds with the fireball afterglow model but quite natural for the presence of a blue kilonova. Our finding reveals a neutron star merger origin of the hybrid GRB 060505 and strongly supports the theoretical speculation that some binary neutron stars can merge ultra-quickly (within $\sim 1$ Myr) after their formation when the surrounding region is still highly star-forming and the metallicity remains low. Gravitational wave and electromagnetic jointed observations are expected to confirm such scenarios in the near future., Comment: 20 pages, 7 figures, 2 tables

We perform a hierarchical Bayesian inference to investigate the population properties of the coalescing compact binaries involving at least one neutron star (NS). With the current observation data, we can not rule out either of the Double Gaussian, Single Gaussian and Uniform NS mass distribution models, although the mass distribution of the Galactic NSs is slightly preferred by the gravitational wave (GW) observations. The mass distribution of black holes (BHs) in the neutron star-black hole (NSBH) population is found to be similar to that for the Galactic X-ray binaries. Additionally, the ratio of the merger rate densities between NSBHs and BNSs is estimated to be about 3 : 7. The spin properties of the binaries, though constrained relatively poor, play nontrivial role in reconstructing the mass distribution of NSs and BHs. We find that a perfectly aligned spin distribution can be ruled out, while a purely isotropic distribution of spin orientation is still allowed.

In a binary system, the gravitational potential of the primary black hole may play an important role in enhancing the fallback accretion onto the lighter compact object newly formed in the second supernova explosion. As a result, the final masses of the binary compact objects would be correlated, as suggested recently by Safarzadeh \& Wysocki. In this work, we analyze the mass distribution of four gravitational-wave events, which are characterized by both a small mass ratio and a low mass ($\leq 5M_\odot$) of the light component, and find tentative evidence for a mass correlation among the objects. To evaluate the feasibility of testing such a hypothesis with upcoming observations, we carry out simulations with a mock population and perform Bayesian hierarchical inference for the mass distribution. We find that with dozens of low mass ratio events, whether there exists correlation in the component mass distributions or not can be robustly tested and the correlation, if it exists, can be well determined., Comment: 8 pages, 3 figures, ApJ published

Recently, the radius of neutron star (NS) PSR J0740+6620 was measured by Neutron Star Interior Composition Explorer (NICER) and an updated measurement of neutron skin thickness of ${}^{208}$Pb ($R_{\rm skin}^{208}$) was reported by the PREX-II experiment. These new measurements can help us better understand the unknown equation of state (EOS) of dense matter. In this work, we adopt a hybrid parameterization method, which incorporates the nuclear empirical parameterization and some widely used phenomenological parameterizations, to analyze the results of nuclear experiments and astrophysical observations. With the joint Bayesian analysis of GW170817, PSR J0030+0451, and PSR J0740+6620, the parameters that characterize the ultradense matter EOS are constrained. We find that the slope parameter $L$ is approximately constrained to $70_{-18}^{+21}$ MeV, which predicts $R_{\rm skin}^{208}=0.204^{+0.030}_{-0.026}\,{\rm fm}$ by using the universal relation between $R_{\rm skin}^{208}$ and $L$. The bulk properties of canonical $1.4\,M_\odot$ NS (e.g., $R_{1.4}$ and $\Lambda_{1.4}$) as well as the pressure ($P_{2\rho_{\rm sat}}$) at two times the nuclear saturation density are well constrained by the data; i.e., $R_{1.4}$, $\Lambda_{1.4}$, and $P_{2\rho_{\rm sat}}$ are approximately constrained to $12.3\pm0.7$ km, $330_{-100}^{+140}$, and $4.1_{-1.2}^{+1.5}\times10^{34}\,{\rm dyn\,cm^{-2}}$, respectively. Besides, we find that the Bayes evidences of the hybrid star and normal NS assumptions are comparable, which indicates that current observation data are compatible with quarkyonic matter existing in the core of massive star. Finally, in the case of normal NS assumption, we obtain a constraint for the maximum mass of nonrotating NS $M_{\rm TOV}=2.30^{+0.30}_{-0.18}$ $M_\odot$. All of the uncertainties reported above are for 68.3% credible levels., Comment: 12 pages, 7 figures, published in PRD

We develop a new method based on Gaussian process to reconstruct the mass distribution of binary black holes (BBHs). Instead of prespecifying the formalisms of mass distribution, we introduce a more flexible and nonparametric model with which the distribution can be mainly determined by the observed data. We first test our method with simulated data, and find that it can well recover the injected distribution. Then we apply this method to analyze the data of BBHs' observations from LIGO-Virgo Gravitational-Wave Transient Catalog 2. By reconstructing the chirp mass distribution, we find that there is a peak or a platform locating at $20-30\,M_{\odot}$ rather than a single-power-law-like decrease from low mass to high mass. Moreover, one or two peaks in the chirp mass range of $\mathcal{M}<20\,M_{\odot}$ may be favored by the data. Assuming a mass-independent mass ratio distribution of $p(q)\propto q^{1.4}$, we further obtain a distribution of primary mass, and find that there is a feature locating in the range of $(30, 40)\,M_{\odot}$, which can be related to \textsc{Broken Power Law} and \textsc{Power Law + Peak} distributions described in The LIGO Scientific Collaboration et al. (2020). Besides, the merger rate of BBHs is estimated to $\mathcal{R} = 26.29^{+14.21}_{-8.96}~{\rm Gpc^{-3}~yr^{-1}}$ supposing there is no redshift evolution.

We analyze the LIGO/Virgo GWTC-2 catalog to study the primary mass distribution of the merging black holes. We perform hierarchical Bayesian analysis, and examine whether the mass distribution has a sharp cutoff for primary black hole masses below $65 M_\odot$, as predicted in pulsational pair instability supernova model. We construct two empirical mass functions. One is a piece-wise function with two power-law segments jointed by a sudden drop. The other consists of a main truncated power-law component, a Gaussian component, and a third very massive component. Both models can reasonably fit the data and a sharp drop of the mass distribution is found at $\sim 50M_\odot$, suggesting that the majority of the observed black holes can be explained by the stellar evolution scenarios in which the pulsational pair-instability process takes place. On the other hand, the very massive sub-population, which accounts for at most several percents of the total, may be formed through hierarchical mergers or other processes., Comment: Accepted for publication in ApJ (This manuscript incorporating the GWTC-2 data is to replace arXiv:2009.03854 that was based on the GWTC-1 data.)

GW190426_152155 was recently reported as one of the 39 candidate gravitational wave (GW) events in \citet{2020arXiv201014527A}, which has an unusual source-frame chirp mass $\sim 2.4M_{\odot}$ and may be the first GW signal from a neutron star-black hole (NSBH) merger. Assuming an astrophysical origin, we reanalyze GW190426_152155 using several waveforms with different characteristics, and consider two different priors for the mass ratio of the binary (Uniform and LogUniform). We find that the results are influenced by the priors of mass ratio, and this candidate could also be from the merger of two low mass black holes (BH). In the case for a binary black hole (BBH) merger, the effective spin is likely negative and the effective precession spin is non-negligible. As for the NSBH merger, supposing the mass of the light object follow the distribution of current neutron stars (NSs) with a reasonably measured/constrained mass, the spin of the low mass BH is so small that is hard to generate bright electromagnetic emission. Finally, we estimate a merger rate of GW190426\_152155-like systems to be $59^{+137}_{-51}~{\rm Gpc}^{-3}~{\rm yr}^{-1}$., Comment: 10 pages, 5 figures

The equation of state (EoS) of the neutron star (NS) matter remains an enigma. In this work we perform the Bayesian parameter inference with the gravitational wave data (GW170817) and mass-radius observations of some NSs (PSR J0030+0451, PSR J0437-4715, and 4U 1702-429) using the phenomenologically constructed EoS models to search for a potential first-order phase transition. Our phenomenological EoS models take the advantages of current widely used parametrizing methods, which are flexible enough to resemble various theoretical EoS models. We find that the current observation data are still not informative enough to support/rule out phase transition, due to the comparable evidences for models with and without phase transition. However, the bulk properties of the canonical $1.4\,M_\odot$ NS and the pressure at around $2\rho_{\rm sat}$ are well constrained by the data, where $\rho_{\rm sat}$ is the nuclear saturation density. Moreover, strong phase transition at low densities is disfavored, and the $1\sigma$ lower bound of transition density is constrained to $1.84\rho_{\rm sat}$., Comment: 9 pages, 5 figures, published in PRD

With the black hole mass function (BHMF; assuming an exponential cutoff at a mass of $\sim 40\,M_\odot$) of coalescing binary black hole systems constructed with the events detected in the O1 run of the advanced LIGO/Virgo network, Liang et al.(2017) predicted that the birth of the lightest intermediate mass black holes (LIMBHs; with a final mass of $\gtrsim 100\,M_\odot$) is very likely to be caught by the advanced LIGO/Virgo detectors in their O3 run. The O1 and O2 observation run data, however, strongly favor a cutoff of the BHMF much sharper than the exponential one. In this work we show that a power-law function followed by a sudden drop at $\sim 40\,M_\odot$ by a factor of $\sim $a few tens and then a new power-law component extending to $\geq 100M_\odot$ are consistent with the O1 and O2 observation run data. With this new BHMF, quite a few LIMBH events can be detected in the O3 observation run of advanced LIGO/Virgo. The first LIMBH born in GW190521, an event detected in the early stage of the O3 run of advanced LIGO/Virgo network, provides additional motivation for our hypothesis., Comment: 9 pages, 4 figures

Low energy cosmic rays are modulated by the solar activity when they propagation in the heliosphere, leading to ambiguities in understanding their acceleration at sources and propagation in the Milky Way. By means of the precise measurements of the $e^-$, $e^+$, $e^-+e^+$, and $e^+/(e^-+e^+)$ spectra by AMS-02 near the Earth, as well as the very low energy measurements of the $e^-+e^+$ fluxes by Voyager-1 far away from the Sun, we derive the local interstellar spectra (LIS) of $e^-$ and $e^+$ components individually. Our method is based on a non-parametric description of the LIS of $e^-$ and $e^+$ and a force-field solar modulation model. We then obtain the evolution of the solar modulation parameters based on the derived LIS and the monthly fluxes of cosmic ray $e^-$ and $e^+$ measured by AMS-02. {\bf To better fit the monthly data, additional renormalization factors for $e^-$ and $e^+$ have been multiplied to the modulated fluxes.} We find that the inferred solar modulation parameters of positrons are in good agreement with that of cosmic ray nuclei, and the time evolutions of the solar modulation parameters of electrons and positrons differ after the reversal of the heliosphere magnetic field polarity, which shows clearly the charge-sign dependent modulation effect., Comment: 12 pages, 8 figures

Gamma-ray bursts (GRBs) are known to be the most violent explosions in the universe, and a variety of correlations between observable GRB properties have been proposed in literature, but none of these correlations is valid for both long GRBs and short GRBs. In this paper we report the discovery of a universal correlation which is suitable for both long and short GRBs using three prompt emission properties of GRBs, i.e. the isotropic peak luminosity $L_{\rm iso}$, the peak energy of the time-integtated prompt emission spectrum $E_{\rm peak}$, and the "high signal" timescale $T_{\rm 0.45}$, $L_{\rm iso} \propto E_{\rm peak}^{1.94} T_{0.45}^{0.37}$. This universal correlation just involves properties of GRB prompt emission and does not require any information of afterglow phase, which can be used as a relatively unbiased redshift estimator. Here we use this correlation to estimate the pseudo-redshifts for short Gamma Ray Bursts and then use Lynden-Bell method to obtain a non-parametric estimate of their luminosity function and formation rate. The luminosity function is $\psi(L_0)\propto{L_0^{-0.63\pm{0.07}}}$ for dim SGRBs and $\psi(L_0)\propto{L_0^{-1.96\pm{0.28}}}$ for bright SGRBs, with the break point $6.95_{-0.76}^{+0.84}\times10^{50} erg/s$. The local formation rate of SGRBs is about 15 events $\rm Gpc^{-3}yr^{-1}$ . This universal correlation may have important implications for GRB physics, implying that the long and short GRBs should share similar radiation processes., Comment: 13 pages, 6 figures, accepted for publication in ApJ

The post-merger gravitational wave (GW) radiation of the remnant formed in the binary neutron star (BNS) coalescence has not been directly measured, yet. We show in this work that the properties of the BNS involved in GW170817, additionally constrained by PSR J0030+0451, the lower limit on the maximum gravitational mass of non-rotating neutron star (NS) and some nuclear data, are in favor of strong post-merger GW radiation. This conclusion applies to the mergers of Galactic BNS systems as well. Significant post-merger GW radiation is also preferred to improve the consistency between the maximum gravitational mass of the non-rotating NS inferred from GW170817/GRB170817A/AT2017gfo and the latest mass measurements of pulsars. The prominent post-merger gravitational radiation of GW170817-like events are expected to be detectable by advanced LIGO/Virgo detectors in the next decade and then shed valuable lights on the properties of the matter in the extremely high density., Comment: 9 pages, 4 figures, Published in ApJ

The discovery of gravitational waves from compact objects coalescence opens a brand-new window to observe the universe. With more events being detected in the future, statistical examinations would be essential to better understand the underlying astrophysical processes. In this work we investigate the prospect of measuring the mass function of black holes that are merging with the neutron stars. Applying Bayesian parameter estimation for hundreds of simulated neutron star$-$black hole (NSBH) mergers, we find that the parameters for most of the injected events can be well recovered. We also take a Bayesian hierarchical model to reconstruct the population properties of the masses of black holes, in the presence of a low mass gap, both the mass gap and power-law index ($\alpha$) of black hole mass function can be well measured, thus we can reveal where the $\alpha$ is different for binary black hole (BBH) and NSBH systems. In the absence of a low mass gap, the gravitational wave data as well as the electromagnetic data can be used to pin down the nature of the merger event and then measure the mass of these very light black holes. However, as a result of the misclassification of BBH into NSBH, the measurement of $\alpha$ is more challenging and further dedicated efforts are needed., Comment: 12 pages, 6 figures, 2 tables, ApJ published

GW190425 is the second neutron star merger event detected by the Advanced LIGO/Virgo detectors. If interpreted as a double neutron star merger, the total gravitational mass is substantially larger than that of the binary systems identified in the Galaxy. In this work we analyze the gravitational-wave data within the neutron star$-$black hole merger scenario. For the black hole, we yield a mass of $2.40^{+0.36}_{-0.32}M_\odot$ and an aligned spin of $0.141^{+0.067}_{-0.064}$. As for the neutron star we find a mass of $1.15^{+0.15}_{-0.13}M_\odot$ and the dimensionless tidal deformability of $1.4^{+3.8}_{-1.2}\times 10^{3}$. These parameter ranges are for 90\% credibility. The inferred masses of the neutron star and the black hole are not in tension with current observations and we suggest that GW190425 is a viable candidate of a neutron star$-$black hole merger event. Benefitting from the continual enhancement of the sensitivities of the advanced gravitational detectors and the increase of the number of the observatories, similar events are anticipated to be much more precisely measured in the future and the presence of black holes below the so-called mass gap will be unambiguously clarified. If confirmed, the mergers of neutron stars with (quickly rotating) low-mass black holes are likely important production sites of the heaviest r-process elements., Comment: 7 pages, 4 figures, ApJL published

Assuming that the differential rotation of the massive neutron star (NS) formed in the double NS (DNS) mergers has been effectively terminated by the magnetic braking and a uniform rotation has been subsequently established (i.e., a supramassive NS is formed), we analytically derive in this work an approximated expression for the critical total gravitational mass ($M_{\rm tot,c}$) to form supramassive NS (SMNS) in the DNS mergers, benefited from some equation of state (EoS) insensitive relationships. The maximum gravitational mass of the nonrotating NSs ($M_{\rm TOV}$) as well as the dimensionless angular momentum of the remnant ($j$) play the dominant roles in modifying $M_{\rm tot,c}$, while the radius and mass differences of the premerger NSs do not. The GW170817/GRB 170817A/AT2017gfo observations have provided so far the best opportunity to quantitatively evaluate $M_{\rm TOV}$. Supposing the central engine for GRB 170817A is a black hole quickly formed in the collapse of an SMNS, we find $M_{\rm TOV}=2.13^{+0.09}_{-0.08}M_\odot$ (68.3% credibility interval, including also the uncertainties of the EoS insensitive relationships), which is consistent with the constraints set by current NS mass measurements., Comment: 13 pages, 11 figures, initially submitted on Sept. 25, 2019

Very recently the NICER collaboration has published the first-ever accurate measurement of mass and radius together for PSR J0030+0451, a nearby isolated quickly-rotating neutron star (NS). In this work we set the joint constraints on the equation of state (EoS) and some bulk properties of NSs with the data of PSR J0030+0451, GW170817 and some nuclear experiments. The piecewise polytropic expansion method and the spectral decomposition method have been adopted to parameterize the EoS. The resulting constraints are consistent with each other. Assuming the maximal gravitational mass of non-rotating NS $M_{\rm TOV}$ lies between $2.04 \rm M_{\odot}$ and $2.4 \rm M_{\odot}$, with the piecewise method the pressure at twice nuclear saturation density is measured to be $3.38^{+2.43}_{-1.50}\times 10^{34}~{\rm dyn~cm^{-2}}$ at the $90\%$ level. For a NS with canonical mass of $1.4 \rm M_\odot$, we have the moment of inertia $I_{1.4} = {1.43}^{+0.28}_{-0.13} \times 10^{38}~{\rm kg \cdot m^2}$, tidal deformability $\Lambda_{1.4} = 390_{-140}^{+320}$, radius $R_{1.4} = 12.2_{-0.9}^{+1.0}~{\rm km}$, and binding energy $BE_{1.4} = {0.16}^{+0.01}_{-0.02} \rm M_{\odot}$ at the $90\%$ level, which are improved in comparison to the constraints with the sole data of GW170817., Comment: 9 pages, 3 figures, 2 tables, ApJ accepted

For some neutron stars (NSs) in the binary systems, the masses have been accurately measured. While for the isolated neutron stars (INSs), no mass measurement has been reported yet. The situation will change soon thanks to the successful performance of the Neutron Star Interior Composition Explorer (NICER), with which the radius and mass of the isolated PSR J0030+0451 can be simultaneously measured. For most INSs, no mass measurements are possible for NICER because of observational limitations. Benefiting from recent significant progress made on constraining the equation of state of NSs, in this work we propose a way to estimate the masses of the INSs with the measured gravitational redshifts. We apply our method to RX J1856.5-3754, RX J0720.4-3125, and RBS 1223, three members of "The Magnificent Seven" (M7), and estimate their masses to be $1.24_{-0.29}^{+0.29}M_{\odot}$, $1.23_{-0.05}^{+0.10}M_{\odot}$, and $1.08_{-0.11}^{+0.20}M_{\odot}$, respectively. These masses are consistent with that of binary NS systems, suggesting no evidence for experiencing significant accretion of these isolated objects., Comment: 9 pages, 7 figures, 3 tables, ApJ published

In this work we parameterize the Equation of State of dense neutron star (NS) matter with four pressure parameters of $\{\hat{p}_1, \hat{p}_2, \hat{p}_3, \hat{p}_4\}$ and then set the combined constraints with the data of GW 170817 and the data of 6 Low Mass X-ray Binaries (LMXBs) with thermonuclear burst or alternatively the symmetry energy of the nuclear interaction. We find that the nuclear data effectively narrow down the possible range of $\hat{p}_1$, the gravitational wave data plays the leading role in bounding $\hat{p}_2$, and the LMXB data as well as the lower bound on maximal gravitational mass of non-rotating NSs govern the constraints on $\hat{p}_3$ and $\hat{p}_4$. Using posterior samples of pressure parameters and some universal relations, we further investigate how the current data sets can advance our understanding of tidal deformability ($\Lambda$), moment of inertia ($I$) and binding energy ($BE$) of NSs. For a canonical mass of $1.4M_\odot$, we have $I_{1.4} = {1.43}^{+0.30}_{-0.13} \times 10^{38}~{\rm kg \cdot m^2}$, $\Lambda_{1.4} = 390_{-210}^{+280}$ , $R_{1.4} = 11.8_{-0.7}^{+1.2}~{\rm km}$ and $BE_{1.4} = {0.16}^{+0.01}_{-0.02} M_{\odot}$ if the constraints from the nuclear data and the gravitational wave data have been jointly applied. For the joint analysis of gravitational wave data and the LMXB data, we have $I_{1.4} = {1.28}^{+0.15}_{-0.08} \times 10^{38}~{\rm kg \cdot m^2}$, $\Lambda_{1.4} = 220_{-90}^{+90}$, $R_{1.4} = 11.1_{-0.6}^{+0.7}~{\rm km}$ and $BE_{1.4} = {0.18}^{+0.01}_{-0.01} M_{\odot}$. These results suggest that the current constraints on $\Lambda$ and $R$ still suffer from significant systematic uncertainties, while $I_{1.4}$ and $BE_{1.4}$ are better constrained., Comment: 12 pages, 5 figures, 6 tables, ApJ published. Corrected the typo error between source frame masses and detector frame masses in section 2.5

The $\gamma$-ray observation of interstellar gas provides a unique way to probe the cosmic rays (CRs) outside the solar system. In this work, we use an updated version of Fermi-LAT data and recent multi-wavelength tracers of interstellar gas to re-analyze a mid-latitude region in the third Galactic quadrant and estimate the local CR proton spectrum. Two $\gamma$-ray production cross section models for $pp$ interaction, the commonly used one from Kamae et al. (2006) and the up-to-date one from Kafexhiu et al. (2014), are adopted separately in the analysis. Both of them can well fit the emissivity and the derived proton spectra roughly resemble the direct measurements from AMS-02 and Voyager 1, but rather different spectral parameters are indicated. A break at $4\pm1~{\rm GeV}\;c^{-1}$ is shown if the cross section model by Kamae et al. (2006) is adopted. The resulting spectrum is $\lesssim 20\%$ larger than the AMS-02 observation above $15~\rm GeV$ and consistent with the de-modulated spectrum within $2\%$. The proton spectrum based on the cross section model of Kafexhiu et al. (2014) is about $1.4-1.8$ times that of AMS-02 at $2-100~\rm GeV$, however the difference decreases to $20\%$ below $10~\rm GeV$ with respect to the de-modulated spectrum. A spectral break at $20\pm11~{\rm GeV}\;c^{-1}$ is required in this model. An extrapolation down to $300~\rm MeV$ is performed to compare with the observation of Voyager 1, and we find a deviation of $\lesssim 2.5\sigma$ for both the models. In general, an approximately consistent CR spectrum can be obtained using $\gamma$-ray observation nowadays, but we still need a better $\gamma$-ray production cross section model to derive the parameters accurately., Comment: 10 pages, 6 figures. Accepted for publication in PRD

For on-axis typical short gamma-ray bursts (sGRBs), the forward shock emission is usually so bright that renders the identification of kilonovae (also known as macronovae) in the early afterglow ($t<0.5$ d) phase rather challenging. This is why previously no thermal-like kilonova component has been identified at such early time except in the off-axis dim GRB 170817A associated with GW170817. Here we report the identification of an unusual optical radiation component in GRB 070809 at $t\sim 0.47$ d, thanks plausibly to the very-weak/subdominant forward shock emission. The optical emission with a very red spectrum is well in excess of the extrapolation of the X-ray emission that is distinguished by an unusually hard spectrum, which is at odds with the forward shock afterglow prediction but can be naturally interpreted as a kilonova. Our finding supports the speculation that kilonovae are ubiquitous , and demonstrates the possibility of revealing the neutron star merger origin with the early afterglow data of some typical sGRBs that take place well beyond the sensitive radius of the advanced gravitational wave detectors and hence the opportunity of organizing dedicated follow-up observations for events of interest., Comment: 20 pages, 5 figures, published in Nature Astronomy

GW170817, the first neutron star merger event detected by advanced LIGO/Virgo detectors, was associated with an underluminous short duration GRB 170817A. In this work we compare the forward shock afterglow emission of GW170817/GRB 170817A to other luminous short GRBs (sGRBs) with both a known redshift and an afterglow emission lasting at least one day after the burst. In the rapid decay phase, the afterglow emission of the bright sGRBs and GW170817/GRB 170817A form a natural and continuous sequence, though separated by an observation time gap. If viewed on-axis, the forward shock afterglow emission of GW170817/GRB 170817A would be among the brightest ones detected so far. This provides a strong evidence for the GW170817-like merger origin of bright sGRBs, and suggests that the detection of the forward shock afterglow emission of most neutron star merger events are more challenging than the case of GW170817 since usually the mergers will be more distant and the viewing angles are plausibly higher., Comment: 7 pages, 3 figures

As revealed recently by the modeling of the multi-wavelength data of the emission following GW170817/GRB 170817A, there was an off-axis energetic relativistic outflow component launched by this historic double neutron star merger event. In this work we use the results of these modeling to examine the energy extraction process of the central engine. We show that the magnetic process (i.e., the Blandford-Znajek mechanism) is favored, while the neutrino process usually requires a too massive accretion disk if the duration of the central engine activity is comparable to the observed $T_{90}$ of GRB 170817A, unless the timescale of the central engine activity is less than $\sim$ 0.2s. We propose that the GRB observations are helpful to constrain the combined tidal parameter $\tilde{\Lambda}$, and by adopting the accretion disk mass distribution estimated in BZ mechanism, the $90\%$ credible interval of $\tilde{\Lambda}$ for the progenitor of GW170817 is inferred as $309-954$., Comment: 10 pages, 5 figures, 2 tables, accepted for publication in ApJ. Comments are welcome

The precise observations of Galactic cosmic ray fluxes of the secondary family, such as Li, Be, B, are expected to have significant implications on our understanding of the cosmic ray origin and propagation. Here we employ the recent very precise measurements of those species by the Alpha Magnetic Spectrometer on the International Space Station, together with their parent species (C and O), as well as the data collected by the Voyager-1 spacecraft outside the heliosphere and the Advanced Composition Explorer, to investigate the propagation of cosmic rays in the Milky Way. We consider the diffusion of cosmic rays plus reacceleration or convection effect during the propagation, and find that the reacceleration model can fit the data significantly better than the convection model. We further find that for the reacceleration model, the spectral hardenings of both the primary and secondary particles can be well described by the injection hardening without including additional propagation hardening. This is due to that the reacceleration effect results in a steeper secondary-to-primary ratio at low energies, and can thus naturally reproduce the fact that the secondary spectra harden more than the primary spectra found by AMS-02., Comment: 12 pages, 3 figures, 1 table; accepted by JCAP

The acceleration of cosmic ray particles and their propagation in the Milky Way and the heliosphere tangle with each other, leading to complexity and degeneracy of the modeling of Galactic cosmic rays (GCRs). The recent measurements of the GCR spectra by Voyager-1 from outside of the heliosphere gave a first direct observation of GCRs in the local interstellar (LIS) environment. Together with the high-precision data near the Earth taken by ACE and AMS-02, we derive the LIS spectra of Helium, Lithium, Beryllium, Boron, Carbon, and Oxygen nuclei from a few MeV/n to TeV/n, using a non-parameterization method. These LIS spectra are helpful in further studying the injection and propagation parameters of GCRs. The nearly 20 years of data recorded by ACE are used to determine the solar modulation parameters over the solar cycles 23 and 24, based on the force-field approximation. We find general agreements of the modulation potential with the results inferred from neutron monitors and other cosmic ray data., Comment: 16 pages, 10 figures, 9 tables, accepted for publication in ApJ

So far, in total 15 double neutron star systems (DNSs) with a reliable measurement of the total gravitational mass ($M_{\rm T}$) have been detected in the Galaxy. In this work we study the distribution of $M_{\rm T}$. The data prefer the double Gaussian distribution over a single Gaussian distribution and the low and high mass populations center at $M_{\rm T}\sim 2.58M_\odot$ and $\sim 2.72M_\odot$, respectively. The progenitor stars of GW170817 have a $M_{\rm T}=2.74^{+0.04}_{-0.01}M_\odot$, falling into the high mass population. With a local neutron star merger rate of $\sim 10^{3}~{\rm Gpc^{-3}~yr^{-1}}$, supposing the $M_{\rm T}$ of those merging neutron stars also follow the double Gaussian distribution, the upcoming runs of the advanced LIGO/Virgo will soon detect some events with a $M_{\rm T}\lesssim 2.6M_\odot$ that can effectively probe the equation of state of the neutron stars and the distribution function is expected to be reliably reconstructed in the next decade., Comment: 6 pages,1 table and 3 figures

The origin of the IceCube neutrinos is still an open question. Upper limits from diffuse gamma-ray observations suggest that the neutrino sources are either distant or hidden from gamma-ray observations. It is possible that the neutrinos are produced in jets that are formed in the core-collapsing massive stars and fail to break out, the so-called choked jets. We study neutrinos from the jets choked in the hydrogen envelopes of red supergiant stars. Fast photo-meson cooling softens the neutrino spectrum, making it difficult to explain the PeV neutrinos observed by IceCube in a one-component scenario, but a two-component model can explain the spectrum. Furthermore, we predict that a newly born jet-driven type-II supernova may be observed to be associated with a neutrino burst detected by IceCube., Comment: 11 pages, 5 figures, 1 table, accepted to ApJ