Source author record

Bing Zhang

Bing Zhang 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

222works
20topics
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

222 published item(s)

preprint2026arXiv

CangLing-KnowFlow: A Unified Knowledge-and-Flow-fused Agent for Comprehensive Remote Sensing Applications

The automated and intelligent processing of massive remote sensing (RS) datasets is critical in Earth observation (EO). Existing automated systems are normally task-specific, lacking a unified framework to manage diverse, end-to-end workflows--from data preprocessing to advanced interpretation--across diverse RS applications. To address this gap, this paper introduces CangLing-KnowFlow, a unified intelligent agent framework that integrates a Procedural Knowledge Base (PKB), Dynamic Workflow Adjustment, and an Evolutionary Memory Module. The PKB, comprising 1,008 expert-validated workflow cases across 162 practical RS tasks, guides planning and substantially reduces hallucinations common in general-purpose agents. During runtime failures, the Dynamic Workflow Adjustment autonomously diagnoses and replans recovery strategies, while the Evolutionary Memory Module continuously learns from these events, iteratively enhancing the agent's knowledge and performance. This synergy enables CangLing-KnowFlow to adapt, learn, and operate reliably across diverse, complex tasks. We evaluated CangLing-KnowFlow on the KnowFlow-Bench, a novel benchmark of 324 workflows inspired by real-world applications, testing its performance across 13 top Large Language Model (LLM) backbones, from open-source to commercial. Across all complex tasks, CangLing-KnowFlow surpassed the Reflexion baseline by at least 4% in Task Success Rate. As the first most comprehensive validation along this emerging field, this research demonstrates the great potential of CangLing-KnowFlow as a robust, efficient, and scalable automated solution for complex EO challenges by leveraging expert knowledge (Knowledge) into adaptive and verifiable procedures (Flow).

preprint2026arXiv

Multimessenger Emission from Very-High-Energy Black Hole-Jet Systems in the Milky Way

Microquasars, compact binary systems with an accreting stellar-mass black hole or neutron star, are promising candidates for high-energy particle acceleration. Recently, the LHAASO collaboration reported on the detection of $>100$ TeV $γ$-ray emission from five microquasars, suggesting that these sources are efficient particle accelerators. In microquasars, high-energy $γ$-rays can be produced in large-scale jets or winds. In this work, we explore the X-ray, $γ$-ray and neutrino emission from SS 433, V4641 Sgr and GRS 1905+105. We consider leptonic and hadronic scenarios to explain the spectra observed by LHAASO and other high-energy $γ$-ray detectors. We estimate the neutrino flux associated with the hadronic component and investigate the detectability of neutrinos from these sources in current and future neutrino telescopes. We find that among the three sources, V4641 Sgr has the best prospects of observation with a combined next-generation neutrino telescopes.

preprint2026arXiv

Runtime-Structured Task Decomposition for Agentic Coding Systems

Agentic coding systems increasingly use large language models (LLMs) for software engineering tasks such as debugging, root cause analysis, and code review. However, many existing systems encode task logic, execution flow, and output generation inside monolithic prompts. This design creates brittle behavior, limited debuggability, and high retry costs because failures often require rerunning the full workflow. We present runtime-structured task decomposition, an architectural approach in which task partitioning and execution flow are managed through executable control logic rather than prompt structure alone. LLMs are used only for focused judgment tasks, and outputs are validated against predefined schemas before downstream execution. We evaluate this approach on two software engineering workloads using three configurations: monolithic execution, static decomposition with fixed subtasks and no runtime branching, and runtime-structured decomposition. Each configuration was evaluated across 10 runs. Our results show that decomposition alone does not necessarily reduce retry cost. In the Kubernetes root cause analysis workload, the static decomposition baseline produced a retry cost of 1,632 +/- 145 tokens versus 904 +/- 17 tokens for the monolithic baseline because failures forced reruns of downstream subtasks. A similar pattern appeared in the multi-file debugging workload, where the static baseline consumed 933 tokens compared to 703 tokens for the monolithic system. The runtime-structured approach reran only failed subtasks, reducing retry costs to 436 +/- 132 tokens for root cause analysis and 460 tokens for debugging. Overall, the approach achieved up to 51.7% lower retry cost than monolithic systems and 73.2% lower retry cost than static decomposition baselines, improving efficiency, debuggability, and operational reliability in agentic coding systems.

preprint2026arXiv

The Advanced X-ray Imaging Satellite (AXIS) Community Science Book

The AXIS Community Science Book represents the collective effort of 592 scientists worldwide to define the transformative science enabled by the Advanced X-ray Imaging Satellite (AXIS), a next-generation X-ray mission selected by NASA's Astrophysics Probe Program for Phase A study. AXIS will advance the legacy of high-angular-resolution X-ray astronomy with ~1.5'' imaging over a wide 24' field of view and an order of magnitude greater collecting area than Chandra in the 0.3-12 keV band. Combining sharp imaging, high throughput, and rapid response capabilities, AXIS will open new windows on virtually every aspect of modern astrophysics, exploring the birth and growth of supermassive black holes, the feedback processes that shape galaxies, the life cycles of stars and exoplanet environments, and the nature of compact stellar remnants, supernova remnants, and explosive transients. This book compiles 138 community-contributed science cases developed by five Science Working Groups focused on AGN and supermassive black holes, galaxy evolution and feedback, compact objects and supernova remnants, stellar physics and exoplanets, and time-domain and multi-messenger astrophysics. Together, these studies establish the scientific foundation for next-generation X-ray exploration in the 2030s and highlight strong synergies with facilities of the 2030s, such as JWST, Roman, Rubin/LSST, SKA, ALMA, ngVLA, and next-generation gravitational-wave and neutrino networks.

preprint2025arXiv

Detection of disk-jet co-precession in a tidal disruption event

Theories and simulations predict that intense spacetime curvature near black holes bends the trajectories of light and matter, driving disk and jet precession under relativistic torques. However, direct observational evidence of disk-jet co-precession remains elusive. Here, we report the most compelling case to date: a tidal disruption event (TDE) exhibiting unprecedented 19.6-day quasi-periodic variations in both X-rays and radio, with X-ray amplitudes exceeding an order of magnitude. The nearly synchronized X-ray and radio variations suggest a shared mechanism regulating the emission regions. We demonstrate that a disk-jet Lense-Thirring precession model successfully reproduces these variations while requiring a low-spin black hole. This study uncovers previously uncharted short-term radio variability in TDEs, highlights the transformative potential of high-cadence radio monitoring, and offers profound insights into disk-jet physics.

preprint2024arXiv

The BOAT GRB 221009A: a Poynting-Flux-Dominated Narrow Jet Surrounded by a Matter-Dominated Structured Jet Wing

We argue that the broad-band observations of the brightest-of-all-time (BOAT) GRB 221009A reveal a physical picture involving two jet components: a narrow ($\sim 0.6$ degree half opening angle) pencil-beam jet that has a Poynting-flux-dominated jet composition, and a broader matter-dominated jet with an angular structure. We discuss various observational evidence that supports such a picture. To treat the problem, we develop an analytical structured jet model for both forward and reverse shock emission from the matter dominated structured jet wing during the deceleration phase. We discuss the physical implications of such a two-component jet configuration for this particular burst and for GRBs in general. We argue that some bright X-ray flares could be similar narrow jets viewed slightly outside the narrow jet cone and that narrow jets may exist in many more GRBs without being detected.

preprint2024arXiv

The Physics of Fast Radio Bursts

Fast radio bursts (FRBs), millisecond-duration bursts prevailing in the radio sky, are the latest big puzzle in the universe and have been a subject of intense observational and theoretical investigations in recent years. The rapid accumulation of the observational data has painted the following sketch about the physical origin of FRBs: They predominantly originate from cosmological distances so that their sources produce the most extreme coherent radio emission in the universe; at least some, probably most, FRBs are repeating sources that do not invoke cataclysmic events; and at least some FRBs are produced by magnetars, neutron stars with the strongest magnetic fields in the universe. Many open questions regarding the physical origin(s) and mechanism(s) of FRBs remain. This article reviews the phenomenology and possible underlying physics of FRBs. Topics include: a summary of the observational data, basic plasma physics, general constraints on FRB models from the data, radiation mechanisms, source and environment models, propagation effects, as well as FRBs as cosmological probes. Current pressing problems and future prospects are also discussed.

preprint2022arXiv

A Channel to Form Fast-spinning Black Hole--Neutron Star Binary Mergers as Multi-messenger Sources

After the successful detection of a gravitational-wave (GW) signal and its associated electromagnetic (EM) counterparts from GW170817, neutron star--black hole (NSBH) mergers have been highly expected to be the next type of multi-messenger source. However, despite the detection of several of NSBH merger candidates during the GW third observation run, no confirmed EM counterparts from these sources have been identified. The most plausible explanation is that these NSBH merger candidates were plunging events mainly because the primary BHs had near-zero projected aligned-spins based on GW observations. In view that NSs can be easily tidally disrupted by BHs with high projected aligned-spins, we study an evolution channel to form NSBH binaries with fast-spinning BHs, the properties of BH mass and spin, and their associated tidal disruption probability. We find that if the NSs are born firstly, the companion helium stars would be tidally spun up efficiently, and would thus finally form fast-spinning BHs. If BHs do not receive significant natal kicks at birth, these NSBH binaries that can merge within the Hubble time would have BHs with the projected aligned-spins $χ_{z}\gtrsim0.8$ and, hence, can certainly allow tidal disruption to happen. Even if significant BH kicks are considered for a small fraction of NSBH binaries, the projected aligned-spins of BHs are $χ_z\gtrsim0.2$. These systems can still be disrupted events unless the NSs are very massive. Thus, NS-first-born NSBH mergers would be promising multi-messenger sources. We discuss various potential EM counterparts associated with these systems and their detectability in the upcoming fourth observation run.

preprint2022arXiv

Damping of long wavelength gravitational waves by the intergalactic medium

The problem of radiation by the charged particles of the intergalactic medium (IGM) when a passing gravitational wave (GW) accelerate them is investigated. The largest acceleration (taking a charge from rest to a maximum speed which remains non-relativistic in the rest frame of the unperturbed spacetime) is found to be limited by the curvature of a propagating spherical gravitational wavefront. Interesting physics arises from the ensuing emission of radiation into the warm hot IGM, which to lowest order is a fully ionized hydrogen plasma with a frozen-in magnetic field $B$. It is found that for a vast majority of propagation directions, the right-handed polarized radiation can penetrate the plasma at frequencies below the plasma frequency $\om_p$, provided $\om<\om_b,$ where $\om_b=eB/m_e$ satisfies $\om_b<\om_p$ for typical IGM conditions. Moreover, the refractive index under such a scenario is $n\gg 1,$ resulting in an enhanced radiative dissipation of GW energy (relative to the vacuum scenario), which is more severe for electrons if both charge species are in thermal equilibrium and accelerated in the same way. The emission by the electrons then prevails, and is further amplified by coherent addition of amplitudes within the size one wavelength. The conversion of GWs of $\lam\gtrsim 5\times 10^{13}$~cm to electromagnetic waves means such GWs can only propagate a distance $\lesssim 1$~Gpc before being significantly damped by an IGM B field of $\sim10^{-8}$ G. The low-frequency GWs \textcolor{black}{targeted by pulsar-timing-arrays} will not survive unless the IGM magnetic field is much lower than expected. The \textcolor{black}{mHz} frequency GW inspirals targeted by future \textcolor{black}{space based} detectors such as the Laser Interferometer Space Antenna remain intact and can be detected.

preprint2022arXiv

Dark Fluxes from Accreting Black Holes and Direct Detections

We discuss the possibility that accreting black hole systems can be sources for dark matter flux through several different mechanisms. We firstly discuss two types of systems: coronal thermal plasmas around supermassive black holes in active galactic nuclei (AGNs), and accretion disks of stellar-mass X-ray black hole binaries (BHBs). We explore how these black hole systems may produce keV light dark matter fluxes and find that the dark fluxes from those sources might be too weak to account for the current XENON1T excess. On the other hand, black holes can be good accelerators to accrete and boost heavy dark matter particles. If considering collisions or dark electromagnetism, those particles can then escape and reach the benchmark speed of 0.1c at the detector. We also extend the black hole mass region to primordial black holes (PBHs) and discuss the possibility of contributing to keV light dark flux via superradiance of PBHs.

preprint2022arXiv

FAST observations of an extremely active episode of FRB 20201124A: II. Energy Distribution

We report the properties of more than 800 bursts detected from the repeating fast radio burst (FRB) source FRB 20201124A with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) during an extremely active episode on UTC September 25-28, 2021 in a series of four papers. In this second paper of the series, we mainly focus on the energy distribution of the detected bursts. The event rate initially increased exponentially but the source activity stopped within 24 hours after the 4th day. The detection of 542 bursts in one hour during the fourth day marked the highest event rate detected from one single FRB source so far. The bursts have complex structures in the time-frequency space. We find a double-peak distribution of the waiting time, which can be modeled with two log-normal functions peaking at 51.22 ms and 10.05 s, respectively. Compared with the emission from a previous active episode of the source detected with FAST, the second distribution peak time is smaller, suggesting that this peak is defined by the activity level of the source. We calculate the isotropic energy of the bursts using both a partial bandwidth and a full bandwidth and find that the energy distribution is not significantly changed. We find that an exponentially connected broken-power-law function can fit the cumulative burst energy distribution well, with the lower and higher-energy indices being $-1.22\pm0.01$ and $-4.27\pm0.23$, respectively. Assuming a radio radiative efficiency of $η_r = 10^{-4}$, the total isotropic energy of the bursts released during the four days when the source was active is already $3.9\times10^{46}$ erg, exceeding $\sim 23\%$ of the available magnetar dipolar magnetic energy. This challenges the magnetar models invoking an inefficient radio emission (e.g. synchrotron maser models).

preprint2022arXiv

FAST observations of an extremely active episode of FRB 20201124A: III. Polarimetry

As the third paper in the multiple-part series, we report the statistical properties of radio bursts detected from the repeating fast radio burst (FRB) source FRB 20201124A with the Five-hundred-meter Aperture Spherical radio telescope (FAST) during an extremely active episode between the 25th and the 28th of September 2021 (UT). We focus on the polarisation properties of 536 bright bursts with $\mathrm{S/N}>50$. We found that the Faraday rotation measures (RMs) monotonically dropped from $-579 \ {\rm rad \ m^{-2}}$ to $-605 \ {\rm rad \ m^{-2}}$ in the 4-day window. The RM values were compatible with the values ($-300$ to $-900\ {\rm rad \ m^{-2}}$ ) reported 4 month ago (Xu et al. 2022). However, the RM evolution rate in the current observation window was at least an order of magnitude smaller than the one ($\sim 500\ {\rm rad \ m^{-2}\, day^{-1}}$) previously reported during the rapid RM-variation phase, but is still higher than the one ($\le 1\ {\rm rad \ m^{-2} day^{-1}}$ ) during the later RM no-evolution phase. The bursts of FRB 20201124A were highly polarised with the total degree of polarisation (circular plus linear) greater than 90% for more than 90\% of all bursts. The distribution of linear polarisation position angles (PAs), degree of linear polarisation ($L/I$), and degree of circular polarisation ($V/I$) can be characterised with unimodal distribution functions. During the observation window, the distributions became wider with time, i.e. with larger scatter, but the centroids of the distribution functions remained nearly constant. For individual bursts, significant PA variations (confidence level 5-$σ$) were observed in 33% of all bursts. The polarisation of single pulses seems to follow certain complex trajectories on the Poincaré sphere, which may shed light on the radiation mechanism at the source or the plasma properties along the path of FRB propagation.

preprint2022arXiv

FAST observations of an extremely active episode of FRB 20201124A: IV. Spin Period Search

We report the properties of more than 800 bursts detected from the repeating fast radio burst (FRB) source FRB 20201124A with the Five-hundred-meter Aperture Spherical radio telescope (FAST) during an extremely active episode on UTC September 25th-28th, 2021 in a series of four papers. In this fourth paper of the series, we present a systematic search of the spin period and linear acceleration of the source object from both 996 individual pulse peaks and the dedispersed time series. No credible spin period was found from this data set. We rule out the presence of significant periodicity in the range between 1 ms to 100 s with a pulse duty cycle $< 0.49\pm0.08$ (when the profile is defined by a von-Mises function, not a boxcar function) and linear acceleration up to $300$ m s$^{-2}$ in each of the four one-hour observing sessions, and up to $0.6$ m s$^{-2}$ in all 4 days. These searches contest theoretical scenarios involving a 1 ms to 100 s isolated magnetar/pulsar with surface magnetic field $<10^{15}$ G and a small duty cycle (such as in a polar-cap emission mode) or a pulsar with a companion star or black hole up to 100 M$_{\rm \odot}$ and $P_b>10$ hours. We also perform a periodicity search of the fine structures and identify 53 unrelated millisecond-timescale "periods" in multi-components with the highest significance of 3.9 $σ$. The "periods" recovered from the fine structures are neither consistent nor harmonically related. Thus they are not likely to come from a spin period. We caution against claiming spin periodicity with significance below $\sim$ 4 $σ$ with multi-components from one-off FRBs. We discuss the implications of our results and the possible connections between FRB multi-components and pulsar micro-structures.

preprint2022arXiv

Frequency-dependent polarization of repeating fast radio bursts-implications for their origin

The polarization of fast radio bursts (FRBs), bright astronomical transients, contains crucial information about their environments. We report polarization measurements of five repeating FRBs, the abundant signals of which enable wide-band observations with two telescopes. A clear trend of lower polarization at lower frequencies was found, which can be well characterized by a single parameter rotation-measure-scatter (σRM) and modeled by multi-path scatter. Sources with higher σRM have higher RM magnitude and scattering timescales. The two sources with the most substantial σRM, FRB 20121102A and FRB 20190520B, are associated with a compact persistent radio source. These properties indicate a complex environment near the repeating FRBs, such as a supernova remnant or a pulsar wind nebula, consistent with their arising from young populations.

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

Limits on the hard X-ray emission from the periodic fast radio burst FRB 180916.J0158+65

FRB 180916.J0158+65 is one of the nearest, periodically repeating, and actively bursting fast radio burst (FRB) which has been localized to the outskirts of a spiral galaxy. In this work we study the FRB with the hard X-ray $14-195$ keV data from the Burst Alert Telescope (BAT) on board The Neil Gehrels Swift Observatory. BAT uses coded mask technology giving a localization of $\lesssim 3$ arc-minute in the hard X-ray band, along with an accurate background estimation. BAT has been observing the source location in survey mode since February 2020. The survey mode observations involves background subtracted spectra, integrated over a time span ranging $300-2000$ seconds, at the source location (from Feb 2020-Jan 2022). We analyzed all the $\sim 230$ survey mode observations from BAT and checked for any signal in any of the observations. We did not detect any signal at $>5σ$ confidence level in any of the observations. We could estimate a $5σ$ upper limit on the $14-195$ keV flux, which ranged between $4.5\times 10^{-10} - 7.6\times 10^{-9}\, \rm erg\, cm^{-2}\, s^{-1}$. At the source distance this relates to a $5σ$ upper limit on luminosity of $5.08\times 10^{44}- 8.5\times 10^{45} \rm erg\, s^{-1}$. With this estimate, we could rule out any persistent X-ray emission, at the source location for these snapshots of BAT observations.

preprint2022arXiv

Luminosity function and event rate density of XMM-Newton-selected supernova shock-breakout candidates

A dozen X-ray supernova shock breakout (SN SBO) candidates were reported recently based on XMM-Newton archival data, which increased the X-ray selected SN SBO sample by an order of magnitude. Assuming they are genuine SN SBOs, we study the luminosity function (LF) by improving upon the method used in our previous work. The light curves and the spectra of the candidates were used to derive the maximum volume within which these objects could be detected with XMM-Newton by simulation. The results show that the SN SBO LF can be described by either a broken power law (BPL) with indices (at the 68$\%$ confidence level) of $0.48 \pm 0.28$ and $2.11 \pm 1.27$ before and after the break luminosity at $\log (L_b/\rm erg\,s^{-1})=$ $45.32 \pm 0.55$ or a single power law (SPL) with index of $0.80 \pm 0.16$. The local event rate densities of SN SBOs above $5\times 10^{42}$ $\rm erg\,s^{-1}$ are consistent for two models, i.e., $4.6 ^{+1.7}_{-1.3} \times 10^4$ and $4.9 ^{+1.9}_{-1.4} \times 10^4$ $\rm Gpc^{-3}\,yr^{-1}$ for BPL and SPL models, respectively. The number of fast X-ray transients of SN SBO origin can be significantly increased by the wide-field X-ray telescopes such as the Einstein Probe.

preprint2022arXiv

Machine learning classification of CHIME fast radio bursts: II. Unsupervised Methods

Fast radio bursts (FRBs) are one of the most mysterious astronomical transients. Observationally, they can be classified into repeaters and apparently non-repeaters. However, due to the lack of continuous observations, some apparently repeaters may have been incorrectly recognized as non-repeaters. In a series of two papers, we intend to solve such problem with machine learning. In this second paper of the series, we focus on an array of unsupervised machine learning methods. We apply multiple unsupervised machine learning algorithms to the first CHIME/FRB catalog to learn their features and classify FRBs into different clusters without any premise about the FRBs being repeaters or non-repeaters. These clusters reveal the differences between repeaters and non-repeaters. Then, by comparing with the identities of the FRBs in the observed classes, we evaluate the performance of various algorithms and analyze the physical meaning behind the results. Finally, we recommend a list of most credible repeater candidates as targets for future observing campaigns to search for repeated bursts in combination of the results presented in Paper I using supervised machine learning methods.

preprint2022arXiv

Magnetar Engines in Fast Blue Optical Transients and Their Connections with SLSNe, SNe Ic-BL, and lGRBs

We fit the multi-band lightcurves of 40 fast blue optical transients (FBOTs) with the magnetar engine model. The mass of the FBOT ejecta, the initial spin period and polar magnetic field of the FBOT magnetars are respectively constrained to $M_{\rm{ej}}=0.18^{+0.52}_{-0.13}\,M_\odot$, $P_{\rm{i}}=9.4^{+8.1}_{-3.9}\,{\rm{ms}}$, and $B_{\rm p}=7^{+16}_{-5}\times10^{14}\,{\rm{G}}$. The wide distribution of the value of $B_{\rm p}$ spreads the parameter ranges of the magnetars from superluminous supernovae (SLSNe) to broad-line Type Ic supernovae (SNe Ic-BL; some are observed to be associated with long-duration gamma-ray bursts), which are also suggested to be driven by magnetars. Combining FBOTs with the other transients, we find a strong universal anti-correlation as $P_{\rm{i}}\propto{M_{\rm{ej}}^{-0.45}}$, indicating them could share a common origin. To be specific, it is suspected that all of these transients originate from collapse of extreme-stripped stars in close binary systems, but with different progenitor masses. As a result, FBOTs distinct themselves by their small ejecta masses with an upper limit of ${\sim}1\,M_\odot$, which leads to an observational separation in the rise time of the lightcurves $\sim12\,{\rm d}$. In addition, the FBOTs together with SLSNe can be separated from SNe Ic-BL by an empirical line in the $M_{\rm peak}-t_{\rm rise}$ plane corresponding to an energy requirement of a mass of $^{56}$Ni of $\sim0.3M_{\rm ej}$, where $M_{\rm peak}$ is the peak absolute magnitude of the transients and $t_{\rm rise}$ is the rise time.

preprint2022arXiv

Neutrino emission from FRB-emitting magnetars

The detection of a bright radio burst (hereafter FRB 200428) in association with a hard X-ray burst from the Galactic magnetar SGR 1935+2154 suggests that magnetars can make FRBs. We study possible neutrino emission from FRB-emitting magnetars by developing a general theoretical framework. We consider three different sites for proton acceleration and neutrino emission, i.e. within the magnetosphere, in the current sheet region beyond the light cylinder, and in relativistic shocks far away from the magnetosphere. Different cooling processes for protons and pions are considered to calculate the neutrino emission suppression factor within each scenario. We find that the flux of the neutrino emission decreases with increasing radius from the magnetar due to the decrease of the target photon number density. We calculate the neutrino flux from FRB 200428 and its associated X-ray burst. The flux of the most optimistic case invoking magnetospheric proton acceleration is still $\sim4$ orders of magnitude below the IceCube sensitivity. We also estimate the diffuse neutrino background from all FRB-emitting magnetars in the universe. The total neutrino flux of magnetars during their FRB emission phases is a negligible fraction of observed diffuse emission even under the most optimistic magnetospheric scenario for neutrino emission. However, if one assumes that many more X-ray bursts without FRB associations can also produce neutrinos with similar mechanisms, magnetars can contribute up to $10^{-8} \ {\rm GeV \ s^{-1} \ sr^{-1} \ cm^{-2}}$ diffuse neutrino background flux in the GeV to multi-TeV range.

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

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

Relativistic oblique shocks with ordered or random magnetic fields: tangential field governs

Relativistic magnetohydrodynamic shocks are efficient particle accelerators, often invoked in the models of gamma-ray bursts (GRBs) and shock-powered fast radio bursts (FRBs). Most theoretical studies assume a perpendicular shock with an ordered magnetic field perpendicular to the shock normal. However, the degree of magnetization $σ$ and the magnetic field geometry in shock-powered GRB/FRB scenarios are still poorly constrained by observations. Analogous to the magnetization $σ$ associated with the total field strength, we define a tangential magnetization $σ_\perp$ associated with the tangential field component. We explore the jump conditions of magnetized relativistic shocks, either with an ordered field of arbitrary inclination angle or with a random field of arbitrary anisotropy. In either case, we find that the jump conditions of relativistic shocks are governed by the tangential magnetization $σ_\perp$ instead of the total magnetization $σ$, insensitive to the inclination angles or the anisotropy of the pre-shock magnetic field. The approximated analytical solution developed in this work could serve as a quick check for numerical simulations and apply to theoretical studies of GRBs/FRBs with a more general field geometry.

preprint2022arXiv

Reverse shock forming condition for magnetized relativistic outflows: reconciling theories and simulations

Reverse shock (RS) emission can be used to probe the properties of the relativistic ejecta, especially the degree of magnetization $σ$, in gamma-ray burst (GRB) afterglows. However, there has been confusion in the literature regarding the physical condition for the RS formation, and the role of magnetic fields in the RS dynamics in the Poynting-flux-dominated regime is not fully understood. Exploiting the shock jump conditions, we characterize the properties of a magnetized RS. We compare the RS dynamics and forming conditions from different theories and numerical simulations, and reconcile the discrepancies among them. The strict RS forming condition is found to be $σ< σ_\mathrm{cr}=(8/3)γ_4^2(n_1/n_4)$, where $n_4$ and $n_1$ are the rest-frame number densities of the ejecta and the ambient medium, respectively, $γ_4$ is the bulk Lorentz factor, and $σ_\mathrm{cr}$ is the critical magnetization. Contrary to previous claims, we prove that this condition agrees with other theoretical and simulated results, which can be further applied to the setup and consistency check of future numerical simulations. Using this condition, we propose a characteristic radius for RS formation, and categorize the magnetized shell into three regimes: 'thick shell' (relativistic RS), 'thin shell' (trans-relativistic RS), and 'no RS' regimes. The critical magnetization $σ_\mathrm{cr}$ is generally below unity for thin shells, but can potentially reaches $\sim 100-1000$ in the 'thick shell' regime. Our results could be applied to the dynamical evolution of Poynting-flux-dominated ejecta, with potential applications to self-consistent lightcurve modelling of magnetized relativistic outflows.

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

SimPO: Simultaneous Prediction and Optimization

Many machine learning (ML) models are integrated within the context of a larger system as part of a key component for decision making processes. Concretely, predictive models are often employed in estimating the parameters for the input values that are utilized for optimization models as isolated processes. Traditionally, the predictive models are built first, then the model outputs are used to generate decision values separately. However, it is often the case that the prediction values that are trained independently of the optimization process produce sub-optimal solutions. In this paper, we propose a formulation for the Simultaneous Prediction and Optimization (SimPO) framework. This framework introduces the use of a joint weighted loss of a decision-driven predictive ML model and an optimization objective function, which is optimized end-to-end directly through gradient-based methods.

preprint2022arXiv

Simultaneous view of the FRB~180301 with FAST and NICER during a bursting phase

FRB180301 is one of the most actively repeating fast radio bursts (FRBs) which has shown polarization angle changes in its radio burst emission, an indication for their likely origin in the magnetosphere of a highly-magnetized neutron star. We carried out a multi-wavelength campaign with the FAST radio telescope and NICER X-ray observatory to investigate any possible X-ray emission temporally coincident with the bright radio bursts. The observations took place on 2021 March 4, 9 and 19. We detected five bright radio bursts with FAST, four of which were strictly simultaneous with the NICER observations. The peak flux-density of the radio bursts ranged between $28-105$ mJy, the burst fluence between $27-170$ mJy-ms, and the burst durations between $1.7-12.3$ ms. The radio bursts from FRB~180301 exhibited complex time domain structure, and sub-pulses were detected in individual bursts, with no significant circular polarisation. The linear degree of polarisation in L-band reduced significantly compared to the 2019 observations. We do not detect any X-ray emission in excess of the background during the 5ms, 10ms, 100ms, 1sec and 100sec time intervals at/around the radio-burst barycenter-corrected arrival times, at a $>5σ$ confidence level. The $5σ$ upper limits on the X-ray a) persistent flux is $<7.64\times 10^{-12}\, \rm erg\, cm^{-2}\, s^{-1}$ , equivalent to $L_{\rm X}<2.50 \times 10^{45} \rm erg\, s^{-1}$ and b) 5 ms fluence is $<2\times 10^{-11} \rm erg\, cm^{-2}$, at the radio burst regions. Using the $5$ ms X-ray fluence upper limit, we can estimate the radio efficiency $η_{R/X} \equiv L_{\rm Radio}/L_{\rm X-ray} \gtrsim 10^{-8}$. The derived upper limit on $η_{R/X}$ is consistent with both magnetospheric models and synchrotron maser models involving relativistic shocks.

preprint2022arXiv

Temporal Scattering, Depolarization, and Persistent Radio Emission from Magnetized Inhomogeneous Environments Near Repeating Fast Radio Burst Sources

Some repeating fast radio burst (FRB) sources exhibit complex polarization behaviors, including frequency-dependent depolarization, variation of rotation measure (RM), and oscillating spectral structures of polarized components. Very recently, Feng et al. (2022) reported that active repeaters exhibit conspicuous frequency-dependent depolarization and a strong correlation between RM scatter ($σ_{\rm RM}$) and the temporal scattering time ($τ_{\rm s}$), $σ_{\rm RM}\proptoτ_{\rm s}^{1.0\pm0.2}$, both of which can be well described by multi-path propagation through a magnetized inhomogeneous plasma screen. This observation strongly suggests that the temporal scattering and RM scatter originate from the same region. Besides, a particular finding of note in Feng et al. (2022) is that the FRBs with compact persistent radio sources (PRS) tend to have extreme $σ_{\rm RM}$. In this work, we focus on some theoretical predictions on the relations among temporal scattering, depolarization by RM scatter, and PRS contributed by the magnetized plasma environment close to a repeating FRB source. The behaviors of the RM scatter imply that the magnetized plasma environment is consistent with a supernova remnant or a pulsar wind nebula, and the predicted $σ_{\rm RM}$-$τ_{\rm s}$ relation is $σ_{\rm RM}\proptoτ_{\rm s}^{(0.54-0.83)}$ for different astrophysical scenarios. We further make a general discussion on PRS that does not depend on specific astrophysical scenarios. We show that the specific luminosity of a PRS should have a positive correlation with the RM contributed by the plasma screen. This is consistent with the observations of FRB 121102 and FRB 190520B.

preprint2022arXiv

Tensor Decompositions for Hyperspectral Data Processing in Remote Sensing: A Comprehensive Review

Owing to the rapid development of sensor technology, hyperspectral (HS) remote sensing (RS) imaging has provided a significant amount of spatial and spectral information for the observation and analysis of the Earth's surface at a distance of data acquisition devices, such as aircraft, spacecraft, and satellite. The recent advancement and even revolution of the HS RS technique offer opportunities to realize the full potential of various applications, while confronting new challenges for efficiently processing and analyzing the enormous HS acquisition data. Due to the maintenance of the 3-D HS inherent structure, tensor decomposition has aroused widespread concern and research in HS data processing tasks over the past decades. In this article, we aim at presenting a comprehensive overview of tensor decomposition, specifically contextualizing the five broad topics in HS data processing, and they are HS restoration, compressed sensing, anomaly detection, super-resolution, and spectral unmixing. For each topic, we elaborate on the remarkable achievements of tensor decomposition models for HS RS with a pivotal description of the existing methodologies and a representative exhibition on the experimental results. As a result, the remaining challenges of the follow-up research directions are outlined and discussed from the perspective of the real HS RS practices and tensor decomposition merged with advanced priors and even with deep neural networks. This article summarizes different tensor decomposition-based HS data processing methods and categorizes them into different classes from simple adoptions to complex combinations with other priors for the algorithm beginners. We also expect this survey can provide new investigations and development trends for the experienced researchers who understand tensor decomposition and HS RS to some extent.

preprint2022arXiv

The CHIME Fast Radio Burst Population Does Not Track the Star Formation History of the Universe

The redshift distribution of fast radio bursts (FRBs) is not well constrained. The association of the Galactic FRB 200428 with the young magnetar SGR 1935+2154 raises the working hypothesis that FRB sources track the star formation history of the universe. The discovery of FRB 20200120E in association with a globular cluster in the nearby galaxy M81, however, casts doubts on such an assumption. We apply the Monte Carlo method developed in a previous work to test different FRB redshift distribution models against the recently released first CHIME FRB catalog in terms of their distributions in specific fluence, external dispersion measure ($\rm DM_E$), and inferred isotropic energy. Our results clearly rule out the hypothesis that all FRBs track the star formation history of the universe. The hypothesis that all FRBs track the accumulated stars throughout history describes the data better but still cannot meet both the $\rm DM_E$ and the energy criteria. The data seem to be better modeled with either a redshift distribution model invoking a significant delay with respect to star formation or a hybrid model invoking both a dominant delayed population and a subdominant star formation population. We discuss the implications of this finding for FRB source models.

preprint2022arXiv

The Future of Gamma-Ray Experiments in the MeV-EeV Range

Gamma-rays, the most energetic photons, carry information from the far reaches of extragalactic space with minimal interaction or loss of information. They bring messages about particle acceleration in environments so extreme they cannot be reproduced on earth for a closer look. Gamma-ray astrophysics is so complementary with collider work that particle physicists and astroparticle physicists are often one in the same. Gamma-ray instruments, especially the Fermi Gamma-ray Space Telescope, have been pivotal in major multi-messenger discoveries over the past decade. There is presently a great deal of interest and scientific expertise available to push forward new technologies, to plan and build space- and ground-based gamma-ray facilities, and to build multi-messenger networks with gamma rays at their core. It is therefore concerning that before the community comes together for planning exercises again, much of that infrastructure could be lost to a lack of long-term planning for support of gamma-ray astrophysics. Gamma-rays with energies from the MeV to the EeV band are therefore central to multiwavelength and multi-messenger studies to everything from astroparticle physics with compact objects, to dark matter studies with diffuse large scale structure. These goals and new discoveries have generated a wave of new gamma-ray facility proposals and programs. This paper highlights new and proposed gamma-ray technologies and facilities that have each been designed to address specific needs in the measurement of extreme astrophysical sources that probe some of the most pressing questions in fundamental physics for the next decade. The proposed instrumentation would also address the priorities laid out in the recent Astro2020 Decadal Survey, a complementary study by the astrophysics community that provides opportunities also relevant to Snowmass.

preprint2021arXiv

"A more probable explanation" is still impossible to explain GN-z11-flash: in response to Steinhardt et al. (arXiv:2101.12738)

In Jiang et al. (2020), we reported a possible bright flash (hereafter GN-z11-flash) from a galaxy GN-z11 at z ~ 11. Recently, Steinhardt et al. (2021; arXiv:2101.12738) found 27 images with transient signals in Keck MOSFIRE archival data and claimed that GN-z11-flash was more likely from a moving object in our Solar system. We show that the Steinhardt et al.'s definition of the chance probability and their methodology of finding GN-z11-flash-like transients are problematic in several aspects. In particular, none of their transients is analogous to GN-z11-flash, and none of them is positionally coincident with a known object in their imaging data. In Jiang et al., we performed a comprehensive analysis of the origin of GN-z11-flash and ruled out, to the best of our knowledge, the possibility of known man-made objects or moving objects in the Solar system, based on all available information and our current understanding of these objects. Steinhardt et al. did not use such information and did not analyse the GN-z11-flash event itself. The majority of their transients are apparently low-Earth orbit satellites or aircrafts. Therefore, their analysis can neither prove nor disprove our results. Finally, we present a method to estimate the chance probability of finding GN-z11-flash-like transients in archival data. Based on this method and the archival data used by Steinhardt et al., we obtain a loose upper limit of the probability that actually support the original results of Jiang et al. (2020).

preprint2021arXiv

Analytical Solution of Magnetically Dominated Astrophysical Jets and Winds: Jet Launching, Acceleration, and Collimation

We present an analytical solution of a highly magnetized jet/wind flow. The left side of the general force-free jet/wind equation (the "pulsar" equation) is separated into a rotating and a nonrotating term. The two equations with either term can be solved analytically, and the two solutions match each other very well. Therefore, we obtain a general approximate solution of a magnetically dominated jet/wind, which covers from the nonrelativistic to relativistic regimes, with the drift velocity well matching the cold plasma velocity. The acceleration of a jet includes three stages. (1) The jet flow is located within the Alfvén critical surface (i.e. the light cylinder), has a nonrelativistic speed, and is dominated by toroidal motion. (2) The jet is beyond the Alfvén critical surface where the flow is dominated by poloidal motion and becomes relativistic. The total velocity in these two stages follows the same law $vΓ=ΩR$. (3) The evolution law is replaced by $vΓ\approx1/\left(θ\sqrt{2-ν}\right)$, where $θ$ is the half-opening angle of the jet and $0\leqν\leq2$ is a free parameter determined by the magnetic field configuration. This is because the earlier efficient acceleration finally breaks the causality connection between different parts in the jet, preventing a global solution. The jet has to carry local charges and currents to support an electromagnetic balance. This approximate solution is consistent with known theoretical results and numerical simulations, and it is more convenient to directly compare with observations. This theory may be used to constrain the spin of black holes in astrophysical jets.

preprint2021arXiv

Binary comb models for FRB 121102

The first repeating fast radio burst source, FRB 121102, is observed to emit bursts periodically. We show that FRB 121102 can be interpreted as an interacting neutron star binary system with an orbital period of 159 days. We develop a binary comb model by introducing an eccentricity in the orbit. Besides the original funnel mode of the binary comb model, which was applied to FRB 180916.J0158+65 by Ioka and Zhang 2020, we also identify two new modes of the binary comb model, i.e. the tau-crossing mode and the inverse funnel mode, and apply them to interpret FRB 121102. These new developments expand the applicable parameter space, allowing the companion star to be a massive star, a massive black hole, or a supermassive black hole, with the latter two having larger parameter spaces. These models are also consistent with other observations, such as the persistent bright radio counterpart associated with the source. We also argue that the observed frequency dependence of the active window does not disfavor the binary comb model, in contrast to recent claims, and propose two possible scenarios to interpret the data.

preprint2021arXiv

Bursts before Burst: A Comparative Study on FRB 200428-associated and FRB-absent X-ray Bursts from SGR J1935+2154

Accompanied by an X-ray burst, the fast radio burst (FRB) FRB 200428 was recently confirmed as originating from the Galactic magnetar soft gamma repeater (SGR) SGR J1935+2154. Just before and after FRB 200428 was detected, the Five-hundred-meter Aperture Spherical radio Telescope (FAST) had been monitoring SGR J1935+2154 for eight hours. From UTC 2020 April 27 23:55:00 to 2020 April 28 00:50:37, FAST detected no pulsed radio emission from SGR J1935+2154, while Fermi/Gamma-ray Burst Monitor registered 34 bursts in the X/soft $γ$-ray band, forming a unique sample of X-ray bursts in the absence of FRBs. After a comprehensive analysis on light curves, time-integrated, and time-resolved spectral properties of these FRB-absent X-ray bursts, we compare this sample with the FRB-associated X-ray burst detected by Insight-HXMT, INTEGRAL, and Konus-Wind. The FRB-associated burst distinguishes itself from other X-ray bursts by its nonthermal spectrum and a higher spectral peak energy, but otherwise is not atypical. We also compare the cumulative energy distribution of our X-ray burst sample with that of first repeating FRB source, FRB 121102, with the calibration of FRB 200428-X-ray burst association. We find a similarity between the two, offering indirect support of the magnetar origin of cosmological FRBs. The event rate density of magnetar bursts is about $\sim 150$ times higher than the FRB event rate density at the energy of FRB 200428. This again suggests that, if all FRBs originate from magnetars, only a small fraction of X-ray bursts are associated with FRBs.

preprint2021arXiv

Characterizing the FRB host galaxy population and its connection to transients in the local and extragalactic Universe

We present the localization and host galaxies of one repeating and two apparently non-repeating Fast Radio Bursts. FRB20180301A was detected and localized with the Karl G. Jansky Very Large Array to a star-forming galaxy at $z=0.3304$. FRB20191228A, and FRB20200906A were detected and localized by the Australian Square Kilometre Array Pathfinder to host galaxies at $z=0.2430$ and $z=0.3688$, respectively. We combine these with 13 other well-localized FRBs in the literature, and analyze the host galaxy properties. We find no significant differences in the host properties of repeating and apparently non-repeating FRBs. FRB hosts are moderately star-forming, with masses slightly offset from the star-forming main-sequence. Star formation and low-ionization nuclear emission-line region (LINER) emission are major sources of ionization in FRB host galaxies, with the former dominant in repeating FRB hosts. FRB hosts do not track stellar mass and star formation as seen in field galaxies (more than 95% confidence). FRBs are rare in massive red galaxies, suggesting that progenitor formation channels are not solely dominated by delayed channels which lag star formation by Gigayears. The global properties of FRB hosts are indistinguishable from core-collapse supernovae (CCSNe) and short gamma-ray bursts (SGRBs) hosts, and the spatial offset (from galaxy centers) of FRBs is mostly inconsistent with that of the Galactic neutron star population (95% confidence). The spatial offsets of FRBs (normalized to the galaxy effective radius) also differ from those of globular clusters (GCs) in late- and early-type galaxies with 95% confidence.

preprint2021arXiv

Coherent inverse Compton scattering by bunches in fast radio bursts

The extremely high brightness temperature of fast radio bursts (FRBs) requires that their emission mechanism must be "coherent", either through concerted particle emission by bunches or through an exponential growth of a plasma wave mode or radiation amplitude via certain maser mechanisms. The bunching mechanism has been mostly discussed within the context of curvature radiation or cyclotron/synchrotron radiation. Here we propose a family of model invoking coherent inverse Compton scattering (ICS) of bunched particles that may operate within or just outside of the magnetosphere of a flaring magnetar. Crustal oscillations during the flaring event may excite low-frequency electromagnetic waves near the magnetar surface. The X-mode of these waves could penetrate through the magnetosphere. Bunched relativistic particles in the charge starved region inside the magnetosphere or in the current sheet outside of the magnetosphere would upscatter these low-frequency waves to produce GHz emission to power FRBs. The ICS mechanism has a much larger emission power for individual electrons than curvature radiation. This greatly reduces the required degree of coherence in bunches, alleviating several criticisms to the bunching mechanism raised in the context of curvature radiation. The emission is $\sim 100\%$ linearly polarized (with the possibility of developing circular polarization) with a constant or varying polarization angle across each burst. The mechanism can account for a narrow-band spectrum and a frequency downdrifting pattern, as commonly observed in repeating FRBs.

preprint2021arXiv

Dissecting the Energy Budget of a Gamma-Ray Burst Fireball

The jet composition and radiative efficiency of GRBs are poorly constrained from the data. If the jet composition is matter-dominated (i.e. a fireball), the GRB prompt emission spectra would include a dominant thermal component originating from the fireball photosphere, and a non-thermal component presumably originating from internal shocks whose radii are greater than the photosphere radius. We propose a method to directly dissect the GRB fireball energy budget into three components and measure their values by combining the prompt emission and early afterglow data. The measured parameters include the initial dimensionless specific enthalpy density ($η$), bulk Lorentz factors at the photosphere radius ($Γ_{\rm ph}$) and before fireball deceleration ($Γ_0$), the amount of mass loading ($M$), as well as the GRB radiative efficiency ($η_γ$). All the parameters can be derived from the data for a GRB with a dominant thermal spectral component, a deceleration bump feature in the early afterglow lightcurve, and a measured redshift. The results only weakly depend on the density $n$ of the interstellar medium when the composition ${\cal Y}$ parameter (typically unity) is specified.

preprint2021arXiv

Magnetospheric curvature radiation by bunches as emission mechanism for repeating fast radio bursts

Coherent curvature radiation as the radiation mechanism for fast radio bursts (FRBs) has been discussed since FRBs were discovered. We study the spectral and polarization properties of repeating FRBs within the framework of coherent curvature radiation by charged bunches in the magnetosphere of a highly magnetized neutron star. The spectra can be generally characterized by multisegment broken power laws, and evolve as bunches move and the line of sight sweeps. Emitted waves are highly linear polarized if the line of sight is confined to the beam within an angle of $1/γ$, while circular polarized degree becomes strong for off-beam cases. The spectro-temporal pulse-to-pulse properties can be a natural consequence due to the magnetospheric geometry. We investigate the relationship between drift rate, central frequency and their temporal duration. The radius-to-frequency mapping is derived and simulated within the assumptions of both dipolar and quadrupolar magnetic configurations. The geometric results show that FRBs are emitted in field lines more curved than open field lines for a dipolar geometry. This suggests that there are most likely existing multipolar magnetic configurations in the emission region.

preprint2021arXiv

MetricUNet: Synergistic Image- and Voxel-Level Learning for Precise CT Prostate Segmentation via Online Sampling

Fully convolutional networks (FCNs), including UNet and VNet, are widely-used network architectures for semantic segmentation in recent studies. However, conventional FCN is typically trained by the cross-entropy or Dice loss, which only calculates the error between predictions and ground-truth labels for pixels individually. This often results in non-smooth neighborhoods in the predicted segmentation. To address this problem, we propose a two-stage framework, with the first stage to quickly localize the prostate region and the second stage to precisely segment the prostate by a multi-task UNet architecture. We introduce a novel online metric learning module through voxel-wise sampling in the multi-task network. Therefore, the proposed network has a dual-branch architecture that tackles two tasks: 1) a segmentation sub-network aiming to generate the prostate segmentation, and 2) a voxel-metric learning sub-network aiming to improve the quality of the learned feature space supervised by a metric loss. Specifically, the voxel-metric learning sub-network samples tuples (including triplets and pairs) in voxel-level through the intermediate feature maps. Unlike conventional deep metric learning methods that generate triplets or pairs in image-level before the training phase, our proposed voxel-wise tuples are sampled in an online manner and operated in an end-to-end fashion via multi-task learning. To evaluate the proposed method, we implement extensive experiments on a real CT image dataset consisting of 339 patients. The ablation studies show that our method can effectively learn more representative voxel-level features compared with the conventional learning methods with cross-entropy or Dice loss. And the comparisons show that the proposed method outperforms the state-of-the-art methods by a reasonable margin.

preprint2021arXiv

SpectralFormer: Rethinking Hyperspectral Image Classification with Transformers

Hyperspectral (HS) images are characterized by approximately contiguous spectral information, enabling the fine identification of materials by capturing subtle spectral discrepancies. Owing to their excellent locally contextual modeling ability, convolutional neural networks (CNNs) have been proven to be a powerful feature extractor in HS image classification. However, CNNs fail to mine and represent the sequence attributes of spectral signatures well due to the limitations of their inherent network backbone. To solve this issue, we rethink HS image classification from a sequential perspective with transformers, and propose a novel backbone network called \ul{SpectralFormer}. Beyond band-wise representations in classic transformers, SpectralFormer is capable of learning spectrally local sequence information from neighboring bands of HS images, yielding group-wise spectral embeddings. More significantly, to reduce the possibility of losing valuable information in the layer-wise propagation process, we devise a cross-layer skip connection to convey memory-like components from shallow to deep layers by adaptively learning to fuse "soft" residuals across layers. It is worth noting that the proposed SpectralFormer is a highly flexible backbone network, which can be applicable to both pixel- and patch-wise inputs. We evaluate the classification performance of the proposed SpectralFormer on three HS datasets by conducting extensive experiments, showing the superiority over classic transformers and achieving a significant improvement in comparison with state-of-the-art backbone networks. The codes of this work will be available at https://github.com/danfenghong/IEEE_TGRS_SpectralFormer for the sake of reproducibility.

preprint2021arXiv

Using Low-rank Representation of Abundance Maps and Nonnegative Tensor Factorization for Hyperspectral Nonlinear Unmixing

Tensor-based methods have been widely studied to attack inverse problems in hyperspectral imaging since a hyperspectral image (HSI) cube can be naturally represented as a third-order tensor, which can perfectly retain the spatial information in the image. In this article, we extend the linear tensor method to the nonlinear tensor method and propose a nonlinear low-rank tensor unmixing algorithm to solve the generalized bilinear model (GBM). Specifically, the linear and nonlinear parts of the GBM can both be expressed as tensors. Furthermore, the low-rank structures of abundance maps and nonlinear interaction abundance maps are exploited by minimizing their nuclear norm, thus taking full advantage of the high spatial correlation in HSIs. Synthetic and real-data experiments show that the low rank of abundance maps and nonlinear interaction abundance maps exploited in our method can improve the performance of the nonlinear unmixing. A MATLAB demo of this work will be available at https://github.com/LinaZhuang for the sake of reproducibility.

preprint2020arXiv

"Slow" Radio Bursts from Galactic Magnetars?

Recently, one fast radio burst, FRB 200428, was detected from the Galactic magnetar SGR J1935+2154 during one X-ray burst. This suggests that magnetars can make FRBs. On the other hand, the majority of X-ray bursts from SGR J1935+2154 are not associated with FRBs. One possible reason for such rarity of FRB-SGR-burst associations is that the FRB emission is much more narrowly beamed than the SGR burst emission. If such an interpretation is correct, one would expect to detect radio bursts with viewing angles somewhat outside the narrow emission beam. These "slow" radio bursts (SRBs) would have broader widths and lower flux densities due to the smaller Doppler factor involved. We derive two "closure relations" to judge whether a long, less luminous radio burst could be an SRB. The 2.2-s, 308 Jy ms, 111 MHz radio burst detected from SGR J1935+2154 by the BSA LPI radio telescope may be such an SRB. The 2-ms, 60 mJy ms faint burst detected by FAST from the same source could be also an SRB if the corresponding FRB has a narrow spectrum. If the FRB beam is narrow, there should be many more SRBs than FRBs from Galactic magnetars. The lack of detection of abundant SRBs from magnetars would disfavor the hypothesis that all SGR bursts are associated with narrow-beam FRBs.

preprint2020arXiv

A Binary Comb Model for Periodic Fast Radio Bursts

We show that the periodic FRB 180916.J0158+65 can be interpreted by invoking an interacting neutron star binary system with an orbital period of $\sim 16$ days. The FRBs are produced by a highly magnetized pulsar, whose magnetic field is ``combed'' by the strong wind from a companion star, either a massive star or a millisecond pulsar. The FRB pulsar wind retains a clear funnel in the companion's wind that is otherwise opaque to induced Compton or Raman scatterings for repeating FRB emission. The 4 day active window corresponds to the time when the funnel points toward Earth. The interaction also perturbs the magnetosphere of the FRB pulsar and may trigger emission of FRBs. We derive the physical constraints on the comb and the FRB pulsar from the observations and estimate the event rate of FRBs. In this scenario, a lower limit on the period of observable FRBs is predicted. We speculate that both the intrinsic factors (strong magnetic field and young age) and the extrinsic factor (interaction) may be needed to generate FRBs in neutron star binary systems.

preprint2020arXiv

A comparative study of host galaxy properties between Fast Radio Bursts and stellar transients

Recent arcsecond localizations of Fast Radio Bursts and identifications of their host galaxies confirmed their extragalactic origin.While FRB 121102 resides in the bright region of a dwarf star forming galaxy, other FRBs reside in more massive galaxies and are related to older stellar populations. We compare the host galaxy properties of {nine} FRBs with those of several types of stellar transients: from young to old population, long duration gamma ray bursts (LGRBs), superluminous supernovae (SLSNe), Type Ib/Ic supernovae (SN Ibc), Type II supernovae (SN II), type Ia supernovae (SN Ia), and short duration gamma ray bursts (SGRBs). We find that as a whole sample, the stellar mass and star formation rate of the FRB host galaxies prefer a medium to old population, and are against a young population similar to LGRBs and SLSNe by a null probability 0.02. Individually, the FRB 121102 host is consistent with that of young population objects; the FRB 180924 environment is similar to that of SGRBs; and the FRB 190523 environment is similar to those of SN Ia. These results are consistent with the magnetar engine model for FRBs, if both magnetars produced from extreme explosions (GRBs/SLSNe) and from regular channels (e.g. those producing Galactic magnetars) can produce FRBs.

preprint2020arXiv

A Comparative Study of Long and Short GRBs. II. A Multi-wavelength Method to distinguish Type II (massive star) and Type I (compact star) GRBs

Gamma Ray Burst (GRBs) are empirically classified as long-duration GRBs (LGRBs, $>$ 2s) and short-duration GRBs (SGRBs, $<$ 2s). Physically they can be grouped into two distinct progenitor categories: those originating from collapse of massive stars (also known as Type II) and those related to mergers of compact stars (also known as Type I). Even though most LGRBs are Type II and most SGRBs are Type I, the duration criterion is not always reliable to determine the physical category of a certain GRB. Based on our previous comprehensive study of the multi-wavelength properties of long and short GRBs, here we utilize the Naive Bayes method to physically classify GRBs as Type I and Type II GRBs based on multi-wavelength criteria. It results in 0.5\% training error rate and 1\% test error rate. Moreover, there is a gap [-1.2, -0.16] in the distribution of the posterior Odds, $\log O({\rm II:I})$, the Type II to Type I probability ratio. Therefore, we propose to use ${\cal O} = \log O({\rm II:I})+0.7$ as the parameter to classify GRBs into Type I ($<0$) or Type II ($>0$). The only confirmed Type I GRB, GRB 170817A, has log $O({\rm II:I})=-10$. According to this criterion, the supernova-less long GRBs 060614 and 060505 belong to Type I, and two controversial short GRBs 090426 and 060121 belong to Type II.

preprint2020arXiv

A Fast Radio Burst discovered in FAST drift scan survey

We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs.

preprint2020arXiv

A serendipitous discovery of GeV gamma-ray emission from supernova 2004dj in a survey of nearby star-forming galaxies with Fermi-LAT

The interaction between a supernova ejecta and the circum-stellar medium drives a strong shock wave which accelerates particles (i.e., electrons and protons). The radio and X-ray emission observed after the supernova explosion constitutes the evidence of the electron acceleration. The accelerated protons are expected to produce GeV-TeV gamma-ray emission via $pp$ collisions, but the flux is usually low since only a small fraction of the supernova kinetic energy is converted into the shock energy at the very early time. The low gamma-ray flux of the nearest supernova explosion, SN 1987A, agrees with this picture. Here we report a serendipitous discovery of a fading GeV gamma-ray source in spatial coincidence with the second nearest supernova--SN 2004dj from our gamma-ray survey of nearby star-forming galaxies with Fermi-LAT. The total gamma-ray energy released by SN 2004dj is about $6\times10^{47}{\rm erg}$. We interpret this gamma-ray emission arising from the supernova ejecta interacting with a surrounding high-density shell, which decelerates the ejecta and converts ~1% of the ejecta's kinetic energy to relativistic protons. In addition, our gamma-ray survey of nearby star-forming galaxies discovers GeV emissions from two star-forming galaxies, i.e., Arp 299 and M33, for the first time.

preprint2020arXiv

A unified picture of Galactic and cosmological fast radio bursts

The discovery of a fast radio burst (FRB) in our galaxy associated with a magnetar (neutron star with strong magnetic field) has provided a critical piece of information to help us finally understand these enigmatic transients. We show that the volumetric rate of Galactic-FRB like events is consistent with the faint end of the cosmological FRB rate, and hence they most likely belong to the same class of transients. The Galactic FRB had an accompanying X-ray burst but many X-ray bursts from the same object had no radio counterpart. Their relative rates suggest that for every FRB there are roughly 10^2 to 10^3 X-ray bursts. The radio lightcurve of the galactic FRB had two spikes separated by 30 ms in the 400-800 MHz frequency band. This is an important clue and highly constraining of the class of models where the radio emission is produced outside the light-cylinder of the magnetar. We suggest that magnetic disturbances close to the magnetar surface propagate to a distance of a few tens of neutron star radii where they damp and produce radio emission. The coincident hard X-ray spikes associated with the two FRB pulses seen in this burst and the flux ratio between the two frequency bands can be understood in this scenario. This model provides a unified picture for faint bursts like the Galactic FRB as well as the bright events seen at cosmological distances.

preprint2020arXiv

Are persistent emission luminosity and rotation measure of fast radio bursts related?

The physical origin of fast radio bursts (FRBs) is still unknown. Multiwavelength and polarization observations of an FRB source would be helpful to diagnose its progenitor and environment. So far only the first repeating source FRB 121102 appears to be spatially coincident with a persistent radio emission. Its bursts also have very large values of the Faraday rotation measure (RM) i.e., $|{\rm RM}|\sim10^5~{\rm rad~m^{-2}}$. We show that theoretically there should be a simple relation between RM and the luminosity of the persistent source of an FRB source if the observed RM mostly arises from the persistent emission region. FRB 121102 follows this relation given that the magnetic field in the persistent emission region is highly ordered and that the number of relativistic electrons powering the persistent emission is comparable to that of non-relativistic electrons that contribute to RM. The non-detections of persistent emission sources from all other localized FRB sources are consistent with their relatively small RMs ($\left|{\rm RM}\right|\lesssim{\rm a~few}\times100~{\rm rad~m^{-2}}$) according to this relation. Based on this picture, the majority of FRBs without a large RM are not supposed to be associated with bright persistent sources.

preprint2020arXiv

Automated Radiological Report Generation For Chest X-Rays With Weakly-Supervised End-to-End Deep Learning

The chest X-Ray (CXR) is the one of the most common clinical exam used to diagnose thoracic diseases and abnormalities. The volume of CXR scans generated daily in hospitals is huge. Therefore, an automated diagnosis system able to save the effort of doctors is of great value. At present, the applications of artificial intelligence in CXR diagnosis usually use pattern recognition to classify the scans. However, such methods rely on labeled databases, which are costly and usually have large error rates. In this work, we built a database containing more than 12,000 CXR scans and radiological reports, and developed a model based on deep convolutional neural network and recurrent network with attention mechanism. The model learns features from the CXR scans and the associated raw radiological reports directly; no additional labeling of the scans are needed. The model provides automated recognition of given scans and generation of reports. The quality of the generated reports was evaluated with both the CIDEr scores and by radiologists as well. The CIDEr scores are found to be around 5.8 on average for the testing dataset. Further blind evaluation suggested a comparable performance against human radiologist.

preprint2020arXiv

Blazar -- IceCube neutrino association revisited

The reported association of high-energy neutrino event IceCube-170922A and blazar TXS 0506+056 has sparkled discussion about blazars as sources of cosmic neutrinos. In this paper, we use publicly released IceCube data and blazar locations from Roma-BzCat to test spacial correlation between neutrino events and blazar locations. We also scrutinize the correlation between $γ$-ray flux and neutrino flux of blazars by applying a temporal filter onto the data based on Fermi monitored source list. We find no compelling evidence to prove blazars as the main source of cosmic neutrinos, as known before the detected IceCube-170922A / TXS 0506+056 association. While we do not rule out the association between IceCube-170922A and TXS 0506+056, the significance level we obtained is not high enough to claim a strong association. If such an association is real, a special physical condition is desired to allow a small fraction of blazars to become bright neutrino sources.

preprint2020arXiv

Charged Compact Binary Coalescence Signal and Electromagnetic Counterpart of Plunging BH-NS mergers

If at least one of the members of a compact binary coalescence is charged, the inspiral of the two members would generate a Poynting flux with an increasing power, giving rise to a brief electromagnetic counterpart temporally associated with the chirp signal of the merger (with possibly a small temporal offset), which we term as the {\em charged Compact Binary Coalescence} (cCBC) signal. We develop a general theory of cCBC for any mass and amount of charge for each member. Neutron stars (NSs), as spinning magnets, are guaranteed to be charged, so the cCBC signal should accompany all neutron star mergers. The cCBC signal is clean in a BH-NS merger with a small mass ratio ($q \equiv m_2/m_1 < 0.2$), in which the NS plunges into the BH as a whole, and its luminosity/energy can reach that of a fast radio burst if the NS is Crab-like. The strength of the cCBC signal in Extreme Mass Ratio Inspiral Systems (EMRIs) is also estimated.

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

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

Coupled Convolutional Neural Network with Adaptive Response Function Learning for Unsupervised Hyperspectral Super-Resolution

Due to the limitations of hyperspectral imaging systems, hyperspectral imagery (HSI) often suffers from poor spatial resolution, thus hampering many applications of the imagery. Hyperspectral super-resolution refers to fusing HSI and MSI to generate an image with both high spatial and high spectral resolutions. Recently, several new methods have been proposed to solve this fusion problem, and most of these methods assume that the prior information of the Point Spread Function (PSF) and Spectral Response Function (SRF) are known. However, in practice, this information is often limited or unavailable. In this work, an unsupervised deep learning-based fusion method - HyCoNet - that can solve the problems in HSI-MSI fusion without the prior PSF and SRF information is proposed. HyCoNet consists of three coupled autoencoder nets in which the HSI and MSI are unmixed into endmembers and abundances based on the linear unmixing model. Two special convolutional layers are designed to act as a bridge that coordinates with the three autoencoder nets, and the PSF and SRF parameters are learned adaptively in the two convolution layers during the training process. Furthermore, driven by the joint loss function, the proposed method is straightforward and easily implemented in an end-to-end training manner. The experiments performed in the study demonstrate that the proposed method performs well and produces robust results for different datasets and arbitrary PSFs and SRFs.

preprint2020arXiv

Fast radio bursts as strong waves interacting with ambient medium

Fast radio bursts (FRBs) are mysterious radio transients whose physical origin is still unknown. Within a few astronomical units near an FRB source, the electric field of the electromagnetic wave is so large that the electron oscillation velocity becomes relativistic, which makes the classical Thomson scattering theory and the linear plasma theory invalid. We discuss FRBs as strong waves interacting with the ambient medium, in terms of both electron motion properties and plasma properties. Several novel features are identified. 1. The cross section of Thomson scattering is significantly enhanced for the scattering photons. 2. On the other hand, because of the nonlinear plasma properties in strong waves, the near-source plasma is more transparent and has a smaller effective dispersion measure (DM) contribution to the observed value. For a repeating FRB source, the brighter bursts would have somewhat smaller DMs contributed by the near-source plasma. 3. The radiation beam undergoes relativistic self-focusing in a dense plasma, the degree of self-focusing (or squeezing) depends on the plasma density. Such a squeezing effect would affect the collimation angle and the true event rate of FRBs. 4. When an FRB propagates in a nearby ambient plasma, a wakefield wave in the plasma will be generated by the ponderomotive force of the FRB and accelerates electrons in the ambient medium. However, such an effect is too weak to be observationally interesting.

preprint2020arXiv

Fast Radio Bursts from Interacting Binary Neutron Star Systems

Recent observations of repeating fast radio bursts (FRBs) suggest that some FRBs reside in an environment consistent with that of binary neutron star (BNS) mergers. The bursting rate for repeaters could be very high and the emission site is likely from a magnetosphere. We discuss a hypothesis of producing abundant repeating FRBs in BNS systems. Decades to centuries before a BNS system coalesces, the magnetospheres of the two neutron stars start to interact relentlessly. Abrupt magnetic reconnection accelerates particles, which emit coherent radio waves in bunches via curvature radiation. FRBs are detected as these bright radiation beams point towards Earth. This model predicts quasi-periodicity of the bursts at the rotation periods of the two merging neutron stars (tens of milliseconds and seconds, respectively) as well as the period of orbital motion (of the order of 100 s). The bursting activities are expected to elevate with time as the two neutron stars get closer. The repeating FRB sources should be gravitational wave (GW) sources for space-borne detectors such as LISA if the sources are close enough, and will be eventually detected by ground-based detectors when the two neutron stars coalesce.

preprint2020arXiv

Finite Interface Dissipation Phase Field Modeling of Ni-Nb under Additive Manufacturing Conditions

During the laser powder bed fusion (L-PBF) process, the built part undergoes multiple rapid heating-cooling cycles, leading to complex microstructures with nonuniform properties. In the present work, a computational framework, which weakly couples a finite element thermal model to a non-equilibrium PF model was developed to investigate the rapid solidification microstructure of a Ni-Nb alloy during L-PBF. The framework is utilized to predict the spatial variation of the morphology and size of cellular segregation structure as well as the microsegregation in single-track melt pool microstructures obtained under different process conditions. A solidification map demonstrating the variation of microstructural features as a function of the temperature gradient and growth rate is presented. A planar to cellular transition is predicted in the majority of keyhole mode melt pools, while a planar interface is predominant in conduction mode melt pools. The predicted morphology and size of the cellular segregation structure agrees well with experimental measurements.

preprint2020arXiv

Growth of stellar mass black holes in dense molecular clouds and GW190521

A stellar mass black hole can grow its mass noticeably through Bondi accretion, if it is embedded in an extremely dense and massive molecular cloud with slow motion with respect to the ambient medium for an extended period of time. This provides a novel, yet challenging channel for the formation of massive stellar-mass black holes. We discuss how this channel may account for the massive binary black hole merger system GW190521 as observed by LIGO/Virgo gravitational wave detectors as well as the claimed massive black hole candidate LB-1.

preprint2020arXiv

HF-UNet: Learning Hierarchically Inter-Task Relevance in Multi-Task U-Net for Accurate Prostate Segmentation

Accurate segmentation of the prostate is a key step in external beam radiation therapy treatments. In this paper, we tackle the challenging task of prostate segmentation in CT images by a two-stage network with 1) the first stage to fast localize, and 2) the second stage to accurately segment the prostate. To precisely segment the prostate in the second stage, we formulate prostate segmentation into a multi-task learning framework, which includes a main task to segment the prostate, and an auxiliary task to delineate the prostate boundary. Here, the second task is applied to provide additional guidance of unclear prostate boundary in CT images. Besides, the conventional multi-task deep networks typically share most of the parameters (i.e., feature representations) across all tasks, which may limit their data fitting ability, as the specificities of different tasks are inevitably ignored. By contrast, we solve them by a hierarchically-fused U-Net structure, namely HF-UNet. The HF-UNet has two complementary branches for two tasks, with the novel proposed attention-based task consistency learning block to communicate at each level between the two decoding branches. Therefore, HF-UNet endows the ability to learn hierarchically the shared representations for different tasks, and preserve the specificities of learned representations for different tasks simultaneously. We did extensive evaluations of the proposed method on a large planning CT image dataset, including images acquired from 339 patients. The experimental results show HF-UNet outperforms the conventional multi-task network architectures and the state-of-the-art methods.

preprint2020arXiv

Is GRB 110715A the progenitor of FRB 171209?

The physical origin of fast radio bursts (FRBs) is unknown. Young magnetars born from gamma-ray bursts (GRBs) have been suggested to be a possible central engine of FRBs. We test such a hypothesis by systematically searching for GRB-FRB spatial associations from 110 FRBs and 1440 GRBs. We find that one FRB event, FRB 171209, discovered by the Parkes telescope is spatially coincident with a historical long-duration GRB 110715A at $z=0.82$. The afterglow of GRB 110715A is consistent with being powered by a millisecond magnetar. The extragalactic dispersion measure of FRB 171209 is in excess of that contributed by the intergalactic medium, which can be interpreted as being contributed by a young supernova remnant associated with the GRB. Overall, the significance of the association is $(2.28 - 2.55) σ$. If the association is indeed physical, our result suggests that the magnetars associated with long GRBs can be the progenitors of at least some FRBs.

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

More Diverse Means Better: Multimodal Deep Learning Meets Remote Sensing Imagery Classification

Classification and identification of the materials lying over or beneath the Earth's surface have long been a fundamental but challenging research topic in geoscience and remote sensing (RS) and have garnered a growing concern owing to the recent advancements of deep learning techniques. Although deep networks have been successfully applied in single-modality-dominated classification tasks, yet their performance inevitably meets the bottleneck in complex scenes that need to be finely classified, due to the limitation of information diversity. In this work, we provide a baseline solution to the aforementioned difficulty by developing a general multimodal deep learning (MDL) framework. In particular, we also investigate a special case of multi-modality learning (MML) -- cross-modality learning (CML) that exists widely in RS image classification applications. By focusing on "what", "where", and "how" to fuse, we show different fusion strategies as well as how to train deep networks and build the network architecture. Specifically, five fusion architectures are introduced and developed, further being unified in our MDL framework. More significantly, our framework is not only limited to pixel-wise classification tasks but also applicable to spatial information modeling with convolutional neural networks (CNNs). To validate the effectiveness and superiority of the MDL framework, extensive experiments related to the settings of MML and CML are conducted on two different multimodal RS datasets. Furthermore, the codes and datasets will be available at https://github.com/danfenghong/IEEE_TGRS_MDL-RS, contributing to the RS community.

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

Nonuniversal interstellar density spectra probed by pulsars

The Galactic interstellar turbulence affects the density distribution and star formation. We introduce a new method of measuring interstellar turbulent density spectra by using the dispersion measures (DMs) of a large sample of pulsars. Without the need of invoking multiple tracers, we obtain nonuniversal density spectra in the multi-phase interstellar medium over different ranges of length scales. By comparing the analytical structure function of DMs with the observationally measured one in different areas of sky, we find a shallow density spectrum arising from the supersonic turbulence in cold interstellar phases, and a Kolmogorov-like density spectrum in the diffuse warm ionized medium (WIM). Both spectra extend up to hundreds of pc. On larger scales, we for the first time identify a steep density spectrum in the diffuse WIM extending up to several kpc. Our results show that the DMs of pulsars can provide unique new information on the interstellar turbulence.

preprint2020arXiv

On the FRB luminosity function -- II. Event rate density

The luminosity function of Fast Radio Bursts (FRBs), defined as the event rate per unit cosmic co-moving volume per unit luminosity, may help to reveal the possible origins of FRBs and design the optimal searching strategy. With the Bayesian modelling, we measure the FRB luminosity function using 46 known FRBs. Our Bayesian framework self-consistently models the selection effects, including the survey sensitivity, the telescope beam response, and the electron distributions from Milky Way / the host galaxy / local environment of FRBs. Different from the previous companion paper, we pay attention to the FRB event rate density and model the event counts of FRB surveys based on the Poisson statistics. Assuming a Schechter luminosity function form, we infer (at the 95% confidence level) that the characteristic FRB event rate density at the upper cut-off luminosity $L^*=2.9_{-1.7}^{+11.9}\times10^{44}\,\rm erg\, s^{-1}$ is $ϕ^*=339_{-313}^{+1074}\,\rm Gpc^{-3}\, yr^{-1}$, the power-law index is $α=-1.79_{-0.35}^{+0.31}$, and the lower cut-off luminosity is $L_0\le9.1\times10^{41}\,\rm erg\, s^{-1}$. The event rate density of FRBs is found to be $3.5_{-2.4}^{+5.7}\times10^4\,\rm Gpc^{-3}\, yr^{-1}$ above $10^{42}\,\rm erg\, s^{-1}$, $5.0_{-2.3}^{+3.2}\times10^3\,\rm Gpc^{-3}\, yr^{-1}$ above $10^{43}\,\rm erg\, s^{-1}$, and $3.7_{-2.0}^{+3.5}\times10^2\,\rm Gpc^{-3}\, yr^{-1}$ above $10^{44}\,\rm erg\, s^{-1}$. As a result, we find that, for searches conducted at 1.4 GHz, the optimal diameter of single-dish radio telescopes to detect FRBs is 30-40 m. The possible astrophysical implications of the measured event rate density are also discussed in the current paper.

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

Peierls-type metal-insulator transition in carbon nanostructures

We report the observation of Peierls-type metal-insulator transition in carbon nanostructures formed by chemical vapor deposition inside the pore network of the ZSM-5 zeolite. The Raman spectrum of this nanocarbon@ZSM-5 indicates a clear signature of the radial breathing mode (RBM) for (3,0) carbon nanotubes that can constitute the carbon network segments. Electrical transport measurements on multiple few-micron-sized nanocarbon@ZSM-5 crystals showed metallic temperature of resistance dependence down to 30 K, at which point the resistance exhibited a sharp upturn that is accompanied by the opening of a quasigap at the Fermi level as indicated by the differential resistance measurements. Further Hall measurements have yielded both the sign of the charge carrier and its density. The latter demonstrated excellent consistency with the quasigap data. We employed first-principles calculations to verify that there can indeed be softening of the phonon modes in the (3,0) carbon nanotubes.

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

Probing the intergalactic turbulence with fast radio bursts

The turbulence in the diffuse intergalactic medium (IGM) plays an important role in various astrophysical processes across cosmic time, but it is very challenging to constrain its statistical properties both observationally and numerically. Via the statistical analysis of turbulence along different sightlines toward a population of fast radio bursts (FRBs), we demonstrate that FRBs provide a unique tool to probe the intergalactic turbulence. We measure the structure function (SF) of dispersion measures (DMs) of FRBs to study the multi-scale electron density fluctuations induced by the intergalactic turbulence. The SF has a large amplitude and a Kolmogorov power-law scaling with angular separations, showing large and correlated DM fluctuations over a range of length scales. Given that the DMs of FRBs are IGM dominated, our result tentatively suggests that the intergalactic turbulence has a Kolmogorov power spectrum and an outer scale on the order of $100$ Mpc.

preprint2020arXiv

Spectral Superresolution of Multispectral Imagery with Joint Sparse and Low-Rank Learning

Extensive attention has been widely paid to enhance the spatial resolution of hyperspectral (HS) images with the aid of multispectral (MS) images in remote sensing. However, the ability in the fusion of HS and MS images remains to be improved, particularly in large-scale scenes, due to the limited acquisition of HS images. Alternatively, we super-resolve MS images in the spectral domain by the means of partially overlapped HS images, yielding a novel and promising topic: spectral superresolution (SSR) of MS imagery. This is challenging and less investigated task due to its high ill-posedness in inverse imaging. To this end, we develop a simple but effective method, called joint sparse and low-rank learning (J-SLoL), to spectrally enhance MS images by jointly learning low-rank HS-MS dictionary pairs from overlapped regions. J-SLoL infers and recovers the unknown hyperspectral signals over a larger coverage by sparse coding on the learned dictionary pair. Furthermore, we validate the SSR performance on three HS-MS datasets (two for classification and one for unmixing) in terms of reconstruction, classification, and unmixing by comparing with several existing state-of-the-art baselines, showing the effectiveness and superiority of the proposed J-SLoL algorithm. Furthermore, the codes and datasets will be available at: https://github.com/danfenghong/IEEE\_TGRS\_J-SLoL, contributing to the RS community.

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

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.

preprint2017arXiv

Cosmological evolution of primordial black holes

The cosmological evolution of primordial black holes (PBHs) is considered. A comprehensive view of the accretion and evaporation histories of PBHs across the entire cosmic history is presented, with focus on the critical mass holes. The critical mass of a PBH for current era evaporation is $M_{cr}\sim 5.1\times10^{14}$ g. Across cosmic time such a black hole will not accrete radiation or matter in sufficient quantity to hasten the inevitable evaporation, if the black hole remains within an average volume of the universe. The accretion rate onto PBHs is most sensitive to the mass of the hole, the sound speed in the cosmological fluid, and the energy density of the accreted components. It is not easy for a PBH to accrete the average cosmological fluid to reach $30M_\odot$ by $z\sim0.1$, the approximate mass and redshift of the merging BHs that were the sources of the gravitational wave events GW150914 and GW151226. A PBH located in an overdense region can undergo enhanced accretion leading to the possibility of growing by many orders of magnitude across cosmic history. Thus, two merging PBHs are a plausible source for the observed gravitational wave events. However, it is difficult for isolated PBHs to grow to supermassive black holes (SMBHs) at high redshift with masses large enough to fit observational constraints.

preprint2016arXiv

A Comparative Study of Long and Short GRBs. I. Overlapping Properties

Gamma ray bursts (GRBs) are classified into long and short categories based on their durations. Broad band studies suggest that these two categories of objects roughly correspond to two different classes of progenitor systems, i.e. compact star mergers (Type I) vs. massive star core collapse (Type II). However, the duration criterion sometimes leads to mis-identification of the progenitor systems. We perform a comprehensive multi-wavelength comparative study between duration-defined long GRBs and short GRBs as well as the so-called "consensus" long GRBs and short GRBs (which are believed to be more closely related to the two types of progenitor systems). The parameters we study include two parts: the prompt emission properties including duration ($T_{90}$), spectral peak energy ($E_{\rm p}$), low energy photon index ($α$), isotropic $γ$-ray energy ($E_{\rm γ, iso}$), isotropic peak luminosity ($L_{\rm p,iso}$), and the amplitude parameters ($f$ and $f_{\rm eff}$); and the host galaxy properties including stellar mass ($M_*$), star formation rate (SFR), metallicity ([X/H]), half light radius ($R_{50}$), angular and physical ($R_{\rm off}$) offset of the afterglow from the center of the host galaxy, the normalized offset ($r_{\rm off}=R_{\rm off}/R_{50}$), and the brightness fraction $F_{\rm light}$. For most parameters, we find interesting overlapping properties between the two populations in both 1D and 2D distribution plots. The three best parameters for the classification purpose are $T_{90}$, $f_{\rm eff}$, and $F_{\rm light}$. However, no single parameter alone is good enough to place a particular burst into the right physical category, suggesting a need of multiple criteria for physical classification.

preprint2016arXiv

A method to constrain mass and spin of GRB black hole within the NDAF model

Black holes (BHs) hide themselves behind various astronomical phenomena, and their properties, i.e., mass and spin, are usually difficult to constrain. One leading candidate for the central engine model of gamma-ray bursts (GRBs) invokes a stellar mass BH and a neutrino-dominated accretion flow (NDAF), with the relativistic jet launched due to neutrino-anti-neutrino annihilations. Such a model gives rise to a matter-dominated fireball, and is suitable to interpret GRBs with a dominant thermal component with a photospheric origin. We propose a method to constrain BH mass and spin within the framework of this model, and apply the method to a thermally-dominant GRB 101219B whose initial jet launching radius $r_0$ is constrained from the data. Using our numerical model of NDAF jets, we estimate the following constraints on the central BH: mass $M_{\rm BH} \sim 5-9~M_\odot$, spin parameter $a_* \gtrsim 0.6$, and disk mass $3~M_\odot \lesssim M_{\rm disk} \lesssim 4~M_\odot$. Our results also suggest that the NDAF model is a competitive candidate for the central engine of GRBs with a strong thermal component.

preprint2016arXiv

A Statistical Study of GRB X-ray Flares: Evidence of Ubiquitous Bulk Acceleration in the Emission Region

When emission in a conical relativistic jet ceases abruptly (or decays sharply), the observed decay light curve is controlled by the high-latitude "curvature effect". Recently, Uhm & Zhang found that the decay slopes of three GRB X-ray flares are steeper than what the standard model predicts. This requires bulk acceleration of the emission region, which is consistent with a Poynting-flux-dominated outflow. In this paper, we systematically analyze a sample of 85 bright X-ray flares detected in 63 Swift GRBs, and investigate the relationship between the temporal decay index $α$ and spectral index $β$ during the steep decay phase of these flares. The $α$ value depends on the choice of the zero time point $t_0$. We adopt two methods. "Method I" takes $t_0^I$ as the first rising data point of each flare, and is the most conservative approach. We find that at 99.9% condifence level 56/85 flares have decay slopes steeper than the simplest curvature effect prediction, and therefore, are in the acceleration regime. "Method II" extrapolates the rising light curve of each flare backwards until the flux density is three orders of magnitude lower than the peak flux density, and defines the corresponding time as the time zero point (t_0^II). We find that 74/85 flares fall into the acceleration regime at 99.9% condifence level. This suggests that bulk acceleration is common, may be even ubiquitous among X-ray flares, pointing towards a Poynting-flux-dominated jet composition for these events.

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

Central Engine Memory of Gamma-Ray Bursts and Soft Gamma-Ray Repeaters

Gamma-ray Bursts (GRBs) are bursts of $γ$-rays generated from relativistic jets launched from catastrophic events such as massive star core collapse or binary compact star coalescence. Previous studies suggested that GRB emission is erratic, with no noticeable memory in the central engine. Here we report a discovery that similar light curve patterns exist within individual bursts for at least some GRBs. Applying the Dynamic Time Warping (DTW) method, we show that similarity of light curve patterns between pulses of a single burst or between the light curves of a GRB and its X-ray flare can be identified. This suggests that the central engine of at least some GRBs carries "memory" of its activities. We also show that the same technique can identify memory-like emission episodes in the flaring emission in Soft Gamma-Ray Repeaters (SGRs), which are believed to be Galactic, highly magnetized neutron stars named magnetars. Such a phenomenon challenges the standard black hole central engine models for GRBs, and suggest a common physical mechanism behind GRBs and SGRs, which points towards a magnetar central engine of GRBs.

preprint2016arXiv

Collision-induced magnetic reconnection and a unified interpretation of polarization properties of GRBs and blazars

The jet composition and energy dissipation mechanism of Gamma-ray bursts (GRBs) and Blazars are fundamental questions which remain not fully understood. One plausible model is to interpret the $γ$-ray emission of GRBs and optical emission of blazars as synchrotron radiation of electrons accelerated from the collision-induced magnetic dissipation regions in Poynting-flux-dominated jets. The polarization observation is an important and independent information to test this model. Based on our recent 3D relativistic MHD simulations of collision-induced magnetic dissipation of magnetically dominated blobs, here we perform calculations of the polarization properties of the emission in the dissipation region and apply the results to model the polarization observational data of GRB prompt emission and blazar optical emission. We show that the same numerical model with different input parameters can reproduce well the observational data of both GRBs and blazars, especially the $90^{\circ}$ polarization angle (PA) change in GRB 100826A and the $180^{\circ}$ PA swing in Blazar 3C279. This supports a unified model for GRB and blazar jets, suggesting that collision-induced magnetic reconnection is a common physical mechanism to power the relativistic jet emission from events with very different black hole masses.

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

Detectable MeV neutrinos from black hole neutrino-dominated accretion flows

Neutrino-dominated accretion flows (NDAFs) around rotating stellar-mass black holes (BHs) have been theorized as the central engine of relativistic jets launched in massive star core collapse events or compact star mergers. In this work, we calculate the electron neutrino/anti-neutrino spectra of NDAFs by fully taking into account the general relativistic effects, and investigate the effects of viewing angle, BH spin, and mass accretion rate on the results. We show that even though a typical NDAF has a neutrino luminosity lower than that of a typical supernova (SN), it can reach $10^{50}-10^{51}~{\rm erg~s^{-1}}$ peaking at $\sim 10$ MeV, making them potentially detectable with the upcoming sensitive MeV neutrino detectors if they are close enough to Earth. Based on the observed GRB event rate in the local universe and requiring that at least 3 neutrinos are detected to claim a detection, we estimate a detection rate up to $\sim$ (0.10-0.25) per century for GRB-related NDAFs by the Hyper-Kamiokande (Hyper-K) detector if one neglects neutrino oscillation. If one assumes that all Type Ib/c SNe have an engine-driven NDAF, the Hyper-K detection rate would be $\sim$ (1-3) per century. By considering neutrino oscillations, the detection rate may decrease by a factor of 2-3. Detecting one such event would establish the observational evidence of NDAFs in the universe.

preprint2016arXiv

Efficient Production of High-energy Nonthermal Particles during Magnetic Reconnection in a Magnetically-dominated Ion-Electron Plasma

Magnetic reconnection is a leading mechanism for dissipating magnetic energy and accelerating nonthermal particles in Poynting-flux dominated flows. In this letter, we investigate nonthermal particle acceleration during magnetic reconnection in a magnetically-dominated ion-electron plasma using fully kinetic simulations. For an ion-electron plasma with the total magnetization $σ_0=B^2/(4πn(m_i+m_e)c^2)$, the magnetization for each species is $σ_i \sim σ_0$ and $σ_e \sim (m_i/m_e) σ_0$, respectively. We have studied the magnetically dominated regime by varying $σ_{e} = 10^3 - 10^5$ with initial ion and electron temperatures $T_i = T_e = 5 - 20 m_ec^2$ and mass ratio $m_i/m_e = 1 - 1836$. The results demonstrate that reconnection quickly establishes power-law energy distributions for both electrons and ions within several ($2-3$) light-crossing times. For the cases with periodic boundary conditions, the power-law index is $1<s<2$ for both electrons and ions. The hard spectra limit the power-law energies for electrons and ions to be $γ_{be} \sim σ_e$ and $γ_{bi} \sim σ_i$, respectively. The main acceleration mechanism is a Fermi-like acceleration through the drift motions of charged particles. When comparing the spectra for electrons and ions in momentum space, the spectral indices $s_p$ are identical as predicted in Fermi acceleration. We also find that the bulk flow can carry a significant amount of energy during the simulations. We discuss the implication of this study in the context of Poynting-flux dominated jets and pulsar winds especially the applications for explaining the nonthermal high-energy emissions.

preprint2016arXiv

Evidence of Bulk Acceleration of the GRB X-ray Flare Emission Region

Applying our recently-developed generalized version of the high-latitude emission theory to the observations of X-ray flares in gamma-ray bursts (GRBs), we present here clear observational evidence that the X-ray flare emission region is undergoing rapid bulk acceleration as the photons are emitted. We show that both the observed X-ray flare light curves and the photon index evolution curves can be simultaneously reproduced within a simple physical model invoking synchrotron radiation in an accelerating emission region far from the GRB central engine. Such an acceleration process demands an additional energy dissipation source other than kinetic energy, which points towards a significant Poynting-flux in the emission region of X-ray flares. As the X-ray flares are believed to share a similar physical mechanism as the GRB prompt emission, our finding here hints that the GRB prompt emission jets may also carry a significant Poynting-flux in their emitting region.

preprint2016arXiv

Extracting host galaxy dispersion measure and constraining cosmological parameters using fast radio burst data

The excessive dispersion measures (DMs) and high Galactic latitudes of fast radio bursts (FRBs) hint toward a cosmological origin of these mysterious transients. Methods of using measured DM and redshift $z$ to study cosmology have been proposed, but one needs to assume a certain amount of DM contribution from the host galaxy (DM$_{\rm HG}$) in order to apply those methods. We introduce a slope parameter $β(z) \equiv d \ln \left< {\rm DM}_{\rm E} \right> / d \ln z$ (where DM$_{\rm E}$ is the observed DM subtracting the Galactic contribution), which can be directly measured when a sample of FRBs have $z$ measured. We show that $\left< {\rm DM_{HG}}\right>$ can be roughly inferred from $β$ and the mean values, $\overline{\rm \left<DM_{\rm E}\right>}$ and $\bar z$, of the sample. Through Monte Carlo simulations, we show that the mean value of local host galaxy DM, $\left<\rm{DM_{HG,loc}}\right>$, along with other cosmological parameters (mass density $Ω_m$ in the $Λ$CDM model, and the IGM portion of the baryon energy density $Ω_b f_{\rm IGM}$) can be independently measured through MCMC fitting to the data.

preprint2016arXiv

Gamma Ray Burst reverse shock emission in early radio afterglows

Reverse shock (RS) emission from Gamma Ray Bursts is an important tool in investigating the nature of the ejecta from the central engine. If the ejecta magnetization is not high enough to suppress the RS, a strong RS emission component, usually peaking in the optical/IR band early on, would give important contribution to early afterglow light curves. In the radio band, synchrotron self-absorption may suppress early RS emission, and also delay the RS peak time. In this paper, we calculate the self-absorbed RS emission in the radio band for different dynamical conditions. In particular, we stress that the RS radio emission is subject to self-absorption in both reverse and forward shocks. We calculate the ratio between the reverse to forward shock flux at the RS peak time for different frequencies, which is a measure of the detectability of the RS emission component. We then constrain the range of physical parameters for a detectable RS, in particular the role of magnetization. We notice that unlike optical RS emission which is enhanced by moderate magnetization, a moderately magnetized ejecta does not necessarily produce a brighter radio RS due to the self-absorption effect. For typical parameters, the RS emission component would not be detectable below 1 GHz unless the medium density is very low (e.g. $n < 10^{-3} ~{\rm cm^{-3}}$ for ISM and $A_* < 5\times 10^{-4}$ for wind). These predictions can be tested with the afterglow observations with current and upcoming radio facilities such as JVLA, LOFAR, FAST, and SKA.

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

Interpretation of the structure function of rotation measure in the interstellar medium

The observed structure function (SF) of rotation measure (RM) varies as a broken power-law function of angular scales. The systematic shallowness of its spectral slope is inconsistent with the standard Kolmogorov scaling. This motivates us to examine the statistical analysis on RM fluctuations. The correlations of RM constructed by Lazarian & Pogosyan (2016) are demonstrated to be adequate in explaining the observed features of RM SFs through a direct comparison between the theoretically obtained and observationally measured SF results. By segregating the density and magnetic field fluctuations and adopting arbitrary indices for their respective power spectra, we find that when the SFs of RM and emission measure have a similar form over the same range of angular scales, the statistics of the RM fluctuations reflect the properties of density fluctuations. RM SFs can be used to evaluate the mean magnetic field along the line of sight, but cannot serve as an informative source on the properties of turbulent magnetic field in the interstellar medium. We identify the spectral break of RM SFs as the inner scale of a shallow spectrum of electron density fluctuations, which characterizes the typical size of discrete electron density structures in the observed region.

preprint2016arXiv

Lorentz factor - Beaming Corrected Energy/Luminosity Correlations and GRB Central Engine Models

We work on a GRB sample whose initial Lorentz factors ($Γ_0$) are constrained with the afterglow onset method and the jet opening angles ($θ_{\rm j}$) are determined by the jet break time. We confirm the $Γ_0$ - $E_{γ,\rm iso}$ correlation by Liang et al. (2010), and the $Γ_0$ - $L_{γ,\rm iso}$ correlation by L{ü} et al. (2012). Furthermore, we find correlations between $Γ_0$ and the beaming corrected $γ$-ray energy ($E_γ$) and mean $γ$-ray luminosity ($L_γ$). By also including the kinetic energy of the afterglow, we find rough correlations (with larger scatter) between $Γ_0$ and the total ($γ$-ray plus kinetic) energy and the total mean luminosity, both for isotropic values and beaming corrected values: these correlations allow us to test the data with GRB central engine models. Limiting our sample to the GRBs that likely have a black hole central engine, we compare the data with theoretical predictions of two types of jet launching mechanisms from BHs, i.e. the non-magnetized $ν\barν$ - annihilation mechanism, and the strongly magnetized Blandford-Znajek (BZ) mechanism. We find that the data are more consistent with the latter mechanism, and discuss the implications of our findings for GRB jet composition.

preprint2016arXiv

Mergers of Charged Black Holes: Gravitational Wave Events, Short Gamma-Ray Bursts, and Fast Radio Bursts

The discoveries of GW 150914, GW 151226, and LVT 151012 suggest that double black hole (BH-BH) mergers are common in the universe. If at least one of the two merging black holes carries certain amount of charge, possibly retained by a rotating magnetosphere, the inspiral of a BH-BH system would drive a global magnetic dipole normal to the orbital plane. The rapidly evolving magnetic moment during the merging process would drive a Poynting flux with an increasing wind power. The magnetospheric activities during the final phase of the merger would make a fast radio burst (FRB) if the BH charge can be as large as a factor of $\hat q \sim (10^{-9}-10^{-8})$ of the critical charge $Q_c$ of the BH. At large radii, dissipation of the Poynting flux energy in the outflow would power a short duration high-energy transient, which would appear as a detectable short-duration gamma-ray burst (GRB) if the charge can be as large as $\hat q \sim (10^{-5}-10^{-4})$. The putative short GRB coincident with GW 150914 recorded by Fermi GBM may be interpreted with this model. Future joint GW/GRB/FRB searches would lead to a measurement or place a constraint on the charges carried by isolate black holes.

preprint2016arXiv

On the afterglow and progenitor of FRB 150418

Keane et al. recently detected a fading radio source following FRB 150418, leading to the identification of a putative host galaxy at $z = 0.492 \pm 0.008$. Assuming that the fading source is the afterglow of FRB 150418, I model the afterglow and constrain the isotropic energy of the explosion to be a few $10^{50}$ erg, comparable to that of a short duration GRB. The outflow may have a jet opening angle of $\sim 0.22$ rad, so that the beaming-corrected energy is below $10^{49}$ erg. The results rule out most FRB progenitor models for this FRB, but may be consistent with either of the following two scenarios. The first scenario invokes a merger of an NS-NS binary, which produced an undetected short GRB and a supra-massive neutron star, which subsequently collapsed into a black hole, probably 100s of seconds after the short GRB. The second scenario invokes a merger of a compact star binary (BH-BH, NS-NS, or BH-NS) system whose pre-merger dynamical magnetospheric activities made the FRB, which is followed by an undetected short GRB-like transient. The gravitational wave (GW) event GW150914 would be a sister of FRB 150418 in this second scenario. In both cases, one expects an exciting prospect of GW/FRB/GRB associations.

preprint2016arXiv

On the origin of the scatter broadening of fast radio burst pulses and astrophysical implications

Fast radio bursts (FRBs) have been identified as extragalactic sources which can make a probe of turbulence in the intergalactic medium (IGM) and their host galaxies. To account for the observed millisecond pulses caused by scatter broadening, we examine a variety of possible models of electron density fluctuations in both the IGM and the host galaxy medium. We find that a short-wave-dominated power-law spectrum of density, which may arise in highly supersonic turbulence with pronounced local dense structures of shock-compressed gas in the host interstellar medium (ISM), can produce the required density enhancements at sufficiently small scales to interpret the scattering timescale of FRBs. It implies that an FRB residing in a galaxy with efficient star formation in action tends to have a broadened pulse. The scaling of the scattering time with dispersion measure (DM) in the host galaxy varies in different turbulence and scattering regimes. The host galaxy can be the major origin of scatter broadening, but contribute to a small fraction of the total DM. We also find that the sheet-like structure of density in the host ISM associated with folded magnetic fields in a viscosity-dominated regime of magnetohydrodynamic (MHD) turbulence cannot give rise to strong scattering. Furthermore, valuable insights into the IGM turbulence concerning the detailed spatial structure of density and magnetic field can be gained from the observed scattering timescale of FRBs. Our results are in favor of the suppression of micro-plasma instabilities and the validity of collisional-MHD description of turbulence properties in the collisionless IGM.

preprint2016arXiv

Radio transient following FRB 150418: afterglow or coincident AGN flare?

Recently, Keane et al. reported the discovery of a fading radio transient following FRB 150418, and interpreted it as the afterglow of the FRB. Williams \& Berger, on the other hand, suggested that the radio transient is analogous to a group of variable radio sources, so that it could be a coincident AGN flare in the observational beam of the FRB. A new observation with VLA showed a re-brightening, which is consistent with the AGN picture. Here, using the radio survey data of Ofek et al., we statistically examine the chance coincidence probability to produce an event like the FRB 150418 transient. We find that the probabilities to produce a variable radio transient with at least the same variability amplitude and signal-to-noise ratio as the FRB 150415 transient, without and with the VLA point, are $P_1 \sim 6 \times 10^{-4}$ and $P_1 \sim 2 \times 10^{-3}$, respectively. In addition, the chance probability to have a fading transient detected following a random time (FRB time) is less than $P_2 \sim 10^{-2.9\pm 1.3}$. Putting these together and assuming that the number of radio sources within one Parkes beam is 16, the final chance coincidence of having an FRB 150418-like radio transient to be unrelated to the FRB is $< 10^{-4.9\pm1.3}$ and $< 10^{-4.4\pm1.3}$, respectively, without and with the VLA point. We conclude that the radio transient following FRB 150418 has a low probability being an unrelated AGN flare, and the possibility of being the afterglow of FRB 150418 is not ruled out.

preprint2016arXiv

Synchrotron heating by a fast radio burst in a self-absorbed synchrotron nebula and its observational signature

Fast radio bursts (FRBs) are mysterious transient sources. If extragalactic, as suggested by their relative large dispersion measures, their brightness temperatures must be extremely high. Some FRB models (e.g. young pulsar model, magnetar giant flare model, or supra-massive neutron star collapse model) suggest that they may be associated with a synchrotron nebula. Here we study a synchrotron-heating process by an FRB in a self-absorbed synchrotron nebula. If the FRB frequency is below the synchrotron self-absorption frequency of the nebula, electrons in the nebula would absorb FRB photons, leading to a harder electron spectrum and enhanced self-absorbed synchrotron emission. In the meantime, the FRB flux is absorbed by the nebula electrons. We calculate the spectra of FRB-heated synchrotron nebulae, and show that the nebula spectra would show a significant hump in several decades near the self-absorption frequency. Identifying such a spectral feature would reveal an embedded FRB in a synchrotron nebula.

preprint2016arXiv

Testing Einstein's weak equivalence principle with a 0.4-nanosecond giant pulse of the Crab pulsar

Einstein's weak equivalence principle (EEP) can be tested through the arrival time delay between photons with different frequencies. Assuming that the arrival time delay is solely caused by the gravitational potential of the Milky Way, we show that a "nano-shot" giant pulse with a time delay between energies corrected for all known effects, e.g. $Δt<0.4~\rm{ns}$, from the Crab pulsar poses a new upper limit on the deviation from EEP, i.e. $Δγ< (0.6-1.8)\times 10^{-15}$. This result provides the hitherto most stringent constraint on the EEP, improving by at least 2 to 3 orders of magnitude from the previous results based on fast radio bursts.

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

Towards an understanding of GRB prompt emission mechanism: I. The origin of spectral lags

Despite decades of investigations, the physical mechanism that powers the bright prompt $γ$-ray emission from gamma-ray bursts (GRBs) is still not identified. One important observational clue that remains not properly interpreted so far is the existence of time lags of broad light curve pulses in different energy bands, named "spectral lags". Here we show that the traditional view invoking the high-latitude emission "curvature effect" of a relativistic jet cannot account for spectral lags. Rather, the observed spectral lags demand the sweep of a spectral peak across the observing energy band in a specific manner. The duration of the broad pulses and inferred typical Lorentz factor of GRBs require that the emission region is in an optically thin emission region far from the GRB central engine. We construct a simple physical model invoking synchrotron radiation from a rapidly expanding outflow. We show that the observed spectral lags appear naturally in our model light-curves given that (1) the gamma-ray photon spectrum is curved (as observed), (2) the magnetic field strength in the emitting region decreases with radius as the region expands in space, and (3) the emission region itself undergoes rapid bulk acceleration as the prompt $γ$-rays are produced. These requirements are consistent with a Poynting-flux-dominated jet abruptly dissipating magnetic energy at a large distance from the engine.

preprint2015arXiv

A Correlated Study of Optical and X-ray Afterglows of GRBs

We study an extensive sample of 87 GRBs for which there are well sampled and simultaneous optical and X-ray light-curves. We extract the cleanest possible signal of the afterglow component, and compare the temporal behaviors of the X-ray light-curve, observed by Swift XRT, and optical data, observed by UVOT and ground-based telescopes for each individual burst. Overall we find 62\% GRBs that are consistent with the standard afterglow model. When more advanced modeling is invoked, up to 91\% of the bursts in our sample may be consistent with the external shock model. A large fraction of these bursts are consistent with occurring in a constant interstellar density medium (ISM) (61\%) while only 39\% of them occur in a wind-like medium. Only 9 cases have afterglow light-curves that exactly match the standard fireball model prediction, having a single power law decay in both energy bands which are observed during their entire duration. In particular, for the bursts with chromatic behavior additional model assumptions must be made over limited segments of the light-curves in order for these bursts to fully agree with the external shock model. Interestingly, for 54\% of the X-ray and 40\% of the optical band observations the end of the shallow decay ($t^{\sim-0.5}$) period coincides with the jet break ($t^{\sim-p}$) time, causing an abrupt change in decay slope. The fraction of the burst that consistent with the external shock model is independent of the observational epochs in the rest frame of GRBs. Moreover, no cases can be explained by the cooling frequency crossing the X-ray or optical band.

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

A Tight $L_{\rm p, iso}-E'_{\rm p}-Γ_0$ Correlation of Gamma-Ray Bursts

We select a sample of 34 gamma-ray bursts (GRBs) whose $Γ_0$ values are derived with the onset peaks observed in the afterglow lightcurves (except for GRB 060218 whose $Γ_0$ is estimated with its radio data), and investigate the correlations among $Γ_0$, the isotropic peak luminosity ($L_{\rm iso}$), and the peak energy ($E_{\rm p,z}$) of the $νf_ν$ spectrum in the cosmological rest frame. An analysis of pair correlations among these observables well confirms the results reported by the previous papers. More interestingly, a tight correlation among $L_{\rm iso}$, $E_{\rm p,z}$, and $Γ_0$ is found from a multiple regression analysis, which takes the form of $L_{\rm iso} \propto E_{\rm p,z}^{1.34\pm 0.14} Γ_0^{1.32\pm 0.19}$ or $E_{\rm p,z} \propto L_{\rm iso}^{0.55\pm 0.06}Γ_0^{-0.50\pm 0.17}$. Nine other GRBs whose $Γ_0$ are derived via the pair production opacity constraint also follow such a correlation. Excluding GRB 060218, the $L_{\rm iso}-E_{\rm p,z}-Γ_0$ correlation is valid, and it even holds in the jet co-moving frame. However, GRB 060218 deviates the $L^{'}_{\rm iso}-E^{'}_{\rm p}$ relation of typical GRBs in the jet co-moving frame with $3σ$. We argue that the $L_{\rm iso} - E_{\rm p, z} - Γ_0$ correlation may be more physical than the $L_{\rm iso} - E_{\rm p,z}$ correlation, since physically the relationship between the observed $L_{\rm iso}$ and $E_{\rm p,z}$ not only depends on radiation physics, but also depends on the bulk motion of the jet. We explore the possible origins of this correlation and discuss its physical implications for understanding GRB jet composition and radiation mechanism.

preprint2015arXiv

Can life survive Gamma-Ray Bursts in the high-redshift universe?

Nearby Gamma-Ray Bursts (GRBs) have been proposed as a possible cause of mass extinctions on Earth. Due to the higher event rate of GRBs at higher redshifts, it has been speculated that life as we know it may not survive above a certain redshift (e.g. $z>0.5$). We examine the duty cycle of lethal (life-threatening) GRBs in the solar neighborhood, in the Sloan Digital Sky Survey (SDSS) galaxies and GRB host galaxies, with the dependence of the long GRB rate on star formation and metallicity properly taken into account. We find that the number of lethal GRBs attacking Earth within the past 500 Myr ($\sim$ epoch of the Ordovician mass extinction) is $0.93$. The number of lethal GRBs hitting a certain planet increases with redshift, thanks to the increasing star formation rate and decreasing metallicity in high-$z$ galaxies. Taking 1 per 500 Myr as a conservative duty cycle for life to survive, as evidenced by our existence, we find that there are still a good fraction of SDSS galaxies beyond $z=0.5$ where the GRB rate at half-mass radius is lower than this value. We derive the fraction of such benign galaxies as a function of redshift through Monte Carlo simulations, and find that the fraction is $\sim 50\%$ at $z\sim 1.5$ and $\sim 10\%$ even at $z \sim 3$. The mass distribution of benign galaxies is dominated by Milky-Way-like ones, thanks to their commonness, relatively large mass, and low star formation rate. GRB host galaxies are among the most dangerous ones.

preprint2015arXiv

Chemical Abundances of Planetary Nebulae in the Substructures of M31

We present deep spectroscopy of planetary nebulae (PNe) that are associated with the substructures of the Andromeda Galaxy (M31). The spectra were obtained with the OSIRIS spectrograph on the 10.4 m GTC. Seven targets were selected for the observations, three in the Northern Spur and four associated with the Giant Stream. The most distant target in our sample, with a rectified galactocentric distance >100 kpc, was the first PN discovered in the outer streams of M31. The [O III] 4363 auroral line was well detected in the spectra of all targets, enabling electron temperature determination. Ionic abundances are derived based on the [O III] temperatures, and elemental abundances of helium, nitrogen, oxygen, neon, sulfur, and argon are estimated. The relatively low N/O and He/H ratios as well as abundance ratios of alpha-elements indicate that our target PNe might belong to populations as old as ~2 Gyr. Our PN sample, including the current seven and the previous three observed by Fang et al., have rather homogeneous oxygen abundances. The study of abundances and the spatial and kinematical properties of our sample leads to the tempting conclusion that their progenitors might belong to the same stellar population, which hints at a possibility that the Northern Spur and the Giant Stream have the same origin. This may be explained by the stellar orbit proposed by Merrett et al. Judging from the position and kinematics, we emphasize that M32 might be responsible for the two substructures. Deep spectroscopy of PNe in M32 will help to assess this hypothesis.

preprint2015arXiv

Evidence for 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. We find that the re-brightening feature in the radio light curve can be naturally explained by the two-component jet model. The possible origin of this structured jet are the Blandford-Znajek and Blandford-Payne mechanisms. We also show that this two-component jet model can interpret the two kinds of quasi-periodic variations in the X-ray light curve: a 200 second quasi-periodic oscillation (QPO) and a 2.7-day quasi-periodic variation. The latter is interpreted by a precessing outer jet launched near the Bardeen-Petterson radius of a warped disk. 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.

preprint2015arXiv

Extra-galactic high-energy transients: event rate densities and luminosity functions

Several types of extra-galactic high-energy transients have been discovered, which include high-luminosity and low-luminosity long-duration gamma-ray bursts (GRBs), short-duration GRBs, supernova shock breakouts (SBOs), and tidal disruption events (TDEs) without or with an associated relativistic jet. In this paper, we apply a unified method to systematically study the redshift-dependent event rate densities and the global luminosity functions (ignoring redshift evolution) of these transients. We introduce some empirical formulae for the redshift-dependent event rate densities for different types of transients, and derive the local specific event rate density, which also represents its global luminosity function. Long GRBs have a large enough sample to reveal features in the global luminosity function, which is best characterized as a triple power law. All the other transients are consistent with having a single power law luminosity function. The total event rate density depends on the minimum luminosity, and we obtain the following values in units of ${\rm Gpc^{-3}~yr^{-1}}$: $0.8^{+0.1}_{-0.1}$ for high-luminosity long GRBs above $ 10^{50}~{\rm erg~s^{-1}}$, $164^{+98}_{-65}$ for low-luminosity long GRBs above $5\times 10^{46}~{\rm erg~s^{-1}}$, $1.3^{+0.4}_{-0.3}$, $1.2^{+0.4}_{-0.3}$, and $3.3^{+1.0}_{-0.8}$ above $ 10^{50}~ {\rm erg~s^{-1}}$ for short GRBs with three different merger delay models (Gaussian, log-normal, and power law), $1.9^{+2.4}_{-1.2}\times 10^4$ above $ 10^{44}~{\rm erg~s^{-1}}$ for SBOs, $ 4.8^{+3.2}_{-2.1}\times10^2 $ for normal TDEs above $10^{44}~ {\rm erg~s^{-1}}$, and $0.03^{+0.04}_{-0.02}$ above $ 10^{48}~ {\rm erg~s^{-1}}$ for TDE jets as discovered by Swift. Intriguingly, the global luminosity functions of different kinds of transients, which cover over 12 orders of magnitude, are consistent with a single power law with an index of -1.6.

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

Metallicity Measurements of Gamma-Ray Burst and Supernova Explosion Sites: Lessons from HII regions in M31

We examine how the small-scale ($<$\,kpc) variation of metallicity within a galaxy, which is found in nearby galaxies, affect the observational estimates of metallicity in the explosion sites of transient events such as core-collapse supernovae (CC SNe) and gamma-ray bursts (GRBs). Assuming the same luminosity, metallicity, and spatial distributions of \ion{H}{ii}\ regions (hereafter HIIRs) as observed in M31, we compute the apparent metallicities that we would obtain when the spectrum of a target region is blended with those of surrounding HIIRs within the length scale of typical spatial resolution. When the spatial resolution of spectroscopy is $\lesssim$ 0.5\,kpc, which is typical for the existing studies of CC SN sites, we find that the apparent metallicities reflect the metallicities of target regions, but with significant systematic uncertainties in some cases. When the spatial resolution is $\gtrsim$ 1.0 kpc, regardless of the target regions (which has a wide range of metallicity that spans $\sim$ 0.6 dex for the M31 HIIRs), we always obtain the apparent metallicities similar to the average metallicity of the M31 HIIRs. Given that the apparent metallicities measured with $\gtrsim$\,kpc scale resolution do not necessarily reflect the immediate environment of the stellar explosions, the current observational estimates of high metallicities for some of the long GRB host galaxies do {\em not} rule out the hypothesis that the long GRBs are exclusively born in a low-metallicity environment.

preprint2015arXiv

On the Curvature Effect of a Relativistic Spherical Shell

We consider a relativistic spherical shell and calculate its spectral flux as received by a distant observer. Using two different methods, we derive a simple analytical expression of the observed spectral flux and show that the well-known relation $\hat α= 2+\hat β$ (between temporal index $\hat α$ and spectral index $\hat β$) of the high-latitude emission is achieved naturally in our derivation but holds only when the shell moves with a constant Lorentz factor $Γ$. Presenting numerical models where the shell is under acceleration or deceleration, we show that the simple $\hat α= 2+\hat β$ relation is indeed deviated as long as $Γ$ is not constant. For the models under acceleration, we find that the light curves produced purely by the high-latitude emission decay initially much steeper than the constant $Γ$ case and gradually resume the $\hat α= 2+\hat β$ relation in about one and half orders of magnitude in observer time. For the models under deceleration, the trend is opposite. The light curves made purely by the high-latitude emission decay initially shallower than the constant $Γ$ case and gradually resume the relation $\hat α= 2+\hat β$ in a similar order of magnitude in observer time. We also show that how fast the Lorentz factor $Γ$ of the shell increases or decreases is the main ingredient determining the initial steepness or shallowness of the light curves.

preprint2015arXiv

On the non-thermal kappa-distributed electrons in planetary nebulae and HII regions: the kappa index and its correlations with other nebular properties

Recently, a suspicion arose that the free electrons in planetary nebulae (PNe) and HII regions might have non-thermal energy distributions. In this scenario, a kappa index is introduced to characterize the electron energy distributions, with smaller kappa values indicating larger deviations from Maxwell-Boltzmann distributions. Assuming that this is the case, we determine the kappa values for a sample of PNe and HII regions by comparing the intensities of [OIII] collisionally excited lines and the hydrogen Balmer jump. We find the average kappa indices of PNe and HII regions to be 27 and 32, respectively. Correlations between the resultant kappa values and various physical properties of the nebulae are examined to explore the potential origin of non-thermal electrons in photoionized gaseous nebulae. However, no positive result is obtained. Thus the current analysis does not lend to support to the idea that kappa-distributed electrons are present in PNe and HII regions.

preprint2015arXiv

On the polarization properties of magnetar giant flare pulsating tails

Three giant flares have been detected so far from soft gamma-ray repeaters, each characterized by an initial short hard spike and a pulsating tail. The observed pulsating tails are characterized by a duration of $\sim100\,\rm{s}$, an isotropic energy of $\sim 10^{44}\,\rm{erg}$, and a pulse period of a few seconds. The pulsating tail emission likely originates from the residual energy after the intense energy release during the initial spike, which forms a trapped fireball composed of a photon-pair plasma in a closed field line region of the magnetars. Observationally the spectra of pulsating tails can be fitted by the superposition of a thermal component and a power-law component, with the thermal component dominating the emission in the early and late stages of the pulsating tail observations. In this paper, assuming that the trapped fireball is from a closed field line region in the magnetosphere, we calculate the atmosphere structure of the optically-thick trapped fireball and the polarization properties of the trapped fireball. By properly treating the photon propagation in a hot, highly magnetized, electron-positron pair plasma, we tally photons in two modes (O mode and E mode) at a certain observational angle through Monte Carlo simulations. Our results suggest that the polarization degree depends on the viewing angle with respect to the magnetic axis of the magnetar, and can be as high as $Π\simeq30\%$ in the $1-30\,\rm{keV}$ band, and $Π\simeq10\%$ in the $30-100\,\rm{keV}$ band, if the line of sight is perpendicular to the magnetic axis.

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

Relativistic MHD simulations of collision-induced magnetic dissipation in Poynting-flux-dominated jets/outflows

We perform 3D relativistic ideal MHD simulations to study the collisions between high-$σ$ (Poynting-flux-dominated) blobs which contain both poloidal and toroidal magnetic field components. This is meant to mimic the interactions inside a highly variable Poynting-flux-dominated jet. We discover a significant electromagnetic field (EMF) energy dissipation with an Alfvénic rate with the efficiency around 35\%. Detailed analyses show that this dissipation is mostly facilitated by the collision-induced magnetic reconnection. Additional resolution and parameter studies show a robust result that the relative EMF energy dissipation efficiency is nearly independent of the numerical resolution or most physical parameters in the relevant parameter range. The reconnection outflows in our simulation can potentially form the multi-orientation relativistic mini-jets as needed for several analytical models. We also find a linear relationship between the $σ$ values before and after the major EMF energy dissipation process. Our results give support to the proposed astrophysical models that invoke significant magnetic energy dissipation in Poynting-flux-dominated jets, such as the internal collision-induced magnetic reconnection and turbulence (ICMART) model for GRBs, and reconnection triggered mini-jets model for AGNs.

preprint2015arXiv

Synchrotron Origin of the Typical GRB Band Function - A Case Study of GRB 130606B

We perform a time-resolved spectral analysis of GRB 130606B within the framework of a fast-cooling synchrotron radiation model with magnetic field strength in the emission region decaying with time, as proposed by Uhm & Zhang. The data from all time intervals can be successfully fit by the model. The same data can be equally well fit by the empirical Band function with typical parameter values. Our results, which involve only minimal physical assumptions, offer one natural solution to the origin of the observed GRB spectra and imply that, at least some, if not all, Band-like GRB spectra with typical Band parameter values can indeed be explained by synchrotron radiation.

preprint2015arXiv

The extension of variability properties in gamma-ray bursts to blazars

Both gamma-ray bursts (GRBs) and blazars have relativistic jets pointing at a small angle from our line of sight. Several recent studies suggested that these two kinds of sources may share similar jet physics. In this work, we explore the variability properties for GRBs and blazars as a whole. We find that the correlation between minimum variability timescale (MTS) and Lorentz factor, $Γ$, as found only in GRBs by Sonbas et al. can be extended to blazars with a joint correlation of $\rm MTS\proptoΓ^{-4.7\pm0.3}$. The same applies to the $\rm MTS\propto \it L_γ^{\rm -1.0\pm0.1}$ correlation as found in GRBs, which can be well extended into blazars as well. These results provide further evidence that the jets in these two kinds of sources are similar despite of the very different mass scale of their central engines. Further investigations of the physical origin of these correlations are needed, which can shed light on the nature of the jet physics.

preprint2015arXiv

The Millisecond Magnetar Central Engine in short GRBs

One favored progenitor model for short duration gamma-ray bursts (SGRBs) is the coalescence of two neutron stars (NS-NS). One possible outcome of such a merger would be a rapidly spinning, strongly magnetized neutron star (known as a millisecond magnetar). These magnetars may be "supra-massive", implying they would collapse to black holes after losing centrifugal support due to magnetic dipole spindown. By systematically analyzing the BAT-XRT light curves of all short GRBs detected by {\em swift}, we test how well the data are consistent with this central engine model of short GRBs. We find that the so-called "extended emission" observed with BAT in some short GRBs are fundamentally the same component as the "internal X-ray plateau" as observed in many short GRBs, which is defined as a plateau in the lightcurve followed by a very rapid drop. Based on how likely a short GRB hosts a magnetar, we characterize the entire {\em Swift} short GRB sample into three categories: the "internal plateau" sample, the "external plateau" sample, and the "no plateau" sample. The derived magnetar surface magnetic field $B_{\rm p}$ and the initial spin period $P_0$ fall into the reasonable range. No GRBs in the internal plateau sample have the total energy exceeding the maximum energy budget of a millisecond magnetar. Assuming that the rapid fall time at the end of the internal plateau is the collapse time of a supra-massive magnetar to a black hole, and applying the measured mass distribution of NS-NS systems in our Galaxy, we constrain the neutron star equation of state (EOS). The data suggest that the NS equation of state is close to the GM1 model, which has a maximum non-rotating NS mass $M_{\rm TOV} \sim 2.37 M_\odot$.

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

A test of the millisecond magnetar central engine model of GRBs with Swift data

A rapidly spinning, strongly magnetized neutron star (magnetar) has been proposed as one possible candidate of the central engine of gamma-ray bursts (GRBs). We systematically analyze the Swift/XRT light curves of long GRBs detected before 2013 August, and characterize them into four categories based on how likely they may harbor a magnetar central engine: Gold, Silver, Aluminum, and Non-magnetar. We also independently analyze the data of short GRBs with a putative magnetar central engine. We then perform a statistical study of various properties of the magnetar samples and the non-magnetar sample, and investigate whether the data are consistent with the hypothesis that there exist two types of central engines. By deriving the physical parameters of the putative magnetars, we find that the observations of the Gold and Silver samples are generally consistent with the predictions of the magnetar model. For a reasonable beaming factor for long GRBs, the derived magnetar surface magnetic field $B_p$ and initial spin period $P_0$ fall into the reasonable range. Magnetar winds in short GRBs, on the other hand, are consistent with being isotropic. No GRB in the magnetar sample has a beam-corrected total energy exceeding the maximum energy budget defined by the initial spin energy of the magnetar, while some non-magnetar GRBs do violate such a limit. With beaming correction, on average the non-magnetar sample is more energetic and luminous than the magnetar samples. Our analysis hints that millisecond magnetars are likely operating in a good fraction, but probably not all, GRBs.

preprint2014arXiv

Cosmological implications of Fast Radio Burst / Gamma-Ray Burst Associations

If a small fraction of Fast Radio Bursts (FRBs) are associated with Gamma-Ray Bursts (GRBs), as recently suggested by Zhang, the combination of redshift measurements of GRBs and dispersion measure (DM) measurements of FRBs opens a new window to study cosmology. At $z<2$ where the universe is essentially fully ionized, detections of FRB/GRB pairs can give an independent measurement of the intergalactic medium portion of the baryon mass fraction, $Ω_b f_{\rm IGM}$, of the universe. If a good sample of FRB/GRB associations are discovered at higher redshifts, the free electron column density history can be mapped, which can be used to probe the reionization history of both hydrogen and helium in the universe. We apply our formulation to GRBs 101011A and 100704A that each might have an associated FRB, and constrained $Ω_b f_{\rm IGM}$ to be consistent with the value derived from other methods. The methodology developed here is also applicable, if the redshifts of FRBs not associated with GRBs can be measured by other means.

preprint2014arXiv

Distributions of Gamma-Ray Bursts and Blazars in the $L_{\rm p}-E_{\rm p}$ Plane and Possible Implications for their Radiation Physics

We present a spectral analysis for a sample of redshift known GRBs observed with {\em Fermi}/GBM. Together with the results derived from our systematical spectral energy distribution modeling with the leptonic models for a {\em Fermi}/LAT blazar sample, we compare the distributions of the GRBs and the blazars by plotting the synchrotron peak luminosity ($L_{\rm s}$) and the corresponding peak photon energy $E_{\rm s}$ of blazars in the $L_{\rm p}-E_{\rm p}$ plane of GRBs, where $L_{\rm p}$ and $E_{\rm p}$ are the peak luminosity and peak photon energy of the GRB time-integrated $νf_ν$ spectrum, respectively. The GRBs are in the high-$L_{\rm p}$, high-$E_{\rm p}$ corner of the plane and a tight $L_{\rm p}-E_{\rm p}$ relation is found, i.e., $L_{\rm p}\propto E_{\rm p}^{2.13^{+0.54}_{-0.46}}$. Both FSRQs and LBLs are clustered in the low-$E_{\rm p}$, low-$L_{\rm p}$ corner. IBLs and HBLs have $E_{\rm s}\sim 2\times 10^{-3} - 10^{2}$ keV and $L_{\rm s} \sim 10^{44} - 10^{47}$ erg s$^{-1}$, but no dependence of $L_{\rm s}$ on $E_{\rm s}$ is found. We show that the tight $L_p-E_p$ relation of GRBs is potentially explained with the synchrotron radiation of fast-cooling electrons in a highly magnetized ejecta, and the weak anti-correlation of $L_{\rm s}-E_{\rm s}$ for FSRQs and LBLs may be attributed to synchrotron radiation of slow-cooling electrons in a moderately magnetized ejecta. The distributions of IBLs and HBLs in the $L_{\rm p}-E_{\rm p}$ plane may be interpreted with synchrotron radiation of fast-cooling electrons in a matter-dominated ejecta. These results may present a unified picture for the radiation physics of relativistic jets in GRBs and blazars within the framework of the leptonic synchrotron radiation models.

preprint2014arXiv

Dynamics and Afterglow Light Curves of GRB Blast Waves Encountering a Density Bump or Void

We investigate the dynamics and afterglow light curves of gamma-ray burst (GRB) blast waves that encounter various density structures (such as bumps, voids, or steps) in the surrounding ambient medium. We present and explain the characteristic response features that each type of density structures in the medium leaves on the forward shock (FS) and reverse shock (RS) dynamics, for blast waves with either a long-lived or short-lived RS. We show that, when the ambient medium density drops, the blast waves exhibit in some cases a period of an actual acceleration (even during their deceleration stage), due to adiabatic cooling of blast waves. Comparing numerical examples that have different shapes of bumps or voids, we propose a number of consistency tests that correct modeling of blast waves needs to satisfy. Our model results successfully pass these tests. Employing a Lagrangian description of blast waves, we perform a sophisticated calculation of afterglow emission. We show that, as a response to density structures in the ambient medium, the RS light curves produce more significant variations than the FS light curves. Some observed features (such as re-brightenings, dips, or slow wiggles) can be more easily explained within the RS model. We also discuss on the origin of these different features imprinted on the FS and RS light curves.

preprint2014arXiv

Fast cooling synchrotron radiation in a decaying magnetic field and $γ$-ray burst emission mechanism

Synchrotron radiation of relativistic electrons is an important radiation mechanism in many astrophysical sources. In the sources where the synchrotron cooling time scale $t_c$ is shorter than the dynamical time scale $t_{dyn}$, electrons are cooled down below the minimum injection energy. It has been believed that such "fast cooling" electrons have an energy distribution $dN_e /dγ_e \propto γ_e^{-2}$, and their synchrotron radiation flux density has a spectral shape $F_ν\propto ν^{-1/2}$. On the other hand, in a transient expanding astrophysical source, such as a gamma-ray burst (GRB), the magnetic field strength in the emission region continuously decreases with radius. Here we study such a system, and find that in a certain parameter regime, the fast cooling electrons can have a harder energy spectrum, and the standard $d N_e / d γ_e \propto γ_e^{-2}$ spectrum is achieved only in the deep fast cooling regime when $t_c \ll t_{dyn}$. We apply this new physical regime to GRBs, and suggest that the GRB prompt emission spectra whose low-energy photon index $α$ has a typical value -1 could be due to synchrotron radiation in this moderately fast cooling regime.

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

Fermi Large Area Telescope Detection of Supernova Remnant RCW 86

Using $5.4$ year Fermi-LAT data, we report the detection of GeV $γ$-ray emission from the shell-type supernova remnant RCW 86 (G315.4-2.3) with a significance of $\sim5.1σ$. The data slightly favors an extended emission of this supernova remnant. The spectral index of RCW 86 is found to be very hard, $Γ\sim1.4$, in the $0.4$ to $300$ GeV range. A one zone leptonic model can well fit the multi-wavelength data from radio to very high energy $γ$-rays. The very hard GeV $γ$-ray spectrum and the inferred low gas density seem to disfavor the hadronic origin of the $γ$-rays. The $γ$-ray behavior of RCW 86 is very similar to several other TeV shell-type supernova remnants, e.g., RX J1713.7-3946, RX J0852.0-4622, SN 1006 and HESS J1731-347.

preprint2014arXiv

Gamma-Ray Burst Prompt Emission

The origin of gamma-ray burst (GRB) prompt emission, bursts of gamma-rays lasting from shorter than one second to thousands of seconds, remains not fully understood after more than 40 years of observations. The uncertainties lie in several open questions in the GRB physics, including jet composition, energy dissipation mechanism, particle acceleration mechanism, and radiation mechanism. Recent broad-band observations of prompt emission with Fermi sharpen the debates in these areas, which stimulated intense theoretical investigations invoking very different ideas. I will review these debates, and argue that the current data suggest the following picture: A quasi-thermal spectral component originating from the photosphere of the relativistic ejecta has been detected in some GRBs. Even though in some cases (e.g. GRB 090902B) this component dominates the spectrum, in most GRBs, this component either forms a sub-dominant "shoulder" spectral component in the low energy spectral regime of the more dominant "Band" component, or is not detectable at all. The main "Band" spectral component likely originates from the optically thin region due to synchrotron radiation. The diverse magnetization in the GRB central engine is likely the origin of the observed diverse prompt emission properties among bursts.

preprint2014arXiv

Gamma-Ray Burst Prompt Emission Light Curves and Power Density Spectra in the ICMART Model

In this paper, we simulate the prompt emission light curves of gamma-ray bursts (GRBs) within the framework of the Internal-Collision-induced MAgnetic Reconnection and Turbulence (ICMART) model. This model applies to GRBs with a moderately-high magnetization parameter $σ$ in the emission region. We show that this model can produce highly variable light curves with both fast and slow components. The rapid variability is caused by many locally Doppler-boosted mini-emitters due to turbulent magnetic reconnection in a moderately-high-$σ$ flow. The run-away growth and subsequent depletion of these mini-emitters as a function time define a broad slow component for each ICMART event. A GRB light curve is usually composed of multiple ICMART events that are fundamentally driven by the erratic GRB central engine activity. Allowing variations of the model parameters, one is able to reproduce a variety of light curves and the power density spectra as observed.

preprint2014arXiv

High Energy Astrophysics

High energy astrophysics is one of the most active branches in contemporary astrophysics. It studies astrophysical objects that emit X-ray and gamma-ray photons, such as accreting super-massive and stellar-size black holes, and various species of neutron stars. With the operations of many space-borne and ground-based observational facilities, high energy astrophysics has enjoyed a rapid development in the past decades. It is foreseen that the field will continue to advance rapidly in the coming decade, with possible ground-breaking discoveries of astrophysical sources in the high-energy neutrino and gravitational wave channels. This Special Issue of Frontiers of Physics is dedicated to a systematic survey of the field of high energy astrophysics as it stands in 2013.

preprint2014arXiv

Internal Energy Dissipation of Gamma-Ray Bursts Observed with Swift: Precursors, Prompt Gamma-rays, Extended emission and Late X-ray Flares

We jointly analyze the gamma-ray burst (GRB) data observed with BAT and XRT on board the Swift mission to present a global view on the internal energy dissipation processes in GRBs, including precursors, prompt gamma-ray emission, extended soft gamma-ray emission, and late X-ray flares. The Bayesian block method is utilized to analyze the BAT lightcurves to identify various emission episodes. Our results suggest that these emission components likely share a same physical origin, which is repeated activation of the GRB central engine. What we observe in the gamma-ray band may be the tip-of-iceberg of more extended underlying activities. The precursor emission, which is detected in about 10% of {\em Swift} GRBs, is preferably detected in those GRBs that have a massive star core-collapse origin. The soft extended emission (EE) tail, on the other hand, is preferably detected in those GRBs that have a compact star merger origin. Bright X-ray emission is detected during the BAT quiescent phases prior to subsequent gamma-ray peaks, implying that X-ray emission may be detectable prior the BAT trigger time. Future GRB alert instruments with soft X-ray capability would be essential to reveal the early stage of GRB central engine activities, sheding light into jet composition and jet launching mechanism in GRBs.

preprint2014arXiv

Low Energy Spectral Index and $E_{p}$ evolution of Quasi-thermal Photosphere Emission of Gamma-Ray Bursts

Recent observations by the Fermi satellite suggest that a photosphere emission component is contributing to the observed spectrum of many GRBs. One important question is whether the photosphere component can interpret the typical ``Band'' function of GRBs with a typical low energy photon spectral index $α\sim -1$. We perform a detailed study of the photosphere emission spectrum by progressively introducing several physical ingredients previously not fully incorporated, including the probability distribution of the location of a dynamically evolving photosphere, superposition of emission from an equal-arrival-time ``volume'' in a continuous wind, the evolution of optical depth of a wind with finite but evolving outer boundary, as well as the effect of different top-hat wind luminosity ($L_w$) profiles. By assuming a co-moving blackbody spectrum emerging from the photosphere, we find that for an outflow with a constant or increasing $L_w$, the low-energy spectrum below the peak energy ($E_{p}$), can be modified to $F_ν\sim ν^{1.5}$ ($α\sim +0.5$). A softer ($-1<α<+0.5$) or flat ($α=-1$) spectrum can be obtained during the $L_w$ decreasing phase or high-latitude-emission-dominated phase. We also study the evolution of $E_{p}$ as a function of wind and photosphere luminosity in this photosphere model. An $E_p-L$ tracking pattern can be reproduced if a certain positive dependence between the dimensionless entropy $η$ and $L_w$ is introduced. However, the hard-to-soft evolution pattern cannot be reproduced unless a contrived condition is invoked. In order to interpret the Band spectrum, a more complicated photosphere model or a different energy dissipation and radiation mechanism are needed.

preprint2014arXiv

Magnetic Field Amplification and Saturation in Turbulence Behind a Relativistic Shock

We have investigated via two-dimensional relativistic MHD simulations the long-term evolution of turbulence created by a relativistic shock propagating through an inhomogeneous medium. In the postshock region, magnetic field is strongly amplified by turbulent motions triggered by preshock density inhomogeneities. Using a long-simulation box we have followed the magnetic-field amplification until it is fully developed and saturated. The turbulent velocity is sub-relativistic even for a strong shock. Magnetic-field amplification is controled by the turbulent motion and saturation occurs when the magnetic energy is comparable to the turbulent kinetic energy. Magnetic-field amplification and saturation depend on the initial strength and direction of the magnetic field in the preshock medium, and on the shock strength. If the initial magnetic field is perpendicular to the shock normal, the magnetic field is first compressed at the shock and then can be amplified by turbulent motion in the postshock region. Saturation occurs when the magnetic energy becomes comparable to the turbulent kinetic energy in the postshock region. If the initial magnetic field in the preshock medium is strong, the postshock region becomes turbulent but significant field amplification does not occur. If the magnetic energy after shock compression is larger than the turbulent kinetic energy in the postshock region, significant field amplification does not occur. We discuss possible applications of our results to gamma-ray bursts and active galactic nuclei.

preprint2014arXiv

Millisecond pulsar interpretation of the Galactic center gamma-ray excess

It was found in the Fermi-LAT data that there is an extended $γ$-ray excess in the Galactic center region. The proposed sources to be responsible for the excess include the dark matter annihilation or an astrophysical alternative from a population of millisecond pulsars (MSPs). Whether or not the MSP scenario can explain the data self-consistently has very important implications for the detection of particle dark matter, which is however, subject to debate in the literature. In this work we study the MSP scenario in detail, based on the detected properties of the MSPs by Fermi-LAT. We build a model of the Milky Way MSPs which can reproduce the $γ$-ray properties of the Fermi-LAT MSPs, and derive the intrinsic luminosity function of the MSPs. The model is then applied to a bulge population of MSPs. We find that the extended $γ$-ray excess can be well explained by the bulge MSPs without violating the detectable flux distribution of MSPs by Fermi-LAT. The spatial distribution of the bulge MSPs as implied by the distribution of low mass X-ray binaries follows a $r^{-2.4}$ profile, which is also consistent with the $γ$-ray excess data. We conclude that the MSP model can explain the Galactic center $γ$-ray excess self-consistently, satisfying all the current observational constraints.

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

Oscillation Driven Magnetospheric Activity In Pulsars

We study the magnetospheric activity in the polar cap region of pulsars under stellar oscillations. The toroidal oscillation of the star propagates into the magnetosphere, which provides additional voltage due to unipolar induction, changes Goldreich-Julian charge density from the traditional value due to rotation, and hence, influences particle acceleration. We present a general solution of the effect of oscillations within the framework of the inner vacuum gap model, and consider three different inner gap modes controlled by curvature radiation, inverse Compton scattering, and two photon annihilation, respectively. With different pulsar parameters and oscillation amplitudes, one of three modes would play a dominant role in defining the gap properties. When the amplitude of oscillation exceeds a critical value, mode changing would occur. Oscillations also lead to change of the size of the polar cap. As applications, we show the inner gap properties under oscillations in both normal pulsars and anomalous X-ray pulsars / soft gamma-ray repeaters (AXPs/SGRs). We interpret the onset of radio emission after glitches/flares in AXPs/SGRs as due to oscilation-driven magnetic activities in these objects, within the framework of both the magnetar model and the solid quark star model. Within the magnetar model, radio activation may be caused by the enlargement of the effective polar cap angle and the radio emission beam due to oscillation; whereas within the solid quark star angle, it may be caused by activation of the pulsar inner gap from below the radio emission death line due to an oscillation-induced voltage enhancement. The model can also explain the glitch-induced radio profile change observed in PSR J1119-6127.

preprint2014arXiv

Photosphere emission in the X-Ray Flares of Swift Gamma-Ray Bursts and Implications for the Fireball Properties

X-ray flares of gamma-ray bursts (GRBs) are usually observed in the soft X-ray range and the spectral coverage is limited. In this paper, we present an analysis of 32 GRB X-ray flares that are simultaneously observed by both BAT and XRT on board the Swift mission, so a joint spectral analysis with a wider spectral coverage is possible. Our results show that the joint spectra of 19 flares are fitted with the absorbed single power-law or the Band function models. More interestingly, the joint spectra of the other 13 X-ray flares are fitted with the absorbed single power-law model plus a black body (BB) component. Phenomenally, the observed spectra of these 13 flares are analogous to several GRBs with a thermal component, but only with a much lower temperature of $kT=1\sim 3$ keV. Assuming that the thermal emission is the photosphere emission of the GRB fireball, we derive the fireball properties of the 13 flares that have redshift measurements, such as the bulk Lorentz factor $Γ_{\rm ph}$ of the outflow. The derived $Γ_{\rm ph}$ range from $50$ to $150$ and a relation of $Γ_{\rm ph}$ to the thermal emission luminosity is found. It is consistent with the $Γ_0-L_{\rm iso}$ relation that are derived for the prompt gamma-ray emission. We discuss the physical implications of these results within the content of jet composition and radiation mechanism of GRBs and X-ray flares.

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

Radio efficiency of pulsars

We investigate radio emission efficiency $ξ$ of pulsars and report a near linear inverse correlation between $ξ$ and the spindown power $\dot E$, as well as a near linear correlation between $ξ$ and pulsar age $τ$. This is a consequence of very weak, if any, dependences of radio luminosity $L$ on pulsar period $P$ and period derivative $\dot{P}$, in contrast to X-ray or $γ$-ray emission luminosities. The analysis of radio fluxes suggests that these correlations are not due to a selection effect, but are intrinsic to the pulsar radio emission physics. We have found that, although with a large variance, the radio luminosity of pulsars is $\left<L\right>\approx 10^{29} \,{\rm erg/s}$, regardless of the position in the $P-\dot P$ diagram. Within such a picture, a model-independent statement can be made that the death line of radio pulsars corresponds to an upper limit in the efficiency of radio emission. If we introduce the maximum value for a radio efficiency into Monte Carlo-based population syntheses we can reproduce the observed sample using the random luminosity model. The Kolmogorov-Smirnov test on a synthetic flux distribution shows high probability of reproducing the observed distribution. Our results suggests that the plasma responsible for generating radio emission is produced under similar conditions regardless of pulsar age, dipolar magnetic field strength, and spin-down rate. The magnetic fields near the pulsar surface are likely dominated by crust-anchored magnetic anomalies, which do not significantly differ among pulsars, leading to similar conditions for generating electron-positron pairs necessary to power radio emission.

preprint2014arXiv

The "amplitude" parameter of Gamma-Ray Bursts and its implications for GRB classification

Traditionally gamma-ray bursts (GRBs) are classified in the $T_{90}$-hardness ratio two-dimensional plane into long/soft and short/hard GRBs. In this paper, we suggest to add the "amplitude" of GRB prompt emission as the third dimension as a complementary criterion to classify GRBs, especially those of short durations. We define three new parameters $f$, $f_{\rm eff}$ and $f_{\rm eff,z}$ as ratios between the measured/simulated peak flux of a GRB/pseudo-GRB and the flux background, and discuss the applications of these parameters to GRB classification. We systematically derive these parameters to find that most short GRBs are likely not "tip-of-iceberg" of long GRBs. However, one needs to be cautious if a short GRB has a relatively small $f$ (e.g. $f<1.5$), since the chance for an intrinsically long GRB to appear as a "disguised" short GRB is higher. Based on avaialble data, we quantify the probability of a disguised short GRB below a certain $f$ value is as $P (<f)\sim 0.78^{+0.71}_{-0.4} f^{-4.33\pm 1.84}$. By progressively "moving" a long GRB to higher redshifts through simulations, we also find that most long GRBs would show up as rest-frame short GRBs above a certain redshift.

preprint2014arXiv

The Physics of Gamma-Ray Bursts and Relativistic Jets

We provide a comprehensive review of major developments in our understanding of gamma-ray bursts, with particular focus on the discoveries made within the last fifteen years when their true nature was uncovered. We describe the observational properties of photons from the radio to multi-GeV bands, both in the prompt emission and the afterglow phases. Mechanisms for the generation of these photons in GRBs are discussed and confronted with observations to shed light on the physical properties of these explosions, their progenitor stars and the surrounding medium. After presenting observational evidence that a powerful, collimated, jet moving at close to the speed of light is produced in these explosions, we describe our current understanding regarding the generation, acceleration, and dissipation of the jet and compare these properties with jets associated with AGNs and pulsars. We discuss mounting observational evidence that long duration GRBs are produced when massive stars die, and that at least some short duration bursts are associated with old, roughly solar mass, compact stars. The question of whether a black-hole or a strongly magnetized, rapidly rotating neutron star is produced in these explosions is also discussed. We provide a brief summary of what we have learned about relativistic collisionless shocks and particle acceleration from GRB afterglow studies, and discuss the current understanding of radiation mechanism during the prompt emission phase. We discuss theoretical predictions of possible high-energy neutrino emission from GRBs and the current observational constraints. Finally, we discuss how these explosions may be used to study cosmology, e.g. star formation, metal enrichment, reionization history, as well as the formation of first stars and galaxies in the universe.

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

A Comprehensive Study of Gamma-Ray Burst Optical Emission: III. Brightness Distributions and Luminosity Functions of Optical Afterglows

We continue our systematic statistical study on optical afterglow data of gamma-ray bursts (GRBs). We present the apparent magnitude distributions of early optical afterglows at different epochs (t= 10^2 s, t = 10^3 s, and 1 hour) for the optical lightcurves of a sample of 93 GRBs (the global sample), and for sub-samples with an afterglow onset bump or a shallow decay segment. For the onset sample and shallow decay sample we also present the brightness distribution at the peak time t_{p} and break time t_{b}, respectively. All the distributions can be fit with Gaussian functions. We further perform Monte Carlo simulations to infer the luminosity function of GRB optical emission at the rest-frame time 10^3 seconds, t_{p}, and t_{b}, respectively. Our results show that a single power-law luminosity function is adequate to model the data, with indices -1.40+/-0.10, -1.06+/- 0.16, and -1.54\+/- 0.22, respectively. Based on the derived rest-frame 10^3 s luminosity function, we generate the intrinsic distribution of the R-band apparent magnitude M_{R} at the observed time 10^{3} seconds post trigger, which peaks at M_{R}=22.5 mag. The fraction of GRBs whose R-band magnitude is fainter than 22 mag, and 25 mag and at the observer time 10^3 seconds are ~63% and ~25%, respectively. The detection probabilities of the optical afterglows with ground-based robotic telescopes and UVOT onboard {Swift} are roughly consistent with that inferred from this intrinsic M_{R} distribution, indicating that the variations of the dark GRB fraction among the samples with different telescopes may be due to the observational selection effect, although the existence of an intrinsically dark GRB population cannot be ruled out.

preprint2013arXiv

A possible connection between Fast Radio Bursts and Gamma-Ray Bursts

The physical nature of Fast Radio Bursts (FRBs), a new type of cosmological transients discovered recently, is not known. It has been suggested that FRBs can be produced when a spinning supra-massive neutron star loses centrifugal support and collapses to a black hole. Here we suggest that such implosions can happen in supra-massive neutron stars shortly (hundreds to thousands of seconds) after their births, and an observational signature of such implosions may have been observed in the X-ray afterglows of some long and short gamma-ray bursts (GRBs). Within this picture, a small fraction of FRBs would be physically connected to GRBs. We discuss possible multi-wavelength electromagnetic signals and gravitational wave signals that might be associated with FRBs, and propose an observational campaign to unveil the physical nature of FRBs. In particular, we strongly encourage a rapid radio follow-up observation of GRBs starting from 100 s after GRB triggers.

preprint2013arXiv

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

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

preprint2013arXiv

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

GRB 121027A: long-lasting, energetic X-ray flares and clues to radiation mechanism and progenitor star

GRB 121027A is un-usual with its extremely long-lasting, energetic X-ray flares. The total energy release in X-ray flares is about one order of magnitude higher than prompt gamma-rays, making it special from most long GRBs. We show that while the prompt gamma-ray emission satisfies the empirical E_{iso}-E_{p} relation of typical long GRBs, the X-ray flares, whose spectra can be fit with a cutoff-power-law model with well-constrained E_p, significantly deviate from such a relation. Nonetheless, a time-resolved spectral analysis of X-ray flares suggest that the X-ray emission is consistent with the L_{iso}-E_{p} relation of long GRBs. We constrain the minimum Lorentz factor of the X-ray flares to be ~14, which is consistent with the Gamma-L_{iso} relation. Our results imply that prompt gamma-ray emission and late X-ray flares share the similar radiation mechanism, but originate from the outflows with different Lorentz factors. We search for similar GRBs from the Swift GRB archives, and find that the z=6.29 GRB 050904 is a carbon copy of GRB 121027A. The long-lasting, energetic X-ray flares in these GRBs demand significant accretion at late times, which point towards a large-radius progenitor star.

preprint2013arXiv

H I Free-Bound Emission of Planetary Nebulae with Large Abundance Discrepancies: Two-Component Models versus Kappa-distributed electrons

The "abundance discrepancy" problem in the study of planetary nebulae (PNe), viz., the problem concerning systematically higher heavy-element abundances derived from optical recombination lines relative to those from collisionally excited lines, has been under discussion for decades, but no consensus on its solution has yet been reached. In this paper we investigate the hydrogen free-bound emission near the Balmer jump region of four PNe that are among those with the largest abundance discrepancies, aiming to examine two recently proposed solutions to this problem: two-component models and Kappa electron energy distributions. We find that the Balmer jump intensities and the spectrum slopes cannot be simultaneously matched by the theoretical calculations based upon single Maxwell-Boltzmann electron-energy distributions, whereas the fitting can be equally improved by introducing Kappa electron energy distributions or an additional Maxwell-Boltzmann component. We show that although H I free-bound emission alone cannot distinguish the two scenarios, it can provide important constraints on the electron energy distributions, especially for cold and low-Kappa plasmas.

preprint2013arXiv

Hyper-accreting black hole as GRB central engine. I: Baryon loading in GRB jets

A hyper-accreting stellar-mass black hole has been long speculated as the best candidate of central engine of gamma-ray bursts (GRBs). Recent rich observations of GRBs by space missions such as Swift and Fermi pose new constraints on GRB central engine models. In this paper, we study the baryon loading processes of a GRB jet launched from a black hole central engine. We consider a relativistic jet powered by $ν\barν$-annihilation or by the Blandford-Znajek (BZ) mechanism. We consider baryon loading from a neutrino-driven wind from a neutrino-cooling-dominated accretion flow. For a magnetically dominated BZ jet, we consider neutron-drifting from the magnetic wall surrounding the jet and subsequent positron capture and proton-neutron inelastic collisions. The minumim baryon loads in both types of jet are calculated. We find that in both cases, a more luminous jet tends to be more baryon poor. A neutrino-driven "fireball" is typically "dirtier" than a magnetically dominated jet, while a magnetically dominated jet can be much cleaner. Both models have the right scaling to interpret the empirical $Γ-L_{\rm iso}$ relation discovered recently. Since some neutrino-driven jets have too much baryon loading as compared with the data, we suggest that at least a good fraction of GRBs should have a magnetically dominated central engine.

preprint2013arXiv

Model-dependent high-energy neutrino flux from Gamma-Ray Bursts

The IceCube Collaboration recently reported a stringent upper limit on the high energy neutrino flux from GRBs, which provides a meaningful constraint on the standard internal shock model. Recent broad band electromagnetic observations of GRBs also challenge the internal shock paradigm for GRBs, and some competing models for gamma-ray prompt emission have been proposed. We describe a general scheme for calculating the GRB neutrino flux, and compare the predicted neutrino flux levels for different models. We point out that the current neutrino flux upper limit starts to constrain the standard internal shock model. The dissipative photosphere models are also challenged if the cosmic ray luminosity from GRBs is at least 10 times larger than the gamma-ray luminosity. If the neutrino flux upper limit continues to go down in the next few years, then it would suggest the following possibilities: 1. the photon-to-proton luminosity ratio in GRBs is anomalously high for shocks, which may be achieved in some dissipative photosphere models and magnetic dissipation models; or 2. the GRB emission site is at a larger radius than the internal shock radius, as expected in some magnetic dissipation models such as the ICMART model.

preprint2013arXiv

On the Non-existence of a Sharp Cooling Break in GRB Afterglow Spectra

Although the widely-used analytical afterglow model of gamma-ray bursts (GRBs) predicts a sharp cooling break $ν_c$ in its afterglow spectrum, the GRB observations so far rarely show clear evidence for a cooling break in their spectra or its corresponding temporal break in their light curves. Employing a Lagrangian description of the blast wave, we conduct a sophisticated calculation of the afterglow emission. We precisely follow the cooling history of non-thermal electrons accelerated into each Lagrangian shell. We show that a detailed calculation of afterglow spectra does not in fact give rise to a sharp cooling break at $ν_c$. Instead, it displays a very mild and smooth transition, which occurs gradually over a few orders of magnitude in energy or frequency. The main source of this slow transition is that different mini-shells have different evolution histories of the comoving magnetic field strength $B$, so that deriving the current value of $ν_c$ of each mini-shell requires an integration of its cooling rate over the time elapsed since its creation. We present the time evolution of optical and X-ray spectral indices to demonstrate the slow transition of spectral regimes, and discuss the implications of our result in interpreting GRB afterglow data.

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.

preprint2013arXiv

Radiation Mechanism and Jet Composition of Gamma-Ray Bursts and GeV-TeV selected Radio Loud Active Galactic Nuclei

Gamma-ray bursts (GRBs) and GeV-TeV selected radio loud Active Galactic Nuclei (AGNs) are compared based on our systematic modeling of the observed spectral energy distributions of a sample of AGNs with a single-zone leptonic model. We show that the correlation between the jet power (P_{jet}) and the prompt gamma-ray luminosity (L_{jet}) of GRBs is consistent, within the uncertainties, with the correlation between jet power and the synchrotron peak luminosity (L_{s, jet}) of flat spectrum radio quasars (FSRQs). Their radiation efficiencies (\varepsilon) are also comparable (>10% for most sources), which increase with the bolometric jet luminosity (L_{bol,jet}) for FSRQs and with the L_{jet} for GRBs with similar power-law indices. BL Lacs do not follow the P_{jet}-L_{s, jet} relation of FSRQs. They have lower \varepsilon and L_{bol, jet} values than FSRQs, and a tentative L_{bol, jet}-\varepsilon relation is also found, with a power-law index being different from that of the FSRQs. The magnetization parameters (σ) of FSRQs are averagely larger than that of BL Lacs. They are anti-correlated with $\varepsilon$ for the FSRQs, but positive correlated with \varepsilon for the BL Lacs. GeV Narrow-line Seyfert 1 galaxies potentially share similar properties with FSRQs. Based on the analogy between GRBs and FSRQs, we suggest that the prompt gamma-ray emission of GRBs is likely produced by synchrotron process in a magnetized jet with high radiation efficiency, similar to FSRQs. The jets of BL Lacs, on the other hand, are less efficient and are likely more matter dominated.

preprint2012arXiv

A Comprehensive Analysis of Fermi Gamma-ray Burst Data: II. $E_{\rm p}$-Evolution Patterns and Implications for the Observed Spectrum-Luminosity Relations

We present a time-resolved spectral analysis of 51 long and 11 short bright GRBs observed with the {\em Femri}/GBM, paying special attention to $E_{\rm p}$ evolution within a same burst. Among 8 single-pulse long GRBs, 5 show hard-to-soft evolution, while 3 show intensity-tracking. The multi-pulse long GRBs have more complicated patterns. Among the GRBs whose time-resolved spectrum is available for the first pulse, almost half (15/32 GRBs) show clear hard-to-soft evolution, and the other half (17/32 GRBs) show clear intensity-tracking. Later pulses typically show the tracking behavior, although a hard-to-soft evolution pattern was identified in the 2nd pulse of 2 GRBs whose pulses are well separated. Statistically, the hard-to-soft evolution pulses tend to be more asymmetric than the intensity-tracking ones, with a steeper rising wing than the falling wing. Short GRBs have $E_{\rm p}$ tracking intensity exclusively with the 16ms time resolution analysis. We performed a simulation analysis, and suggest that at least for some bursts, the late intensity-tracking pulses could be a consequence of overlapping hard-to-soft pulses. However, the fact that the intensity-tracking pattern exists in the first pulse of multi-pulse long GRBs and some single-pulse GRBs suggest that intensity tracking is an independent component, which may operate in some late pulses as well. For the GRBs with measured redshifts, we present a time-resolved $E_{\rm p}-L_{γ, \rm iso}$ correlation analysis and show that the scatter of the correlation is comparable to that of the global Amati/Yonetoku relation. We discuss the predictions of various radiation models regarding $E_{\rm p}$ evolution, as well as the possibility of a precession jet in GRBs. It seems that the data pose great challenge to all these models, and hold the key to unveil the physics of GRB prompt emission.

preprint2012arXiv

A Comprehensive Analysis of Fermi Gamma-ray Burst Data: III. Energy-Dependent T90 Distributions of GBM GRBs and Instrumental Selection Effect on Duration Classification

The durations (T90) of 315 GRBs detected with Fermi/GBM (8-1000 keV) by 2011 September are calculated using the Bayesian Block method. We compare the T90 distributions between this sample and those derived from previous/current GRB missions. We show that the T90 distribution of this GRB sample is bimodal, with a statistical significance level being comparable to those derived from the BeppoSAX/GRBM sample and the Swift/BAT sample, but lower than that derived from the CGRO/BATSE sample. The short-to-long GRB number ratio is also much lower than that derived from the BATSE sample, i.e., 1:6.5 vs 1:3. We measure T90 in several bands, i.e., 8-15, 15-25, 25-50, 50-100, 100-350, and 350-1000 keV, to investigate the energy-dependence effect of the bimodal T90 distribution. It is found that the bimodal feature is well observed in the 50-100 and 100-350 keV bands, but is only marginally acceptable in the 25-50 keV and 350-1000 keV bands. The hypothesis of the bimodality is confidently rejected in the 8-15 and 15-25 keV bands. The T90 distributions in these bands are roughly consistent with those observed by missions with similar energy bands. The parameter T90 as a function of energy follows \bar T90 \propto E^{-0.20\pm 0.02} for long GRBs. Considering the erratic X-ray and optical flares, the duration of a burst would be even much longer for most GRBs. Our results, together with the observed extended emission of some short GRBs, indicate that the central engine activity time scale would be much longer than T90} for both long and short GRBs and the observed bimodal T90 distribution may be due to an instrumental selection effect.

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

Dynamics and Afterglow Light Curves of GRB Blast Waves with a Long-lived Reverse Shock

We perform a detailed study on the dynamics of a relativistic blast wave with the presence of a long-lived reverse shock (RS). Although a short-lived RS has been widely considered, the RS is believed to be long-lived as a consequence of a stratification expected on the ejecta Lorentz factors. The existence of a long-lived RS makes the forward shock (FS) dynamics to deviate from a self-similar Blandford-McKee solution. Employing the "mechanical model" that correctly incorporates the energy conservation, we present an accurate solution for both the FS and RS dynamics. We conduct a sophisticated calculation of the afterglow emission. Adopting a Lagrangian description of the blast wave, we keep track of an adiabatic evolution of numerous shells between the FS and RS. An evolution of the electron spectrum is also followed individually for every shell. We then find the FS and RS light curves by integrating over the entire FS and RS shocked regions, respectively. In particular, we make use of an analytic expression for observed spectral flux, which we derive here in terms of an observed frequency and observer time. Exploring a total of 20 different ejecta stratifications, we explain in detail how a stratified ejecta affects its blast wave dynamics and afterglow light curves. We show that, while the FS light curves are not sensitive to the ejecta stratifications, the RS light curves exhibit much richer features, including steep declines, plateaus, bumps, re-brightenings, and a variety of temporal decay indices. These distinctive RS features may be observable if the RS has higher values of the micophysics parameters than the FS. We discuss possible applications of our results in understanding the GRB afterglow data.

preprint2012arXiv

Episodic Jets as the Central Engine of Gamma-Ray Bursts

Most Gamma-ray bursts (GRBs) have erratic light curves, which demand that the GRB central engine launches an episodic outflow. Recent Fermi observations of some GRBs indicate a lack of the thermal photosphere component as predicted by the baryonic fireball model, which suggests a magnetic origin of GRBs. In view that powerful episodic jets have been observed along with continuous jets in other astrophysical black hole systems, here we propose an intrinsically episodic, magnetically-dominated jet model for GRB central engine. Accumulation and eruption of free magnetic energy in the corona of a differentially-rotating, turbulent accretion flow around a hyperaccreting black hole lead to ejections of episodic, magnetically dominated plasma blobs. These blobs are accelerated magnetically, collide with each other at large radii, trigger rapid magnetic reconnection and turbulence, efficient particle acceleration and radiation, and power the observed episodic prompt gamma-ray emission from GRBs.

preprint2012arXiv

Faint high-energy gamma-ray photon emission of GRB 081006A from Fermi observations

Since the launch of the Fermi gamma - ray Space Telescope on June 11, 2008, the LAT instrument has solidly detected more than 20 GRBs with high energy photon emission above 100 MeV. Using the matched filter technique, 3 more GRBs have also shown evidence of correlation with high energy photon emission as demonstrated by Akerlof et al. In this paper, we present another GRB unambiguously detected by the matched filter technique, GRB 081006A. This event is associated with more than 13 high energy photons above 100 MeV. The likelihood analysis code provided by the $Fermi$ Science Support Center (FSSC) generated an independent verification of this detection by comparison of the Test Statistics (TS) value with similar calculations for random LAT data fields. We have performed detailed temporal and spectral analysis of photons from 8 keV up to 0.8 GeV from the GBM and the LAT. The properties of GRB 081006A can be compared to the other two long duration GRBs detected at similar significance, GRB 080825C and GRB 090217A. We find that GRB 081006A is more similar to GRB 080825C with comparable appearances of late high energy photon emission. As demonstrated previously, there appears to be a surprising dearth of faint LAT GRBs, with only one additional GRB identified in a sample of 74 GRBs. In this unique period when both $Swift$ and $Fermi$ are operational, there is some urgency to explore this aspect of GRBs as fully as possible.

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

GRB 110709A, 111117A and 120107A: Faint high-energy gamma-ray photon emission from Fermi/LAT observations and demographic implications

Launched on June 11, 2008, the LAT instrument onboard the $Fermi$ Gamma-ray Space Telescope has provided a rare opportunity to study high energy photon emission from gamma-ray bursts. Although the majority of such events (27) have been iden tified by the Fermi LAT Collaboration, four were uncovered by using more sensiti ve statistical techniques (Akerlof et al 2010, Akerlof et al 2011, Zheng et al 2 012). In this paper, we continue our earlier work by finding three more GRBs ass ociated with high energy photon emission, GRB 110709A, 111117A and 120107A. To s ystematize our matched filter approach, a pipeline has been developed to identif y these objects in near real time. GRB 120107A is the first product of this anal ysis procedure. Despite the reduced threshold for identification, the number of GRB events has not increased significantly. This relative dearth of events with low photon number prompted a study of the apparent photon number distribution. W e find an extremely good fit to a simple power-law with an exponent of -1.8 $\pm $ 0.3 for the differential distribution. As might be expected, there is a substa ntial correlation between the number of lower energy photons detected by the GBM and the number observed by the LAT. Thus, high energy photon emission is associ ated with some but not all of the brighter GBM events. Deeper studies of the pro perties of the small population of high energy emitting bursts may eventually yi eld a better understanding of these entire phenomena.

preprint2012arXiv

GRB 110721A: photosphere "death line" and the physical origin of the GRB "Band" function

The prompt emission spectra of gamma-ray bursts (GRBs) usually have a dominant component that is well described by a phenomenological "Band" function. The physical origin of this spectral component is debated. Although the traditional interpretation is synchrotron radiation of non-thermal electrons accelerated in internal shocks or magnetic dissipation regions, a growing trend in the community is to interpret this component as modified thermal emission from a dissipative photosphere of a GRB fireball. We analyze the time dependent spectrum of GRB 110721A detected by {\em Fermi} GBM and LAT, and pay special attention to the rapid evolution of the peak energy $E_p$. We define a "death line" of thermally-dominated dissipative photospheric emission in the $E_p - L$ plane, and show that $E_p$ of GRB 110721A at the earliest epoch has a very high $E_p \sim 15$ MeV that is beyond the "death line". Together with the finding that an additional "shoulder" component exists in this burst that is consistent with a photospheric origin, we suggest that at least for some bursts, the "Band" component is not from a dissipative photosphere, but must invoke a non-thermal origin (e.g. synchrotron or inverse Compton) in the optically thin region of a GRB outflow. We also suggest that the rapid "hard-to-soft" spectral evolution is consistent with the quick discharge of magnetic energy in a magnetically-dominated outflow in the optically thin region.

preprint2012arXiv

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

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

preprint2012arXiv

GRB Progenitors and Observational Criteria

Phenomenologically, two classes of GRBs (long/soft vs. short/hard) are identified based on their gamma-ray properties. The boundary between the two classes is vague. Multi-wavelength observations lead to identification of two types of GRB progenitor: one related to massive stars (Type II), and another related to compact stars (Type I). Evidence suggests that the majority of long GRBs belong to Type II, while at least the majority of nearby short GRBs belong to Type I. Nonetheless, counter examples do exist. Both long-duration Type I and short-duration Type II GRBs have been observed. In this talk, I review the complications in GRB classification and efforts in diagnosing GRB progenitor based on multiple observational criteria. In particular, I raise the caution to readily accept that all short/hard GRBs detected by BATSE are due to compact star mergers. Finally, I propose to introduce "amplitude" as the third dimension (besides "duration" and "hardness") to quantify burst properties, and point out that the "tip-of-iceberg" effect may introduce confusion in defining the physical category of GRBs, especially for low-amplitude, high-redshift GRBs.

preprint2012arXiv

Hadronic-Origin TeV gamma-Rays and Ultra-High Energy Cosmic Rays from Centaurus A

Centaurus A (Cen A) is the nearest radio-loud AGN and is detected from radio to very high energy gamma-rays. Its nuclear spectral energy distribution (SED) shows a double-peak feature, which is well explained by the leptonic synchrotron + synchrotron self-Compton model. This model however cannot account for the observed high energy photons in the TeV range, which display a distinct component. Here we show that ~ TeV photons can be well interpreted as the neutral pion decay products from p-gamma interactions of Fermi accelerated high energy protons in the jet with the seed photons around the second SED peak at ~170 keV. Extrapolating the inferred proton spectrum to high energies, we find that this same model is consistent with the detection of 2 ultra-high-energy cosmic ray events detected by Pierre Auger Observatory from the direction of Cen A. We also estimate the GeV neutrino flux from the same process, and find that it is too faint to be detected by current high-energy neutrino detectors.

preprint2012arXiv

Lorentz Factor-Isotropic Luminosity/Energy Correlations of GRBS and Their Interpretation

The bulk Lorentz factor of the gamma-ray burst (GRB) ejecta (Gamma_0) is a key parameter to understand the GRB physics. Liang et al. have discovered a correlation between Gamma_0 and isotropic gamma-ray energy: Gamma_0\proptoE_{gamma,iso,52}^{0.25}. By including more GRBs with updated data and more methods to derive Gamma_0, we confirm this correlation and obtain Gamma_0~91E_{gamma,iso,52}^{0.29}. Evaluating the mean isotropic gamma-ray luminosities L_{gamma,iso} of the GRBs in the same sample, we discover an even tighter correlation Gamma_0~249 L_{gamma,iso,52}^{0.30}. We propose an interpretation to this later correlation. Invoking a neutrino-cooled hyperaccretion disk around a stellar mass black hole as the central engine of GRBs, we derive jet luminosity powered by neutrino annihilation and baryon loading from a neutrino-driven wind. Applying beaming correction, we finally derive Gamma_0\proptoL_{gamma,iso}^{0.22}, which is well consistent with the data. This suggests that the central engine of long GRBs is likely a stellar mass black hole surrounded by a hyper-accreting disk.

preprint2012arXiv

Observational constraints on the external shock prior emission hypothesis of GRBs

An intriguing hypothesis, that there exists a decaying X-ray emission component before the GRB trigger, has been suggested in order to explain the shallow decay phase of the X-ray afterglow detected in many Swift GRBs. If this "prior emission" is from an external shock, one would expect a corresponding optical emission component during the GRB prompt emission phase. In this paper we apply the available prompt optical emission data (both detections and upper limits) to constrain such a scenario. We fit the shallow and normal decay segments of the XRT light curves in our sample with a T_Δ-shifted single power law, and extrapolate the X-ray flux back to the time of the early optical observations. We then use the synchrotron spectrum predicted by the standard external shock model to extrapolate from the X-ray flux to the optical band, and obtain the possible range of the predicted optical flux. Finally, we compare the predictions with the observations. In the cases where later optical data are available, we also compare the shapes of the optical light curves to the predicted optical light curves from the external shock prior emission model. We find that for a good fraction of GRBs (4 out of 8; up to December 2006), the available data already impose severe constraints on the hypothesis. In particular, the expected optical flux from the prior external shock model is higher than what the data allow. We conclude that if the shallow-decay X-ray component were from a prior emission component, it would have to be of an internal origin with optical flux suppressed.

preprint2012arXiv

Spectral and temporal analysis of the joint Swift/BAT-Fermi/GBM GRB sample

Using the gamma-ray bursts simultaneously detected by Swift/BAT and Fermi/GBM we performed a joint spectral and temporal analysis of the prompt emission data and confirm the rough correlation between the BAT-band photon index Gamma_BAT and the peak spectral energy Epeak. With the redshift known sub-sample, we derived the isotropic gamma-ray energy E_gamma,iso and also confirm the E_gamma,iso - Epeak,rest relation, with a larger scatter than the Amati sample but consistent with GBM team analyses. We also compare the T_90 values derived in the GBM band with those derived in the BAT band and find that for long GRBs the BAT T_90 is usually longer than the GBM T_90, while for short GRBs the trend reverses. This is consistent with the soft/hard nature of long/short GRBs and suggests the importance of an energy-dependent temporal analysis of GRBs.

preprint2012arXiv

Statistical Properties of Multiple Optical Emission Components in Gamma-Ray Bursts and Implications

Well-sampled optical lightcurves of 146 gamma-ray bursts (GRBs) are complied from the literature. Multiple optical emission components are extracted with power-law function fits to these lightcurves. We present a systematical analysis for statistical properties and their relations to prompt gamma-ray emission and X-ray afterglow for each component. We show that peak luminosity in the prompt and late flares are correlated and the evolution of the peak luminosity may signal the evolution of the accretion rate. No tight correlation between the shallow decay phase/plateau and prompt gamma-ray emission is found. Assuming that they are due to a long-lasting wind injected by a compact object, we show that the injected behavior favors the scenarios of a long-lasting wind after the main burst episode. The peak luminosity of the afterglow onset is tightly correlated with Eiso, and it is dimmer as peaking later. Assuming that the onset bump is due to the fireball deceleration by the external medium, we examine the Gamma_0-Eiso relation and find that it is confirmed with the current sample. Optical re-brightening is observed in 30 GRBs in our sample. It shares the same relation between the width and the peak time as found in the onset bump, but no clear correlation between the peak luminosity and Eiso as observed in the onset bumps is found. Although its peak luminosity also decays with time, the slope is much shallower than that of the onset peak. We get L t^{-1}_{p}$, being consistent with off-axis observations to an expanding external fireball in a wind-like circum medium. The late re-brightening may signal another jet component. Mixing of different emission components may be the reason for the observed chromatic breaks in different energy bands.

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.

preprint2012arXiv

Unusual Central Engine Activity in the Double Burst GRB 110709B

The double burst, GRB 110709B, triggered Swift/BAT twice at 21:32:39 UT and 21:43:45 UT, respectively, on 9 July 2011. This is the first time we observed a GRB with two BAT triggers. In this paper, we present simultaneous Swift and Konus-WIND observations of this unusual GRB and its afterglow. If the two events originated from the same physical progenitor, their different time-dependent spectral evolution suggests they must belong to different episodes of the central engine, which may be a magnetar-to-BH accretion system.

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

A statistical model for the $γ$-ray variability of the Crab nebula

A statistical scenario is proposed to explain the $γ$-ray variability and flares of the Crab nebula, which were observed recently by the Fermi/LAT. In this scenario electrons are accelerated in a series of knots, whose sizes follow a power-law distribution. These knots presumably move outwards from the pulsar and have a distribution in the Doppler boost factor. The maximal electron energy is assumed to be proportional to the size of the knot. Fluctuations at the highest energy end of the overall electron distribution will result in variable $γ$-ray emission via the synchrotron process in the $\sim 100$ MeV range. Since highly boosted larger knots are rarer than smaller knots, the model predicts that the variability of the synchrotron emission increases with the photon energy. We realize such a scenario with a Monte-Carlo simulation and find that the model can reproduce both the two $γ$-ray flares over a period of $\sim$ year and the monthly scale $γ$-ray flux fluctuations as observed by the Fermi/LAT. The observed $γ$-ray spectra in both the steady and flaring states are also well reproduced.

preprint2011arXiv

A statistical model to explain the gamma-ray variability and flares of the Crab nebula

Recently the AGILE and Fermi/LAT detectors uncovered giant $γ$-ray flares from the Crab nebula. The duration of these flares is a few days. The Fermi/LAT data with monthly time binning further showed significant variability of the synchrotron tail of the emission, while the inverse Compton component was stable. The simultaneous or follow-up observations in X-ray, optical, infrared and radio bands did not find significant flux variation. Based on these observations, we propose that the $γ$-ray variability and flares are due to statistical fluctuations of knots that can accelerate electrons to $\sim$PeV energies. The maximum achievable energy of electrons is adopted to be proportional to the size of the knot, which is assumed to follow a power-law distribution. Thus the low energy electron flux will be stable due to the large number of small knots, while the high energy electron flux may experience large fluctuations. Monte Carlo realization of such a picture can reproduce the observational data quite well given proper model parameters.

preprint2011arXiv

Black hole Spin in Sw J1644+57 and Sw J2058+05

Recently a hard X-ray transient event, Sw J1644+57, was discovered by the Swift satellite. It likely marks the onset of a relativistic jet from a supermassive black hole, possibly triggered by a tidal disruption event. Another candidate in the same category, Sw J2058+05, was also reported. The low event rate suggests that only a small fraction of TDEs launch relativistic jets. A common speculation is that these rare events are related to rapidly spinning black holes. We attribute jet launching to the Blandford-Znajek mechanism, and use the available data to constrain the black hole spin parameter for the two events. It is found that the two black holes indeed carry a moderate to high spin, suggesting that black hole spin is likely the crucial factor behind the Sw J1644+57 - like events.

preprint2011arXiv

Constraint on dark matter annihilation with dark star formation using Fermi extragalactic diffuse gamma-ray background data

It has been proposed that during the formation of the first generation stars there might be a "dark star" phase in which the power of the star comes from dark matter annihilation. The adiabatic contraction process to form the dark star would result in a highly concentrated density profile of the host halo at the same time, which may give enhanced indirect detection signals of dark matter. In this work we investigate the extragalactic $γ$-ray background from dark matter annihilation with such a dark star formation scenario, and employ the isotropic $γ$-ray data from Fermi-LAT to constrain the model parameters of dark matter. The results suffer from large uncertainties of both the formation rate of the first generation stars and the subsequent evolution effects of the host halos of the dark stars. We find, in the most optimistic case for $γ$-ray production via dark matter annihilation, the expected extragalactic $γ$-ray flux will be enhanced by 1-2 orders of magnitude. In such a case, the annihilation cross section of the supersymmetric dark matter can be constrained to the thermal production level, and the leptonic dark matter model which is proposed to explain the positron/electron excesses can be well excluded. Conversely, if the positron/electron excesses are of a dark matter annihilation origin, then the early Universe environment is such that no dark star can form.

preprint2011arXiv

Cosmic ray spectral hardening due to dispersion in the source injection spectra

Recent cosmic ray (CR) experiments discovered that the CR spectra experience a remarkable hardening for rigidity above several hundred GV. We propose that this is caused by the superposition of the CR energy spectra of many sources that have a dispersion in the injection spectral indices. Adopting similar parameters as those of supernova remnants derived from the Fermi $γ$-ray observations, we can reproduce the observational CR spectra of different species well. This may be interpreted as evidence to support the supernova remnant origin of CRs below the knee. We further propose that the same mechanism may explain the "ankle" of the ultra high energy CR spectrum.

preprint2011arXiv

Cosmic ray spectral hardening due to dispersion of source injection spectra

The cosmic ray (CR) energy spectra measured with ATIC, CREAM and PAMELA showed that there is remarkable hardening for rigidity of several hundred GV. We propose that this hardening is due to the superposition of spectra from a population of sources, e.g., supernova remnants (SNRs), whose injection spectral indices have a dispersion. Adopting proper model parameters the observational data can be well explained. It is interesting that the injection source parameters are similar with that derived from gamma-ray observations of SNRs, which may support the SNR-origin of CRs. Furthermore this mechanism provides an alternative explanation of the "ankle-cutoff" structure of the ultra high energy CR spectra.

preprint2011arXiv

Gamma-ray burst rate: high-redshift excess and its possible origins

Prompted by various analyses of long (Type II) GRB rates and their relationship to the cosmic star-formation history, metallicity and luminosity function evolution, we systematically analyze these effects with a Monte Carlo code. We test various cosmic star-formation history models including analytical and empirical models as well as those derived from cosmological simulations. We also explore expressions for metallicity enhancement of the GRB rate with redshift, as presented in the literature, and discuss improvements to these analytic expressions from the point of view of galactic evolution. These are also compared to cosmological simulations on metal enrichment. Additionally we explore possible evolutionary effects of the GRB rate and luminosity function with redshift. The simulated results are tested with the observed Swift sample including the L, z, and peak flux (log N-log P) distributions. The observational data imply that an increase in the GRB rate is necessary to account for the observations at high redshift, although the form of this enhancement is unclear. A rate increase due to lower metallicity at higher redshift may not be the singular cause and is subject to a variety of uncertainties. Alternatively, evolution of the GRB luminosity function break with redshift shows promise as a possible alternative.

preprint2011arXiv

Open Questions in GRB Physics

Open questions in GRB physics are summarized as of 2011, including classification, progenitor, central engine, ejecta composition, energy dissipation and particle acceleration mechanism, radiation mechanism, long term engine activity, external shock afterglow physics, origin of high energy emission, and cosmological setting. Prospects of addressing some of these problems with the upcoming Chinese-French GRB mission, SVOM, are outlined.

preprint2010arXiv

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

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

preprint2010arXiv

A New Classification Method for Gamma-Ray Bursts

Recent Swift observations suggest that the traditional long vs. short GRB classification scheme does not always associate GRBs to the two physically motivated model types, i.e. Type II (massive star origin) vs. Type I (compact star origin). We propose a new phenomenological classification method of GRBs by introducing a new parameter epsilon=E_{gamma, iso,52}/E^{5/3}_{p,z,2}, where E_{γ,iso} is the isotropic gamma-ray energy (in units of 10^{52} erg), and E_{p,z} is the cosmic rest frame spectral peak energy (in units of 100 keV). For those short GRBs with "extended emission", both quantities are defined for the short/hard spike only. With the current complete sample of GRBs with redshift and E_p measurements, the epsilon parameter shows a clear bimodal distribution with a separation at epsilon ~ 0.03. The high-epsilon region encloses the typical long GRBs with high-luminosity, some high-z "rest-frame-short" GRBs (such as GRB 090423 and GRB 080913), as well as some high-z short GRBs (such as GRB 090426). All these GRBs have been claimed to be of the Type II origin based on other observational properties in the literature. All the GRBs that are argued to be of the Type I origin are found to be clustered in the low-epsilon region. They can be separated from some nearby low-luminosity long GRBs (in 3sigma) by an additional T_{90} criterion, i.e. T_{90,z}<~ 5 s in the Swift/BAT band. We suggest that this new classification scheme can better match the physically-motivated Type II/I classification scheme.

preprint2010arXiv

Are all short-hard gamma-ray bursts produced from mergers of compact stellar objects?

The origin and progenitors of short-hard gamma-ray bursts remain a puzzle and a highly debated topic. Recent Swift observations suggest that these GRBs may be related to catastrophic explosions in degenerate compact stars, denoted as "Type I" GRBs. The most popular models include the merger of two compact stellar objects (NS-NS or NS-BH). We utilize a Monte Carlo approach to determine whether a merger progenitor model can self-consistently account for all the observations of short-hard GRBs, including a sample with redshift measurements in the Swift era (z-known sample) and the CGRO/BATSE sample. We apply various merger time delay distributions invoked in compact star merger models to derive the redshift distributions of these Type I GRBs, and then constrain the unknown luminosity function of Type I GRBs using the observed luminosity-redshift (L - z) distributions of the z-known sample. The best luminosity function model, together with the adopted merger delay model, are then applied to confront the peak flux distribution (log N - log P distribution) of the BATSE and Swift samples. We find that for all the merger models invoking a range of merger delay time scales (including those invoking a large fraction of "prompt mergers"), it is difficult to reconcile the models with all the data. The data are instead statistically consistent with the following two possible scenarios. First, that short/hard GRBs are a superposition of compact-star-merger origin (Type I) GRBs and a population of GRBs that track the star formation history, which are probably related to the deaths of massive stars (Type II GRBs). Sec- ond, the entire short/hard GRB population is consistent with a typical delay of 2 Gyr with respect to the star formation history with modest scatter. This may point towards a different Type I progenitor than the traditional compact star merger models.

preprint2010arXiv

Effect of Resonant Neutrino Oscillation on TeV Neutrino Flavor Ratio from Choked GRBs

In the collapsar scenario of the long duration Gamma-Ray Bursts (GRBs), multi-TeV neutrino emission is predicted as the jet makes its way through the stellar envelope. Such a neutrino signal is also expected for more general ``failed'' GRBs in which a putative jet is ``choked'' by a heavy envelope. If the ν_e \rightarrow ν_μneutrino oscillation parameters are in the atmospheric neutrino oscillation range, we show that the resonant oscillation of ν_e\leftrightarrowν_{μ,τ} can take place within the inner high density region of the choked jet progenitor with a heavy envelope, altering the neutrino flavor ratio on its surface to Φ^s_{ν_e}:Φ^s_{ν_μ}:Φ^s_{ν_τ}=5:11:2. Considering vacuum oscillation of these neutrinos on their way to Earth, the final flavor ratio detected on Earth is further modified to either $1:1.095:1.095$ for the large mixing angle solution to the solar neutrino data, or 1:1.3:1.3 for maximal mixing among the muon and tau neutrinos in vacuum.

preprint2010arXiv

GRB 090417B and its Host Galaxy: A Step Towards an Understanding of Optically-Dark Gamma-Ray Bursts

GRB 090417B was an unusually long burst with a T_90 duration of at least 2130 s and a multi-peaked light curve at energies of 15-150 keV. It was optically dark and has been associated with a bright star-forming galaxy at a redshift of 0.345 that is broadly similar to the Milky Way. This is one of the few cases where a host galaxy has been clearly identified for a dark gamma-ray burst and thus an ideal candidate for studying the origin of dark bursts. We find that the dark nature of GRB 090417B cannot be explained by high redshift, incomplete observations, or unusual physics in the production of the afterglow. Assuming the standard relativistic fireball model for the afterglow we find that the optical flux is at least 2.5 mag fainter than predicted by the X-ray flux. The Swift/XRT X -ray data are consistent with the afterglow being obscured by a dense, localized sheet of dust approximately 30-80 pc from the burst along the line of sight. Our results suggest that this dust sheet imparts an extinction of A_V >~ 12 mag, which is sufficient to explain the missing optical flux. GRB 090417B is an example of a gamma-ray burst that is dark due to the localized dust structure in its host galaxy.

preprint2010arXiv

Is GeV Emission from Gamma-Ray Bursts of External Shock Origin?

Recent observations of Gamma-Ray Bursts (GRBs) by the Fermi Large Area Telescope (LAT) revealed a power law decay feature of the high energy emission (above 100 MeV), which led to the suggestion that it originates from a (probably radiative) external shock. We analyze four GRBs (080916C, 090510, 090902B and 090926A) jointly detected by Fermi LAT and Gamma-ray Burst Monitor (GBM), which have high quality lightcurves in both instrument energy bands. Using the MeV prompt emission (GBM) data, we can record the energy output from the central engine as a function of time. Assuming a constant radiative efficiency, we are able to track energy accumulation in the external shock using our internal/external shell model code. By solving for the early evolution of both an adiabatic and a radiative blastwave, we calculate the high energy emission lightcurve in the LAT band and compare it with the observed one for each burst. The late time LAT light curves after T90 can be well fit by the model. However, due to continuous energy injection into the blastwave during the prompt emission phase, the early external shock emission cannot account for the observed GeV flux level. The high energy emission during the prompt phase (before T90) is most likely a superposition of a gradually enhancing external shock component and a dominant emission component that is of an internal origin.

preprint2010arXiv

Late-time detections of the X-ray afterglow of GRB 060729 with Chandra - the latest detections ever of an X-ray afterglow

We report on 5 Chandra observations of the X-ray afterglow of the Gamma-Ray Burst GRB 060729 performed between 2007 March and 2008 May. In all five observations the afterglow is clearly detected. The last Chandra pointing was performed on 2008-May-04, 642 days after the burst - the latest detection of a GRB X-ray afterglow ever. A reanalysis of the Swift XRT light curve together with the three detections by Chandra in 2007 reveals a break at about 1.0 Ms after the burst with a slight steepening of the decay slope from alpha = 1.32 to 1.61. This break coincides with a significant hardening of the X-ray spectrum, consistent with a cooling break in the wind medium scenario, in which the cooling frequency of the afterglow crosses the X-ray band. The last two Chandra observations in 2007 December and 2008 May provide evidence for another break at about one year after the burst. If interpreted as a jet break, this late-time break implies a jet half opening angle of about 14 degrees for a wind medium. Alternatively, this final break may have a spectral origin, in which case no jet break has been observed and the half-opening angle of the jet of GRB 060729 must be larger than about 15 degrees for a wind medium. We compare the X-ray afterglow of GRB 060729 in a wind environment with other bright X-ray afterglows, in particular GRBs 061121 and 080319B, and discuss why the X-ray afterglow of GRB 060729 is such an exceptionally long-lasting event.

preprint2010arXiv

Luminosity Distribution of Gamma-Ray Burst Host Galaxies at redshift z=1 in Cosmological Smoothed Particle Hydrodynamic Simulations: Implications for the Metallicity Dependence of GRBs

We study the relationship between the metallicity of gamma-ray burst (GRB) progenitors and the probability distribution function (PDF) of GRB host galaxies as a function of luminosity using cosmological hydrodynamic simulations of galaxy formation. We impose a maximum limit to the gas metallicity in which GRBs can occur, and examine how the predicted luminosity PDF of GRB host galaxies changes in the simulation. We perform the Kolmogorov-Smirnov test, and show that the result from our simulation agrees with the observed luminosity PDF of core-collapse supernovae (SNe) host galaxies when we assume that the core-collapse SNe trace star formation. When we assume that GRBs occur only in a low-metallicity environment with $Z\lesssim 0.1 \Zsun$, GRBs occur in lower luminosity galaxies, and the simulated luminosity PDF becomes quantitatively consistent with the observed luminosity PDF. The observational bias against the host galaxies of optically dark GRBs owing to dust extinction may be another reason for the lower luminosities of GRB host galaxies, but the observed luminosity PDF of GRB host galaxies cannot be reproduced solely by the dust bias in our simulation.

preprint2010arXiv

Magnetic Field Amplification by Turbulence in A Relativistic Shock Propagating through An Inhomogeneous Medium

We perform two-dimensional relativistic magnetohydrodynamic simulations of a mildly relativistic shock propagating through an inhomogeneous medium. We show that the postshock region becomes turbulent owing to preshock density inhomogeneity, and the magnetic field is strongly amplified due to the stretching and folding of field lines in the turbulent velocity field. The amplified magnetic field evolves into a filamentary structure in two-dimensional simulations. The magnetic energy spectrum is flatter than the Kolmogorov spectrum and indicates that a so-called small-scale dynamo is occurring in the postshock region. We also find that the amount of magnetic-field amplification depends on the direction of the mean preshock magnetic field, and the time scale of magnetic-field growth depends on the shock strength.

preprint2010arXiv

Probing the Nature of High-z Short GRB 090426 with Its Early Optical and X-ray Afterglows

GRB 090426 is a short duration burst detected by Swift ($T_{90}\sim 1.28$ s in the observer frame, and $T_{90}\sim 0.33$ s in the burst frame at $z=2.609$). Its host galaxy properties and some $γ$-ray related correlations are analogous to those seen in long duration GRBs, which are believed to be of a massive-star origin (so-called Type II GRBs). We present the results of its early optical observations with the 0.8-m TNT telescope at Xinglong observatory, and the 1-m LOAO telescope at Mt. Lemmon Optical Astronomy Observatory in Arizona. Our well-sampled optical afterglow lightcurve covers from $\sim 90$ seconds to $\sim 10^4$ seconds post the GRB trigger. It shows two shallow decay episodes that are likely due to energy injection, which end at $\sim 230$ seconds and $\sim 7100$ seconds, respectively. The decay slopes post the injection phases are consistent with each other ($α\simeq 1.22$). The X-ray afterglow lightcurve appears to trace the optical, although the second energy injection phase was missed due to visibility constraints introduced by the {\em Swift} orbit. The X-ray spectral index is $β_X\sim 1.0$ without temporal evolution. Its decay slope is consistent with the prediction of the forward shock model. Both X-ray and optical emission is consistent with being in the same spectral regime above the cooling frequency ($ν_c$). The fact that $ν_c$ is below the optical band from the very early epoch of the observation provides a constraint on the burst environment, which is similar to that seen in classical long duration GRBs. We therefore suggest that death of a massive star is the possible progenitor of this short burst.

preprint2010arXiv

The Internal-Collision-Induced Magnetic Reconnection and Turbulence (ICMART) Model of Gamma-Ray Bursts

The recent Fermi observation of GRB 080916C shows that the bright photosphere emission associated with a putative fireball is missing, which suggests a Poynting-flux-dominated outflow. We propose a model of gamma-ray burst (GRB) prompt emission in the Poynting-flux-dominated regime, namely, the Internal-Collision-induced MAgnetic Reconnection and Turbulence (ICMART) model. It is envisaged that the GRB central engine launches an intermittent, magnetically-dominated wind, and that in the GRB emission region, the ejecta is still moderately magnetized. Similar to the internal shock (IS) model, the mini-shells interact internally at the traditional internal shock radius. Most of these early collision have little energy dissipation, but serve to distort the ordered magnetic field lines. At a certain point, the distortion of magnetic field configuration reaches the critical condition to allow fast reconnection seeds to occur, which induce relativistic MHD turbulence in the interaction regions. The turbulence further distorts field lines easing additional magnetic reconnections, resulting in a runway release of the stored magnetic field energy (an ICMART event). Particles accelerated in the ICMART region radiate synchrotron photons that power the observed gamma-rays. Each ICMART event corresponds to a broad pulse in the GRB lightcurve, and a GRB is composed of multiple ICMART events. This model retains the merits of the IS and other models, but may overcome several difficulties/issues faced by the IS model (e.g. low efficiency, fast cooling, electron number excess, Amati/Yonetoku relation inconsistency, and missing bright photosphere). It predicts two-component variability time scales, and a decreasing Ep and polarization degree during each ICMART event. The model may be applied to most Fermi LAT GRBs that have time-resolved, featureless Band-function spectra (abridged).

preprint2010arXiv

Very High Energy Cosmic Rays from Centaurus A

Centaurus A is the nearest radio-loud AGN and is detected from radio to very high energy gamma-rays. Its nuclear spectral energy distribution shows two peaks, one in the far-infrared band and another at about 150 keV. By assuming the second peak is due to the electron synchrotron emission and the power index for the differential spectrum of the very high energy cosmic ray proton to be 2.7 we show that only pp interaction is responsible for the observed GeV-TeV emission from Centaurus A. We also found that indeed many very high energy cosmic ray protons from Centaurus A can arrive on Earth thus supporting the recent observation of two events by Pierre Auger Observatory.

preprint2010arXiv

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

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

preprint2009arXiv

A Comprehensive Analysis of Swift/XRT Data: IV. Single Power-Law Decaying Lightcurves vs. Canonical Lightcurves and Implications for a Unified Origin of X-rays

By analyzing the Swift/XRT lightcurves detected before 2009 July, we find 19 cases that monotonously decay as a single power law (SPL) with an index of 1 ~ 1.7 from tens (or hundreds) to ~ 10^5 seconds post the GRB trigger, apparently different from the canonical lightcurves characterized by a shallow-to-normal decay transition. No statistical difference is found in their prompt gamma-rays, and the X-ray properties of the two samples are also similar, although the SPL sample tend to have a slightly lower NH value of the host galaxies and larger energy release compared with the canonical sample. The SPL XRT lightcurves in the burst frame gradually merge into a conflux.The normal decay segment for the canonical sample has the same feature. Similar to the normal decay segment, the SPL lightcurves satisfy the closure relations of external shock models. If the X-rays are the afterglow of the GRB fireball, our results indicate that the shallow decay would be due to energy injection and the total energy budget after injection for both samples of GRBs is comparable. More intriguing, we find that a prior X-ray emission model proposed by Yamazaki is more straightforward to interpret the data. We show that the zero times (T0) of the X-rays atisfy a log-normal distribution, and the negligible T0's of the SPL sample are consistent with being the tail of T0 distributions at low end. Referenced to T0, the canonical XRT lightcurves well trace the SPL lightcurves. The T0's of the canonical lightcurves in our analysis are usually much larger than the offsets of the known precursors from the main GRBs, indicating the X-rays would be an emission component from external shocks prior to GRB trigger. The lack of detection of a jet-like break in most XRT lightcurves implies that the opening angle of its jet would be usually large.

preprint2009arXiv

Constraining GRB Initial Lorentz Factor with the Afterglow Onset Feature and Discovery of a Tight Gamma_0-E_iso Correlation

The onset of GRB afterglow is characterized by a smooth bump in the early afterglow lightcurve. We make an extensive search for such a feature. Twenty optically selected GRBs and 12 X-ray selected GRBs are found, among which 17 optically selected GRBs and 2 X-ray-selected GRBs have redshift measurements. We fit the lightcurves with a smooth broken power-law and measure the temporal characteristic timescales of the bumps at FWHM. Strong mutual correlations among these timescales are found, and a dimmer and broader bump tends to peak at a later peak time. The ratio of rising to decaying timescales is almost universal among bursts, but the ratio of the rising time to the peak time varies from 0.3~1. The E_iso is tightly correlated with the peak luminosity and the peak time of the bump in the burst frame. Assuming that the bumps signal the deceleration of the GRB fireballs in a constant density medium, we calculate the initial Lorentz factor (Gamma_0) and the deceleration radius (R_dec) of the GRBs in the optical-selected sample. It is found that Gamma_0 are typically a few hundreds, and the typical deceleration radius is R_dec~10^{17} cm. More intriguingly, a tight correlation between the Gamma_0 and E_iso is found, namely Gamma_0 ~ 195 E_iso, 52}^{0.27} (satisfied for both the optical and X-ray z-known samples). It is helpful to understand GRB physics, and may serve as an indicator of Gamma_0. We find that the early bright X-rays are usually dominated by a different component from the external shock emission, but occasionally (for one case) an achromatic deceleration feature is observed. Components in X-rays would contribute to the diversity of the observed X-ray lightcurves (abridge).

preprint2009arXiv

Evidence of an initially magnetically dominated outflow in GRB 080916C

The composition of gamma-ray burst (GRB) ejecta is still a mystery. The standard model invokes an initially hot ``fireball'' composed of baryonic matter. Here we analyze the broad band spectra of GRB 080916C detected by the Fermi satellite. The featureless Band-spectrum of all five epochs as well as the detections of >~10 GeV photons in this burst place a strong constraint on the prompt emission radius, which is typically >~ 10^{15} cm, independent on the details of the emission process. The lack of detection of a thermal component as predicted by the baryonic models strongly suggests that a significant fraction of the outflow energy is initially not in the ``fireball'' form, but is likely in a Poynting flux entrained with the baryonic matter. The ratio between the Poynting and baryonic fluxes is at least ~(15-20) at the photosphere radius, if the Poynting flux is not directly converted to kinetic energy below the photosphere.

preprint2009arXiv

Left-right asymmetry for pion and kaon production in the semi-inclusive deep inelastic scattering process

We analyze the left-right asymmetry in the semi-inclusive deep inelastic scattering (SIDIS) process without introducing any weighting functions. With the current theoretical understanding, we find that the Sivers effect plays a key role in our analysis. We use the latest parametrization of the Sivers and fragmentation functions to reanalyze the $π^\pm$ production process and find that the results are sensitive to the parametrization. We also extend our calculation on the $K^{\pm}$ production, which can help us know more about the Sivers distribution of the sea quarks and the unfavored fragmentation processes. HERMES kinematics with a proton target, COMPASS kinematics with a proton, deuteron, and neutron target (the information on the neutron target can be effectively extracted from the $^3$He target), and JLab kinematics (both 6 GeV and 12 GeV) with a proton and neutron target are considered in our paper.

preprint2009arXiv

Modeling Gamma-Ray Burst X-Ray Flares within the Internal Shock Model

X-ray afterglow light curves have been collected for over 400 Swift gamma-ray bursts with nearly half of them having X-ray flares superimposed on the regular afterglow decay. Evidence suggests that gamma-ray prompt emission and X-ray flares share a common origin and that at least some flares can only be explained by long-lasting central engine activity. We have developed a shell model code to address the question of how X-ray flares are produced within the framework of the internal shock model. The shell model creates randomized GRB explosions from a central engine with multiple shells and follows those shells as they collide, merge and spread, producing prompt emission and X-ray flares. We pay special attention to the time history of central engine activity, internal shocks, and observed flares, but do not calculate the shock dynamics and radiation processes in detail. Using the empirical E_p - E_iso (Amati) relation with an assumed Band function spectrum for each collision and an empirical flare temporal profile, we calculate the gamma-ray (Swift/BAT band) and X-ray (Swift/XRT band) lightcurves for arbitrary central engine activity and compare the model results with the observational data. We show that the observed X-ray flare phenomenology can be explained within the internal shock model. The number, width and occurring time of flares are then used to diagnose the central engine activity, putting constraints on the energy, ejection time, width and number of ejected shells. We find that the observed X-ray flare time history generally reflects the time history of the central engine, which reactivates multiple times after the prompt emission phase with progressively reduced energy...

preprint2009arXiv

Possible Effects of Pair Echoes on Gamma-Ray Burst Afterglow Emission

High-energy emission from gamma-ray bursts (GRBs) is widely expected but had been sparsely observed until recently when the Fermi satellite was launched. If >TeV gamma rays are produced in GRBs and can escape from the emission region, they are attenuated by the cosmic infrared background photons, leading to regeneration of GeV-TeV secondary photons via inverse-Compton scattering. This secondary emission can last for a longer time than the duration of GRBs, and it is called a pair echo. We investigate how this pair echo emission affects spectra and light curves of high energy afterglows, considering not only prompt emission but also afterglow as the primary emission. Detection of pair echoes is possible as long as the intergalactic magnetic field (IGMF) in voids is weak. We find (1) that the pair echo from the primary afterglow emission can affect the observed high-energy emission in the afterglow phase after the jet break, and (2) that the pair echo from the primary prompt emission can also be relevant, but only when significant energy is emitted in the TeV range, typically E_{gamma, >0.1 TeV} > (Y/(1+Y)) epsilon_e E_k. Even non-detections of the pair echoes could place interesting constraints on the strength of IGMF. The more favorable targets to detect pair echoes may be the "naked" GRBs without conventional afterglow emission, although energetic naked GRBs would be rare. If the IGMF is weak enough, it is predicted that the GeV emission extends to >30-300 s.

preprint2009arXiv

Probing the birth of fast rotating magnetars through high-energy neutrinos

We investigate the high-energy neutrino emission expected from newly born magnetars surrounded by their stellar ejecta. Protons might be accelerated up to 0.1-100 EeV energies possibly by, e.g., the wave dissipation in the winds, leading to hadronic interactions in the stellar ejecta. The resulting PeV-EeV neutrinos can be detected by IceCube/KM3Net with a typical peak time scale of a few days after the birth of magnetars, making the characteristic soft-hard-soft behavior. Detections would be important as a clue to the formation mechanism of magnetars, although there are ambiguities coming from uncertainties of several parameters such as velocity of the ejecta. Non-detections would also lead to useful constraints on the scenario.

preprint2009arXiv

Prompt optical emission and synchrotron self-absorption constraints on emission site of GRBs

We constrain the distance of the Gamma-Ray Burst (GRB) prompt emission site from the explosion centre, R, by determining the location of the electron's self absorption frequency in the GRB prompt optical-to-X/gamma-ray spectral energy distribution, assuming that the optical and the gamma-ray emissions are among the same synchrotron radiation continuum of a group of hot electrons. All possible spectral regimes are considered in our analysis. The method has only two assumed parameters, namely, the bulk Lorentz factor of the emitting source Gamma, and the magnetic field strength B in the emission region (with a weak dependence). We identify a small sample of 4 bursts that satisfy the following three criteria: (1) they all have simultaneous optical and gamma-ray detections in multiple observational time intervals; (2) they all show temporal correlations between the optical and gamma-ray light curves; and (3) the optical emission is consistent with belonging to the same spectral component as the gamma-ray emission. For all the time intervals of these 4 bursts, it is inferred that R \geq 10^{14} (Gamma/300)^{3/4} (B/10^5 Gauss)^{1/4} cm. For a small fraction of the sample, the constraint can be pinned down to R \approx 10^{14} - 10^{15} cm for Gamma ~ 300. For a second sample of bursts with prompt optical non-detections, only upper limits on R can be obtained. We find no inconsistency between the R-constraints for this non-detection sample and those for the detection sample.

preprint2009arXiv

Statistical Properties of Gamma-Ray Burst Polarization

The emission mechanism and the origin and structure of magnetic fields in gamma-ray burst (GRB) jets are among the most important open questions concerning the nature of the central engine of GRBs. In spite of extensive observational efforts, these questions remain to be answered and are difficult or even impossible to infer with the spectral and lightcurve information currently collected. Polarization measurements will lead to unambiguous answers to several of these questions. Recent developments in X-ray and gamma-ray polarimetry techniques have demonstrated a significant increase in sensitivity enabling several new mission concepts, e.g. POET (Polarimeters for Energetic Transients), providing wide field of view and broadband polarimetry measurements. If launched, missions of this kind would finally provide definitive measurements of GRB polarizations. We perform Monte Carlo simulations to derive the distribution of GRB polarizations in three emission models; the synchrotron model with a globally ordered magnetic field (SO model), the synchrotron model with a locally random magnetic field (SR model), and the Compton drag model (CD model). The results show that POET, or other polarimeters with similar capabilities, can constrain the GRB emission models by using the statistical properties of GRB polarizations. In particular, the ratio of the number of GRBs for which the polarization degrees can be measured to the number of GRBs that are detected (N_m/N_d) and the distributions of the polarization degrees (Pi) can be used as the criteria. If N_m/N_d > 30% and Pi is clustered between 0.2 and 0.7, the SO model will be favored. If instead N_m/N_d < 15%, then the SR or CD model will be favored. If several events with Pi > 0.8 are observed, then the CD model will be favored.

preprint2008arXiv

Incidence Rate of GRB-host-DLAs at High Redshift

We study the incidence rate of damped Ly-a systems associated with the host galaxies of gamma-ray bursts (GRB-host-DLAs) as functions of neutral hydrogen column density (N_HI) and projected star formation rate (SFR) using cosmological SPH simulations. Assuming that the occurrence of GRBs is correlated with the local SFR, we find that the median N_HI of GRB-host-DLAs progressively shifts to lower N_HI values with increasing redshift, and the incidence rate of GRB-host-DLAs with log N_HI > 21.0 decreases rapidly at z>=6. Our results suggest that the likelihood of observing the signature of IGM attenuation in GRB afterglows increases towards higher redshift, because it will not be blocked by the red damping wing of DLAs in the GRB host galaxies. This enhances the prospects of using high-redshift GRBs to probe the reionization history of the Universe. The overall incidence rate of GRB-host-DLAs decreases monotonically with increasing redshift, whereas that of QSO-DLAs increases up to z=6. A measurement of the difference between the two incidence rates would enable an estimation of the value of η_grb, which is the mass fraction of stars that become GRBs for a given amount of star formation. Our predictions can be tested by upcoming high-z GRB missions, including JANUS (Joint Astrophysics Nascent Universe Scout) and SVOM (Space multi-band Variable Object Monitor).

preprint2007arXiv

Constraining Galactic $pγ$ Interactions with Cosmic Ray Electron and Positron Spectra

High energy protons produced by various sources of cosmic rays, {\it e.g.} supernova remnants, pulsar wind nebulae, active galactic nuclei and gamma-ray bursts participate in $pγ$ and $pp$ interactions. Although $pp$ interactions may be the dominant mechanism in our Galaxy, it is unclear how important $pγ$ process is. We show that the upper bound on the fraction of protons participating in $pγ$ interactions inside all Galactic astrophysical sources of cosmic rays is 10%.

preprint2007arXiv

Diagnosing GRB Prompt Emission Site with Spectral Cut-Off Energy

The site and mechanism of gamma-ray burst (GRB) prompt emission is still unknown. Although internal shocks have been widely discussed as the emission site of GRBs, evidence supporting other emission sites have been also suggested recently, including the closer-in photosphere where the fireball becomes transparent and further-out radii near the fireball deceleration radius where magnetic dissipation may be important. With the successful operation of the GLAST experiment, prompt high energy emission spectra from many GRBs would be detected in the near future. We suggest that the cut-off energy of the prompt emission spectrum from a GRB depends on both the fireball bulk Lorentz factor and the unknown emission radius from the central engine. If the bulk Lorentz factor could be independently measured (e.g. from early afterglow observations), the observed spectral cutoff energy can be used to diagnose the emission site of gamma-rays. This would provide valuable information to understand the physical origin of the GRB promp emission.

preprint2007arXiv

Gamma-Ray Bursts in the Swift Era

Since the successful launch of NASA's dedicated gamma-ray burst (GRB) mission, Swift, the study of cosmological GRBs has entered a new era. Here I review the rapid observational and theoretical progress in this dynamical research field during the first two-year of the Swift mission, focusing on how observational breakthroughs have revolutionized our understanding of the physical origins of GRBs. Besides summarizing how Swift helps to solve some pre-Swift mysteries, I also list some outstanding problems raised by the Swift observations. An outlook of GRB science in the future, especially in the GLAST era, is briefly discussed.

preprint2007arXiv

Low Luminosity Gamma-Ray Bursts as a Unique Population: Luminosity Function, Local Rate, and Beaming Factor

Swift/BAT has detected ~ 200 long-duration GRBs, with redshift measurements for ~50 of them. We derive the luminosity function (Phi^{HL}) and the local event rate (rho_0^{HL}) of the conventional high luminosity (HL) GRBs by using the z-known Swift GRBs. Our results are generally consistent with that derived from the CGRO/BATSE data. However, the fact that Swift detected a low luminosity (LL) GRB, GRB 060218, at z=0.033 within ~ 2 year of operation, together with the previous detection of the nearby GRB 980425, suggests a much higher local rate for these LL-GRBs. We explore the possibility that LL-GRBs as a distinct GRB population from the HL-GRBs. We find that rho_0^{LL} is 325_{-177}^{+352} Gpc^{-3} yr^{-1}, which is much higher than rho_0^{HL}(1.12_{-0.20}^{+0.43} Gpc^{-3} yr^{-1}). This rate is ~ 0.7% of the local Type Ib/c SNe. Our results, together with the finding that less than 10% of Type Ib/c SNe are associated with off-beam GRBs, suggest that LL-GRBs have a beaming factor typically less than 14, or a jet angle typically wider than 31^{o}. The high local GRB rate, the small beaming factor, and low luminosity make the LL-GRBs distinct from the HL-GRBs. Although the current data could not fully rule out the possibility that both HL- and LL-GRBs are the same population, our results suggest that LL-GRBs are likely a unique GRB population and the observed low redshift GRB sample is dominated by the LL-GRBs.

preprint2007arXiv

Swift and XMM-Newton Observations of the Extraordinary GRB 060729: An afterglow with a more than 100 days X-ray light curve

We report the results of the Swift and XMM observations of the Swift-discovered long Gamma-Ray Burst GRB 060729 ($T_{90}$=115s). The afterglow of this burst was exceptionally bright in X-rays as well as at UV/Optical wavelengths showing an unusually long slow decay phase ($α$=0.14\plm0.02) suggesting a larger energy injection phase at early times than in other bursts. The X-ray light curve displays a break at about 60 ks after the burst. The X-ray decay slope after the break is $α$=1.29\plm0.03. Up to 125 days after the burst we do not detect a jet break, suggesting that the jet opening angle is larger than 28 degrees. In the first 2 minutes after the burst (rest frame) the X-ray spectrum of the burst changed dramatically from a hard X-ray spectrum to a very soft one. We find that the X-ray spectra at this early phase can all be fitted by an absorbed single power law model or alternatively by a blackbody plus power law model. The power law fits show that the X-ray spectrum becomes steeper while the absorption column density decreases. In Swift's UV/Optical telescope the afterglow was clearly detected up to 9 days after the burst in all 6 filters and even longer in some of the UV filters with the latest detection in the UVW1 31 days after the burst. A break at about 50 ks is clearly detected in all 6 UVOT filters from a shallow decay slope of about 0.3 and a steeper decay slope of 1.3. In addition to the \swift observations we also present and discuss the results from a 61 ks ToO observation by XMM. (Abriviated)

preprint2006arXiv

Contribution of GRB Emission to the GeV Extragalactic Diffuse Gamma-Ray Flux

TeV gamma rays emitted by GRBs are converted into electron-positron pairs via interactions with the extragalactic infrared radiation fields. In turn the pairs produced, whose trajectories are randomized by magnetic fields, will inverse Compton scatter off the cosmic microwave background photons. The beamed TeV gamma ray flux from GRBs is thus transformed into a GeV isotropic gamma ray flux, which contributes to the total extragalactic gamma-ray background emission. Assuming a model for the extragalactic radiation fields, for the GRB redshift distribution and for the GRB luminosity function, we evaluate the contribution of the GRB prompt and scattered emissions to the measured extragalactic gamma-ray flux. To estimate this contribution we optimistically require that the energy flux at TeV energies is about 10 times stronger than the energy flux at MeV energies. The resulting gamma-ray diffuse background is only a small fraction of what is observed, allowing blazars and other sources to give the dominant contribution.

preprint2006arXiv

Swift Observations of GRB 050603: An afterglow with a steep late time decay slope

We report the results of Swift observations of the Gamma Ray Burst GRB 050603. With a V magnitude V=18.2 about 10 hours after the burst the optical afterglow was the brightest so far detected by Swift and one of the brightest optical afterglows ever seen. The Burst Alert Telescope (BAT) light curves show three fast-rise-exponential-decay spikes with $T_{90}$=12s and a fluence of 7.6$\times 10^{-6}$ ergs cm$^{-2}$ in the 15-150 keV band. With an $E_{\rm γ, iso} = 1.26 \times 10^{54}$ ergs it was also one of the most energetic bursts of all times. The Swift spacecraft began observing of the afterglow with the narrow-field instruments about 10 hours after the detection of the burst. The burst was bright enough to be detected by the Swift UV/Optical telescope (UVOT) for almost 3 days and by the X-ray Telescope (XRT) for a week after the burst. The X-ray light curve shows a rapidly fading afterglow with a decay index $α$=1.76$^{+0.15}_{-0.07}$. The X-ray energy spectral index was $β_{\rm X}$=0.71\plm0.10 with the column density in agreement with the Galactic value. The spectral analysis does not show an obvious change in the X-ray spectral slope over time. The optical UVOT light curve decays with a slope of $α$=1.8\plm0.2. The steepness and the similarity of the optical and X-ray decay rates suggest that the afterglow was observed after the jet break. We estimate a jet opening angle of about 1-2$^{\circ}$

preprint2006arXiv

The late time evolution of Gamma-Ray Bursts: ending hyperaccretion and producing flares

We consider the properties of a hyperaccretion model for gamma-ray bursts (GRBs) at the late time when the mass supply rate is expected to decrease with time. We point out that the region in the vicinity of the accretor and the accretor itself can play an important role in determining the rate and time behavior of the accretion and ultimately the energy output. Motivated by numerical simulations and theoretical results, we conjecture that the energy release can be repeatedly stopped and then restarted by the magnetic flux accumulated around the accretor. We propose that the episode or episodes when the accretion resumes correspond to X-ray flares discovered recently in a number of GRBs.

preprint2004arXiv

GeV and higher energy photon interactions in gamma-ray burst fireballs and surroundings

We have calculated the opacities and secondary production mechanisms of high energy photons arising in gamma-ray burst internal shocks, using exact cross-sections for the relevant processes. We find that for reasonable choices of parameters, photons in the range of 10's to 100's of GeV may be emitted in the prompt phase. Photons above this range are subject to electron-positron pair production with fireball photons and would be absent from the spectrum escaping the gamma-ray burst. We find that, in such cases, the fireball becomes optically thin again at ultra-high energies ($\gtrsim$ PeV). On the other hand, for sufficiently large fireball bulk Lorentz factors, the fireball is optically thin at all energies. Both for $γγ$ self-absorbed and optically thin cases, the escaping high energy photons can interact with infra-red and microwave background photons to produce delayed secondary photons in the GeV-TeV range. These may be observable with GLAST, or at low redshifts with ground-based air Cherenkov telescopes. Detection of the primary prompt spectrum constrains the bulk Lorentz factor, while detection of delayed secondary gamma-rays would provide a consistency check for the primary spectrum and the bulk Lorentz factor as well as constraints on the intergalactic magnetic field strength.

preprint2003arXiv

Gamma-ray burst early optical afterglows: implications for the initial Lorentz factor and the central engine

Early optical afterglows have been observed from GRB 990123, GRB 021004, and GRB 021211, which reveal rich emission features attributed to reverse shocks. It is expected that Swift will discover many more early afterglows. Here we investigate in a unified manner both the forward and the reverse external shock emission components, and introduce a straightforward recipe for directly constraining the initial Lorentz factor of the fireball using early optical afterglow data. The scheme is largely independent of the shock microphysics. We identify two types of combinations of the reverse and forward shock emission, and explore their parameter regimes. We also discuss a possible diagnostic for magnetized ejecta. There is evidence that the central engine of GRB 990123 is strongly magnetized.

preprint2002arXiv

Magnetars and pulsars: a missing link

There is growing evidence that soft gamma-ray repeaters (SGRs) and anomalous X-ray pulsars (AXPs) are isolated neutron stars with superstrong magnetic fields, i.e., magnetars, marking them a distinguished species from the conventional species of spindown-powered isolated neutron stars, i.e., radio pulsars. The current arguments in favor of the magnetar interpretation of SGR/AXP phenomenology will be outlined, and the two energy sources in magnetars, i.e. a magnetic dissipation energy and a spindown energy, will be reviewed. I will then discuss a missing link between magnetars and pulsars, i.e., lack of the observational evidence of the spindown-powered behaviors in known magnetars. Some recent theoretical efforts in studying such behaviors will be reviewed along with some predictions testable in the near future.

preprint2002arXiv

Regimes of Pulsar Pair Formation and Particle Energetics

We investigate the conditions required for the production of electron-positron pairs above a pulsar polar cap (PC) and the influence of pair production on the energetics of the primary particle acceleration. Assuming space-charge limited flow acceleration including the inertial frame-dragging effect, we allow both one-photon and two-photon pair production by either curvature radiation (CR) photons or photons resulting from inverse-Compton scattering of thermal photons from the PC by primary electrons. We find that, while only the younger pulsars can produce pairs through CR, nearly all known radio pulsars are capable of producing pairs through non-resonant inverse-Compton scatterings. The effect of the neutron star equations of state on the pair death lines is explored. We show that pair production is facilitated in more compact stars and more massive stars. Therefore accretion of mass by pulsars in binary systems may allow pair production in most of the millisecond pulsar population. We also find that two-photon pair production may be important in millisecond pulsars if their surface temperatures are above about three million degrees K. Pulsars that produce pairs through CR will have their primary acceleration limited by the effect of screening of the electric field. In this regime, the high-energy luminosity should follow a $L_{HE} \propto \dot E_{rot}^{1/2}$ dependence. The acceleration voltage drop in pulsars that produce pairs only through inverse-Compton emission will not be limited by electric field screening. In this regime, the high-energy luminosity should follow a $L_{HE} \propto \dot E_{rot}$ dependence. Thus, older pulsars will have significantly lower gamma-ray luminosity.

preprint2000arXiv

What if pulsars are born as strange stars?

The possibility and the implications of the idea, that pulsars are born as strange stars, are explored. Strange stars are very likely to have atmospheres with typical mass of $\sim 5\times 10^{-15}M_\odot$ but bare polar caps almost throughout their lifetimes, if they are produced during supernova explosions. A direct consequence of the bare polar cap is that the binding energies of both positively and negatively charged particles at the bare quark surface are nearly infinity, so that the vacuum polar gap sparking scenario as proposed by Ruderman & Sutherland should operate above the cap, regardless of the sense of the magnetic pole with respect to the rotational pole. Heat can not accumulate on the polar cap region due to the large thermal conductivity on the bare quark surface. We test this ``bare polar cap strange star'' (BPCSS) idea with the present broad band emission data of pulsars, and propose several possible criteria to distinguish BPCSSs from neutron stars.