Source author record

He Gao

He Gao 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

63works
11topics
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

63 published item(s)

preprint2026arXiv

Multimessenger Prospects for Low-Luminosity Gamma-Ray Bursts: Joint Neutrino and X-Ray Observations

Low--luminosity gamma-ray bursts (LLGRBs) are promising candidates for high-energy neutrinos, yet no coincident neutrino events have been detected so far. Recent advances in X-ray time-domain astronomy, together with the development of next-generation neutrino telescopes, open new opportunities for joint X-ray and neutrino observations of these transients. We calculate the jet dynamical evolution and the associated neutrino production for both non-magnetized and magnetized outflows. For individual events, joint X-ray and neutrino detection is generally limited to nearby LLGRBs or sources with high luminosities. Thus, we consider a next-generation neutrino telescope with an effective area enhanced by a factor of $\sim30$ relative to IceCube. In the non-magnetized scenario, joint detection of individual events is enabled for sources with typical isotropic luminosities of $L_{\mathrm{iso}}\sim10^{47}\,\mathrm{erg\,s^{-1}}$ out to luminosity distances of $D_L\sim1.6\times10^{2}\,\mathrm{Mpc}$, corresponding to an expected detection rate of order $1$ per year. In contrast, for the magnetized scenario at the same luminosity, the accessible distance is significantly reduced, with joint observations confined to sources within $D_L\sim6.5\times10^{1}\,\mathrm{Mpc}$ and an expected detection rate of order $0.5$ per year. For stacked samples of $\sim100$ magnetized LLGRBs, stacking substantially enlarges the accessible distance range, enabling joint observations for sources with representative luminosities of $L_{\mathrm{iso}}\sim1\times10^{47}\,\mathrm{erg\,s^{-1}}$ out to $D_L\lesssim7.0\times10^{2}\,\mathrm{Mpc}$ and corresponding to an expected detection rate of order $2$ per year. These results demonstrate that joint X-ray and next-generation neutrino observations enable a practical multimessenger probe of LLGRBs.

preprint2025arXiv

Constraining the Properties of GRB Accreting Magnetar with $R/I$ Evolutionary Effects Using \emph{Swift}/XRT Data

A newly born millisecond magnetar has been proposed as one possible central engine of some long gamma-ray bursts (LGRBs) with X-ray plateau. In this work, we used a universal correlation between initial spin period ($P_0$) and surface magnetic field ($B_p$) of newborn magnetar based on an LGRB sample in \cite{Lan2025} to explore the propeller properties of accreting magnetars with $R/I$ evolutionary effects. We found that $B_p-P_0$ relation is approximately consistent with $B_p\propto P_{\rm eq}^{7/6}$. Here $P_{\rm eq}$ is equilibrium spin period in magnetic propeller model. The $B_p-P_0$ relation indicates that $P_0$ may not be true initial spin period of newborn magnetar but had reached an equilibrium spin period via fallback accretion in propeller model. The magnetar accretion rate in our LGRBs is in range of $\dot{M}\sim10^{-5}-10^{-2} M_{\odot} \rm s^{-1}$ by incorporating $R/I$ evolutionary effects and using the transition relation between gravitational mass $M_g$ and baryonic mass $M_b$ in different equations of state. Such accretion rates ensure that the accreting magnetars in our sample survive until reaching the equilibrium spin period, and the accretion rate is one order of magnitude lower compared to the statistical results in \cite{Stratta2018} and \cite{Linweili2020}, which used constant $R/I/M_g$ scenario. We suggested that adopting a constant $R/I/M_g$ scenario for modeling propeller regime in accreting magnetar results in a higher mass accretion rate, which may impair our understanding of the physical nature and its surroundings of accreting magnetar, and low-metallicity progenitors can provide enough material to satisfy the accretion requirements of newborn accreting magnetar in LGRBs.

preprint2025arXiv

FRB cosmology with the RM-PRS Luminosity Correlation

Fast Radio Bursts (FRBs) have emerged as a powerful tool for cosmological studies, particularly through the dispersion measure-redshift ($\mathrm{DM}-z$) relation. This work proposes a novel calibration method for FRBs using the Yang-Li-Zhang (YLZ) empirical relation, which links the rotation measure (RM) of FRBs to the luminosity of their associated persistent radio sources (PRS). We demonstrate that this approach provides independent constraints on cosmological parameters, bypassing limitations inherent to traditional $\mathrm{DM}-z$ method. Utilizing the current sample of four YLZ-calibrated FRBs, we derive a Hubble constant measurement of $H_0 = 86.18_{-14.99}^{+18.03}\ \mathrm{km\ s^{-1}\ Mpc^{-1}}$ (68\% CL). Monte Carlo simulations indicate that a future catalog of 400 FRB-PSR systems could reduce the relative uncertainty of $H_0$ to 4.5\%. Combining YLZ-calibrated FRBs with $\mathrm{DM}-z$ sample reveals critical synergies: joint analysis of equalized samples ($N=100$ for both methods) reduces the relative uncertainty of $H_0$ to 2.9\%, mainly because the incorporation of PRS observations substantially mitigates the degeneracy between the parameters such as IGM baryon mass fraction ($f_{\rm IGM}$) and other cosmological parameters inherent to the $\mathrm{DM}-z$ relation.

preprint2022arXiv

A Search for Millilensing Gamma-Ray Bursts in the Observations of Fermi GBM

Millilensing of Gamma-Ray Bursts (GRBs) is expected to manifest as multiple emission episodes in a single triggered GRB with similar light-curve patterns and similar spectrum properties. Identifying such lensed GRBs could help improve constraints on the abundance of compact dark matter. Here we present a systemic search for millilensing among 3000 GRBs observed by the \textit{Fermi} GBM up to 2021 April. Eventually we find 4 interesting candidates by performing auto-correlation test, hardness test, and time-integrated/resolved spectrum test to the whole sample. GRB 081126A and GRB 090717A are ranked as the first class candidate based on their excellent performance both in temporal and spectrum analysis. GRB 081122A and GRB 110517B are ranked as the second class candidates (suspected candidates), mainly because their two emission episodes show clear deviations in part of the time-resolved spectrum or in the time-integrated spectrum. Considering a point mass model for the gravitational lens, our results suggest that the density parameter of lens objects with mass $M_{\rm L}\sim10^{6} M_{\odot}$ is larger than $1.5\times10^{-3}$.

preprint2022arXiv

Constraints on the abundance of primordial black holes with different mass distributions from lensing of fast radio bursts

Primordial black holes (PBHs) has been considered to form a part of dark matter for a long time but the possibility has been poorly constrained over a wide mass range, including the stellar mass range ($1-100~M_{\odot}$). However, due to the discovery of merger events of black hole binaries by LIGO-Virgo gravitational wave observatories, the interest for PBHs in the stellar mass window has been aroused again. Fast radio bursts (FRBs) are bright radio transients with millisecond duration and very high all-sky occurrence rate. Lensing effect of these bursts has been proposed as one of the optimal probes for constraining the abundance of PBHs in the stellar mass range. In this paper, we first investigate constraints on the abundance of PBHs from the latest $593$ FRB observations for both the monochromatic mass distribution and three other popular extended mass distributions related to different formation mechanisms of PBHs. It is found that constraints from currently public FRB observations are relatively weaker than those from existing gravitational wave detections. Furthermore, we forecast constraining power of future FRB observations on the abundance of PBHs with different mass distributions of PBHs and different redshift distributions of FRBs taken into account. Finally, We find that constraints of parameter space on extended mass distributions from $\sim10^5$ FRBs with $\overline{Δt}\leq1 ~\rm ms$ would be comparable with what can be constrained from gravitational wave events. It is foreseen that upcoming complementary multi-messenger observations will yield considerable constraints on the possibilities of PBHs in this intriguing mass window.

preprint2022arXiv

Gamma-ray Burst Prompt Emission Spectrum and $E_p$ Evolution Patterns in the ICMART Model

In this paper, we simulate the gamma-ray bursts (GRBs) prompt emission light curve, spectrum and $E_p$ evolution patterns within the framework of the Internal-Collision-induced MAgnetic Reconnection and Turbulence (ICMART) model. We show that this model can produce a Band shape spectrum, whose parameters ($E_p$, $α$, $β$) could distribute in the typical distribution from GRB observations, as long as the magnetic field and the electron acceleration process in the emission region are under appropriate conditions. On the other hand, we show that for one ICMART event, $E_p$ evolution is always a hard-to-soft pattern. However, a GRB light curve is usually composed of multiple ICMART events that are fundamentally driven by the erratic GRB central engine activity. In this case, we find that if one individual broad pulse in the GRB light curve is composed of multiple ICMART events, the overall $E_p$ evolution could be disguised as the intense-tracking pattern. Therefore, mixed $E_p$ evolution patterns can coexist in the same burst, with a variety of combined patterns. Our results support the ICMART model to be a competitive model to explain the main properties of GRB prompt emission. The possible challenges faced by the ICMART model are also discussed in details.

preprint2022arXiv

GRB 211211A: a Prolonged Central Engine under a Strong Magnetic Field Environment

Recently, a kilonova-associated gamma-ray burst (GRB 211211A), whose light curve consists of a precursor ($\sim$0.2 s), a hard spiky emission ($\sim$10 s), and a soft long extended emission ($\sim$40 s), has attracted great attention. Kilonova association could prove its merger origin, while the detection of the precursor can be used to infer at least one highly magnetized neutron star (NS) being involved in the merger. In this case, a strong magnetic flux $Φ$ is expected to surround the central engine of GRB 211211A. Here we suggest that when $Φ$ is large enough, the accretion flow could be halted far from the innermost stable radius, which will significantly prolong the lifetime of the accretion process, and so the GRB duration. For example, we show that as long as the central black hole (BH) is surrounded by a strong magnetic flux $Φ\sim 10^{29}\rm cm^{2}G$, an accretion flow with $\dot{M}_{\rm ini} \simeq 0.1 M_\odot s^{-1}$ could be halted at 40 times the gravitational radius and be slowly transferred into the black hole on the order of $\sim$10 s, which naturally explains the duration of hard spiky emission. After most of the disk mass has been accreted onto the BH, the inflow rate will be reduced, so a long and soft extended emission is expected when a new balance between the magnetic field and the accretion current is reconstructed at a larger radius. Our results further support that the special behavior of GRB 211211A is mainly due to the strong magnetic field of its progenitor stars.

preprint2022arXiv

Kilonova and Optical Afterglow from Binary Neutron Star Mergers. I. Luminosity Function and Color Evolution

In the first work of this series, we adopt a GW170817-like, viewing-angle-dependent kilonova model and the standard afterglow model with the lightcurve distribution based on the properties of cosmological short gamma-ray bursts afterglows to simulate the luminosity functions and color evolution of both kilonovae and optical afterglow emissions from binary neutron star (BNS) mergers. {We find that $\sim10\%$ afterglows are brighter than the associated kilonovae at the peak time, most of which are on-axis or nearly-on-axis. These kilonovae would be significantly polluted by the associated afterglow emission. Only at large viewing angles with $\sinθ_{\rm v}\gtrsim 0.20$, the EM signals of most BNS mergers would be kilonova-dominated and some off-axis afterglows may emerge at $\sim5-10$\,day after the mergers.} At brightness dimmer than $\sim23-24$\,mag, according to their luminosity functions, the number of afterglows is much larger than that of kilonovae. Since the search depth of the present survey projects is $<22$\,mag, the number of afterglow events detected via serendipitous observations would be much higher than that of kilonova events, consistent with the current observations. {For the foreseeable survey projects (e.g., Mephisto, WFST, LSST) whose search depths can reach $\gtrsim23-24$\,mag, the detection rate of kilonovae could have the same order of magnitude as that of afterglows.} We also find that it may be difficult to use the fading rate in a single band to directly identify kilonovae and afterglows among various fast-evolving transients by serendipitous surveys. However, the color evolution between optical and infrared bands can identify them, since their color evolution patterns are unique compared with those of other fast-evolving transients.

preprint2022arXiv

Population Properties of Gravitational-Wave Neutron Star--Black Hole Mergers

Over the course of the third observing run of LIGO-Virgo-KAGRA Collaboration, several gravitational-wave (GW) neutron star--black hole (NSBH) candidates have been announced. By assuming these candidates are real signals and of astrophysical origins, we analyze the population properties of the mass and spin distributions for GW NSBH mergers. We find that the primary BH mass distribution of NSBH systems, whose shape is consistent with that inferred from the GW binary BH (BBH) primaries, can be well described as a power-law with an index of $α= 4.8^{+4.5}_{-2.8}$ plus a high-mass Gaussian component peaking at $\sim33^{+14}_{-9}\,M_\odot$. The NS mass spectrum could be shaped as a near flat distribution between $\sim1.0-2.1\,M_\odot$. The constrained NS maximum mass agrees with that inferred from NSs in our Galaxy. If GW190814 and GW200210 are NSBH mergers, the posterior results of the NS maximum mass would be always larger than $\sim2.5\,M_\odot$ and significantly deviate from that inferred in the Galactic NSs. The effective inspiral spin and effective precession spin of GW NSBH mergers are measured to potentially have near-zero distributions. The negligible spins for GW NSBH mergers imply that most events in the universe should be plunging events, which supports the standard isolated formation channel of NSBH binaries. More NSBH mergers to be discovered in the fourth observing run would help to more precisely model the population properties of cosmological NSBH mergers.

preprint2022arXiv

Prospects of strongly lensed fast radio bursts: Simultaneous measurement of post-Newtonian parameter and Hubble constant

Strong gravitational lensing effect is a powerful tool to probe cosmological models and gravity theories. Recently, the time-delay cosmography from strong lensing and the stellar kinematics of the deflector, which encode the Hubble constant and the post-Newtonian parameter via two distance ratios reflecting the lensing mass and dynamical mass respectively, have been proposed to investigate these two parameters simultaneously. Among strong lensing systems with different sources, strongly lensed fast radio bursts (FRBs) have been proposed as precision probes of the universe since the time delay $\sim$ 10 days between images could be measured extremely precisely because of their short duration of a few milliseconds. In this work, we investigate the ability of strongly lensed FRBs on simultaneously estimating these two parameters via simulations. Take the expected FRB detection rate of upcoming facilities and lensing probability into consideration, it is likely to accumulate 10 lensed FRBs in several years and we find that $H_0$ could be determined to a $\sim1.5\%$ precision and $γ_{\rm PPN}$ could be constrained to a $\sim8.7\%$ precision simultaneously from them. These simultaneous estimations will be helpful for properly reflecting the possible correlation between these two fundamental parameters.

preprint2022arXiv

Quasi-periodic oscillations of the X-ray burst from the magnetar SGR J1935+2154 and associated with the fast radio burst FRB 200428

The origin(s) and mechanism(s) of fast radio bursts (FRBs), which are short radio pulses from cosmological distances, have remained a major puzzle since their discovery. We report a strong Quasi-Periodic Oscillation(QPO) of 40 Hz in the X-ray burst from the magnetar SGR J1935+2154 and associated with FRB 200428, significantly detected with the Hard X-ray Modulation Telescope (Insight-HXMT) and also hinted by the Konus-Wind data. QPOs from magnetar bursts have only been rarely detected; our 3.4 sigma (p-value is 2.9e-4) detection of the QPO reported here reveals the strongest QPO signal observed from magnetars (except in some very rare giant flares), making this X-ray burst unique among magnetar bursts. The two X-ray spikes coinciding with the two FRB pulses are also among the peaks of the QPO. Our results suggest that at least some FRBs are related to strong oscillation processes of neutron stars. We also show that we may overestimate the significance of the QPO signal and underestimate the errors of QPO parameters if QPO exists only in a fraction of the time series of a X-ray burst which we use to calculate the Leahy-normalized periodogram.

preprint2022arXiv

Search for lensing signatures from the latest fast radio burst observations and constraints on the abundance of primordial black holes

The possibility that primordial black holes (PBHs) form a part of dark matter has been considered for a long time but poorly constrained over a wide mass range. Fast radio bursts (FRBs) are bright radio transients with millisecond duration. Lensing effect of them has been proposed as one of the cleanest probes for constraining the presence of PBHs in the stellar mass window. In this paper, we first apply the normalised cross-correlation algorithm to search and identify candidates of lensed FRBs in the latest public FRB observations, i.e. $593$ FRBs which mainly consist of the first Canadian Hydrogen Intensity Mapping Experiment FRB catalog, and then derive constraints on the abundance of PBHs from the null search result of lensing signature. For a monochromatic mass distribution, the fraction of dark matter made up of PBHs could be constrained to $\leq87\%$ for $\geq500~M_{\odot}$ at 95\% confidence level by assuming flux ratio thresholds dependent signal-to-noise ratio for each FRB and that apparently one-off events are intrinsic single bursts. This result would be improved by a three times factor when a conventional constant flux ratio threshold is considered. Moreover, we derive constraints on PBHs with a log-normal mass function naturally predicted by some popular inflation models and often investigated with gravitational wave detections. We find that, in this mass distribution scenario, the constraint from currently public FRB observations is relatively weaker than the one from gravitational wave detections. It is foreseen that upcoming complementary multi-messenger observations will yield considerable constraints on the possibilities of PBHs in this intriguing mass window.

preprint2022arXiv

The Number of Possible CETIs within Our Galaxy and the Communication Probability among These CETIs

As the only known intelligent civilization, human beings are always curious about the existence of other communicating extraterrestrial intelligent civilizations (CETIs). Based on the latest astrophysical information, we carry out Monte Carlo simulations to estimate the number of possible CETIs within our Galaxy and the communication probability among them. Two poorly known parameters have a great impact on the results. One is the probability of life appearing on terrestrial planets and eventually evolving a into CETI ($f_c$), and the other determines at what stage of their host star's evolution CETIs would be born ($F$). In order to ensure the completeness of the simulation, we consider a variety of combinations of $f_c$ and $F$. Our results indicate that for optimistic situations (e.g. $F=25\%$ and $f_c=0.1\%$), there could be $42777_{-369}^{+267}$ CETIs and they need to survive for $3_{-2}^{+17}$ yr ($2000_{-1400}^{+2000}$ yr) to achieve one-way communication (two-way communication). In this case, human beings need to survive $0.3_{-0.298}^{+0.6}$ Myr to receive one alien signal. For pessimistic situations (e.g. $F=75\%$ and $f_c=0.001\%$), only $111_{-17}^{+28}$ CETIs could be born and they need to survive for $0.8_{-0.796}^{+1.2}$ Myr ($0.9_{-0.88}^{+4.1}$ Myr) to achieve one-way communication (two-way communication). In this case, human beings need to survive $50_{-49.6}^{+250}$ Myr to receive one signal from other CETIs. Our results may quantitatively explain why we have not detected any alien signals so far. The uncertainty of the results has been discussed in detail and would be alleviated with the further improvement of our astronomical observation ability in the future.

preprint2021arXiv

Searching for gravitationally lensed gamma-ray bursts with their afterglows

Gamma-ray bursts (GRBs) at high redshifts are expected to be gravitationally lensed by objects of different mass scales. Besides a single recent claim, no lensed GRB has been detected so far by using the gamma-ray data only. In this paper, we suggest that the multi-band afterglow data might be an efficient way to search for lensed GRB events. Using the standard afterglow model we calculate the characteristics of the lensed afterglow lightcurves under the assumption of two popular analytic lens models: point mass and Singular Isothermal Sphere (SIS) model. In particular, when different lensed images cannot be resolved, their signals would be superimposed together with a given time delay. In this case, the X-ray afterglows are likely to contain several X-ray flares of similar width in linear scale and similar spectrum, and the optical afterglow lightcurve will show rebrightening signatures. Since the lightcurves from the image arriving later would be compressed and deformed in the logarithmic time scale, the larger time delay (i.e. the larger mass of the lens), the easier is to identify the lensing effect. We analyzed the archival data of optical afterglows and found one potential candidate of the lensed GRB (130831A) with time delay $\sim$500 s, however, observations of this event in gamma-ray and X-ray band seem not to support the lensing hypothesis. In the future, with the cooperation of the all-sky monitoring gamma-ray detectors and multi-band sky survey projects, our method proposed in this paper would be more efficient in searching for strongly lensed GRBs.

preprint2020arXiv

Constraining the long-lived magnetar remnants in short gamma-ray bursts from late-time radio observations

The joint detection of GW 170817 and GRB 170817A indicated that at least a fraction of short gamma ray bursts (SGRBs) originate from binary neutron star (BNS) mergers. One possible remnant of a BNS merger is a rapidly rotating, strongly magnetized neutron star, which has been discussed as one possible central engine for GRBs. For a rapidly rotating magnetar central engine, the deposition of the rotation energy into the ejecta launched from the merger could lead to bright radio emission. The brightness of radio emission years after an SGRB would provide an estimation of the kinetic energy of ejecta and hence, a possible constraint on the BNS merger product. We perform a more detailed calculation on the brightness of radio emission from the interaction between the merger ejecta and circumburst medium in the magnetar scenario, invoking several important physical processes such as generic hydrodynamics, relativistic effects, and the deep Newtonian phase. We use the model to constrain the allowed parameter space for 15 SGRBs that have late radio observations. Our results show that an injected energy $E_{\rm inj} \sim 10^{52}$ erg is allowed for all the cases, which suggests that the possibility of a supra-massive or hyper-massive neutron star remnant is not disfavored by the available radio data.

preprint2020arXiv

Constraints on compact dark matter with fast radio burst observations

Fast Radio Bursts (FRBs) are bright radio transients with millisecond duration at cosmological distances. Since compact dark matter/objects (COs) could act as lenses and cause split of this kind of very short duration signals, Mu$\rm{\tilde{n}}$oz et al. (2016) has proposed a novel method to probe COs with lensing of FRBs. In this Letter, we for the first time apply this method to real data and give constraints of the nature of COs with currently available FRB observations. We emphasize the information from dynamic spectra of FRBs is quite necessary for identifying any lensed signals and find no echoes in the existing data. The null search gives a constraint comparable to that from galactic wide binaries, though the methods of redshift inference from dispersion measure would impact a little. Furthermore, we make an improved forecast basing on the distributions of real data for the ongoing and upcoming telescopes. Finally, we discuss the situation where one or more lensed signals will be detected. In such a case, the parameter space of COs can be pinned down very well since the lens mass can be directly determined through the observed flux ratio and time delay between split images.

preprint2020arXiv

Cosmology-insensitive estimate of IGM baryon mass fraction from five localized fast radio bursts

Five fast radio bursts (FRBs), including three apparently non-repeating ones FRB 180924, FRB 181112, and FRB 190523, and two repeaters, FRB 121102 and FRB 180916.J0158+65, have already been localized so far. We apply a method developed recently by us (Li et al. 2019) to these five localized FRBs to give a cosmology-insensitive estimate of the fraction of baryon mass in the IGM, $f_{\rm IGM}$. Using the measured dispersion measure (DM) and luminosity distance $d_{\rm L}$ data (inferred from the FRB redshifts and $d_{\rm L}$ of type Ia supernovae at the same redshifts) of the five FRBs, we constrain the local $f_{\rm IGM} = 0.84^{+0.16}_{-0.22}$ with no evidence of redshift dependence. This cosmology-insensitive estimate of $f_{\rm IGM}$ from FRB observations is in excellent agreement with previous constraints using other probes. Moreover, using the three apparently non-repeating FRBs only we get a little looser but consistent result $f_{\rm IGM} = 0.74^{+0.24}_{-0.18}$. In these two cases, reasonable estimations for the host galaxy DM contribution (${\rm DM_{host}}$) can be achieved by modelling it as a function of star formation rate. The constraints on both $f_{\rm IGM}$ and ${\rm DM_{host}}$ are expected to be significantly improved with the rapid progress in localizing FRBs.

preprint2020arXiv

Discovery of oscillations above 200 keV in a black hole X-ray binary with Insight-HXMT

Low-frequency quasi-periodic oscillations (LFQPOs) are commonly found in black hole X-ray binaries, and their origin is still under debate. The properties of LFQPOs at high energies (above 30 keV) are closely related to the nature of the accretion flow in the innermost regions, and thus play a crucial role in critically testing various theoretical models. The Hard X-ray Modulation Telescope (Insight-HXMT) is capable of detecting emissions above 30 keV, and is therefore an ideal instrument to do so. Here we report the discovery of LFQPOs above 200 keV in the new black hole MAXI J1820+070 in the X-ray hard state, which allows us to understand the behaviours of LFQPOs at hundreds of kiloelectronvolts. The phase lag of the LFQPO is constant around zero below 30 keV, and becomes a soft lag (that is, the high-energy photons arrive first) above 30 keV. The soft lag gradually increases with energy and reaches ~0.9s in the 150-200 keV band. The detection at energies above 200 keV, the large soft lag and the energy-related behaviors of the LFQPO pose a great challenge for most currently existing models, but suggest that the LFQPO probably originates from the precession of a small-scale jet.

preprint2020arXiv

Giant X-ray and optical Bump in GRBs: evidence for fall-back accretion model

The successful operation of dedicated detectors has brought us valuable information for understanding the central engine and the progenitor of gamma-ray bursts (GRBs). For instance, the giant X-ray and optical bumps found in some long-duration GRBs (e.g. GRBs 121027A and 111209A) imply that some extended central engine activities, such as the late X-ray flares, are likely due to the fall-back of progenitor envelope materials. Here we systemically search for long GRBs that consist of a giant X-ray or optical bump from the Swift GRB sample, and eventually we find 19 new possible candidates. The fall-back accretion model could well interpret the X-ray and optical bump for all candidates within a reasonable parameter space. Six candidates showing simultaneous bump signatures in both X-ray and optical observations, which could be well fitted at the same time when scaling down the X-ray flux into optical by one order of magnitude, are consistent with the standard $F_ν\proptoν^{1/3}$ synchrotron spectrum. The typical fall-back radius is distributed around $10^{10}\rm-10^{12}$ cm, which is consistent with the typical radius of a Wolf-Rayet star. The peak fall-back accretion rate is in the range of $\sim 10^{-11}-10^{-4}M_{\odot} \ \text{s} ^{-1}$ at time $\sim10^{2}- 10^{5}~\rm s$, which is relatively easy to fulfill as long as the progenitor's metallicity is not too high. Combined with the sample we found, future studies of the mass supply rate for the progenitors with different mass, metallicity, and angular momentum distribution would help us to better constrain the progenitor properties of long GRBs.

preprint2020arXiv

Kilonova Emission From Black Hole-Neutron Star Mergers. I. Viewing-Angle-Dependent Lightcurves

In this paper, we present a numerical method to study the predicted lightcurves as a function of viewing angle. We extrapolate the fitting formulae for the mass and velocity of tidal dynamical ejecta across a wide range of mass ratio validated with 66 simulations and use them in the calculations of the kilonova lightcurves. The calculated peak luminosity of a BH-NS merger kilonova is typically about a few times $10^{41}\ {\rm erg\ s^{-1}}$, which is always $\lesssim4.5\times10^{41}\ {\rm erg\ s^{-1}}$. This corresponds to the AB absolute magnitudes fainter than $\sim -15\ {\rm mag}$ in optical and $\sim -16\ {\rm mag}$ in infrared. Since the projected photosphere area of the dynamical ejecta is much larger than that of the disk wind outflows, the dynamical ejecta usually contribute to the majority of the kilonova emission from BH-NS mergers. The fitted blackbody temperature and the shape of the observed multi-band lightcurves are insensitive to the line of sight. The peak time of the observed multi-band lightcurves, affected by the light propagation effect, is related to the relative motion direction between the dynamical ejecta and the observer. The observed luminosity varies with the projected photosphere area determined by the viewing angles. However, the predicted peak luminosity only varies by a factor of $\sim (2 - 3)$ (or by $\sim1\ {\rm mag}$) for different viewing angles. When the short-duration gamma-ray burst afterglow is taken into account, for an on-axis geometry, the kilonova emission is usually outshone by the afterglow emission and can be only observed in the redder bands, especially in the $K$-band at late times. Compared with GW170817/AT2017gfo, the BH-NS merger kilonovae are optically dim but possibly infrared bright. At the same epoch after the merger, the blackbody fitting temperature of the BH-NS merger kilonovae is lower than that of GW170817/AT2017gfo.

preprint2020arXiv

Lorentz Invariance Violation Limits from the Spectral Lag Transition of GRB~190114C

The spectral lags of gamma-ray bursts (GRBs) have been viewed as the most promising probes of the possible violations of Lorentz invariance (LIV). However, these constraints usually depend on the assumption of the unknown intrinsic time lag in different energy bands and the use of a single highest-energy photon. A new approach to test the LIV effects has been proposed by directly fitting the spectral lag behavior of a GRB with a well-defined transition from positive lags to negative lags. This method simultaneously provides a reasonable formulation of the intrinsic time lag and robust lower limits on the quantum-gravity energy scales ($E_{\rm QG}$). In this work, we perform a global fitting to the spectral lag data of GRB~190114C by considering the possible LIV effects based on a Bayesian approach. We then derive limits on $E_{\rm QG}$ and the coefficients of the Standard Model Extension. The Bayes factors output in our analysis shows a very strong evidence for the spectral-lag transition in GRB~190114C. Our constraints on a variety of isotropic and anisotropic coefficients for LIV are somewhat weaker than existing bounds, but they can be viewed as comparatively robust and have the promise to complement existing LIV constraints. The observations of GRBs with higher-energy emissions and higher temporal resolutions will contribute to a better formulation of the intrinsic time lag and more rigorous LIV constraints in the dispersive photon sector.

preprint2020arXiv

Neutron Star Mergers in AGN Accretion Disks: Cocoon and Ejecta Shock Breakouts

Neutron star mergers are believed to occur in accretion disks around supermassive black holes. Here we show that a putative jet launched from the merger of a binary neutron star (BNS) or a neutron star--black hole (NSBH) merger occurring at the migration trap in an active galactic nucleus (AGN) disk would be choked. The jet energy is deposited within the disk materials to power a hot cocoon. The cocoon is energetic enough to break out from the AGN disk and produce a bright X-ray shock breakout transient peaking at $\sim0.15\,{\rm d}$ after the merger. The peak luminosity is estimated as $\sim 10^{46}\,{\rm erg}\,{\rm s}^{-1}$, which can be discovered by Einstein Probe from $z\lesssim 0.5$. Later on, the non-relativistic ejecta launched from the merger would break out the disk, powering an X-ray/UV flare peaking at $\sim 0.5\,{\rm d}$ after the merger. This second shock breakout signal may be detected by UV transient searches. The cocoon cooling emission and kilonova emission are outshone by the disk emission and difficult to be detected. Future joint observations of gravitational waves from BNS/NSBH mergers and associated two shock breakout signatures can provide a strong support for the compact binary coalescence formation channel in AGN disks.

preprint2020arXiv

On the true fractions of repeating and non-repeating FRB sources

Observationally, fast radio bursts (FRBs) can be divided into repeating and apparently non-repeating (one-off) ones. It is unclear whether all FRBs repeat and whether there are genuine non-repeating FRBs. We attempt to address these questions using Monte Carlo simulations. We define a parameter $T_c$ at which the accumulated number of non-repeating sources becomes comparable to the total number of the repeating sources, which is a good proxy to denote the intrinsic repeater fraction among FRBs. Assuming that both types of sources exist and that their burst energies follow power law distributions, we investigate how the {\em observed} repeater fraction evolves with time for different parameters. If the lifetime of repeaters is sufficiently long so that the evolutionary effect can be neglected within the observational time span, unless $T_c \rightarrow \infty$ (i.e. there is no genuine non-repeating FRB source) the observed repeater fraction should increase with time first, reach a peak, and then decline. The peak time $T_p$ and the peak fraction $F_{\rm r,obs,p}$ depend on $T_c$ and other repeating rate parameters. With the current data, we pose a lower limit $T_c > 0.1$ d for reasonable parameter values. We predict that future continuous monitoring of FRBs with CHIME or similar wide-field radio telescopes would obtain an $F_{\rm r,obs}$ less than $0.04$. The detection of a smaller peak value $F_{\rm r,obs,p}<0.04$ in the near future would disfavor the ansatz that "all FRB sources repeat".

preprint2020arXiv

Physical Implications of the Sub-threshold GRB GBM-190816 and its Associated Sub-threshold Gravitational Wave Event

The LIGO-Virgo and Fermi collaborations recently reported a possible joint detection of a sub-threshold gravitational wave (GW) event and a sub-threshold gamma-ray burst (GRB), GBM-190816, that occurred 1.57 s after the merger. We perform an independent analysis of the publicly available data and investigate the physical implications of this potential association. By carefully studying the following properties of GBM-190816 using Fermi/GBM data, including signal-to-noise ratio, duration, f-parameter, spectral properties, energetic properties, and its compliance with some GRB statistical correlations, we confirm that this event is likely a typical short GRB. Assuming its association with the sub-threshold GW event, the inferred luminosity is $1.47_{-1.04}^{+3.40} \times 10^{49}$ erg s$^{-1}$. Based on the available information of the sub-threshold GW event, we infer the mass ratio q of the compact binary as $q=2.26_{-1.43}^{+2.75}$ according to the reported range of luminosity distance. If the heavier compact object has a mass > 3 solar masses, q can be further constrained to $q=2.26_{-0.12}^{+2.75}$. The leading physical scenario invokes an NS-BH merger system with the NS tidally disrupted. Within this scenario, we constrain the physical properties of such a system to produce a GRB. The GW data may also allow an NS-BH system with no tidal disruption of the NS or a BH-BH merger. We apply the charged compact binary coalescence (cCBC) theory (for both a constant charge and an increasing charge for the merging members) to derive the model parameters to account for GBM-190816 and found that the required parameters are extreme. Finally, we argue that the fact that the observed GW-GRB delay time scale is comparable to that of GW170817/GRB 170817A suggests that the GW-GRB time delay of these two cases is mainly defined by the time scale for the jet to propagate to the GRB emission site.

preprint2020arXiv

Test the hypothesis of compact-binary-coalescence origin of fast radio bursts through a multi-messenger approach

In the literature, compact binary coalescences (CBCs) have been proposed as one of the main scenarios to explain the origin of some non-repeating fast radio bursts (FRBs). The large discrepancy between the FRB and CBC event rate densities suggest their associations, if any, should only apply at most for a small fraction of FRBs. Through a Bayesian estimation method, we show how a statistical analysis of the coincident associations of FRBs with CBC gravitational wave (GW) events may test the hypothesis of these associations. We show that during the operation period of advanced LIGO, the detection of $\sim100$ ($\sim1000$) GW-less FRBs with dispersion measure (DM) values smaller than 500 ${\rm pc~cm^{-3}}$ could reach the constraint that less than $10\%$ (or $1\%$) FRBs are related to binary black hole (BBH) mergers. The same number of FRBs with DM values smaller than 100 ${\rm pc~cm^{-3}}$ is required to reach the same constraint for binary neutron star (BNS) mergers. With the upgrade of GW detectors, the same constraints for BBH and BNS mergers can be reached with less FRBs or looser requirements for the DM values. It is also possible to pose constraints on the fraction of each type of CBCs that are able to produce observable FRBs based on the event density of FRBs and CBCs. This would further constrain the dimensionless charge of black holes in binary BH systems.

preprint2020arXiv

The evolution of a newborn millisecond magnetar with a propeller-recycling disk

A rapidly rotating and highly magnetized neutron star (NS) could be formed from the explosive phenomena such as superluminous supernovae and gamma-ray bursts. This newborn NS can substantially influence the emission of these explosive transients through its spin-down. The spin-down evolution of the NS can sometimes be affected by fallback accretion, although it is usually regulated by the magnetic dipole radiation and gravitational wave radiation of the NS. Under appropriate conditions, the accreting material can be firstly ejected and subsequently recycled back, so that the accretion disk can keep in a quasi-steady state for a long time. Here we describe the interaction of the NS with such a propeller-recycling disk and their co-evolution. Our result shows that, the spin-down of the NS can be initially dominated by the propeller, which prevents the disk material from falling onto the NS until hundreds or thousands of seconds later. It is suggested that the abrupt fall of the disk material onto the NS could significantly suppress the magnetic dipole radiation and then convert the NS from a normal magnetar to a low-field magnetar. This evolution behavior of the newborn NS can help to understand the very different influence of the NS on the early GRB afterglows and the late supernova/kilonova emission.

preprint2020arXiv

The Second Plateau in X-ray Afterglow Providing Additional Evidence for Rapidly Spinning Magnetars as the GRB Central Engine

Evidence for the central engine of gamma-ray bursts (GRBs) has been collected in the Neil Gehrels \emph{Swift} data. For instance, some GRBs show an internal X-ray plateau followed by very steep decay, which is difficult to be interpreted within the framework of a black hole (BH) central engine, but are consistent within a rapidly spinning magnetar engine picture. The very steep decay at the end of the plateau suggests a sudden cessation of the central engine, which is explained as the collapse of a supra-massive magnetar into a black hole when it spins down. Here we propose that some additional evidence, such as a second X-ray plateau feature would show up, if the fall-back accretion could activate the newborn BH and sufficient energy could be transferred from the newborn BH to the GRB blast wave. With a systematic data analysis for all long GRBs, we find three candidates in \emph{Swift} sample, i.e., GRBs 070802, 090111, and 120213A, whose X-ray afterglow light curves contain two plateaus, with the first one being an internal plateau. We find that in a fairly loose and reasonable parameter space, the second X-ray plateau data for all 3 GRBs could be well interpreted with our proposed model. Future observations are likely to discover more similar events, which could offer more information of the properties of the magnetar as well as the newborn BH.

preprint2020arXiv

What binary systems are the most likely sources for periodically repeating FRBs?

The newly discovered 16.35 days period for repeating FRB 180916.J0158+65 provides an essential clue for understanding the sources and emission mechanism of repeating FRBs. Many models propose that the periodically repeating FRBs might be related to binary star systems that contain at least one neutron star (NSC-FRB system). It has been suggested that the NS "combed" by the strong wind from a companion star might provide a solution. Following the "Binary Comb" model, we use the population synthesis method to study in detail the properties of the companion stars and the nature of NSC-FRB systems. Our main findings are: 1) the companion star is most likely to be a B-type star; 2) the period of 16 days of FRB 180916 happens to fall in the most probable period range, which may explain why FRB 180916 was the first detected periodically repeating FRB, and we expect to observe more periodically repeating FRBs with periods around 10-30 days; 3) the birth rate for NSC-FRB system is large enough to fulfill the event rate requirement set by the observation of FRB 180916, which supports the proposal that NSC-FRB can provide one significant channel for producing periodically repeating FRBs.

preprint2020arXiv

What constraints on the neutron star maximum mass can one pose from GW170817 observations?

The post-merger product of the first binary neutron star merger event detected in gravitational waves, GW170817, depends on neutron star equation of state (EoS) and is not well determined. We generally discuss the constraints one may pose on the maximum mass of a non-spinning neutron star, $M_{\rm TOV}$, based on the observations and some EoS-independent universal relations of rapidly-spinning neutron stars. If the merger product is a black hole after a brief hypermassive neutron star (HMNS) phase, we derive $M_{\rm TOV} < 2.09^{+0.11}_{-0.09}(^{+0.06}_{-0.04}) M_{\odot}$ at the 2$σ$ (1$σ$) level. The cases for a massive neutron star (MNS), either a supra-massive neutron star (SMNS) or even a stable neutron star (SNS), are also allowed by the data. We derive $2.09^{+0.11}_{-0.09}(^{+0.06}_{-0.04} M_{\odot}) \leq M_{\rm TOV}< 2.43^{+0.10}_{-0.08}(^{+0.06}_{-0.04}) M_{\odot}$ for the SMNS case and $M_{\rm TOV} \geq 2.43^{+0.10}_{-0.08}(^{+0.06}_{-0.04})M_{\odot}$ for the SNS case, at the $2σ$ ($1σ$) confidence level. In the MNS cases, we also discuss the constraints on the neutron star parameters (the dipolar magnetic field strength at the surface $B_p$ and the ellipticity $ε$) that affect the spindown history, by considering different MNS survival times, e.g. 300 s, 1 d, and 155 d after the merger, as suggested by various observational arguments. We find that once an SMNS is formed, without violating the EM observational constraints, there always exist a set of ($B_p, ε$) parameters that allow the SMNS to survive for 300s, 1 d, 155 d, or even longer.

preprint2019arXiv

Overview to the Hard X-ray Modulation Telescope (Insight-HXMT) Satellite

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

preprint2016arXiv

Black hole central engine for ultra-long gamma-ray burst 111209A and its associated supernova 2011kl

Recently, the first association between an ultra-long gamma-ray burst (GRB) and a supernova is reported, i.e., GRB 111209A/SN 2011kl, which enables us to investigate the physics of central engines or even progenitors for ultra-long GRBs. In this paper, we inspect the broad-band data of GRB 111209A/SN 2011kl. The late-time X-ray lightcurve exhibits a GRB 121027A-like fall-back bump, suggesting a black hole central engine. We thus propose a collapsar model with fall-back accretion for GRB 111209A/SN 2011kl. The required model parameters, such as the total mass and radius of the progenitor star, suggest that the progenitor of GRB 111209A is more likely a Wolf-Rayet star instead of blue supergiant, and the central engine of this ultra-long burst is a black hole. The implications of our results is discussed.

preprint2016arXiv

Catching jetted tidal disruption events early in millimetre

Relativistic jets can form from at least some tidal disruption events (TDEs) of (sub-)stellar objects around supermassive black holes. We detect the millimeter (MM) emission of IGR J12580+0134 --- the nearest TDE known in the galaxy NGC 4845 at the distance of only 17 Mpc, based on Planck all-sky survey data. The data show significant flux jumps after the event, followed by substantial declines, in all six high frequency Planck bands from 100 GHz to 857 GHz. We further show that the evolution of the MM flux densities are well consistent with our model prediction from an off-axis jet, as was initially suggested from radio and X-ray observations. This detection represents the second TDE with MM detections; the other is Sw J1644+57, an on-axis jetted TDE at redshift of 0.35. Using the on- and off-axis jet models developed for these two TDEs as templates, we estimate the detection potential of similar events with the Large Millimeter Telescope (LMT) and the Atacama Large Millimeter/submillimeter Array (ALMA). Assuming an exposure of one hour, we find that the LMT (ALMA) can detect jetted TDEs up to redshifts $z\sim1$ (2), for a typical disrupted star mass of $\sim1$ M$_\odot$. The detection rates of on- and off-axis TDEs can be as high as $\sim0.6$ (13) and 10 (220) per year, respectively, for the LMT (ALMA). We briefly discuss how such observations, together with follow-up radio monitoring, may lead to major advances in understanding the jetted TDEs themselves and the ambient environment of the CNM.

preprint2016arXiv

Constraints on binary neutron star merger product from short GRB observations

Binary neutron star mergers are strong gravitational wave (GW) sources and the leading candidates to interpret short duration gamma-ray bursts (SGRBs). Under the assumptions that SGRBs are produced by double neutron star mergers and that the X-ray plateau followed by a steep decay as observed in SGRB X-ray light curves marks the collapse of a supra-massive neutron star to a black hole (BH), we use the statistical observational properties of {\em Swift} SGRBs and the mass distribution of Galactic double neutron star systems to place constraints on the neutron star equation of state (EoS) and the properties of the post-merger product. We show that current observations already put following interesting constraints: 1) A neutron star EoS with a maximum mass close to a parameterization of $M_{\rm max} = 2.37\,M_\odot (1+1.58\times10^{-10} P^{-2.84})$ is favored; 2) The fractions for the several outcomes of NS-NS mergers are as follows: $\sim40\%$ prompt BHs, $\sim30\%$ supra-massive NSs that collapse to BHs in a range of delay time scales, and $\sim30\%$ stable NSs that never collapse; 3) The initial spin of the newly born supra-massive NSs should be near the breakup limit ($P_i\sim1 {\rm ms}$), which is consistent with the merger scenario; 4) The surface magnetic field of the merger products is typically $\sim 10^{15}$ G; 5) The ellipticity of the supra-massive NSs is $ε\sim (0.004 - 0.007)$, so that strong GW radiation is released post the merger; 6) Even though the initial spin energy of the merger product is similar, the final energy output of the merger product that goes into the electromagnetic channel varies in a wide range from several $10^{49}$ erg to several $10^{52}$ erg, since a good fraction of spin energy is either released in the form of GW or falls into the black hole as the supra-massive NS collapses.

preprint2016arXiv

Constraints on the Photon Mass with Fast Radio Bursts

Fast radio bursts (FRBs) are radio bursts characterized by millisecond durations, high Galactic latitude positions, and high dispersion measures. Very recently, the cosmological origin of FRB 150418 has been confirmed by \cite{kea16}, and FRBs are now strong competitors as cosmological probes. The simple sharp feature of the FRB signal is ideal for them to probe some of the fundamental laws of physics. Here we show that by analyzing the delay time between different frequencies, the FRB data can place stringent upper limits on the rest mass of the photon. For FRB 150418 at $z=0.492$, one can potentially reach $m_γ \leq 5.2 \times 10^{-47}$ g, which is $10^{20}$ times smaller than the rest mass of electron, and is about $10^3$ times smaller than that obtained using other astrophysical sources with the same method.

preprint2016arXiv

Internal X-ray plateau in short GRBs: Signature of supramassive fast-rotating quark stars?

A supramassive, strongly-magnetized millisecond neutron star (NS) has been proposed to be the candidate central engine of at least some short gamma-ray bursts (SGRBs), based on the "internal plateau" commonly observed in the early X-ray afterglow. While a previous analysis shows a qualitative consistency between this suggestion and the Swift SGRB data, the distribution of observed break time $t_b$ is much narrower than the distribution of the collapse time of supramassive NSs for the several NS equations-of-state (EoSs) investigated. In this paper, we study four recently-constructed "unified" NS EoSs, as well as three developed strange quark star (QS) EoSs within the new confinement density-dependent mass model. All the EoSs chosen here satisfy the recent observational constraints of the two massive pulsars whose masses are precisely measured. We construct sequences of rigidly rotating NS/QS configurations with increasing spinning frequency $f$, from non-rotating ($f = 0$) to the Keplerian frequency ($f = f_{\rm K}$), and provide convenient analytical parametrizations of the results. Assuming that the cosmological NS-NS merger systems have the same mass distribution as the Galactic NS-NS systems, we demonstrate that all except the BCPM NS EoS can reproduce the current $22\%$ supramassive NS/QS fraction constraint as derived from the SGRB data. We simultaneously simulate the observed quantities (the break time $t_b$, the break time luminosity $L_b$ and the total energy in the electromagnetic channel $E_{\rm total}$) of SGRBs, and find that while equally well reproducing other observational constraints, QS EoSs predict a much narrower $t_b$ distribution than that of the NS EoSs, better matching the data. We therefore suggest that the post-merger product of NS-NS mergers might be fast-rotating supramassive QSs rather than NSs.

preprint2016arXiv

Limits on the Neutrino Velocity, Lorentz Invariance, and the Weak Equivalence Principle with TeV Neutrinos from Gamma-Ray Bursts

Five TeV neutrino events weakly correlated with five gamma-ray bursts (GRBs) were detected recently by IceCube. This work is an attempt to show that if the GRB identifications are verified, the observed time delays between the TeV neutrinos and gamma-ray photons from GRBs provide attractive candidates for testing fundamental physics with high accuracy. Based on the assumed associations between the TeV neutrinos and GRBs, we find that the limiting velocity of the neutrinos is equal to that of photons to an accuracy of $\sim1.9\times10^{-15}-2.5\times10^{-18}$, which is about $10^{4}-10^{7}$ times better than the constraint obtained with the neutrino possibly from a blazar flare. In addition, we set the most stringent limits up to date on the energy scale of quantum gravity for both the linear and quadratic violations of Lorentz invariance, namely $E_{\rm QG, 1}>6.3\times10^{18}-1.5\times10^{21}$ GeV and $E_{\rm QG, 2}>2.0\times10^{11}-4.2\times10^{12}$ GeV, which are essentially as good as or are an improvement of one order of magnitude over the results previously obtained by the GeV photons of GRB 090510 and the PeV neutrino from a blazar flare. Assuming that the Shapiro time delay is caused by the gravitational potential of the Laniakea supercluster of galaxies, we also place the tightest limits to date on Einstein's weak equivalence principle through the relative differential variations of the parameterized post-Newtonian parameter $γ$ values for two different species of particles (i.e., neutrinos and photons), yielding $Δγ\sim 10^{-11}-10^{-13}$. However, it should be emphasized again that these limits here obtained are at best forecast of what could be achieved if the GRB/neutrino correlations would be finally confirmed.

preprint2016arXiv

Testing Einstein's Weak Equivalence Principle With Gravitational Waves

A conservative constraint on the Einstein Weak Equivalence Principle (WEP) can be obtained under the assumption that the observed time delay between correlated particles from astronomical sources is dominated by the gravitational fields through which they move. Current limits on the WEP are mainly based on the observed time delays of photons with different energies. It is highly desirable to develop more accurate tests that include the gravitational wave (GW) sector. The detection by the advanced LIGO/VIRGO systems of gravitational waves will provide attractive candidates for constraining the WEP, extending the tests to gravitational interactions, with potentially higher accuracy. Considering the capabilities of the advanced LIGO/VIRGO network and the source direction uncertainty, we show that the joint detection of GWs and electromagnetic signals could probe the WEP to an accuracy down to $10^{-10}$, which is one order of magnitude tighter than previous limits, and seven orders of magnitude tighter than the multi-messenger (photons and neutrinos) results by supernova 1987A.

preprint2016arXiv

Tests of the Einstein Equivalence Principle using TeV Blazars

The observed time delays between different energy bands from TeV blazars provide a new interesting way of testing the Einstein Equivalence Principle (EEP). If the whole time delay is assumed to be dominated by the gravitational field of the Milky Way, the conservative upper limit on the EEP can be estimated. Here we show that the strict limits on the differences of the parameterized post-Newtonian parameter $γ$ values are $γ_{\rm TeV}-γ_{\rm keV}<3.86\times10^{-3}$ for Mrk 421 and $γ_{\rm TeV}-γ_{\rm keV}<4.43\times10^{-3}$ for Mrk 501, while expanding the scope of the tested EEP energy range out to the TeV--keV range for the first time. With the small time lag from the 0.2--0.8 TeV and $>0.8$ TeV light curves of PKS 2155-304, a much more severe constraint on $γ$ differences of $\sim10^{-6}$ can be achieved, although the energy difference is of order of $\sim$ TeV. Furthermore, we can combine these limits on the energy dependence of $γ$ with the bound on the absolute $γ$ value $γ-1\sim0.3\%$ from light deflection measurements at the optical (eV) bands, and conclude that this absolute bound on $γ$ can be extended from optical to TeV energies.

preprint2015arXiv

A morphological analysis of gamma-ray burst early optical afterglows

Within the framework of the external shock model of gamma-ray bursts (GRBs) afterglows, we perform a morphological analysis of the early optical lightcurves to directly constrain model parameters. We define four morphological types, i.e. the reverse shock dominated cases with/without the emergence of the forward shock peak (Type I/ Type II), and the forward shock dominated cases without/with $ν_m$ crossing the band (Type III/IV). We systematically investigate all the Swift GRBs that have optical detection earlier than 500 s and find 3/63 Type I bursts (4.8%), 12/63 Type II bursts (19.0%), 30/63 Type III bursts (47.6%), 8/63 Type IV bursts (12.7%) and 10/63 Type III/IV bursts (15.9%). We perform Monte Carlo simulations to constrain model parameters in order to reproduce the observations. We find that the favored value of the magnetic equipartition parameter in the forward shock ($ε_B^f$) ranges from $10^{-6}$ to $10^{-2}$, and the reverse-to-forward ratio of $ε_B$ ($R_B$) is about 100. The preferred electron equipartition parameter $ε_e^{r,f}$ value is 0.01, which is smaller than the commonly assumed value, e.g., 0.1. This could mitigate the so- called "efficiency problem" for the internal shock model, if $ε_e$ during the prompt emission phase (in the internal shocks) is large (say, $\sim 0.1$). The preferred $R_B$ value is in agreement with the results in previous works that indicates a moderately magnetized baryonic jet for GRBs.

preprint2015arXiv

Constraints on Lorentz Invariance Violation with gamma-ray bursts via a Markov Chain Monte Carlo approach

In quantum theory of gravity, we expect the Lorentz Invariance Violation (LIV) and the modification of the dispersion relation between energy and momentum for photons. The effect of the energy-dependent velocity due to the modified dispersion relation for photons was studied in the standard cosmological context by using a sample of Gamma Ray Bursts (GRBs). In this paper we mainly discuss the possible LIV effect by using different cosmological models for the accelerating universe. Due to the degeneracies among model parameters, the GRBs' time delay data are combined with the cosmic microwave background data from the Planck first year release, the baryon acoustic oscillation data at six different redshifts, as well as Union2 type Ia supernovae data, to constrain both the model parameters and the LIV effect. We find no evidence of LIV.

preprint2015arXiv

Cosmic Transients Test Einstein's Equivalence Principle out to GeV Energies

The Einstein Equivalence Principle (EEP) can be probed with astrophysical sources emitting simultaneously different types of neutral particles, or particles with varying energies, by testing their time of flight through the same gravitational field. Here we use the time delays between correlated photons from cosmological transients to constrain the accuracy of the EEP. We take data from two gamma-ray bursts as an example, and use, as a lower limit to the theoretical time delays between different energies, delays arising from only the gravitational field of our own galaxy. We then show that the parameterized post-Newtonian parameter $γ$ is the same for photons over energy ranges between eV and MeV and between MeV and GeV to a part in $10^{-7}$, which is at least one order of magnitude better than previous limits. Combining this bound on the wavelength dependence of $γ$ with the absolute bound $|γ-1|<0.3\%$ from light-deflection measurements at optical (eV) wavelengths, we thus extend this absolute bound on $γ$ to GeV energies.

preprint2015arXiv

Gamma Ray Bursts in the HAWC Era

Gamma-Ray Bursts are the most energetic explosions in the Universe, and are among the most promising for detecting multiple non-electromagnetic signals, including cosmic rays, high energy neutrinos and gravitational waves. The multi-GeV to TeV gamma-ray range of GRB could have significant contributions from hadronic interactions, mixed with more conventional leptonic contributions. This energy range is important for probing the source physics, including overall energetics, the shock parameters and the Lorentz factor. We discuss some of the latest observational and theoretical developments in the field.

preprint2015arXiv

GRB 080503 late afterglow re-brightening: signature of a magnetar powered merger-nova

GRB 080503 is a short gamma-ray burst (GRB) detected by \emph{Swift} and has been classified as a compact-star-merger-origin GRB. The soft extended emission and the simultaneous late re-brightening in both the X-ray and optical afterglow lightcurves raise interesting questions regarding its physical origin. We show that the broad-band data of GRB 080503 can be well explained within the framework of the double neutron star merger model, provided that the merger remnant is a rapidly-rotating massive neutron star with an extremely high magnetic field (i.e. a millisecond magnetar). We show that the late optical re-brightening is consistent with the emission from a magnetar-powered "merger-nova". This adds one more case to the growing sample of merger-novae associated with short GRBs. The soft extended emission and the late X-ray excess emission are well connected through a magnetar dipole spin-down luminosity evolution function, suggesting that direct magnetic dissipation is the mechanism to produce these X-rays. The X-ray emission initially leaks from a hole in the merger ejecta pierced by the short GRB jet. The hole subsequently closes after the magnetar spins down and the magnetic pressure drops below ram pressure. The X-ray photons are then trapped behind the mergernova ejecta until the ejecta becomes optically thin at a later time. This explains the essentially simultaneous re-brightening in both the optical and X-ray lightcurves. Within this model, future gravitational wave sources could be associated with a bright X-ray counterpart along with the mergernova, even if the short GRB jet beams away from Earth.

preprint2015arXiv

How Bad/Good Are the External Forward Shock Afterglow Models of Gamma-Ray Bursts?

The external forward shock (EFS) models have been the standard paradigm to interpret the broad-band afterglow data of gamma-ray bursts (GRBs). One prediction of the models is that some afterglow temporal breaks at different energy bands should be achromatic. Observations in the Swift era have revealed chromatic afterglow behaviors at least in some GRBs, casting doubts on the EFS origin of GRB afterglows. In this paper, we perform a systematic study to address the question: how bad/good are the external forward shock models? Our sample includes 85 GRBs well-monitored X-ray and optical lightcurves. Based on how well the data abide by the EFS models, we categorize them as: Gold sample: (Grade I and II) include 45/85 GRBs. They show evidence of, or are consistent with having, an achromatic break. The temporal/spectral behaviors in each afterglow segment are consistent with the predictions (closure relations) of the EFS models. Silver sample: (Grade III and IV) include 37/85 GRBs. They are also consistent with having an achromatic break, even though one or more afterglow segments do not comply with the closure relations. Bad sample: (Grade V), 3/85 shows direct evidence of chromatic behaviors, suggesting that the EFS models are inconsistent with the data. These are included in the Bad sample. We further perform statistical analyses of various observational properties ($α$, $β$, $t_b$ and model parameters (energy injection index q, p, $θ_j$, $η_γ$, etc) of the GRBs in the Gold Sample, and derive constraints on the magnetization parameter $ε_B$ in the EFS. Overall, we conclude that the simplest EFS models can account for the multi-wavelength afterglow data of at least half of the GRBs. When more advanced modeling (e.g., long-lasting reverse shock, structured jets) is invoked, up to $>90 \%$ of the afterglows may be interpreted within the framework of the ESF models.

preprint2015arXiv

Photosphere emission from a hybrid relativistic outflow with arbitrary dimensionless entropy and magnetization in GRBs

In view of the recent Fermi observations of GRB prompt emission spectra, we develop a theory of photosphere emission of a hybrid relativistic outflow with a hot fireball component (defined by dimensionless entropy $η$) and a cold Poynting-flux component (defined by magnetization $σ_0$ at the central engine). We consider the scenarios both without and with sub-photospheric magnetic dissipations. Based on a simplified toy model of jet dynamics, we develop two approaches: a "bottom-up" approach to predict the temperature (for a non-dissipative photosphere) and luminosity of the photosphere emission and its relative brightness for a given pair of $(η,σ_0)$; and a "top-down" approach to diagnose central engine parameters ($η$ and $σ_0$) based on the observed quasi-thermal photosphere emission properties. We show that a variety of observed GRB prompt emission spectra with different degrees of photosphere thermal emission can be reproduced by varying $η$ and $σ_0$ within the non-dissipative photosphere scenario. In order to reproduce the observed spectra, the outflows of most GRBs need to have a significant $σ$, both at the central engine, and at the photosphere. The $σ$ value at $10^{15}$ cm from the central engine (a possible non-thermal emission site) is usually also greater than unity, so that internal-collision-induced magnetic reconnection and turbulence (ICMART) may be the mechanism to power the non-thermal emission. We apply our top-down approach to GRB 110721A, and find that the temporal evolution behavior of its blackbody component can be well interpreted with a time-varying $(η,σ_0)$ at the central engine, instead of invoking a varying engine base size $r_0$ as proposed by previous authors.

preprint2015arXiv

Relation between the intrinsic and observed central engine activity time: implications for ultra-long GRBs

The GRB central engine intrinsic activity time $T_{\rm ce}$ is usually described through either the $γ$-ray duration $T_{90}$ or through a generalized burst duration $t_{\rm burst}$ which includes both the $γ$-ray emission and (when present) an extended flaring X-ray plateau. Here, we define a more specific operational description of $T_{\rm ce}$, and within the framework of the internal-external shock model, we develop a numerical code to study the relationship between $T_{90}$ and $T_{\rm ce}$, as well as between $t_{\rm burst}$ and $T_{\rm ce}$, for different initial conditions. We find that when $T_{\rm ce}\lesssim 10^4$ s, late internal collisions or refreshed external collisions result in values of $T_{\rm 90}$ and $t_{\rm burst}$ larger than $T_{\rm ce}$, usually by factors of $2-3$. For $T_{\rm ce}\gtrsim 10^4$ s, the $t_{\rm burst}$ is always a good estimator for $T_{\rm ce}$, while $T_{90}$ can underpredict $T_{\rm ce}$ when the late central engine activity is moderate. We find a clear bimodal distribution for $T_{\rm ce}$, based on our simulations as well as on the observational data for $T_{90}$ and $t_{\rm burst}$. We suggest that $t_{\rm burst}$ is a reliable measure for defining "ultra-long" GRBs. Bursts with $T_{90}$ of order $10^3$ s need not belong to a special population, while bursts with $t_{\rm burst} > 10^4$ s, where the late central engine activity is more moderate and shows up in X-rays, may represent a new population. These conclusions are insensitive to the initial conditions assumed in the models.

preprint2015arXiv

Reverse Shock Emission in Gamma-ray Bursts Revisited

A generic synchrotron external shock model is the widely preferred paradigm used to interpret the broad-band afterglow data of gamma-ray bursts (GRBs), including predicted observable signatures from a reverse shock which have been confirmed by observations. Investigations of the nature of the reverse shock emission can provide valuable insights into the intrinsic properties of the GRB ejecta. Here we briefly review the standard and the extended models of the reverse shock emission, discussing the connection between the theory and observations, including the implications of the latest observational advances.

preprint2014arXiv

A double neutron star merger origin for the cosmological relativistic fading source PTF11agg?

The Palomar Transient Factory (PTF) team recently reported the discovery of a rapidly fading optical transient source, PTF11agg. A long-lived scintillating radio counterpart was identified, but the search for a high energy counterpart showed negative results. The PTF team speculated that PTF11agg may represent a new class of relativistic outbursts. Here we suggest that a neutron star (NS)-NS merger system with a supra-massive magnetar central engine could be a possible source to power such a transient, if our line of sight is not on the jet axis direction of the system. These systems are also top candidates for gravitational wave sources to be detected in the advanced LIGO/Virgo era. We find that the PTF11agg data could be explained well with such a model, suggesting that at least some gravitational wave bursts due to NS-NS mergers may be associated with such a bright electromagnetic counterpart without a γ-ray trigger.

preprint2014arXiv

Fast Radio Burst/Gamma-Ray Burst Cosmography

Recently, both theoretical arguments and observational evidence suggested that a small fraction of fast radio bursts (FRBs) could be associated with gamma-ray bursts (GRBs). If such FRB/GRB association systems are commonly detected in the future, the combination of dispersion measure (DM) derived from FRBs and redshifts derived from GRBs makes these systems a plausible tool to conduct cosmography. We quantify uncertainties in deriving the redshift-dependent DM_{IGM} as a function of z, and test how well dark energy models can be constrained with Monte Carlo simulations. We show that with potentially several 10s of FRB/GRB systems detected in a decade or so, one may reach reasonable constraints on wCDM models. When combined with SN Ia data, unprecedented constraints on dark energy equation of state may be achieved, thanks to the prospects of detecting FRB/GRB systems at relatively high redshifts. The ratio between the mean value < DM_{IGM} (z)> and luminosity distance (D_{L} (z)) is insensitive to dark energy models. This gives the prospects of applying SN Ia data to calibrate < DM_{IGM} (z)> using a relatively small sample of FRB/GRB systems, allowing a reliable constraint on the baryon inhomogeneity distribution as a function of redshift. The methodology developed in this paper can also be applied, if the FRB redshifts can be measured by other means. Some caveats of putting this method into practice are also discussed.

preprint2014arXiv

Multi-wavelength afterglows of fast radio bursts

The physical nature of fast radio bursts (FRBs) is not identified. Detecting electromagnetic counterparts in other wavelengths is essential to measure their distances and to settle down their physical nature. Assuming that at least some of them are of a cosmological origin, we calculate their afterglow lightcurves in multi-wavelengths (X-rays, optical and radio) by assuming a range of their total kinetic energies and redshifts. We focus on forward shock emission, but also consider the possibility that some of them might have bright reverse shock emission. In general, the FRB afterglows are too faint to be detected by current detectors. Only if an FRB has a very low radiative efficiency in radio (hence, a very large kinetic energy), and when it is close enough, can its afterglow be detected in the optical and radio bands. We discuss observational strategies to detect these faint afterglows using future telescopes such as LSST and EVLA.

preprint2014arXiv

Quasi-Periodic Variations in X-ray Emission and Long-Term Radio Observations: Evidence for a Two-Component Jet in Sw J1644+57

The continued observations of Sw J1644+57 in X-ray and radio bands accumulated a rich data set to study the relativistic jet launched in this tidal disruption event. The X-ray light curve of Sw J1644+57 from 5-30 days presents two kinds of quasi-periodic variations: a 200 second quasi-periodic oscillation (QPO) and a 2.7-day quasi-periodic variation. The latter has been interpreted by a precessing jet launched near the Bardeen-Petterson radius of a warped disk. Here we suggest that the $\sim$ 200s QPO could be associated with a second, narrower jet sweeping the observer line-of-sight periodically, which is launched from a spinning black hole in the misaligned direction with respect to the black hole's angular momentum. In addition, we show that this two-component jet model can interpret the radio light curve of the event, especially the re-brightening feature starting $\sim 100$ days after the trigger. From the data we infer that inner jet may have a Lorentz factor of $Γ_{\rm j} \sim 5.5$ and a kinetic energy of $E_{\rm k,iso} \sim 3.0 \times 10^{52} {\rm erg}$, while the outer jet may have a Lorentz factor of $Γ_{\rm j} \sim 2.5$ and a kinetic energy of $E_{\rm k,iso} \sim 3.0 \times 10^{53} {\rm erg}$.

preprint2014arXiv

Variability of the giant X-ray bump in GRB 121027A and possible origin

The particular giant X-ray bump of GRB 121027A triggered by \emph{Swift} is quite different from the typical X-ray flares in gamma-ray bursts. There exhibit four parts of the observed structural variabilities in the rise and decay phase of the bump. Considering the quality of four parts of the data, we can only analyze the data from about 5300 s to about 6100 s in the bump using the stepwise filter correlation method (Gao et al. 2012), and find that the $86^{+5.9}_{-9.4}~\rm s$ periodic oscillation may exist, which is confirmed by the Lomb-Scargle method (Scargle 1982). Furthermore, a jet precession model (Liu et al. 2010) is proposed to account for such a variability.

preprint2013arXiv

A Complete Reference of the Analytical Synchrotron External Shock Models of Gamma-Ray Bursts

Gamma-ray bursts are most luminous explosions in the universe. Their ejecta are believed to move towards Earth with a relativistic speed. The interaction between this "relativistic jet" and a circum burst medium drives a pair of (forward and reverse) shocks. The electrons accelerated in these shocks radiate synchrotron emission to power the broad-band afterglow of GRBs. The external shock theory is an elegant theory, since it invokes a limit number of model parameters, and has well predicted spectral and temporal properties. On the other hand, depending on many factors (e.g. the energy content, ambient density profile, collimation of the ejecta, forward vs. reverse shock dynamics, and synchrotron spectral regimes), there is a wide variety of the models. These models have distinct predictions on the afterglow decaying indices, the spectral indices, and the relations between them (the so-called "closure relations"), which have been widely used to interpret the rich multi-wavelength afterglow observations. This review article provides a complete reference of all the analytical synchrotron external shock afterglow models by deriving the temporal and spectral indices of all the models in all spectral regimes, including some regimes that have not been published before. The review article is designated to serve as a useful tool for afterglow observers to quickly identify relevant models to interpret their data. The limitations of the analytical models are reviewed, with a list of situations summarized when numerical treatments are needed.

preprint2013arXiv

A Comprehensive Study of Gamma-Ray Burst Optical Emission: II. Afterglow Onset and Late Re-Brightening Components

We continue our systematic statistical study of various components in gamma-ray burst (GRB) optical lightcurves. We decompose the early onset bump and the late re-brightening bump with empirical fits and analyze their statistical properties. Among the 146 GRBs that have well-sampled optical lightcurves, the onset and re-brightening bumps are observed in 38 and 26 GRBs, respectively. It is found that the typical rising and decaying slopes for both the onset and re-brightening bumps are ~1.5 and -1.15, respectively. No early onset bumps in the X-ray band are detected to be associated with the optical onset bumps, while an X-ray re-brightening bump is detected for half of the re-brightening optical bumps. The peak luminosity is anti-correlated with the peak time, L_p\propto t_{p}^{-1.81+/-0.32} and L_p\propto t_{p}^{-0.83+/-0.17} for the onset and re-brightening bumps, respectively. Both L_p and the isotropic energy release of the onset bumps are correlated with E_{γ, iso}, whereas no similar correlation is found for the re-brightening bumps. Taking the onset bumps as probes for the properties of the fireballs and their ambient medium, we find that the typical power-law index of the relativistic electrons is 2.5 and the medium density profile behaves as n\propto r^{-1} within the framework of the synchrotron external shock models. With the medium density profile obtained from our analysis, we also confirm the correlation between initial Lorentz factor (Γ_0) and E_{γ, iso} in our previous work. The jet component that produces the re-brightening bump seems to be on-axis and independent of the prompt emission jet component. Its typical kinetic energy budget would be about one order of magnitude larger than the prompt emission component, but with a lower Γ_0, typically several tens.

preprint2013arXiv

Bright "merger-nova" from the remnant of a neutron star binary merger: A signature of a newly born, massive, millisecond magnetar

A massive millisecond magnetar may survive a merger of a neutron star (NS) binary, which would continuously power the merger ejecta. We develop a generic dynamic model for the merger ejecta with energy injection from the central magnetar. The ejecta emission (the "merger-nova") powered by the magnetar peaks in the UV band and the peak of lightcurve progressively shifts to an earlier epoch with increasing frequency. A magnetar-powered mergernova could have an optical peak brightness comparable to a supernova, which is a few tens or hundreds times brighter than the radioactive-powered merger-novae (the so-called macro-nova or kilo-nova). On the other hand, such a merger-nova would peak earlier and have a significantly shorter duration than that of a supernova. An early collapse of the magnetar could suppress the brightness of the optical emission and shorten its duration. Such millisecond-magnetar-powered merger-novae may be detected from NS-NS merger events without an observed short gamma-ray burst, and could be a bright electromagnetic counterpart for gravitational wave bursts due to NS-NS mergers. If detected, it suggests that the merger leaves behind a massive NS, which has important implications for the equation-of-state of nuclear matter.

preprint2013arXiv

Bright broad-band afterglows of gravitational wave bursts from mergers of binary neutron stars

If double neutron star mergers leave behind a massive magnetar rather than a black hole, a bright early afterglow can follow the gravitational wave burst (GWB) even if there is no short gamma-ray burst (SGRB) - GWB association or there is an association but the SGRB does not beam towards earth. Besides directly dissipating the proto-magnetar wind as suggested by Zhang, we here suggest that the magnetar wind could push the ejecta launched during the merger process, and under certain conditions, would reach a relativistic speed. Such a magnetar-powered ejecta, when interacting with the ambient medium, would develop a bright broad-band afterglow due to synchrotron radiation. We study this physical scenario in detail, and present the predicted X-ray, optical and radio light curves for a range of magnetar and ejecta parameters. We show that the X-ray and optical lightcurves usually peak around the magnetar spindown time scale (10^3-10^5s), reaching brightness readily detectable by wide-field X-ray and optical telescopes, and remain detectable for an extended period. The radio afterglow peaks later, but is much brighter than the case without a magnetar energy injection. Therefore, such bright broad-band afterglows, if detected and combined with GWBs in the future, would be a probe of massive millisecond magnetars and stiff equation-of-state for nuclear matter.

preprint2013arXiv

Compton Scattering of Self-Absorbed Synchrotron Emission

Synchrotron self-Compton (SSC) scattering is an important emission mechanism in many astronomical sources, such as gamma-ray bursts (GRBs) and active galactic nuclei (AGNs). We give a complete presentation of the analytical approximations for the Compton scattering of synchrotron emission with both weak and strong synchrotron self-absorption. All possible orders of the characteristic synchrotron spectral breaks ($ν_{\rm a}$, $ν_{\rm m}$, and $ν_{\rm c}$) are studied. In the weak self-absorption regime, i.e., $ν_{\rm a} < ν_c$, the electron energy distribution is not modified by the self-absorption process. The shape of the SSC component broadly resembles that of synchrotron, but with the following features: The SSC flux increases linearly with frequency up to the SSC break frequency corresponding to the self-absorption frequency $ν_{\rm a}$; and the presence of a logarithmic term in the high-frequency range of the SSC spectra makes it harder than the power-law approximation. In the strong absorption regime, i.e. $ν_{\rm a} > ν_{\rm c}$, heating of low energy electrons due to synchrotron absorption leads to pile-up of electrons, and form a thermal component besides the broken power-law component. This leads to two-component (thermal + non-thermal) spectra for both the synchrotron and SSC spectral components. For $ν_{\rm c} < ν_{\rm a} < ν_{\rm m}$, the spectrum is thermal (non-thermal) -dominated if $ν_a > \sqrt{ν_m ν_c}$ ($ν_a < \sqrt{ν_m ν_c}$). Similar to the weak-absorption regime, the SSC spectral component is broader than the simple broken power law approximation. We derive the critical condition for strong absorption (electron pile-up), and discuss a case of GRB reverse shock emission in a wind medium, which invokes $ν_{\rm a} > {\rm max} (ν_{\rm m}, ν_{\rm c})$.

preprint2013arXiv

Possible High-Energy Neutrino and Photon Signals from Gravitational Wave Bursts due to Double Neutron Star Mergers

As the technology of gravitational-wave and neutrino detectors becomes increasingly mature, a multi-messenger era of astronomy is ushered in. Advanced gravitational wave detectors are close to making a ground-breaking discovery of gravitational wave bursts (GWBs) associated with mergers of double neutron stars (NS-NS). It is essential to study the possible electromagnetic (EM) and neutrino emission counterparts of these GWBs. Recent observations and numerical simulations suggest that at least a fraction of NS-NS mergers may leave behind a massive millisecond magnetar as the merger product. Here we show that protons accelerated in the forward shock powered by a magnetar wind pushing the ejecta launched during the merger process would interact with photons generated in the dissipating magnetar wind and emit high energy neutrinos and photons. We estimate the typical energy and fluence of the neutrinos from such a scenario. We find that $\sim$PeV neutrinos could be emitted from the shock front as long as the ejecta could be accelerated to a relativistic speed. The diffuse neutrino flux from these events, even under the most optimistic scenarios, is too low to account for the two events announced by the IceCube Collaboration, but it is only slightly lower than the diffuse flux of GRBs, making it an important candidate for the diffuse background of $\sim$PeV neutrinos. The neutron-pion decay of these events make them a moderate contributor to the sub-TeV gamma-ray diffuse background.

preprint2012arXiv

A Comprehensive Study of Gamma-Ray Burst Optical Emission: I. Flares and Early Shallow Decay Component

Well-sampled optical lightcurves of 146 gamma-ray bursts (GRBs) are compiled from the literature. By empirical fitting we identify eight possible emission components and summarize the results in a "synthetic" lightcurve. Both optical flare and early shallow-decay components are likely related to long-term central engine activities. We focus on their statistical properties in this paper. Twenty-four optical flares are obtained from 19 GRBs. The isotropic R-band energy is smaller than 1% of $E_{γ, \rm iso}$. The relation between isotropic luminosities of the flares and gamma-rays follows $L^{\rm F}_{\rm R, iso}\propto L_{γ, \rm iso}^{1.11\pm 0.27}$. Later flares tend to be wider and dimmer, i.e., $w^{\rm F}\sim t^{\rm F}_{\rm p}/2$ and $L^{\rm F}_{\rm R, iso}\propto [t^{\rm F}_{\rm p}/(1+z)]^{-1.15\pm0.15}$. The detection probability of the optical flares is much smaller than that of X-ray flares. An optical shallow decay segment is observed in 39 GRBs. The relation between the break time and break luminosity is a power-law, with an index of $-0.78\pm 0.08$, similar to that derived from X-ray flares. The X-ray and optical breaks are usually chromatic, but a tentative correlation is found. We suggest that similar to the prompt optical emission that tracks $γ$-rays, the optical flares are also related to the erratic behavior of the central engine. The shallow decay component is likely related to a long-lasting spinning-down central engine or piling up of flare materials onto the blastwave. Mixing of different emission components may be the reason of the diverse chromatic afterglow behaviors.

preprint2012arXiv

Frame-dragging, disk warping, jet precessing, and dipped X-ray lightcurve of Sw J1644+57

The X-ray transient source Sw J1644+57 recently discovered by Swift is believed to be triggered by tidal disruption of a star by a rapidly spinning supermassive black hole (SMBH). For such events, the outer disk is very likely misaligned with respect to the equatorial plane of the spinning SMBH, since the incoming star before disruption most likely has an inclined orbital plane. The tilted disk is subject to the Lense-Thirring torque, which tends to twist and warp due to the Bardeen-Petterson effect. The inner disk tends to align with the SMBH spin, while the outer region tends to remain in the stellar orbital plane, with a transition zone around the Bardeen-Petterson radius. The relativistic jet launched from the spinning SMBH would undergo precession. The 5-30 day X-ray lightcurve of Sw J1644+57 shows a quasi-periodic (2.7-day) variation with noticeable narrow dips. We numerically solve a warped disk and propose a jet-precessing model by invoking a Blandford-Znajek jet collimated by a wind launched near the Bardeen-Petterson radius. Through simulations, we show that the narrow dips in the X-ray lightcurve can be reproduced for a range of geometric configurations. From data we infer that the inclination angle of the initial stellar orbit is in the range of $10^{\circ}-20^{\circ}$ from the SMBH equatorial plane, that the jet should have a moderately high Lorentz factor, and that the inclination angle, jet opening angle, and observer's viewing angle are such that the duty cycle of the line-of-sight sweeping the jet cone is somewhat less than 0.5.

preprint2012arXiv

Stepwise Filter Correlation Method and Evidence of Superposed Variability Components in GRB Prompt Emission Lightcurves

Gamma-ray bursts (GRBs) have variable lightcurves. Although most models attribute the observed variability to one physical origin (e.g. central engine activity, clumpy circumburst medium, relativistic turbulence), some models invoke two physically distinct variability components. We develop a method, namely, the stepwise filter correlation (SFC) method, to decompose the variability components in a GRB lightcurve. Based on a low-pass filter technique, we progressively filter the high frequency signals from the lightcurve, and then perform a correlation analysis between each adjunct pair of filtered lightcurves. Our simulations suggest that if a mock lightcurve contains a slow variability component superposed on a rapidly varying time sequence, the correlation coefficient as a function of the filter frequency would display a prominent dip feature around the frequency of the slow component. Through simulations, we demonstrate that this method can identify significant clustering structures of a lightcurve in the frequency domain, and proved that it can catch superposed signals that are otherwise not easy to retrieve based on other methods (e.g. the power density spectrum analysis method). We apply this method to 266 BATSE bright GRBs. We find that the majority of the bursts have clear evidence of such a superposition effect. We perform a statistical analysis of the identified variability components, and discuss the implications for GRB physics.

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

Calibration of GRB Luminosity Relations with Cosmography

For the use of Gamma-Ray Bursts (GRBs) to probe cosmology in a cosmology-independent way, a new method has been proposed to obtain luminosity distances of GRBs by interpolating directly from the Hubble diagram of SNe Ia, and then calibrating GRB relations at high redshift. In this paper, following the basic assumption in the interpolation method that objects at the same redshift should have the same luminosity distance, we propose another approach to calibrate GRB luminosity relations with cosmographic fitting directly from SN Ia data. In cosmography, there is a well-known fitting formula which can reflect the Hubble relation between luminosity distance and redshift with cosmographic parameters which can be fitted from observation data. Using the Cosmographic fitting results from the Union set of SNe Ia, we calibrate five GRB relations using GRB sample at $z\leq1.4$ and deduce distance moduli of GRBs at $1.4< z \leq 6.6$ by generalizing above calibrated relations at high redshift. Finally, we constrain the dark energy parameterization models of the Chevallier-Polarski-Linder (CPL) model, the Jassal-Bagla-Padmanabhan (JBP) model and the Alam model with GRB data at high redshift, as well as with the Cosmic Microwave Background radiation (CMB) and the baryonic acoustic oscillation (BAO) observations, and we find the $Λ$CDM model is consistent with the current data in 1-$σ$ confidence region.