Source author record

Yi-Zhong Fan

Yi-Zhong Fan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

80works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

80 published item(s)

preprint2026arXiv

Current Status of Inert Higgs Dark Matter with Dark Fermions

The precision measurements of the muon magnetic moment and the $W$ boson mass have sparked interest in the potential deviations from standard model (SM) predictions. While it may be premature to attribute any excesses in these precision measurements to new physics, they do offer a valuable indication of potential directions for physics beyond the SM. Additionally, the particle nature of dark matter (DM) remains a crucial enigma. Despite the absence of any definitive DM signal in direct detection and collider experiments, the Galactic Center GeV $γ$-ray excess and the AMS-02 antiproton ($\overline{p}$) excess could potentially offer hints related to the evidence of DM. Motivated by these observations, we propose a simple DM model that addresses all these issues. This model extends the SM by incorporating singlet and doublet Dirac fermion fields, along with a doublet complex scalar field. For the viable parameter regions in this model, we find that future upgrades of the Large Hadron Collider and DM direct detection experiments can only partially probe them, while future high-energy muon colliders hold promise for exploring the unexplored parameter space.

preprint2025arXiv

Ultrahigh-Energy Gamma-ray Emission Associated with Black Hole-Jet Systems

Black holes (BH), one of the most intriguing objects in the universe, can manifest themselves through electromagnetic radiation initiated by the accretion flow. Some stellar-mass BHs drive relativistic jets when accreting matter from their companion stars, forming microquasars. Non-thermal emission from the radio to tera-electronvolt (TeV) gamma-ray band has been observed from microquasars, indicating the acceleration of relativistic particles. Here we report detection of four microquasars (SS 433, V4641 Sgr, GRS 1915+105, MAXI J1820+070) of spectrum extending to the ultrahigh-energy (UHE; photon energy $E>100$ TeV) band and one microquasar (Cygnus X-1) of spectrum approaching 100 TeV, using the Large High Altitude Air Shower Observatory (LHAASO). Notably, the total emission associated with SS 433 cannot be interpreted with a single leptonic component. In the UHE band, its emission is in spatial coincidence with a giant atomic cloud, which is consistent with a hadronic origin. An elongated source is discovered from V4641 Sgr with the spectrum continuing up to 800 TeV. The detection of UHE gamma rays demonstrates that accreting BHs and their environments can operate as extremely efficient accelerators of particles out of 1 peta-electronvolt (PeV), suggesting microquasars to be important contributors to Galactic cosmic rays especially around the `knee' region.

preprint2025arXiv

Verification of the Black Hole Area Law with GW230814

We present an observational confirmation of Hawking's black-hole area theorem using the newly released gravitational-wave data from the GWTC-4.0. We analyze the high signal-to-noise ratio binary black hole (BBH) merger GW230814 and measure the (total) horizon area of the black holes before and after the merger. For preferred (and reasonable) choices of the post-truncation start time, the horizon area of the remnant black hole is found to be greater than the total horizon area of the two pre-merger black holes at a high possibility (at least $\gtrsim 99.5\%$). Importantly, our analysis accounts for sky-location uncertainty. These results provide a stringent observational confirmation of the black-hole area law, further bolstering the validity of classical general relativity in the dynamical, strong-field regime.

preprint2023arXiv

Potential Subpopulations and Assembling Tendency of the Merging Black Holes

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}μ_{\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.

preprint2022arXiv

Divergence in mass ratio distributions between low-mass and high-mass coalescing binary black holes

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.

preprint2022arXiv

Explanation of nearby SNRs for primary electron excess and proton spectral bump

Several groups have reported a possible excess of primary electrons at high energies with the joint fit of the positron fraction and total electron/positron spectra. With the latest release of high-precision electron/positron spectra measured by AMS-02, we further confirm this excess by fitting $ΔΦ$ $\rm(i.e., Φ_{e^-}-Φ_{e^+})$ data in this work. Then we investigate the contribution of a single nearby supernova remnant to the primary electron excess and find that Monogem can reasonably account for this excess. Moreover, we predict that the electron spectrum may harden again at a few TeVs due to Vela's contribution. DAMPE, which can accurately measure electrons at TeV scale, is expected to provide the robust test of this new spectral feature in the near future. Finally, we fit the proton spectrum data of DAMPE with Monogem or Loop I. We find that both the primary electron excess and the proton spectral bump could be mainly generated by Monogem.

preprint2022arXiv

GB6 J2113+1121: A multi-wavelength flaring gamma-ray blazar temporally and spatially coincident with the neutrino event IceCube-191001A

A radio-emitting tidal disruption event (AT2019dsg) is proposed as a likely counterpart of the IceCube neutrino event IC-191001A. In this work we have revisited the {\it Fermi}-LAT data in the direction of the neutrino and confirmed no signal at the site of AT2019dsg. Instead, at the edge of the 90\% confidential level error region of this neutrino there is a $γ$-ray transient source associated with a blazar GB6 J2113+1121. In May 2019, GB6 J2113+1121 was undergoing an unprecedented $γ$-ray flare since the start of the {\it Fermi}-LAT operation, with a variability amplitude about 20-fold. Similar violent flares of GB6 J2113+1121, unobserved before, have been also detected observed in optical bands. Moreover, the blazar remained in a high flux state in the infrared bands when IceCube-191001A arrived, though its $γ$-ray and optical activities has temporally ceased. Motivated by the spatial and temporal coincidence, we suggest that GB6 J2113+1121 is a candidate of the counterpart of IC-191001A. The jet properties of GB6 J2113+1121 are investigated, which are found to be comparable with that of the neutrino-emitting blazars (candidates). A specific analysis of archival IceCube data in this direction and future observations would put a further constraint on the origin of the neutrino.

preprint2022arXiv

Inert Higgs Dark Matter for CDF-II W-boson Mass and Detection Prospects

The $W$-boson mass, which was recently measured at FermiLab with an unprecedented precision, suggests the presence of new multiplets beyond the Standard Model (SM). One of the minimal extensions of the SM is to introduce an additional scalar doublet, in which the non-SM scalars can enhance $W$-boson mass via the loop corrections. On the other hand, with a proper discrete symmetry, the lightest new scalar in the doublet can be stable and play the role of dark matter particle. We show that the inert two Higgs doublet model can naturally handle the new $W$-boson mass without violating other constraints, and the preferred dark matter mass is between $54$ and $74$ GeV. We identify three feasible parameter regions for the thermal relic density: the $SA$ co-annihilation, the Higgs resonance, and the $SS \to WW^*$ annihilation. We find that the first region can be fully tested by the HL-LHC, the second region will be tightly constrained by direct detection experiments, and the third region could yield detectable GeV gamma-ray and antiproton signals in the Galaxy that may have been observed by Fermi-LAT and AMS-02.

preprint2022arXiv

Muonphilic Dark Matter explanation of gamma-ray galactic center excess: a comprehensive analysis

The Galactic center gamma-ray excess (GCE) is a long-standing unsolved problem. One of candidate solutions, the dark matter (DM) annihilation, has been recently tested with other astrophysical observations, such as AMS-02 electron-positron spectra, Fermi Dwarf spheroidal galaxies gamma-ray data, and so on. By assuming that the DM particles annihilate purely into a normal charged fermion pair, Di Mauro and Winkle (2021) claimed that only a muon-pair is compatible with the null detection of all the corresponding astrophysical measurements and can explain GCE simultaneously. On the other hand, a muonphilic DM model may also lead to a signal in the recent Fermilab muon $g-2$ measurement or be constrained by the latest PandaX-4T limit. In this work, we comprehensively study interactions between DM and muon, including various combinations of DM and mediator spins. In agreement with GCE (not only $2μ$ but also $4μ$ final states), we test these interactions against all the thermal DM constraints. Our results show that only the parameter space near the resonance region of mediator can explain GCE and relic density simultaneously, and larger parameter spaces are still allowed if other poorly-known systematic uncertainties are included. Regardless of the DM spin, only the interactions with the spin-0 mediator can explain the recent muon $g-2$ excess on top of GCE, relic density, and other DM and mediator constraints.

preprint2022arXiv

Prospects of calibrating afterglow modeling of short GRBs with gravitational wave inclination angle measurements and resolving the Hubble tension with a GW-GRB association event

In the numerical modeling of the GRB afterglow, some approximations have been made for simplicity, and different groups developed their codes. A robust test of these models/approaches is challenging because of the lack of directly measured physical parameters. Fortunately, the viewing angle inferred from the afterglow modeling is widely anticipated to be the same as the inclination angle of the binary neutron star (BNS) mergers that can be evaluated with the gravitational wave (GW) data. Therefore in the future, it is possible to calibrate the afterglow modeling with the GW inclination angle measurements. We take three methods, including both analytical estimations and direct simulations, to project the uncertainties of the inclination angle measurements. For some BNS mergers accompanied with electromagnetic counterparts detected in the O4/O5 runs of LIGO/Virgo/KAGRA/LIGO-India detectors, we show that the inclination angle can be determined within an uncertainty of $\leq 0.1$ rad, supposing that the Hubble constant is known with an accuracy of $\leq 3\%$ and the uncertainty of Hubble flow velocity is within $\sim 1 \%$. The off-axis GRB outflow will give rise to afterglow emission, and the most energetic ones may be detectable at the distance of $\sim$ 100-200 Mpc even for a viewing angle of $\geq 0.3$ rad. Such events can thus serve as a robust test of the afterglow modeling approach. We have also evaluated the prospect of resolving the so-called Hubble tension with a single GW/GRB association event. We find out that a $\sim 3\%$ precision Hubble constant is obtainable if the uncertainty of the viewing angle can be constrained to be within $\sim 0.1$ rad, which is expected to be the case for some nearby ($\leq 250$ Mpc) bright/on-axis GRBs with a well-behaved afterglow light curve displaying a clear achromatic break at early times.

preprint2021arXiv

A systematic study of $γ$-ray flares from the Crab Nebula with Fermi-LAT: I. flare detection

Significant flares of GeV $γ$-ray emission from the Crab Nebula have been found by AGILE and Fermi-LAT years ago, indicating that extreme particle acceleration and radiation occurs in young pulsar wind nebulae. To enlarge the flare sample and to investigate their statistical properties will be very useful in understanding the nature of the $γ$-ray flares. In this paper, we investigate the flaring emission from the Crab Nebula with eleven year observations of the Fermi-LAT. We identify 17 significant flares in the light curve of the low-energy (synchrotron) component of the $γ$-ray emission. The flare rate is about 1.5 per year, without any significant change or clustering during the 11 years of the observation. We detect a special flare with an extremely long duration of nearly one month, occurred in October, 2018, with synchrotron photons up to energies of about 1GeV. The synchrotron component could be fitted by a steady power-law background and a variable flare component with an expotentially cutoff power-law spectrum, not only for individual flare but also for the combined data, which may favor a similar emission mechanism for all flares. However, we do not find a universal relation between the cutoff energy and the energy fluxes of the flares, which may reflect the complicated acceleration and/or cooling processes of the involved particles.

preprint2021arXiv

Constraining $U(1)_{L_μ-L_τ}$ charged dark matter model for muon $g-2$ anomaly with AMS-02 electron and positron data

Very recently, the Fermi-Lab reported the new experimental combined results on the magnetic momentum of muon with a 4.2$σ$ discrepancy compared with the expectation of the Standard Model \cite{Fermi_Lab}. A new light gauge boson $X$ in the $L_μ-L_τ$ model provides a good explanation for the $g-2$ anomaly. A Dirac fermion dark matter with a large $L_μ-L_τ$ charge can explain both the $g-2$ anomaly and the dark matter relic density \cite{Asai_2021}. In this work, we focus on the case that the mass of the dark matter is larger than the mass of muon (i.e. $m_Ψ > m_μ$) for which the channel $ΨΨ\rightarrow μ^- μ^+$ opens. Although the cross section $(σv)_{μ^{-}μ^{+}}$ is smaller by a factor of $1/q_Ψ^2$ ($q_Ψ$ represents the $L_μ-L_τ$ charge of the dark matter) compared with the channel $ΨΨ\rightarrow XX \rightarrow νν\barν\barν$, the resulting secondary electrons and positrons could imprint on their spectra above GeV energies due to the reacceleration effect of cosmic ray propagation. We use the AMS-02 measurements of electrons and positrons to constrain the annihilation cross section of the channel $ΨΨ\rightarrow μ^{-}μ^{+}$, which rules out part of the parameter space of the large $L_μ-L_τ$ charged dark matter model to account for the muon $g-2$ anomaly.

preprint2020arXiv

A kilonova associated with GRB 070809

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.

preprint2020arXiv

Awakening of two $γ$-ray high redshift flat-spectrum radio quasars in the southern hemisphere

High-redshift blazars are valuable tools to study the early Universe. So far only a handful of $γ$-ray blazars have been found at redshifts above 3. $γ$-ray signals are detected in the direction of PMN J2219-2719 ($z=3.63$) and PMN J2321-0827 ($z=3.16$) by analyzing the 10-year $Fermi$-LAT Pass 8 data. PMN J2219-2719 is not distinguished from the background in the global analysis. During the 5-month epoch, the TS value is 47.8 and the flux is more than 10 times of the 10-year averaged flux. In addition, the angular distance between the $γ$-ray position and the radio position of PMN J2219-2719 is only ${0.04}^{\circ}$. Moreover, the $γ$-ray and infrared light curves of long time scale are very similar, which support the association between the $γ$-ray source and PMN J2219-2719. The global analysis of PMN J2321-0827 suggest a new $γ$-ray source, during the flare phase, the TS value is 61.4 and the $γ$-ray flux increased significantly. The association probability suggests that PMN J2321-0827 may be the counterpart of the new $γ$-ray source. In the future, the number of high-redshift $γ$-ray sources will increase by combining $Fermi$-LAT and the upcoming Large Synoptic Survey Telescope.

preprint2020arXiv

Black Hole Mass Function of Coalescing Neutron Star-Black Hole Binary Systems: The Prospect of Reconstruction with the Gravitational Wave Observations

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 ($α$) of black hole mass function can be well measured, thus we can reveal where the $α$ 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 $α$ is more challenging and further dedicated efforts are needed.

preprint2020arXiv

Detections of simultaneous brightening of gamma-ray and optical emissions of a distant blazar GB 1508+5714 at redshift 4.3

GB 1508+5714 is a high-redshift blazar ($z$ = 4.3) and a spectrally soft $γ$-ray source has been detected in its direction. By analyzing 11.4-yr {\it Fermi}-LAT data, significant long-term variability of the $γ$-ray source is confirmed. More importantly, a $γ$-ray emission flare appeared in an epoch of several tens of days in year 2018, when the flux is about four times of the value from the global fit. Meanwhile, optical flares displayed in both $r$- and $i$-bands from the {\it Zwicky} Transient Facility light curves. Detections of the simultaneous $γ$-ray and optical brightening provide a decisive evidence to pin down the association between the $γ$-ray source and GB 1508+5714, which makes it the {\it first} identified $γ$-ray blazar beyond redshift 4. A broadband spectral energy distribution in the high flux state is constructed and the origin of the multiwavelength brightening is also briefly discussed. Upcoming wide-deep-fast optical telescopes together with the $γ$-ray sky surveyors will shed lights on the role that the AGN jets play in the early cosmic time.

preprint2020arXiv

Estimating the maximum gravitational mass of nonrotating neutron stars from the GW170817/GRB 170817A/AT2017gfo observation

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.

preprint2020arXiv

GW190521 and the GWTC-1 Events: Implication on the Black Hole Mass Function of Coalescing Binary Black Hole Systems

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.

preprint2020arXiv

Interpretation of the cosmic ray positron and electron excesses with an annihilating-decaying dark matter scenario

The precise measurements of energy spectra of cosmic ray positrons and/or electrons by recent experiments show clear excesses above 10 GeV. Moreover, a potential sharp spectral feature was suggested by the Dark Matter Particle Explorer (DAMPE) data. These results inspire quite a number of discussions on the connection with either the annihilation/decay of dark matter (DM) or the astrophysical origins. Here we discuss a DM scenario in which DM particles could annihilate and decay into standard model particle pairs simultaneously. In this model, the peak structure is due to the DM annihilation in a nearby subhalo and the broad positron/electron excesses are due to the decay of DM in the Milky Way. This model can reasonably explain the DAMPE and AMS-02 data of the total $e^+e^-$ spectra and the positron fraction, with model parameters being consistent with existing constraints. A simple realization of such a DM model is the spin-1 vector DM model.

preprint2020arXiv

Is GW190425 consistent with being a neutron star$-$black hole merger?

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.

preprint2020arXiv

Maximum mass cutoff in the neutron star mass distribution and the prospect of forming supramassive objects in the double neutron star mergers

The sample of neutron stars with a measured mass is growing quickly. With the latest sample, we adopt both a flexible Gaussian mixture model and a Gaussian plus Cauchy-Lorentz component model to infer the mass distribution of neutron stars and use the Bayesian model selection to explore evidence for multimodality and a sharp cutoff in the mass distribution. The two models yield rather similar results. Consistent with previous studies, we find evidence for a bimodal distribution together with a cutoff at a mass of $M_{\rm max}=2.26_{-0.05}^{+0.12}M_\odot$ (68% credible interval; for the Gaussian mixture model). If such a cutoff is interpreted as the maximum gravitational mass of nonrotating cold neutron stars, the prospect of forming supramassive remnants is found to be quite promising for the double neutron star mergers with a total gravitational mass less than or equal to 2.7$M_\odot$ unless the thermal pions could substantially soften the equation of state for the very hot neutron star matter. These supramassive remnants have a typical kinetic rotational energy of approximately $1-2\times 10^{53}$ ergs. Together with a high neutron star merger rate approximately $10^{3}~{\rm Gpc^{-3}~yr^{-3}}$, the neutron star mergers are expected to be significant sources of EeV($10^{18}$eV) cosmic-ray protons.

preprint2020arXiv

Plasma dark matter and electronic recoil events in XENON1T

Dark matter might be in the form of a dark plasma in the Milky Way halo. Specifically, we consider here a hidden sector consisting of a light `dark electron' and a much heavier `dark proton', each charged under an unbroken $U(1)'$ gauge symmetry. These self-interacting dark sector particles can also interact with ordinary matter via the kinetic mixing interaction, and lead to a signal in dark matter direct detection experiments. Indeed, keV electron recoils can arise quite naturally in such models from dark electron scattering off loosely bound atomic electrons. Here we examine the recently reported XENON1T excess in the context of such a plasma dark matter model. We find that the observed excess can be explained if kinetic mixing is in the approximate range: $10^{-12} \lesssim ε\lesssim 10^{-10}$. The allowed parameter space is consistent with astrophysical and cosmological constraints and consistent also with other direct detection experiments.

preprint2020arXiv

PSR J0030+0451, GW170817 and the nuclear data: joint constraints on equation of state and bulk properties of neutron stars

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 $Λ_{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.

preprint2020arXiv

The Masses of Isolated Neutron Stars Inferred from the Gravitational Redshift Measurements

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.

preprint2020arXiv

The roles of individual advanced LIGO detectors on revealing the neutron star properties

In this work we re-analyze the data of GW170817, the first binary neutron star (BNS) merger event, in two ways, including the parameterized equation of state (EoS) method and gravitational wave (GW) parameter estimation analysis. Three sets of data, including the combined LIGO/Virgo (HLV) data, the LIGO-Hanford (H1) data and the LIGO-Livingston (L1) data, have been analyzed, respectively. We confirm the bimodal probability distribution of the tidal deformability of GW170817 with the HLV data and resolve the different contributions of H1 and L1 detectors. Our simulation reveals a tendency of ``overestimating" the tidal parameter with the decreasing signal-to-noise ratio (SNR). Such a tendency is likely caused by the dependence of the result on the prior (which is not well understood and has been assumed to be widely-distributed) in the low SNR case. Though this effect is interesting, the slightly higher SNR of H1 than L1 at frequencies above $\sim 100$ Hz is most-likely not the main reason for the lower tidal parameter region favored by H1. In view of the large fluctuation of the expected tide parameter in the case of low SNR, the different probability distributions favored by the L1 and H1 data of GW170817 are reasonable. We have also explored the dependence of expected difference of SNRs of H1 and L1 detectors on the source locations in the O3 and full sensitivity runs. Such maps are helpful in evaluating the relative powers of individual detectors on measuring the tidal deformabilities for the new double neutron star merger events.

preprint2019arXiv

Circinus Galaxy Revisted with 10 Years of Fermi-LAT Data

Circinus galaxy is a nearby composite starburst/AGN system. In this work we re-analyze the GeV emission from Circinus with 10 years of {\it Fermi}-LAT Pass 8 data. In the energy range of 1-500 GeV, the spectrum can be well fitted by a power-law model with a photon index of $Γ$ = $2.20\pm0.14$, and its photon flux is $(5.90\pm1.04) \times 10^{-10}$ photons cm$^{-2}$ s$^{-1}$. Our 0.1-500 GeV flux is several times lower than that reported in the previous literature, which is roughly in compliance with the empirical relation for star-forming and local group galaxies and might be reproduced by the interaction between cosmic rays and the interstellar medium. The ratio between the $γ$-ray luminosity and the total infrared luminosity is near the proton calorimetric limit, indicating that Circinus may be a proton calorimeter. However, marginal evidence for variability of the $γ$-ray emission is found in the timing analysis, which may indicate the activity of AGN jet. More {\it Fermi}-LAT data and future observation of CTA are required to fully reveal the origin of its $γ$-ray emission.

preprint2019arXiv

Searching for the possible signal of the photon-axionlike particle oscillation in the combined GeV and TeV spectra of supernova remnants

The conversion between photons and axionlike particles (ALPs) in the Milky Way magnetic field could result in the detectable oscillation phenomena in $γ$-ray spectra of Galactic sources. In this work, the GeV (Fermi-LAT) and TeV (MAGIC/VERITAS/H.E.S.S.) data of three bright supernova remnants (SNRs, ie. IC443, W51C and W49B) have been adopted together to search such the oscillation effect. Different from our previous analysis of the sole Fermi-LAT data of IC443, we do not find any reliable signal for the photon-ALP oscillation in the joint broadband spectrum of each SNR. The reason for the inconsistence is that in this work we use the latest revision (P8R3) of Fermi-LAT data, updated diffuse emission templates and the new version of the source catalog (4FGL), which lead to some modification of the GeV spectrum of IC443. Then we set constraints on ALP parameters based on the combined analysis of all the three sources. Though these constraints are somewhat weaker than limits from the CAST experiment and globular clusters, they are supportive of and complementary to these other results.

preprint2016arXiv

A Monte Carlo Bayesian Search for the Plausible Source of the Telescope Array Hotspot

The Telescope Array (TA) collaboration has reported a hotspot of 19 ultrahigh-energy cosmic rays (UHECRs). Using a universal model with one source and energy-dependent magnetic deflections, we show that the distribution of the TA hotspot events is consistent with a single source hypothesis, although multiple sources cannot be ruled out. The chance probability of this distribution arising from a homogeneous distribution is $0.2\%$. We describe a Monte Carlo Bayesian (MCB) inference approach, which can be used to derive parameters of the magnetic fields as well as the source coordinates, and we apply this method to the TA hotspot data, inferring the location of the likely source. We discuss possible applications of the same approach to future data.

preprint2016arXiv

Discovery of $γ$-ray emission from the radio-intermediate quasar III Zw 2: violent jet activity with intraday $γ$-ray variability

III Zw 2 is the prototype of radio-intermediate quasars. Although there is the evidence of possessing strong jet, significant $γ$-ray emission has not been reported before. In this work, we carry out a detailed analysis of the latest {\it Fermi}-LAT {\it Pass} 8 data. No significant $γ$-ray signal has been detected in the time-averaged 7-year {\it Fermi}-LAT data of III Zw 2, we however have identified two distinct $γ$-ray flares with isotropic luminosities of $\sim 10^{45}$ erg $\rm s^{-1}$. Multiwavelength data analysis (including also the optical photometric observations from Yunnan Observatories) are presented and the main finding is the simultaneous optical and $γ$-ray flares of III Zw 2 appearing in Nov. 2009. Violent $γ$-ray variability with doubling timescale of 2.5 hours was detected in another $γ$-ray flare in May 2010, for which the 3-hour $γ$-ray peak flux is $\sim 250$ times of the average flux in 7 years. Rather similar behaviors are observed in blazars and the blazar model can reasonably reproduce the spectral energy distribution of III Zw 2 in a wide energy range, strongly suggesting that its central engine resembles that of blazars. In view of its core which shares radio similarities with young radio sources together with the weak extended radio lobe emission, we suggest that III Zw 2 harbors a recurrent activity core and serves as a valuable target for investigating the fuelling and triggering of the activity in radio loud AGNs.

preprint2016arXiv

GRB/GW association: Long-short GRB candidates, time-lag, measuring gravitational wave velocity and testing Einstein's equivalence principle

Short-duration gamma-ray bursts (SGRBs) are widely believed to be powered by the mergers of compact binaries, such as binary neutron stars or possibly neutron star-black hole binaries. Though the prospect of detecting SGRBs with gravitational wave (GW) signals by the advanced Laser Interferometer Gravitational-Wave Observatory (LIGO)/VIRGO network is promising, no known SGRB has been found within the expected advanced LIGO/VIRGO sensitivity range for binary neutron star systems. We find, however, that the two long-short GRBs (GRB 060505 and GRB 060614) may be within the horizon of advanced GW detectors. In the upcoming era of GW astronomy, the merger origin of some long-short GRBs, as favored by the macronova signature displayed in GRB 060614, can be unambiguously tested. The model-dependent time lags between the merger and the onset of the prompt emission of the GRB are estimated. The comparison of such time lags between model predictions and the real data expected in the era of the GW astronomy would be helpful in revealing the physical processes taking place at the central engine (including the launch of the relativistic outflow, the emergence of the outflow from the dense material ejected during the merger, and the radiation of gamma rays). We also show that the speed of GWs, with or without a simultaneous test of Einstein's equivalence principle, can be directly measured to an accuracy of $\sim 3\times 10^{-8}~{\rm cm~s^{-1}}$ or even better in the advanced LIGO/VIRGO era. The Astrophysical Journal, Volume

preprint2016arXiv

Implications of the Tentative Association between GW150914 and a {\it Fermi}-GBM Transient

The merger-driven Gamma-ray Bursts (GRBs) and their associated gravitational wave (GW) radiation, if both successfully detected, have some far-reaching implications, including for instance: (i) The statistical comparison of the physical properties of the short/long-short GRBs with and without GW detection can test the general origin model; (ii) Revealing the physical processes taking place at the central engine; (iii) Measuring the velocity of the Gravitational wave directly/accurately. In this work we discuss these implications in the case of possible association of GW150914/ GBM transient 150914. We compared GBM transient 150914 with other SGRBs and found that such an event {may be} a distinct outlier in some statistical diagrams, possibly due to its specific binary-black-hole merger origin. However, the presence of a "new" group of SGRBs with "unusual" physical parameters is also possible. If the outflow of GBM transient 150914 was launched by the accretion onto the nascent black hole, the magnetic activity rather than the neutrino process is likely responsible for the energy extraction and the accretion disk mass is estimated to be $\sim 10^{-5}~M_\odot$. The GW150914/GBM transient 150914 association, {if confirmed, would} provide the first opportunity to directly measure the GW velocity and its departure from the speed of the light {should be within} a factor of $\sim 10^{-17}$.

preprint2016arXiv

Interpretations of the possible 42.7 GeV $γ$-ray line

Recently Liang et al. found a tentative line signal at about 43 GeV in the directions of sixteen nearby Galaxy Clusters. If arising from dark matter annihilation, the mass of the dark matter particles should be $m_χ\sim 43$ GeV and the annihilation cross section $\langle σv \rangle \sim 5\times10^{-28}(\overline{\rm BF}/10^{3})^{-1}~{\rm cm^3~s^{-1}}$ is needed, where $\overline{\rm BF}$ is the averaged boost factor of the annihilation signal of these Galaxy Clusters. In this work we discuss several models which could interpret these features and estimate the model parameters. Usually sizeable coupling parameter is needed. The $2σ$ upper limits on the cross section of the dark matter annihilation into various final states such as $b\bar{b}$, $μ^+μ^-$ and so on are also presented.

preprint2016arXiv

Search for a gamma-ray line feature from a group of nearby Galaxy clusters with Fermi LAT Pass 8 data

Galaxy clusters are the largest gravitationally bound objects in the universe and may be suitable targets for indirect dark matter searches. With 85 months of Fermi-LAT Pass 8 publicly available data, we analyze the gamma-ray emission in the directions of 16 nearby Galaxy Clusters with an unbinned likelihood analysis. No globally statistically-significant $γ-$ray line feature is identified and a tentative line signal may be present at $\sim 43$ GeV. The 95\% confidence level upper limits on the velocity-averaged cross section of dark matter particles annihilating into double $γ-$rays (i.e., $\langle σv \rangle_{χχ\rightarrow γγ}$) are derived. Unless very optimistic boost factors of dark matter annihilation in these Galaxy Clusters have been assumed, such constraints are much weaker than the bounds set by the Galactic $γ-$ray data.

preprint2016arXiv

Search for gamma-ray emission from eight dwarf spheroidal galaxy candidates discovered in Year Two of Dark Energy Survey with Fermi-LAT data

Very recently the Dark Energy Survey (DES) Collaboration has released their second group of Dwarf spheroidal (dSph) galaxy candidates. With the publicly-available Pass 8 data of Fermi-LAT we search for $γ-$ray emissions from the directions of these eight newly discovered dSph galaxy candidates. No statistically significant $γ-$ray signal has been found in the combined analysis of these sources. With the empirically estimated J-factors of these sources, the constraint on the annihilation channel of $χχ\rightarrow τ^{+}τ^{-}$ is comparable to the bound set by the joint analysis of fifteen previously known dSphs with kinematically constrained J-factors for the dark matter mass $m_χ>250$ GeV. In the direction of Tucana III (DES J2356-5935), one of the nearest dSph galaxy candidates that is $\sim 25$ kpc away, there is a weak $γ-$ray signal and its peak test statistic (TS) value for the dark matter annihilation channel $χχ\rightarrow τ^{+}τ^{-1}$ is $\approx 6.7$ at $m_χ\sim 15$ GeV. The significance of the possible signal likely increases with time. More data is highly needed to pin down the physical origin of such a GeV excess.

preprint2016arXiv

Search for gamma-ray line features from Milky Way satellites with Fermi-LAT Pass 8 data

With 91 months of the publicly available Fermi-LAT Pass 8 data, we analyze the gamma-ray emission from the Milky Way satellites to search for potential line signals due to the annihilation of dark matter particles into double photons. The searched targets include a sample of dwarf spheroidal galaxies, the Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC). No significant line emission has been found neither in the stacked dwarf galaxy sample nor in the direction of LMC/SMC. The corresponding upper limits on the cross section of DM annihilation into two photons are derived. Compared with results of previous gamma-ray line searches with the Pass 7 data, the current constraints on the line emission from dwarf spheroidal galaxies has been significantly improved in a wide energy range. With the rapid increase of the sample of dwarf spheroidal galaxies (candidates), we expect that the sensitivity of gamma ray line searches will be significantly improved in the near future.

preprint2016arXiv

Testing the dark matter subhalo hypothesis of the gamma-ray source 3FGL J2212.5+0703

N-body simulations predict that galaxies at the Milky Way scale host a large number of dark matter (DM) subhalos. Some of these subhalos, if they are massive enough or close enough to the Earth, might be detectable in $γ$ rays due to the DM annihilation. 3FGL J2212.5+0703, an unidentified gamma-ray source, has been suggested to be the counterpart candidate of a DM subhalo by Bertoni et al. (2015, 2016). In this work we analyze the Fermi-LAT Pass 8 data of 3FGL J2212.5+0703 to independently test the DM subhalo hypothesis of this source. In order to suppress the possible contamination from two nearby very-bright blazars, we just take into account the front-converting gamma-rays which have better angular resolutions than that of the back-converting photons. In addition to the spatial distribution analysis, we have extended the spectrum analysis down to the energies of $\sim 100$ MeV, and thoroughly examined the variability of the emission during the past 8 years. We confirm that 3FGL J2212.5+0703 is a steady and spatially-extended gamma-ray emitter at a high confidence level. The spectrum is well consistent with that expected from DM annihilation into $b\bar{b}$. The introduction of a phenomenological LogParabola spectrum just improves the fit slightly. All these results suggest that 3FGL J2212.5+0703 could be indicative of a DM subhalo.

preprint2016arXiv

The Macronova in GRB 050709 and the GRB/macronova connection

GRB 050709 was the first short Gamma-ray Burst (sGRB) with an identified optical afterglow. In this work we report a re-analysis of the publicly available data of this event and the discovery of a Li-Paczynski macronova/kilonova that dominates the optical/IR signal at t $>$ 2.5 days. Such a signal would arise from ~0.05 M$_\odot$ r-process material launched by a compact binary merger. The implied mass ejection supports the suggestion that compact binary mergers are significant and possibly main sites of heavy r-process nucleosynthesis. We have re-analyzed all afterglow data from nearby short and hybrid GRBs. A statistical study of sGRB/macronova connection reveals that macronova may have taken place in all these GRBs though the fraction as low as 0.18 cannot be ruled out. The identification of two of the three macronova candidates in the I-band implies a more promising detection prospect for ground-based surveys.

preprint2015arXiv

A possible Macronova in the late afterglow of the `long-short' burst GRB 060614

Long-duration ($>2$ s) $γ$-ray bursts that are believed to originate from the death of massive stars are expected to be accompanied by supernovae. GRB 060614, that lasted 102 s, lacks a supernova-like emission down to very stringent limits and its physical origin is still debated. Here we report the discovery of near-infrared bump that is significantly above the regular decaying afterglow. This red bump is inconsistent with even the weakest known supernova. However, it can arise from a Li-Paczyński macronova $-$ the radioactive decay of debris following a compact binary merger. If this interpretation is correct GRB 060614 arose from a compact binary merger rather than from the death of a massive star and it was a site of a significant production of heavy r-process elements. The significant ejected mass favors a black hole-neutron star merger but a double neutron star merger cannot be ruled out.

preprint2015arXiv

An $r-$process macronova/kilonova in GRB 060614: evidence for the merger of a neutron star-black hole binary

After the jet break at $t\sim 1.4$ days, the optical afterglow emission of the long-short burst GRB 060614 can be divided into two components. One is the power-law decaying forward shock afterglow emission. The other is an excess of flux in several multi-band photometric observations, which emerges at $\sim$4 days after the burst, significantly earlier than that observed for a supernova associated with a long-duration GRB. At $t>13.6$ days, the F814W-band flux drops faster than $t^{-3.2}$. Moreover, the spectrum of the excess component is very soft and the luminosity is extremely low. These observed signals are incompatible with those from weak supernovae, but the ejection of $\sim 0.1~M_\odot$ of $r-$process material from a black hole-neutron star merger, as recently found in some numerical simulations, can produce it. If this interpretation is correct, it represents the first time that a multi-epoch/band lightcurve of a Li-Paczynski macronova (also known as kilonova) has been obtained and black hole-neutron star mergers are sites of significant production of $r-$process elements.

preprint2015arXiv

Discovery of gamma-ray emission from a strongly lobe-dominated quasar 3C 275.1

We systematically analyze the 6-year {\it Fermi}/LAT data of the lobe-dominated quasars (LDQs) in the complete LDQ sample from 3CRR survey and report the discovery of high-energy $γ$-ray emission from 3C 275.1. The $γ$-ray emission of 3C 207 is confirmed and significant variability of the lightcurve is identified. We do not find statistically significant $γ$-ray emission from other LDQs. 3C 275.1 is the known $γ$-ray quasar with the lowest core dominance parameter (i.e., $R=0.11$). We also show that both the northern radio hotspot and parsec jet models can reasonably reproduce the $γ$-ray data. The parsec jet model, however, is favored by the potential $γ$-ray variability at the timescale of months. We suggest that some dimmer $γ$-ray LDQs will be detected in the future and LDQs could contribute non-negligibly to the extragalactic $γ$-ray background.

preprint2015arXiv

GeV excess in the Milky Way: The Role of Diffuse Galactic gamma ray Emission template

Several groups have analyzed the publicly-available Fermi-LAT data and reported a spatially extended $γ-$ray excess of around $1-3$ GeV from the region surrounding the Galactic Center that might originate from annihilation of dark matter particles with a rest mass $m_χ\sim 30-40$ GeV. In this work we examine the role of the diffuse Galactic gamma ray emission (DGE) templates played in suppressing the GeV excess. For such a purpose, we adopt in total 128 background templates that have been generated by Ackermann et al. \cite{FermiLAT:2012aa} in the study of the {Fermi-LAT} observations of the diffuse gamma ray emission considering the effects of cosmic rays and the interstellar medium. The possible GeV excess, assumed to follow the spatial distribution of the prompt gamma-rays produced in the annihilation of dark matter particles taking a generalized NFW profile with an inner slope $α=1.2$, has been analyzed in some regions of interest. The introduction of such an additional component centered at the Galactic center is found to have improved the goodness of fit to the data significantly in all background template models regardless of whether the excess spectrum is fixed or not. Our results thus suggest that the presence of a statistically significant GeV excess in the inner Galaxy is robust thought its spectrum depends on the DGE model adopted in the analysis. The possible physical origin of the GeV excess component is discussed and in the dark matter model the annihilation cross section of such particles is evaluated.

preprint2015arXiv

The `excess' of primary cosmic ray electrons

With the accurate cosmic ray (CR) electron and positron spectra (denoted as $Φ_{\rm e^{-}}$ and $Φ_{\rm e^{+}}$, respectively) measured by AMS-02 collaboration, the difference between the electron and positron fluxes (i.e., $ΔΦ=Φ_{\rm e^{-}}-Φ_{\rm e^{+}}$), dominated by the propagated primary electrons, can be reliably inferred. In the standard model, the spectrum of propagated primary CR electrons at energies $\geq 30$ GeV softens with the increase of energy. The absence of any evidence for such a continuous spectral softening in $ΔΦ$ strongly suggests a significant `excess' of primary CR electrons and at energies of $100-400$ GeV the identified excess component has a flux comparable to that of the observed positron excess. Middle-age but `nearby' supernova remnants (e.g., Monogem and Geminga) are favored sources for such an excess.

preprint2015arXiv

The lightcurve of the macronova associated with the long-short burst GRB 060614

The {\it Swift}-detected GRB 060614 was a unique burst that straddles an imaginary divide between long- and short-duration gamma-ray bursts (GRBs), and its physical origin has been heavily debated over the years. Recently, a distinct very-soft F814W-band excess at $t\sim 13.6$ days after the burst was identified in a joint-analysis of VLT and HST optical afterglow data of GRB~060614, which has been interpreted as evidence for an accompanying Li-Paczynski macronova (also called a kilonova). Under the assumption that the afterglow data in the time interval of $1.7-3.0$ days after the burst are due to external forward shock emission, when this assumption is extrapolated to later times it is found that there is an excess of flux in several multi-band photometric observations. This component emerges at $\sim$4 days after the burst, and it may represent the first time that a multi-epoch/band lightcurve of a macronova has been obtained. The macronova associated with GRB 060614 peaked at $t\lesssim 4$ days after the burst, which is significantly earlier than that observed for a supernova associated with a long-duration GRB. Due to the limited data, no strong evidence for a temperature evolution is found. We derive a conservative estimate of the macronova rate of $\sim 16.3^{+16.3}_{-8.2}~{\rm Gpc^{-3}}{\rm yr^{-1}}$, implying a promising prospect for detecting the gravitational wave radiation from compact object mergers by upcoming Advanced LIGO/VIRGO/KAGRA detectors (i.e., the rate is ${\cal R}_{\rm GW} \sim 0.5^{+0.5}_{-0.25}(D/200~{\rm Mpc})^{3}~{\rm yr^{-1}}$).

preprint2014arXiv

Fast radio bursts as a cosmic probe?

We discuss the possibility of using fast radio bursts (FRBs), if cosmological, as a viable cosmic probe. We find out that the contribution of the host galaxies to the detected dispersion measures can be inapparent for the FRBs not from galaxy centers or star forming regions. The inhomogeneity of the intergalactic medium (IGM), however, causes significant deviation of the dispersion measure from that predicted in the simplified homogeneous IGM model for individual event. Fortunately, with sufficient FRBs along different sightlines but within a very narrow redshift interval (e.g., $Δz \sim 0.05$ or $Δz \sim 0.05(1+z)$), the mean from averaging observed dispersion measures does not suffer such a problem and hence may be used as a cosmic probe. We show that in the optimistic case (e.g., tens FRBs in each $Δz$ have been measured; the most distant FRBs were at redshift $\geq 3$; the host galaxies and the FRB sources contribute little to the detected dispersion measures) and with all the uncertainties (i.e. the inhomogeneity of the IGM, the contribution and uncertainty of host galaxies as well as the evolution and error of $f_{\rm IGM}$) considered, FRBs could help constrain the equation of state of dark energy.

preprint2014arXiv

GRB 131231A: Implications of the GeV emission

GRB 131231A was detected by the Large Area Telescope onboard Fermi Space Gamma-ray Telescope. The high energy gamma-ray ($> 100$ MeV) afterglow emission spectrum is $F_ν\propto ν^{-0.54\pm0.15}$ in the first $\sim 1300$ s after the trigger and the most energetic photon has an energy $\sim 62$ GeV arriving at $t\sim 520$ s. With reasonable parameters of the GRB outflow as well as the density of the circum-burst medium, the synchrotron radiation of electrons or protons accelerated at an external forward shock have difficulty accounting for the data. The synchrotron self-Compton radiation of the forward shock-accelerated electrons, instead, can account for both the spectrum and temporal behavior of the GeV afterglow emission. We also show that the prospect for detecting GRB 131231A$-$like GRBs with Cherenkov Telescope Array (CTA) is promising.

preprint2014arXiv

Model-dependent estimate on the connection between fast radio bursts and ultra-high energy cosmic rays

The existence of fast radio bursts (FRBs), a new type of extragalatic transients, has been established recently and quite a few models have been proposed. In this work we discuss the possible connection between the FRB sources and ultra-high energy ($>10^{18}$ eV) cosmic rays. We show that in the blitzar model and the model of merging binary neutron stars, the huge energy release of each FRB central engine together with the rather high rate of FRBs, the accelerated EeV cosmic rays may contribute significantly to the observed ones. In other FRB models including for example the merger of double white dwarfs and the energetic magnetar radio flares, no significant EeV cosmic ray is expected. We also suggest that the mergers of double neutron stars, even if they are irrelevant to FRBs, may play a non-ignorable role in producing EeV cosmic ray protons if supramassive neutron stars were formed in a good fraction of mergers and the merger rate is $\gtrsim 10^{3}~{\rm yr^{-1}~ Gpc^{-3}}$. Such a possibility will be unambiguously tested in the era of gravitational wave astronomy.

preprint2013arXiv

A supra-massive magnetar central engine for short GRB 130603B

We show that the peculiar early optical and in particular X-ray afterglow emission of the short duration burst GRB 130603B can be explained by continuous energy injection into the blastwave from a supra-massive magnetar central engine. The observed energetics and temporal/spectral properties of the late infrared bump (i.e., the "kilonova") are also found consistent with emission from the ejecta launched during an NS-NS merger and powered by a magnetar central engine. The isotropic-equivalent kinetic energies of both the GRB blastwave and the kilonova are about $E_{\rm k}\sim 10^{51}$ erg, consistent with being powered by a near-isotropic magnetar wind. However, this relatively small value demands that most of the initial rotational energy of the magnetar $(\sim {\rm a~ few \times 10^{52}~ erg})$ is carried away by gravitational wave radiation. Our results suggest that (i) the progenitor of GRB 130603B would be a NS-NS binary system, whose merger product would be a supra-massive neutron star that lasted for about $\sim 1000$ seconds; (ii) the equation-of-state of nuclear matter would be stiff enough to allow survival of a long-lived supra-massive neutron star, so that it is promising to detect bright electromagnetic counterparts of gravitational wave triggers without short GRB associations in the upcoming Advanced LIGO/Virgo era.

preprint2013arXiv

AMS-02 positron excess: new bounds on dark matter models and hint for primary electron spectrum hardening

The data collected by ATIC, CREAM and PAMELA all display remarkable cosmic-ray-nuclei spectrum hardening above the magnetic rigidity $\sim$ 240 GV. One natural speculation is that the primary electron spectrum also gets hardened (possibly at $\sim 80$ GV) and the hardening partly accounts for the electron/positron total spectrum excess discovered by ATIC, HESS and Fermi-LAT. If it is the case, the increasing behavior of the subsequent positron-to-electron ratio will get flattened and the spectrum hardening should be taken into account in the joint fit of the electron/psoitron data otherwise the inferred parameters will be biased. Our joint fits to the latest AMS-02 positron fraction data together with the PAMELA/Fermi-LAT electron/positron spectrum data suggest that the primary electron spectrum hardening is needed in most though not all modelings. The bounds on dark matter models have also been investigated. In the presence of spectrum hardening of primary electrons, the amount of dark-matter-originated electron/positron pairs needed in the modeling is smaller. Even with such a modification, the annihilation channel $χχ\rightarrow μ^{+}μ^{-}$ has been tightly constrained by the Fermi-LAT Galactic diffuse emission data. The decay channel $χ\rightarrow μ^{+}μ^{-}$ is found to be viable.

preprint2013arXiv

Constraining the distribution of dark matter in inner galaxy with indirect detection signal: The case of tentative 130 gev γ-ray line

The dark matter distribution in the very inner region of our Galaxy is still in debate. In the N-body simulations a cuspy dark matter halo density profile is favored. Several dissipative baryonic processes however are found to be able to significantly flatten dark matter distribution and a cored dark matter halo density profile is possible. The baryons dominate the gravitational potential in the inner Galaxy, hence a direct constrain on the abundance of the dark matter particles is rather challenging. Recently, a few groups have identified a tentative 130 GeV line signal in the Galactic center, which could be interpreted as the signal of the dark matter annihilation. With current 130 GeV line data and adopting the generalized Navarro-Frenk-White profile of the dark matter halo, for local dark matter density ρ_0=0.4 GeV cm^{-3} and r_s=20 kpc we obtain a 95% confidence level lower (upper) limit on the inner slope of dark matter density distribution α= 1.06 (the cross section of dark matter annihilation into gamma-rays <σv >_{χχ-> γγ}= 1.3\times 10^{-27} cm^3 s^{-1}). Such a slope is consistent with the results of some N-body simulations, and if the signal is due to dark matter, suggests that baryonic processes may be unimportant.

preprint2013arXiv

Cosmological Time Dilation in Durations of Swift Long Gamma-Ray Bursts

Cosmological time dilation is a fundamental phenomenon in an expanding universe, which stresses that both the duration and wavelength of the emitted light from a distant object at the redshift $z$ will be dilated by a factor of $1+z$ at the observer. By using a sample of 139 \emph{Swift} long GRBs with known redshift ($z\leq8.2$), we measure the observed duration ($T_{90}$) in the observed energy range between $140/(1+z)$ keV and $350/(1+z)$ keV, corresponding to a fixed energy range of 140-350 keV in the rest frame. We obtain a significant correlation between the duration and the factor $1+z$, i.e., $T_{\rm{90}}=10.5(1+z)^{0.94\pm0.26}$, which is well consistent with that expected from cosmological time dilation effect.

preprint2013arXiv

Diffuse PeV neutrino emission from Ultra-Luminous Infrared Galaxies

Ultra-luminous infrared galaxies (ULIRGs) are the most luminous and intense starburst galaxies in the Universe. Both their star-formation rate (SFR) and gas surface mass density are very high, implying a high supernovae rate and an efficient energy conversion of energetic protons. A small fraction of these supernovae is the so-called hypernovae with a typical kinetic energy ~1e52 erg and a shock velocity >=1e9 cm/s. The strong shocks driven by hypernovae are able to accelerate cosmic ray protons up to 1e17 eV. These energetic protons lose a good fraction of their energy through proton-proton collision when ejected into very dense interstellar medium, and as a result, produce high energy neutrinos (<=5 PeV). Recent deep infrared surveys provide solid constraints on the number density of ULIRGs across a wide redshift range 0<z<2.3, allowing us to derive the flux of diffuse neutrinos from hypernovae. We find that at PeV energies, the diffuse neutrinos contributed by ULIRGs are comparable with the atmosphere neutrinos with the flux of 2e-9GeV cm^-2/s/sr, by assuming the injected cosmic ray power law spectrum with an index of -2.

preprint2013arXiv

High energy emission of GRB 130427A: evidence for inverse Compton radiation

A nearby super-luminous burst GRB 130427A was simultaneously detected by six $γ$-ray space telescopes ({\it Swift}, Fermi-GBM/LAT, Konus-Wind, SPI-ACS/INTEGRAL, AGILE and RHESSI) and by three RAPTOR full-sky persistent monitors. The isotropic $γ-$ray energy release is of $\sim 10^{54}$ erg, rendering it the most powerful explosion among the GRBs with a redshift $z\leq 0.5$. The emission above 100 MeV lasted about one day and four photons are at energies greater than 40 GeV. We show that the count rate of 100 MeV-100 GeV emission may be mainly accounted for by the forward shock synchrotron radiation and the inverse Compton radiation likely dominates at GeV-TeV energies. In particular, an inverse Compton radiation origin is established for the $\sim (95.3,~47.3,~41.4,~38.5,~32)$ GeV photons arriving at $t\sim (243,~256.3,~610.6,~3409.8,~34366.2)$ s after the trigger of Fermi-GBM. Interestingly, the external-inverse-Compton-scattering of the prompt emission (the second episode, i.e., $t\sim 120-260$ s) by the forward-shock-accelerated electrons is expected to produce a few $γ-$rays at energies above 10 GeV, while five were detected in the same time interval. A possible unified model for the prompt soft $γ-$ray, optical and GeV emission of GRB 130427A, GRB 080319B and GRB 090902B is outlined. Implication of the null detection of $>1$ TeV neutrinos from GRB 130427A by IceCube is discussed.

preprint2013arXiv

High energy emission of GRB 130821A: constraining the density profile of the circum-burst medium as well as the initial Lorentz factor of the outflow

GRB 130821A was detected by Fermi-GBM/LAT, Konus-Wind, SPI-ACS/INTEGRAL, RHESSI and Mars Odyssey-HEND. Although the data of GRB 130821A are very limited, we show in this work that the high energy gamma-ray emission (i.e., above 100 MeV) alone imposes tight constraint on the density profile of the circum-burst medium as well as the initial Lorentz factor of the outflow. The temporal behavior of the high energy gamma-ray emission is consistent with the forward shock synchrotron radiation model and the circum-burst medium likely has a constant-density profile. The Lorentz factor is about a few hundred, similar to other bright GRBs.

preprint2013arXiv

Is the late near-infrared bump in short-hard GRB 130603B due to the Li-Paczynski kilonova?

Short-hard gamma-ray bursts (GRBs) are widely believed to be produced by the merger of two binary compact objects, specifically by two neutron stars or by a neutron star orbiting a black hole. According to the Li-Paczynski kilonova model, the merger would launch sub-relativistic ejecta and a near-infrared/optical transient would then occur, lasting up to days, which is powered by the radioactive decay of heavy elements synthesized in the ejecta. The detection of a late bump using the {\em Hubble Space Telescope} ({\em HST}) in the near-infrared afterglow light curve of the short-hard GRB 130603B is indeed consistent with such a model. However, as shown in this Letter, the limited {\em HST} near-infrared lightcurve behavior can also be interpreted as the synchrotron radiation of the external shock driven by a wide mildly relativistic outflow. In such a scenario, the radio emission is expected to peak with a flux of $\sim 100 μ$Jy, which is detectable for current radio arrays. Hence, the radio afterglow data can provide complementary evidence on the nature of the bump in GRB 130603B. It is worth noting that good spectroscopy during the bump phase in short-hard bursts can test validity of either model above, analogous to spectroscopy of broad-lined Type Ic supernova in long-soft GRBs.

preprint2013arXiv

PAMELA/Fermi-LAT electron cosmic ray spectrum at $\sim$100 GeV: implication for dark matter annihilation signal in accordance with the 130 GeV $γ$-ray line

Recently, a tentative 130 GeV $γ$-ray line signal was identified by quite a few groups. If correct it would constitute a ``smoking gun'' for dark matter annihilations. Interestingly, the spectra of the cosmic ray electrons detected by PAMELA and Fermi-LAT both show tiny wiggle-like structure at $\sim 100$ GeV, which might indicate a weak signal of the annihilation of $\sim 130$ GeV dark matter particles into electrons/positrons with a velocity-weighted cross section $\langleσv\rangle_{\rm χχ\rightarrow e^{+}e^{-}} \sim 4\times10^{-26}~{\rm cm^{3}~s^{-1}}$. The prospect of identifying such a potential weak dark-matter-annihilation electron and/or positron component by AMS-02, a mission to measure the high energy cosmic ray spectra with unprecedented accuracy, is investigated.

preprint2013arXiv

Signature of gravitational radiation in afterglow of short Gamma-Ray Bursts?

Short Gamma-Ray Bursts (GRBs), brief intense emission of $γ-$rays characterized by a duration shorter than 2 seconds that are plausibly powered by the coalescence of binary neutron stars, are believed to be strong gravitational wave radiation (GWR) sources. The test of such a speculation has been thought to be impossible until the performance of the detectors like advanced LIGO. Recently there has been growing evidence for the formation of highly-magnetized neutron star (i.e., magnetar) in the double neutron star mergers. In this work we re-examine the interpretation of the X-ray plateau followed by an abrupt decline detected in some short GRB afterglows within the supramassive magnetar model and find that the maximum gravitational mass of the non-rotating neutron stars is $\sim 2.3M_{\odot}$ and the observed duration of some X-ray plateaus are significantly shorter than that expected in the magnetic dipole radiation scenario, suggesting that the collapse of the supramassive magnetars has been considerably enhanced by the energy loss via GWR. Such a result demonstrates that the signature of GWR may have already existed in current electromagnetic data of short GRBs.

preprint2013arXiv

TeV-PeV neutrinos over the atmospheric background: originating from two groups of sources?

In addition to the two ~1 PeV neutrinos, the IceCube Collaboration recently reported a detection of 26 neutrino candidates at energies from 30 TeV to 250 TeV, implying a confidence level of 4.3σover the atmospheric background. We suggest that these TeV-PeV non-atmospheric neutrinos may originate from two groups of sources, motivated by the non-detection of neutrinos in the energy range 250 TeV- 1 PeV in current data. If intrinsic, the non-detection of 250 TeV-1 PeV neutrinos disfavors the single power-law spectrum model for the TeV-PeV non-atmospheric neutrinos at a confidence level of ~ 2σ. We then interpret the current neutrino data with a two-component spectrum model. One has a flat spectrum with a cutoff at the energy ~ 250 TeV and the other has a sharp peak at ~1 PeV. The former is likely via pp collision while the latter may be generated by the photomeson interaction.

preprint2012arXiv

Constraining Asymmetric Bosonic Non-interacting Dark Matter with Neutron Stars

The Hawking evaporation of small black holes formed by the collapse of dark matter at the center of neutron stars plays a key role in loosing the constraint on the mass of asymmetric bosonic non-interacting dark matter particles. Different from previous works we show that such a kind of dark matter is viable in the mass range from 3.3 GeV to ~ 10 TeV, which covers the most attractive regions, including the preferred asymmetric dark matter mass ~ 5.7 GeV as well as the 5-15 GeV range favored by DAMA and CoGeNT.

preprint2012arXiv

Fermi Large Area Telescope observations of GRB 110625A

Gamma-ray bursts (GRBs) that emit photons at GeV energies form a small but significant population of GRBs. However, the number of GRBs whose GeV-emitting period is simultaneously observed in X-rays remains small. We report gamma-ray observations of GRB 110625A using Fermi's Large Area Telescope (LAT) in the energy range 100 MeV to 20 GeV. Gamma-ray emission at these energies was clearly detected using data taken between 180s and 580s after the burst, an epoch after the prompt emission phase. The GeV light curve differs from a simple power-law decay, and probably consists of two emission periods. Simultaneous Swift/XRT observations did not show flaring behaviors as in the case of GRB 100728A. We discuss the possibility that the GeV emission is the synchrotron self-Compton radiation of underlying ultraviolet flares.

preprint2012arXiv

GRB 120422A: A Low-luminosity Gamma-ray Burst Driven by Central Engine

GRB 120422A is a low-luminosity Gamma-ray burst (GRB) associated with a bright supernova, which distinguishes itself by its relatively short T90 ~ 5 s and an energetic X-ray tail. We analyze the Swift BAT and XRT data and discuss the physical implications. We show that the early steep decline in the X-ray light curve can be interpreted as the curvature tail of a late emission episode around 58-86 s, with a curved instantaneous spectrum at the end of the emission episode. Together with the main activity in the first ~ 20 s and the weak emission from 40 s to 60 s, the prompt emission is variable, which points towards a central engine origin, in contrast to the shock breakout origin as invoked to interpret some other nearby low-luminosity supernova GRBs. The curvature effect interpretation and interpreting the early shallow decay as the coasting external forward shock emission in a wind medium both give a constraint on the bulk Lorentz factor Γto be around several. Comparing the properties of GRB 120422A and other supernova GRBs, we found that the main criterion to distinguish engine-driven GRBs from the shock breakout GRBs is the time-averaged luminosity, with a separation line of about ~ 10^48 erg s^-1.

preprint2012arXiv

Statistical interpretation of the spatial distribution of current 130 GeV gamma-ray line signal within the dark matter annihilation scenario

Recently, several groups identified a tentative $γ$-ray line signal with energy $\sim 130$ GeV in the central Galaxy from the Fermi-LAT data. %The morphology study shows that the signal is consistent with dark matter %annihilation, but with an offset $\sim 220$ pc ($1.5^{\circ}$) of the %center from the Galactic center Sgr A$^{\star}$, Such a $γ-$ray line can be interpreted as the signal of dark matter annihilation. However, the offset $\sim 220$ pc ($1.5^{\circ}$) of the center of the most prominent signal region from the Galactic center Sgr A$^{\star}$ has been thought to challenge the dark matter annihilation interpretation. Considering the fact that such a 130 GeV $γ$-ray line signal consists of only $\sim14$ photons, we suggest that the "imperfect" consistency of these photons with the expected dark matter distribution is due to the limited statistics. The offset will be smaller as more signal photons have been collected in the near future. Our Monte Carlo simulation supports the above speculation.

preprint2012arXiv

The photospheric radiation model for the prompt emission of Gamma-ray Bursts: Interpreting four observed correlations

We show that the experiential $E_{\rm p}-L$, $Γ-L$, $E_{\rm p}-Γ$ and ${\rm \barη_γ}-E_{\rm p}$ correlations (where $L$ is the time-averaged luminosity of the prompt emission, $E_{\rm p}$ is the spectral peak energy, $Γ$ is the bulk Lorentz factor and $\barη_γ$ is the emission efficiency of Gamma-ray bursts) are well consistent with the relations between the resembling parameters predicted in the photospheric radiation model of the prompt emission of Gamma-ray bursts. The time-resolved thermal radiation of GRB 090902B does follow the $E_{\rm p}-L$ and $Γ-L$ correlations. A reliable interpretation of the four correlations in alternative models is still lacking. These facts may point towards a photospheric origin of prompt emission of some Gamma-ray Bursts.

preprint2011arXiv

A Comprehensive Analysis of Fermi Gamma-Ray Burst Data. I. Spectral Components and Their Possible Physical Origins of LAT/GBM GRBs

We present a systematic analysis of the spectral and temporal properties of 17 GRBs co-detected by GBM and LAT on board the Fermi satellite by May 2010. We performed a time-resolved spectral analysis of all the bursts with the finest temporal resolution allowed by statistics, in order to avoid temporal smearing of different spectral components. We found that the time-resolved spectra of 14 out of 17 GRBs are best modeled with the Band function over the entire Fermi spectral range, which may suggest a common origin for emissions detected by LAT and GBM. GRB 090902B and GRB 090510 require the superposition between an MeV component and an extra power law component, with the former having a sharp cutoff above E_p. For GRB 090902B, this MeV component becomes progressively narrower as the time bin gets smaller, and can be fit with a Planck function as the time bin becomes small enough. In general, we speculate that phenomenologically there may be three elemental spectral components : (I) a Band-function component (e.g. in GRB 080916C) that extends in a wide energy range and does not narrow with reducing time bins, which may be of the non-thermal origin; (II) a quasi-thermal component (e.g. in GRB 090902B) with the spectra progressively narrowing with reducing time bins; and (III) another non-thermal power law component extending to high energies. The spectra of different bursts may be decomposed into one or more of these elemental components. We compare this sample with the BATSE sample and investigate some correlations among spectral parameters. We discuss the physical implications of the data analysis results for GRB prompt emission, including jet compositions (matter-dominated vs. Poynting-flux-dominated outflow), emission sites (internal shock, external shock or photosphere), as well as radiation mechanisms (synchrotron, synchrotron self-Compton, or thermal Compton upscattering).

preprint2011arXiv

Dark Matter Mini-halo around the Compact Objects: the Formation, Evolution and Possible Contribution to the Cosmic Ray Electrons/Positrons

Dark matter particles may be captured by a star and then thermalized in the star's core. At the end of its life a massive star collapses suddenly and a compact object is formed. The dark matter particles redistribute accordingly. In the inelastic dark matter model, an extended dense dark matter mini-halo surrounding the neutron star may be formed. Such mini-halos may be common in the Galaxy. The electron/positron flux resulting in the annihilation of dark matter particles, however, is unable to give rise to observable signal unless a nascent mini-halo is within a distance \sim a few 0.1 pc from the Earth.

preprint2011arXiv

GRB 081029: A Gamma-Ray Burst with a Multi-Component Afterglow

We present an analysis of the unusual optical light curve of the gamma-ray burst GRB 081029, a long-soft burst with a redshift of z = 3.8479. We combine X-ray and optical observations from the Swift X-Ray Telescope and the Swift UltraViolet/Optical Telescope with ground-based optical and infrared data obtained using the REM, ROTSE, and CTIO 1.3-m telescopes to construct a detailed data set extending from 86 s to approximately 100,000 s after the BAT trigger. Our data cover a wide energy range, from 10 keV to 0.77 eV (1.24 to 16,000 Angstrom). The X-ray afterglow shows a shallow initial decay followed by a rapid decay starting at about 18,000 s. The optical and infrared afterglow, however, shows an uncharacteristic rise at about 3000 s that does not correspond to any feature in the X-ray light curve. Our data are not consistent with synchrotron radiation from a jet interacting with an external medium, a two-component jet, or continuous energy injection from the central engine. We find that the optical light curves can be broadly explained by a collision between two ejecta shells within a two-component jet. A growing number of gamma-ray burst afterglows are consistent with complex jets, which suggests that some (or all) gamma-ray burst jets are complex and will require detailed modelling to fully understand them.

preprint2011arXiv

Implications of Understanding Short Gamma-Ray Bursts Detected by {\it Swift}

In an effort to understand the puzzle of classifying gamma-ray bursts (GRBs), we perform a systematic study of {\it Swift} GRBs and investigate several short GRB issues. Though short GRBs have a short ($\lesssim2$ s) prompt duration as monitored by the Burst Alert Telescope, the composite light curves including both the prompt and afterglow emission suggest that most of the short GRBs have a similar radiative feature to long GRBs. Further, some well-studied short GRBs might also have an intrinsically long prompt duration, which renders them as a type of short GRB imposters. Genuine short GRBs detected by {\it Swift} might be rare that discriminating the observed short GRBs is, not surprisingly, troublesome. In particular, the observational biases in the host identification and redshift measurement of GRBs should be taken with great caution. The redshift distribution which has been found to be different for long and short GRBs might have been strongly affected by the measurement methods. We find that the redshifts measured from the presumed host galaxies of long and short GRBs appear to have a similar distribution.

preprint2011arXiv

Short Gamma-ray Bursts: the mass of the accretion disk and the initial radius of the outflow

In this work we estimate the accretion-disk mass in the specific scenario of binary-neutron-star-merger with current observational data. Assuming that the outflows of short Gamma-ray Bursts (GRBs) are driven via neutrino-antineutrino annihilation we estimate the disk mass of about half of short bursts in the sample to be $\sim 0.01-0.1 M_\odot$, in agreement with that obtained in the numerical simulations. Massive disks ($\sim {\rm several} 0.1 M_\odot$) found in some other short GRBs may point to the more efficient magnetic process of energy extraction or the neutron star and black hole binary progenitor. Our results suggest that some short bursts may be really due to the coalescence of double neutron stars and are promising gravitational wave radiation sources. For future short GRBs with simultaneous gravitational-wave detections, the disk mass may be reliably inferred and the validity of our approach will be tested. We also propose a method to constrain the initial radius of a baryonic outflow where it is launched ($R_0$) without the need of identifying an ideal thermal spectrum. We then apply it to GRB 090510 and get that $R_0 \lesssim 6.5\times 10^{6}(Γ_{\rm ph}/2000)^{-4}$ cm, suggesting that the central engine is a black hole with a mass $< 22 M_\odot(Γ_{\rm ph}/2000)^{-4}$, where $Γ_{\rm ph}$ is the bulk Lorentz factor of the outflow at the photospheric radius.

preprint2011arXiv

Sw 1644+57/GRB 110328A: the physical origin and the composition of the relativistic outflow

Sw 1644+57/GRB 110328A is a remarkable cosmological X-ray outburst detected by the {\it Swift} satellite. Its early-time ($t\lesssim 0.1$ days since the trigger) X-ray emission resembles some gamma-ray bursts (GRBs), e.g., GRB 090417B. But the late-time flaring X-ray plateau lasting $> 40$ days renders it unique. We examine the possibilities that the outburst is a super-long GRB powered either by the fallback accretion onto a nascent black hole or by a millisecond pulsar, and find out that these two scenarios can address some but not all of the main observational features. We then focus on the model of tidal disruption of a (giant) star by a massive black hole. The mass of the tidal-disrupted star is estimated to be $\gtrsim$ a few solar masses. A simple/straightforward argument for a magnetic origin of the relativistic outflow is presented.

preprint2011arXiv

Very Old Isolated Compact Objects as Dark Matter Probes

Very old isolated neutron stars and white dwarfs have been suggested to be probes of dark matter. To play such a role, two requests should be fulfilled, i.e., the annihilation luminosity of the captured dark matter particles is above the thermal emission of the cooling compact objects (request-I) and also dominate over the energy output due to the accretion of normal matter onto the compact objects (request-II). Request-I calls for very dense dark matter medium and the critical density sensitively depends on the residual surface temperature of the very old compact objects. The accretion of interstellar/intracluster medium onto the compact objects is governed by the physical properties of the medium and by the magnetization and rotation of the stars and may outshine the signal of dark matter annihilation. Only in a few specific scenarios both requests are satisfied and the compact objects are dark matter burners. The observational challenges are discussed and a possible way to identify the dark matter burners is outlined.

preprint2010arXiv

$e^\pm$ Excesses in the Cosmic Ray Spectrum and Possible Interpretations

The data collected by ATIC, PPB-BETS, FERMI-LAT and HESS all indicate that there is an electron/positron excess in the cosmic ray energy spectrum above $\sim$ 100 GeV, although different instrumental teams do not agree on the detailed spectral shape. PAMELA also reported a clear excess feature of the positron fraction above several GeV, but no excess in anti-protons. Here we review the observational status and theoretical models of this interesting observational feature. We pay special attention to various physical interpretations proposed in the literature, including modified supernova remnant models for the $e^\pm$ background, new astrophysical sources, and new physics (the dark matter models). We suggest that although most models can make a case to interpret the data, with the current observational constraints the dark matter interpretations, especially those invoking annihilation, require much more exotic assumptions than some astrophysical interpretations. Future observations may present some ``smoking-gun'' observational tests to differentiate among different models and to identify the correct interpretation to the phenomenon.

preprint2010arXiv

A Revised Limit of the Lorentz Factors of GRBs with Two Emitting Regions

Fermi observations of GeV emission from GRBs have suggested that the Lorentz factor of some GRBs is around a thousand or even higher. At the same time the same Fermi observations have shown an extended GeV emission indicating that this higher energy emission might be a part of the afterglow and it does not come from the same region as the lower energy prompt emission. If this interpretation is correct than we should reconsider the opacity limits on the Loretnz factor of the emitting regions which are based on a one-zone model. We describe here a two-zone model in which the GeV photons are emitted in a larger radius than the MeV photons and we calculate the optical depth for pair creation of a GeV photon passing the lower energy photons shell. We find that, as expected, the new two-zone limits on the Lorentz factor are significantly lower. The corresponding limits for the Fermi bursts are lower by a factor of five compared to the one-zone model and it is possible that both the MeV and GeV regions have relatively modest Lorentz factors (~200 - 400).

preprint2010arXiv

Interpreting the high energy emission of Fermi GRBs

The high energy emission from Gamma-ray Bursts has some interesting features, including the absence of the GeV excess in the prompt spectrum, the delayed onset of the GeV emission, and the longer duration of the GeV emission than the prompt soft gamma-ray emission. We suggest that the non-detection of a GeV excess in most GRB spectrum may favor the magnetic fireball model and the early prompt emission may be dominated by the photosphere radiation of the breakout material and is thus very soft. The synchrotron radiation in GeV band can be the dominant component of the high energy afterglow emission, as speculated in GRB 080319B and then confirmed in GRB 080916C and GRB 090510. A simple estimate of the thermal radiation of the breakout material has been presented.

preprint2010arXiv

The bulk Lorentz factor of outflow powering X-ray flare in GRB afterglow

We develop two methods to estimate the bulk Lorentz factor of X-ray flare outflow. In the first method the outflow is assumed to be baryonic and is accelerated by the thermal pressure, for which the final bulk Lorentz factor is limited by the outflow luminosity as well as the initial radius of the outflow getting accelerated. Such a method may be applied to a considerable fraction of flares. The second method, based on the curvature effect interpretation of the quick decline of the flare, can give a tightly constrained estimate of the bulk Lorentz factor but can only be applied to a few giant flares. The results obtained in these two different ways are consistent with each other. The obtained bulk Lorentz factor (or just upper limit) of the X-ray flare outflows, ranging from ten to a few hundred, is generally smaller than that of the Gamma-ray Burst (GRB) outflows.

preprint2010arXiv

Underlying global features of the x-ray light curves of {\it swift} gamma-ray bursts

With its rapid response, {\it Swift} has revealed plenty of unexpected properties of gamma-ray bursts (GRBs). With an abundance of observations, our current understanding is only limited by complexity of early X-ray light curves. In this work, based on the public {\it Swift} data of 150 well-monitored GRBs with measured redshifts, we find some interesting global features in the rest-frame X-ray light curves. The distinct spectral behaviors between the prompt emission and the afterglow emission implies dissimilar radiation scenarios. Interestingly, an unforeseen plateau is exhibited in the prompt X-ray light curves despite the presence of complex spikes, which might indicate the presence of a steady central engine. In particular, the seemingly continuous evolution with a single power law from the prompt to the afterglow of most GRBs might place strong constraints on the theoretical models.

preprint2010arXiv

XRF 100316D/SN 2010bh: clue to the diverse origin of nearby supernova-associated GRBs

X-ray Flash (XRF) 100316D, a nearby super-long under-luminous burst with a peak energy E_p \sim 20 keV, was detected by Swift and was found to be associated with an energetic supernova SN 2010bh. Both the spectral and the temporal behavior of this burst are rather similar to that of XRF 060218, except that the latter was associated with a "less energetic" SN 2006aj and had a prominent soft thermal emission component in the spectrum. We analyze the spectral and temporal properties of this burst, and interpret the prompt gamma-ray emission and the early X-ray plateau emission as synchrotron emission from a dissipating Poynting-flux-dominated outflow, probably powered by a magnetar with a spin period of $P \sim 10$ ms and the polar cap magnetic field $B_{\rm p} \sim 3\times 10^{15}$ G. The energetic supernova SN 2010bh associated with this burst is, however, difficult to interpret within the slow magnetar model, which implies that the nascent magnetar may spin much faster with an initial rotation period $\sim 1$ ms, and thus suggests a delay between the core collapse and the emergence of the relativistic magnetar wind from the star. The diverse behaviors of low-luminosity GRBs and their associated SNe may be understood within a unified picture that invokes different initial powers of the central engine and different delay times between the core collapse and the emergence of the relativistic jet from the star.

preprint2009arXiv

GRB 080916C and GRB 090510: the high energy emission and the afterglow

We constrain the physical composition of the outflows of GRBs 080916C and 090510 with the prompt emission data and find that the former is likely magnetic while the latter may be baryonic. The X-ray and optical afterglow emission of both GRBs can be reasonably fitted using the standard external shock model but the density profiles of the circum-burst medium are different. We also propose a simple method to estimate the number of the seed photons suppose the GeV afterglow photons are due to the inverse Compton radiation of external forward shock electrons. The seed photons needed in the modeling are too many to be realistic for both events. The synchrotron radiation of the forward shock seems able to account for the GeV afterglow data.

preprint2009arXiv

The spectrum of Gamma-ray Burst: a clue

In this work we numerically calculate the thermal radiation efficiency of the baryonic outflow. The possible outflow acceleration in the transparent stage, which lowers thermal radiation efficiency, has been taken into account. In the standard internal shock model for the prompt emission, the fast shells should move with a typical Lorentz factor $\gtrsim 5 Γ_{\rm i}$ otherwise the GRB efficiency will be in disagreement with the observations, where $Γ_{\rm i}$ is the bulk Lorentz factor of the shocked/emitting region. The photosphere radius of these fast shells is small and the thermal radiation is too strong to be effectively outshone by the internal shock emission. This is particularly the case for some extremely bright events having $Γ_{\rm i} \sim 10^{3}$, like GRBs 080319B and 080916C. The absence of a distinct thermal component in the spectrum of most GRBs challenges the standard internal shock model and may suggest a non-baryonic (magnetic) outflow component. Though the magnetic outflow model seems favored by more and more data, it can hardly reproduce the typical GRB spectrum. In the photosphere-gradual magnetic dissipation scenario, the spectrum cuts off at $\sim 1$ GeV, too low to account for the observations of GRBs 080916C. In the sudden magnetic energy dissipation model, the low energy spectrum is expected to be $F_ν\propto ν^{-1/2}$, too soft to be consistent with the data $F_ν\propto ν^{0}$. We speculate that the low energy spectrum puzzle could be unveiled by the mechanism that particles, in the magnetic dissipation process, are repeatedly accelerated.

preprint2007arXiv

The synchrotron-self-Compton radiation accompanying shallow decaying X-ray afterglow: the case of GRB 940217

High energy emission (> tens MeV) of Gamma-Ray Bursts (GRBs) provides an important clue to understand the physical processes involved in GRBs, which may be correlated with the GRB early afterglow. A shallow decline phase has been well detected in about half {\it Swift} Gamma-ray Burst X-ray afterglows. The widely considered interpretation involves a significant energy injection and possibly time-evolving shock parameter(s). This work we calculate the synchrotron-self-Compton (SSC) radiation of such an external forward shock and show that it could explain the well-known long term high energy (i.e., tens MeV to GeV) afterglow of GRB 940217. We propose that the cooperation of Swift and GLAST will help to reveal the nature of GRBs.

preprint2006arXiv

The interpretation of the Swift GRB X-ray afterglows

We discuss the current interpretations of the {\it Swift} GRB X-ray afterglows, mainly focusing on the sharp decline at the prompt tail emission, and the shallow decay afterward, which is then followed by the conventional pre-\emph{Swift} decay behavior, and the possible X-ray flares during the latter two stages. We emphasize the role of the central engine in interpreting the GRB afterglows.