Source author record

Kunihito Ioka

Kunihito Ioka 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

71works
9topics
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

71 published item(s)

preprint2026arXiv

One-dimensional PIC Simulation of Induced Compton Scattering in Magnetized Electron-Positron Pair Plasma

We investigate induced Compton scattering of a circularly polarized Alfvén wave propagating in a magnetized electron-positron pair plasma using one-dimensional Particle-in-Cell (PIC) simulations. In this system, two distinct modes of density fluctuations, referred to as the charged mode and the neutral mode, are theoretically expected to arise through parametric instabilities. Our simulations confirm these predictions: in the charged mode, the electron and positron densities fluctuate oppositely (Langmuir-like), while in the neutral mode, the charge is Debye-screened and both species fluctuate in phase (acoustic-like). The linear growth rates obtained from the simulations are in good agreement with analytical estimates for both modes. We also find that, in some cases, the linear growth saturates before full scattering occurs, allowing the incident wave to propagate without significant attenuation. Our results allow us to determine whether induced Compton scattering grows linearly in magnetized pair plasmas, offering a foundation for studies of fast radio bursts and laser-plasma experiments.

preprint2023arXiv

Expanding Fireball in Magnetar Bursts and Fast Radio Bursts

A fireball of radiation plasma created near the surface of a neutron star (NS) expands under its own pressure along magnetic field lines, and produces photon emission and relativistic matter outflow. We comprehensively classify the expanding fireball evolution into five cases and obtain the photospheric luminosity and the kinetic energy of the outflow, taking into account key processes; lateral diffusion of photons escaping from a magnetic flux tube, effects of strong magnetic field, baryon loading from the NS surface, and radiative acceleration via cyclotron resonant scattering, some of which have not been considered in the context of gamma-ray bursts. Applying our model to magnetar bursts with fast radio bursts (FRBs), in particular the X-ray short bursts from SGR 1935+2154 associated with the Galactic FRB 20200428A, we show that the burst radiation can accelerate the outflow to high Lorentz factor with sufficient energy to power FRBs.

preprint2022arXiv

Cocoon breakout and escape from the ejecta of neutron star mergers

The cocoon is an inevitable product of a jet propagating through ambient matter, and takes a fair fraction of the jet energy. In short gamma-ray bursts, the ambient matter is the ejecta from the merger of neutron stars, expanding with a high velocity $\sim 0.2 c$, in contrast to the static stellar envelope in collapsars. Using 2D relativistic hydrodynamic simulations with the ejecta density profile as $ρ\propto r^{-2}$, we find that the expansion makes a big difference; only 0.5--5\% of the cocoon mass escapes from (faster than) the ejecta, with an opening angle $20^{\circ}$--$30^{\circ}$, while it is $\sim 100\%$ and spherical in collapsars. We also analytically obtain the shares of mass and energies for the escaped and trapped cocoons. Because the mass of the escaped cocoon is small and the trapped cocoon is concealed by the ejecta and the escaped cocoon, we conclude that it is unlikely that the cocoon emission was observed as a counterpart to the gravitational wave event GW170817.

preprint2022arXiv

Probing cosmology and gastrophysics with fast radio bursts: Cross-correlations of dark matter haloes and cosmic dispersion measures

For future surveys of fast radio bursts (FRBs), we clarify information available from cosmic dispersion measures (DMs) through cross-correlation analyses of foreground dark matter haloes (hosting galaxies and galaxy clusters) with their known redshifts. With a halo-model approach, we predict that the cross-correlation with cluster-sized haloes is less affected by the details of gastrophysics, providing robust cosmological information. For less massive haloes, the cross-correlation at angular scales of $<10\, \mathrm{arcmin}$ is sensitive to gas expelled from the halo centre due to galactic feedback. Assuming $20000$ FRBs over $20000 \, {\rm deg}^2$ with a localisation error being 3 arcmin, we expect that the cross-correlation signal at halo masses of $10^{12}-10^{14}\, M_\odot$ can be measured with a level of $\sim 1\%$ precision in a redshift range of $0<z<1$. Such precise measurements enable one to put a 1.5\% level constraint on $σ_8\, (Ω_\mathrm{M}/0.3)^{0.5}$ and a 3\% level constraint on $(Ω_\mathrm{b}/0.049)(h/0.67)(f_\mathrm{e}/0.95)$ ($σ_8$, $Ω_\mathrm{M}$, $Ω_\mathrm{b}$, $h$ and $f_\mathrm{e}$ are the linear mass variance smoothed at $8\, h^{-1}\mathrm{Mpc}$, mean mass density, mean baryon density, the present-day Hubble parameter and fraction of free electrons in cosmic baryons today), whereas the gas-to-halo mass relation in galaxies and clusters can be constrained with a level of $10\%-20\%$. Furthermore the cross-correlation analyses can break the degeneracy among $Ω_\mathrm{b}$, $h$ and $f_\mathrm{e}$, inherent in the DM-redshift relation. Our proposal opens new possibilities for FRB cosmology, while it requires extensive galaxy redshift catalogues and further improvement of the halo model.

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

Diverse Jet Structures Consistent with the Off-axis Afterglow of GRB 170817A

The jet structure of short gamma-ray bursts (GRBs) has been controversial after the detection of GRB 170817A as the electromagnetic counterparts to the gravitational wave event GW170817. Different authors use different jet structures for calculating the afterglow light curves. We formulated a method to inversely reconstruct a jet structure from a given off-axis GRB afterglow, without assuming any functional form of the structure. By systematically applying our inversion method, we find that more diverse jet structures are consistent with the observed afterglow of GRB 170817A within errors: such as hollow-cone, spindle, Gaussian, and power-law jet structures. In addition, the total energy of the reconstructed jet is arbitrary, proportional to the ambient density $n_0$, with keeping the same jet shape if the parameters satisfy the degeneracy combination $n_0 \varepsilon_\mathrm{B}^{(p+1)/(p+5)} \varepsilon_\mathrm{e}^{4(p-1)/(p+5)} = \mathrm{const.}$. Observational accuracy less than $\sim 6$ per cent is necessary to distinguish the different shapes, while the degeneracy of the energy scaling would be broken by observing the spectral breaks and viewing angle. Future events in denser environment with brighter afterglows and observable spectral breaks are ideal for our inversion method to pin down the jet structure, providing the key to the jet formation and propagation.

preprint2021arXiv

Statistical modelling of the cosmological dispersion measure

We have investigated the basic statistics of the cosmological dispersion measure (DM) -- such as its mean, variance, probability distribution, angular power spectrum and correlation function -- using the state-of-the-art hydrodynamic simulations, IllustrisTNG300, for the fast radio burst (FRB) cosmology. To model the DM statistics, we first measured the free-electron abundance and the power spectrum of its spatial fluctuations. The free-electron power spectrum turns out to be consistent with the dark matter power spectrum at large scales, but it is strongly damped at small scales ($\lesssim 1$Mpc) owing to the stellar and active galactic nucleus feedback. The free-electron power spectrum is well modelled using a scale-dependent bias factor (the ratio of its fluctuation amplitude to that of the dark matter). We provide analytical fitting functions for the free-electron abundance and its bias factor. We next constructed mock sky maps of the DM by performing standard ray-tracing simulations with the TNG300 data. The DM statistics are calculated analytically from the fitting functions of the free-electron distribution, which agree well with the simulation results measured from the mock maps. We have also obtained the probability distribution of source redshift for a given DM, which helps in identifying the host galaxies of FRBs from the measured DMs. The angular two-point correlation function of the DM is described by a simple power law, $ξ(θ) \approx 2400 (θ/{\rm deg})^{-1} \, {\rm pc}^2 \, {\rm cm}^{-6}$, which we anticipate will be confirmed by future observations when thousands of FRBs are available.

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

Analytic properties of electromagnetic field of binary compact stars and electromagnetic precursors to gravitational waves

We analytically study the properties of the electromagnetic field in vacuum around close binary compact stars containing at least one neutron star. We show that the orbital motion of the neutron star induces high multipole modes of the electromagnetic field just before the merger. These modes are superimposed to form a spiral arm configuration and its edge is found to be a likely site for magnetic reconnection. These modes also enhance the total Poynting flux from neutron star binaries by a factor of 2--4. We also indicate that the electric field induced by the orbital motion lead to a magnetosphere around binaries and estimate its plasma density, which has a different parameter dependence than Goldreich-Julian density. With these properties, we discuss possible electromagnetic counterparts to gravitational wave events, and identify radio precursors, such as fast radio bursts, as the most promising observational targets.

preprint2020arXiv

Current status of space gravitational wave antenna DECIGO and B-DECIGO

Deci-hertz Interferometer Gravitational Wave Observatory (DECIGO) is the future Japanese space mission with a frequency band of 0.1 Hz to 10 Hz. DECIGO aims at the detection of primordial gravitational waves, which could be produced during the inflationary period right after the birth of the universe. There are many other scientific objectives of DECIGO, including the direct measurement of the acceleration of the expansion of the universe, and reliable and accurate predictions of the timing and locations of neutron star/black hole binary coalescences. DECIGO consists of four clusters of observatories placed in the heliocentric orbit. Each cluster consists of three spacecraft, which form three Fabry-Perot Michelson interferometers with an arm length of 1,000 km. Three clusters of DECIGO will be placed far from each other, and the fourth cluster will be placed in the same position as one of the three clusters to obtain the correlation signals for the detection of the primordial gravitational waves. We plan to launch B-DECIGO, which is a scientific pathfinder of DECIGO, before DECIGO in the 2030s to demonstrate the technologies required for DECIGO, as well as to obtain fruitful scientific results to further expand the multi-messenger astronomy.

preprint2020arXiv

Inverse reconstruction of jet structure from off-axis gamma-ray burst afterglows

The gravitational wave event GW170817 and the slowly-rising afterglows of short gamma-ray burst GRB 170817A clearly suggest that the GRB jet has an angular structure. However the actual jet structure remains unclear as different authors give different structures. We formulate a novel method to inversely reconstruct the jet structure from off-axis GRB afterglows, without assuming any functional form of the structure in contrast to the previous studies. The jet structure is uniquely determined from the rising part of a light curve for a given parameter set by integrating an ordinary differential equation, which is derived from the standard theory of GRB afterglows. Applying to GRB 170817A, we discover that a non-trivial hollow-cone jet is consistent with the observed afterglows, as well as Gaussian and power-law jets within errors, which implies the Blandford-Znajek mechanism or an ejecta-jet interaction. The current observations only constrain the jet core, not in principle the outer jet structure around the line of sight. More precise and high-cadence observations with our inversion method will fix the jet structure, providing a clue to the jet formation and propagation.

preprint2019arXiv

Jet Propagation in Neutron Star Mergers and GW170817

The gravitational wave event from the binary neutron star (BNS) merger GW170817 and the following multi-messenger observations present strong evidence for i) merger ejecta expanding with substantial velocities and ii) a relativistic jet which had to propagate through the merger ejecta. The ejecta's expansion velocity is not negligible for the jet head motion, which is a fundamental difference from the other systems like collapsars and active galactic nuclei. Here we present an analytic model of the jet propagation in an expanding medium. In particular, we notice a new term in the expression of the breakout time and velocity. In parallel, we perform a series of over a hundred 2D numerical simulations of jet propagation. The BNS merger ejecta is prepared based on numerical relativity simulations of a BNS merger with the highest-resolution to date. We show that our analytic results agree with numerical simulations over a wide parameter space. Then we apply our analytic model to GW170817, and obtain two solid constraints on: i) the central engine luminosity as $L_{iso,0} \sim 3\times10^{49}-2.5\times10^{52}$ erg s$^{-1}$, and on ii) the delay time between the merger and engine activation $t_0-t_m < 1.3$ s. The engine power implies that the apparently-faint \textit{short} gamma-ray burst (\textit{s}GRB) \textit{s}GRB 170817A is similar to typical \textit{s}GRBs if observed on-axis.

preprint2016arXiv

Are Ultra-Long Gamma-Ray Bursts Caused by Blue Supergiant Collapsars, Newborn Magnetars, or White Dwarf Tidal Disruption Events?

Ultra-long gamma-ray bursts (ulGRBs) are a new population of GRBs with extreme durations of $\sim 10^{4}$ s. Leading candidates for their origin are blue supergiant collapsars, magnetars, and white dwarf tidal disruption events (WD-TDEs) caused by massive black holes (BHs). Recent observations of supernova-like (SN-like) bumps associated with ulGRBs challenged both the WD-TDE and the blue supergiant models because of the detection of SNe and the absence of hydrogen lines, respectively. We propose that WD-TDEs can accommodate the observed SN-like bumps if the fallback WD matter releases energy into the unbound WD ejecta. The observed ejecta energy, luminosity, and velocity are explained by the gravitational energy, Eddington luminosity, and escape velocity of the formed accretion disk, respectively. We also show that the observed X-rays can ionize the ejecta, eliminating lines. The SN-like light curves (SN 2011kl) for the ulGRB 111209A are consistent with all three models, although a magnetar model is unnatural because the spin-down time required to power the SN-like bump is a hundred times longer than the GRB. Our results imply that TDEs are a possible energy source for SN-like events in general and for ulGRBs in particular.

preprint2016arXiv

Cosmic-Ray Hardenings in the Light of AMS-02

Recent precise observations of cosmic rays (CRs) by AMS-02 experiment clearly show (1) harder spectra of helium and carbon compared to protons by $\propto R^{0.08}$, and (2) concave breaks in proton and helium spectra at a rigidity $R \sim 300$ GV. In particular the helium and carbon spectra are exactly similar, pointing to the same acceleration site. We examine possible interpretations of these features and identify a chemically enriched region, that is, superbubbles as the most probable origin of Galactic CRs in high rigidity $R>30$ GV. The similar spectra of CR carbon and helium further suggest that the CRs with $R>30$ GV originate from the supernova ejecta in the superbubble core, mixed with comparable or less amount of interstellar medium. We predict similar spectra for heavy nuclei.

preprint2016arXiv

Jet-powered supernovae of $\sim 10^5\,M_{\odot}$ population III stars are observable by $Euclid$, $WFIRST$, $WISH$, and $JWST$

Supermassive black holes observed at high redshift $z\gtrsim6$ could grow from direct collapse black holes (DCBHs) with mass $\sim10^5\,M_{\odot}$, which result from the collapse of supermassive stars (SMSs). If a relativistic jet is launched from a DCBH, it can break out of the collapsing SMS and produce a gamma-ray burst (GRB). Although most of the GRB jets are off-axis from our line of sight, we show that the energy injected from the jet into a cocoon is huge $\sim10^{55-56}\,{\rm{erg}}$, so that the cocoon fireball is observed as ultra-luminous supernovae of $\sim10^{45-46}\rm{\,erg\,s^{-1}}$ for $\sim5000 [(1+z)/16] \rm{\,days}$. They are detectable by the future telescopes with near infrared bands, such as, $Euclid$, $WFIRST$, $WISH$, and $JWST$ up to $z\sim20$ and $\sim 100$ events per year, providing a direct evidence of the DCBH scenario.

preprint2016arXiv

The unreasonable weakness of r-process cosmic rays in the neutron-star-merger nucleosynthesis scenario

We reach the robust conclusion that, by combining the observed cosmic rays of r-process elements with the fact that the velocity of the neutron-star-merger ejecta is much higher than that of the supernova ejecta, either (1) the reverse shock in the neutron-star-merger ejecta is a very inefficient accelerator that converts less than 0.003% of the ejecta kinetic energy to the cosmic-ray energy or (2) the neutron star merger is not the origin of the Galactic r-process elements. We also find that the acceleration efficiency should be less than 0.1% for the reverse shock of the supernova ejecta with the observed cosmic rays lighter than the iron.

preprint2015arXiv

130GeV gamma-ray line through axion conversion

We apply the axion-photon conversion mechanism to the 130 GeV $γ$-ray line observed by the Fermi satellite. Near the Galactic center, some astrophysical sources and/or particle dark matter can produce energetic axions (or axionlike particles), and the axions convert to $γ$ rays in Galactic magnetic fields along their flight to the Earth. Since continuum $γ$-ray and antiproton productions are sufficiently suppressed in axion production, the scenario fits the 130 GeV $γ$-ray line without conflicting with cosmic ray measurements. We derive the axion production cross section and the decay rate of dark matter to fit the $γ$-ray excess as functions of axion parameters. In the scenario, the $γ$-ray spatial distributions depend on both the dark matter profile and the magnetic field configuration, which will be tested by future $γ$-ray observations, e.g., H.E.S.S. II, CTA, and GAMMA-400. As an illustrative example, we study realistic supersymmetric axion models, and show the favored parameters that nicely fit the $γ$-ray excess.

preprint2015arXiv

Can direct collapse black holes launch gamma-ray bursts and grow to supermassive black holes?

The existence of black holes (BHs) of mass ~ 10^{9} M_sun at z > 6 is a big puzzle in astrophysics because even optimistic estimates of the accretion time are insufficient for stellar mass BHs of ~ 10 M_sun to grow into such supermassive BHs. A resolution of this puzzle might be the direct collapse of supermassive stars with mass M ~ 10^{5} M_sun into massive seed BHs. We find that if a jet is launched from the accretion disk around the central BH, the jet can break out the star because of the structure of the radiation pressure-dominated envelope. Such ultra-long gamma-ray bursts with duration of ~ 10^{4} - 10^{6} s and flux of 10^{-11} - 10^{-8} erg s^{-1} cm^{-2} could be detectable by Swift. We estimate an event rate of < 1yr^{-1}. The total explosion energy is > 10^{55} - 10^{56} erg. The resulting negative feedback delays the growth of the remnant BH by about 70 Myr or evacuates the host galaxy completely.

preprint2015arXiv

Can we explain AMS-02 antiproton and positron excesses simultaneously by nearby supernovae without pulsars nor dark matter?

We explain the excess of the antiproton fraction recently reported by the AMS-02 experiment by considering collisions between cosmic-ray protons accelerated by a local supernova remnant (SNR) and the surrounding dense cloud. The same "$pp$-collisions" provide the right ratio of daughter particles to fit the observed positron excess simultaneously in the natural model parameters. The supernova happened in relatively lower metalicity than the major cosmic-ray sources. The cutoff energy of electrons marks the supernova age of $\sim 10^{5}$ years, while the antiproton excess may extend to higher energy. Both antiproton and positron fluxes are completely consistent with our predictions in Fujita, Kohri, Yamasaki and Ioka (2009).

preprint2015arXiv

Dynamical mass ejection from black hole-neutron star binaries

We investigate properties of material ejected dynamically in the merger of black hole-neutron star binaries by numerical-relativity simulations. We systematically study the dependence of ejecta properties on the mass ratio of the binary, spin of the black hole, and equation of state of the neutron-star matter. Dynamical mass ejection is driven primarily by tidal torque, and the ejecta is much more anisotropic than that from binary neutron star mergers. In particular, the dynamical ejecta is concentrated around the orbital plane with a half opening angle of 10--20deg and often sweeps out only a half of the plane. The ejecta mass can be as large as ~0.1M_sun, and the velocity is subrelativistic with ~0.2--0.3c for typical cases. The ratio of the ejecta mass to the bound mass (disk and fallback components) is larger, and the ejecta velocity is larger, for larger values of the binary mass ratio, i.e., for larger values of the black-hole mass. The remnant black hole-disk system receives a kick velocity of O(100)km/s due to the ejecta linear momentum, and this easily dominates the kick velocity due to gravitational radiation. Structures of postmerger material, velocity distribution of the dynamical ejecta, fallback rates, and gravitational waves are also investigated. We also discuss the effect of ejecta anisotropy on electromagnetic counterparts, specifically a macronova/kilonova and synchrotron radio emission, developing analytic models.

preprint2015arXiv

Energy Sources and Light Curves of Macronovae

A macronova (kilonova) was discovered with a short gamma-ray burst, GRB 130603B, which is widely believed to be powered by the radioactivity of $r$-process elements synthesized in the ejecta of a neutron star binary merger. As an alternative, we propose that macronovae are energized by the central engine, i.e., a black hole or neutron star, and the injected energy is emitted after the adiabatic expansion of ejecta. This engine model is motivated by extended emission of short GRBs. In order to compare the theoretical models with observations, we develop analytical formulae for the light curves of macronovae. The engine model allows a wider parameter range, especially smaller ejecta mass, and better fit to observations than the $r$-process model. Future observations of electromagnetic counterparts of gravitational waves should distinguish energy sources and constrain the activity of central engine and the $r$-process nucleosynthesis.

preprint2015arXiv

Isotropic Detectable X-ray Counterparts to Gravitational Waves from Neutron Star Binary Mergers

Neutron star binary mergers are strong sources of gravitational waves (GWs). Promising electromagnetic counterparts are short gamma-ray bursts (GRBs) but the emission is highly collimated. We propose that the scattering of the long-lasting plateau emission in short GRBs by the merger ejecta produces nearly isotropic emission for $\sim 10^4$ s with flux $10^{-13}-10^{-10}$ erg cm$^{-2}$ s$^{-1}$ at 100 Mpc in X-ray. This is detectable by Swift XRT and wide field X-ray detectors such as ISS-Lobster, Einstein Probe, eROSITA and WF-MAXI, which are desired by the infrared and optical follow-ups to localize and measure the distance to the host galaxy. The scattered X-rays obtain linear polarization, which correlates with the jet direction, X-ray luminosity and GW polarizations. The activity of plateau emission is also a natural energy source of a macronova (or kilonova) detected in short GRB 130603B without the $r$-process radioactivity.

preprint2015arXiv

Long-Lasting Black-Hole Jets in Short Gamma-Ray Bursts

Whether a short gamma-ray burst (GRB) is caused by a black hole (BH) or a neutron star (NS) after the merger of a NS binary is a crucial problem. We propose a BH model that explains short GRBs with long-lasting activities such as extended emission and plateau emission up to $\sim10000$ s. To extract the BH rotational energy, the topological evolution of the magnetic field should accompany the mass ejection, mass fallback, and magnetic field reconnection. The observations suggest the magnetic field decay from $\sim10^{14}$ G to $\sim10^{13} - 10^{11}$ G at the BH, bounded below by the pre-merger strength and kept constant while the luminosity is constant, and the fallback mass of $\sim10^{-4} - 10^{-2} M_{\odot}$, comparable to the ejecta mass implied by the macronova (or kilonova) in GRB 130603B. The BH model has implications for gravitational waves and the equation of state of NS matter.

preprint2015arXiv

Neutrino Flavor Ratios Modified by Cosmic Ray Secondary Acceleration

Acceleration of $π$'s and $μ$'s modifies the flavor ratio at Earth (at astrophysical sources) of neutrinos produced by $π$ decay, $ν_e:ν_μ:ν_τ$, from $1:1:1$ ($1:2:0$) to $1:1.8:1.8$ ($0:1:0$) at high energy, because $π$'s decay more than $μ$'s during secondary-acceleration. The neutrino spectrum accompanies a flat excess, differently from the case of energy losses. With the flavor spectra, we can probe timescales of cosmic-ray acceleration and shock dynamics. We obtain general solutions of convection-diffusion equations and apply to gamma-ray bursts, which may have the flavor modification at around PeV -- EeV detectable by IceCube and next-generation experiments.

preprint2015arXiv

Testing the millisecond pulsar scenario of the Galactic center gamma-ray excess with very high energy gamma-rays

The recent analyses of the Fermi Large Area Telescope data show an extended GeV $γ$-ray excess on top of the expected diffuse background in the Galactic center region, which can be explained with annihilating dark matter or a population of millisecond pulsars (MSPs). We propose to observe the very high energy $γ$-rays for distinguishing the MSP scenario from the dark matter scenario. The GeV $γ$-ray MSPs should release most energy to the relativistic $e^{\pm}$ wind, which will diffuse in the Galaxy and radiate TeV $γ$-rays through inverse Compton scattering and bremsstrahlung processes. By calculating the spectrum and spatial distribution, we show that such emission is detectable with the next generation very high energy $γ$-ray observatory, the Cherenkov Telescope Array (CTA), under reasonable model parameters. It is essential to search for the multi-wavelength counterparts to the GeV $γ$-ray excess for solving this mystery in the high energy universe.

preprint2015arXiv

X-ray-powered macronovae

A macronova (or kilonova) was observed as an infrared excess several days after short gamma-ray burst, GRB 130603B. Although the $r$-process radioactivity is widely discussed as an energy source, it requires huge mass of ejecta from a neutron star (NS) binary merger. We propose a new model that the X-ray excess gives rise to the simultaneously observed infrared excess via thermal re-emission and explore what constraints this would place on the mass and velocity of the ejecta. This X-ray-powered model explains both the X-ray and infrared excesses with a single energy source by the central engine like a black hole, and allows for broader parameter region, in particular smaller ejecta mass $\sim10^{-3}-10^{-2}M_{\odot}$ with iron mixed as suggested by general relativistic simulations for typical NS-NS mergers, than the previous models. We also discuss the other macronova candidates in GRB 060614 and GRB 080503, and implications for the search of electromagnetic counterparts to gravitational waves.

preprint2015arXiv

X-raying extended emission and rapid decay of short gamma-ray bursts

Extended emission is a mystery in short gamma-ray bursts (SGRBs). By making time resolved spectral analyses of brightest nine events observed by ${\it Swift}$ XRT, we obviously classify the early X-ray emission of SGRBs into two types. One is the extended emission with exponentially rapid decay, which shows significant spectral softening during hundreds seconds since the SGRB trigger and is also detected by ${\it Swift}$-BAT. The other is a dim afterglow only showing power-law decay over $10^4$ s. The correlations between the temporal decay and spectral indices of the extended emissions are inconsistent with the $α$-$β$ correlation expected for the high-latitude curvature emission from a uniform jet. The observed too-rapid decay suggests the emission from a photosphere or a patchy surface, and manifests the stopping central engine via such as magnetic reconnection at the black hole.

preprint2014arXiv

Dust formation in macronovae

We examine dust formation in macronovae (as known as kilonovae), which are the bright ejecta of neutron star binary mergers and one of the leading sites of r-process nucleosynthesis. We find that dust grains of r-process elements are difficult to form because of the low number density of the r-process atoms, while carbon or elements lighter than irons can condense into dust if they are abundant, in light of the first macronova candidate associated with GRB 130603B. Dust grains absorb emission from ejecta with opacity even greater than that of the r-process elements, and re-emit photons at infrared wavelengths. Such dust emission can potentially account for the macronova without r-process nucleosynthesis as an alternative model. This dust scenario predicts a more featureless spectrum than the r-process model and day-scale optical-to-ultraviolet emission.

preprint2014arXiv

IceCube PeV-EeV Neutrinos and Secret Interactions of Neutrinos

We show that the PeV neutrinos detected by IceCube put unique constraints on "secret" interactions of neutrinos with the cosmic neutrino background (C$ν$B). The coupling must be $g <0.03$ for the mediating boson mass $m_{X} \lesssim 2$ MeV, $g/m_{X} < 5$ GeV$^{-1}$ for $m_{X} \gtrsim 20$ MeV, and $g/m_{X} < 0.07$ GeV$^{-1}$ in between. We also investigate the possibility that neutrino cascades degrade high-energy neutrinos to PeV energies by upgrading C$ν$B where the energy flux of PeV neutrinos can coincide with the Waxman-Bahcall bound or the cosmogenic neutrino flux for protons, thanks to energy conservation. However a large coupling is required, which is disfavored by laboratory decay constraints. The suppression of PeV-EeV neutrinos is a testable prediction for the Askaryan Radio Array.

preprint2014arXiv

Jet collimation in the ejecta of double neutron star merger: new canonical picture of short gamma-ray bursts

The observations of jet breaks in the afterglows of short gamma-ray bursts (SGRBs) indicate that the jet has a small opening angle of < 10°. The collimation mechanism of the jet is a longstanding theoretical problem. We numerically analyze the jet propagation in the material ejected by double neutron star merger, and demonstrate that if the ejecta mass is > 10^{-2} M_{sun}, the jet is well confined by the cocoon and emerges from the ejecta with the required collimation angle. Our results also suggest that there are some populations of choked (failed) SGRBs or low-luminous new types of event. By constructing a model for SGRB 130603B, which is associated with the first kilonova/macronova can- didate, we infer that the equation-of-state of neutron stars would be soft enough to provide sufficient ejecta to collimate the jet, if this event was associated with a double neutron star merger.

preprint2013arXiv

Cosmological Fast Radio Bursts from Binary White Dwarf Mergers

Recently, Thornton et al. reported the detection of four fast radio bursts (FRBs). The dispersion measures indicate that the sources of these FRBs are at cosmological distance. Given the large full sky event rate ~ 10^4 sky^-1 day^-1, the FRBs are a promising target for multi-messenger astronomy. Here we propose double degenerate, binary white-dwarf (WD) mergers as the source of FRBs, which are produced by coherent emission from the polar region of a rapidly rotating, magnetized massive WD formed after the merger. The basic characteristics of the FRBs, such as the energetics, emission duration and event rate, can be consistently explained in this scenario. As a result, we predict that some FRBs can accompany type Ia supernovae (SNe Ia) or X-ray debris disks. Simultaneous detection could test our scenario and probe the progenitors of SNe Ia, and moreover would provide a novel constraint on the cosmological parameters. We strongly encourage future SN and X-ray surveys that follow up FRBs.

preprint2013arXiv

High-Energy Radiation from Remnants of Neutron Star Binary Mergers

We study high-energy emission from the mergers of neutron star binaries as electromagnetic counterparts to gravitational waves aside from short gamma-ray bursts. The mergers entail significant mass ejection, which interacts with the surrounding medium to produce similar but brighter remnants than supernova remnants in a few years. We show that electrons accelerated in the remnants can produce synchrotron radiation in X-rays detectable at $\sim 100$ Mpc by current generation telescopes and inverse Compton emission in gamma rays detectable by the \emph{Fermi} Large Area Telescopes and the Cherenkov Telescope Array under favorable conditions. The remnants may have already appeared in high-energy surveys such as the Monitor of All-sky X-ray Image and the \emph{Fermi} Large Area Telescope as unidentified sources. We also suggest that the merger remnants could be the origin of ultra-high-energy cosmic rays beyond the knee energy, $\sim 10^{15}$ eV, in the cosmic-ray spectrum.

preprint2013arXiv

Opening Angles of Collapsar Jets

We investigate the jet propagation and breakout from the stellar progenitor for gamma-ray burst (GRB) collapsars by performing two-dimensional relativistic hydrodynamic simulations and analytical modeling. We find that the jet opening angle is given by $θ_j \sim 1/5 Γ_{0}$, and infer the initial Lorentz factor of the jet at the central engine, $Γ_0$, is a few for existing observations of $θ_j$. The jet keeps the Lorentz factor low inside the star by converging cylindrically via collimation shocks under the cocoon pressure, and accelerates at jet breakout before the free expansion to a hollow-cone structure. In this new picture the GRB duration is determined by the sound crossing time of the cocoon, after which the opening angle widens, reducing the apparent luminosity. Some bursts violating the maximum opening angle $θ_{j,\max}\sim 1/5 \sim 12^{\circ}$ imply the existence of a baryon-rich sheath or a long-acting jet. We can explain the slopes in both Amati and Yonetoku spectral relations using an off-centered photosphere model, if we make only one assumption that the total jet luminosity is proportional to the initial Lorentz factor of the jet. We also numerically calibrate the pre-breakout model (Bromberg et al.) for later use.

preprint2013arXiv

TeV-PeV Neutrinos from Low-Power Gamma-Ray Burst Jets inside Stars

We study high-energy neutrino production in collimated jets inside progenitors of gamma-ray bursts (GRBs) and supernovae, considering both collimation and internal shocks. We obtain simple, useful constraints, using the often overlooked point that shock acceleration of particles is ineffective at radiation-mediated shocks. Classical GRBs may be too powerful to produce high-energy neutrinos inside stars, which is consistent with IceCube nondetections. We find that ultralong GRBs avoid such constraints and detecting the TeV signal will support giant progenitors. Predictions for low-power GRB classes including low-luminosity GRBs can be consistent with the astrophysical neutrino background IceCube may detect, with a spectral steepening around PeV. The models can be tested with future GRB monitors.

preprint2012arXiv

Escape of cosmic-ray electrons from supernova remnants

We investigate escape of cosmic ray (CR) electrons from a supernova remnant (SNR) to interstellar space. We show that CR electrons escape in order from high energies to low energies like CR nuclei, while the escape starts later than the beginning of the Sedov phase at an SNR age of 10^3 - 7*10^3 yrs and the maximum energy of runaway CR electrons is below the knee about 0.3 - 50 TeV because unlike CR nuclei, CR electrons lose their energy due to synchrotron radiation. Highest energy CR electrons will be directly probed by AMS-02, CALET, CTA and LHAASO experiments, or have been already detected by H.E.S.S. and MAGIC as a cutoff in the CR electron spectrum. Furthermore, we also calculate the spatial distribution of runaway CR electrons and their radiation spectra around SNRs. Contrary to common belief, maximum-energy photons of synchrotron radiation around 1 keV are emitted by runaway CR electrons which have been caught up by the shock. Inverse Compton scattering by runaway CR electrons can dominate the gamma-ray emission from runaway CR nuclei via pion decay, and both are detectable by CTA and LHAASO as clues to the CR origin and the amplification of magnetic fluctuations around the SNR. We also discuss middle-aged and/or old SNRs as unidentified very-high-energy gamma-ray sources.

preprint2012arXiv

Gamma-ray flare and absorption in Crab Nebula: Lovely TeV-PeV astrophysics

We spectrally fit the GeV gamma-ray flares recently-observed in the Crab Nebula by considering a small blob Lorentz-boosted towards us. We point out that the corresponding inverse-Compton flare at TeV--PeV region is more enhanced than synchrotron by a Lorentz factor square \sim Γ^2, which is already excluding Γ\gtrsim 200 and will be detected by future TeV - PeV observatories, CTA, Tibet AS + MD and LHAASO for Γ\gtrsim 30. We also show that PeV photons emitted from the Crab Nebula are absorbed by Cosmic Microwave Background radiation through electron-positron pair creation.

preprint2012arXiv

Population III Gamma-Ray Bursts and Breakout Criteria for Accretion-Powered Jets

We investigate the propagation of accretion-powered jets in various types of massive stars such as Wolf-Rayet stars, light Population III (Pop III) stars, and massive Pop III stars, all of which are the progenitor candidates of Gamma-Ray Bursts (GRBs). We perform two dimensional axisymmetric simulations of relativistic hydrodynamics taking into account both the envelope collapse and the jet propagation (i.e., the negative feedback of the jet on the accretion). Based on our hydrodynamic simulations, we show for the first time that the accretion-powered jet can potentially break out relativistically from the outer layers of Pop III progenitors. In our simulations, the accretion rate is estimated by the mass flux going through the inner boundary, and the jet is injected with a fixed accretion-to-jet conversion efficiency $η$. By varying the efficiency $η$ and opening angle $θ_{op}$ for more than 40 models, we find that the jet can make a relativistic breakout from all types of progenitors for GRBs if a simple condition $η\gtrsim 10^{-4} (θ_{op}/8^{\circ})^2$ is satisfied, which is consistent with analytical estimates. Otherwise no explosion or some failed spherical explosions occur.

preprint2012arXiv

Prospects for Detecting Gamma-Ray Bursts at Very High Energies with the Cherenkov Telescope Array

We discuss the prospects for the detection of gamma-ray bursts (GRBs) by the Cherenkov Telescope Array (CTA), the next generation, ground-based facility of imaging atmospheric Cherenkov telescopes (IACTs) operating above a few tens of GeV. By virtue of its fast slewing capabilities, the lower energy threshold compared to current IACTs, and the much larger effective area compared to satellite instruments, CTA can measure the spectra and variability of GRBs with excellent photon statistics at multi-GeV energies. Employing a model of the GRB population whose properties are broadly consistent with observations by the Gamma-ray Burst Monitor (GBM) and Large Area Telescope (LAT) onboard Fermi, we simulate follow-up observations of GRBs with the Large Size Telescopes (LSTs), the component of CTA with the fastest slew speed and the best sensitivity at energies below a few hundred GeV. For our fiducial assumptions, we foresee that the LSTs can detect ~0.1 GRBs per year during the prompt phase and ~0.5 per year in the afterglow phase, considering only one array site and both GBM and the Space-based multi-band astronomical Variable Object Monitor (SVOM) as the alert instruments. The detection rates can be enhanced by a factor of about 5 and 6 for the prompt emission and the afterglow, respectively, assuming two array sites with the same sensitivity and that the GBM localization error can be reduced to less than 1 deg. The expected distribution of redshift and photon counts are presented, showing that despite the modest event rate, hundreds or more multi-GeV photons can be anticipated from a single burst once they are detected. We also study how the detection rate depends on the intrinsic GRB properties and the delay time between the burst trigger and the follow-up observation.

preprint2012arXiv

The survival of nuclei in jets associated with core-collapse supernovae

Heavy nuclei such as nickel-56 are synthesized in a wide range of core-collapse supernovae (CCSN), including energetic supernovae associated with gamma-ray bursts (GRBs). Recent studies suggest that jet-like outflows are a common feature of CCSN. These outflows may entrain synthesized nuclei at launch or during propagation, and provide interesting multi-messenger signals including heavy ultra-high energy cosmic rays. Here, we investigate the destruction processes of nuclei during crossing from the stellar material into the jet material via a cocoon, and during propagation after being successfully loaded into the jet. We find that nuclei can survive for a range of jet parameters because collisional cooling is faster than spallation. While canonical high-luminosity GRB jets may contain nuclei, magnetic dominated models or low-luminosity jets with small bulk Lorentz factors are more favorable for having a more significant heavy nuclei component.

preprint2012arXiv

Upper Limit on the Cosmological Gamma-ray Background

We show that the current extragalactic gamma-ray background (EGB) measurement below 100 GeV sets an upper limit on EGB itself at very high energy (VHE) above 100 GeV. The limit is conservative for the electromagnetic cascade emission from VHE EGB interacting with the cosmic microwave-to-optical background radiation not to exceed the current EGB measurement. The cascade component fits the measured VHE EGB spectrum rather well. However, once we add the contribution from known source classes, the Fermi VHE EGB observation exceeds or even violates the limit, which is approximated as E^2dN/dE < 4.5x10^-5 (E/100 GeV)^-0.7 MeV/cm^2/s/sr. The upper limit above 100 GeV is useful in the future to probe the EGB origin and the new physics like axion-like particles and Lorentz-invariance violation.

preprint2011arXiv

Afterglow of binary neutron star merger

The merger of two neutron stars results often in a rapidly and differentially rotating hypermassive neutron star (HMNS). We show by numerical-relativity simulation that the magnetic-field profile around such HMNS is dynamically varied during its subsequent evolution, and as a result, electromagnetic radiation with a large luminosity ~ 0.1 B^2 R^3 Omega is emitted with baryon (B, R, and Omega are poloidal magnetic-field strength at stellar surface, stellar radius, and angular velocity of a HMNS). The predicted luminosity of electromagnetic radiation, which is primarily emitted along the magnetic-dipole direction, is ~ 10^{47} (B/10^{13} G)^2(R/10 km)^3(Omega/10^4 rad/s) ergs/s, that is comparable to the luminosity of quasars.

preprint2011arXiv

Cosmic Ray Helium Hardening

Recent observations by CREAM and ATIC-2 experiments suggest that (1) the spectrum of cosmic ray (CR) helium is harder than that of CR proton below the knee 10^15 eV and (2) all CR spectra become hard at > 10^11 eV/n. We propose a new picture that higher energy CRs are generated in more helium-rich region to explain the hardening (1) without introducing different sources for CR helium. The helium to proton ratio at ~100 TeV exceeds the Big Bang abundance Y=0.25 by several times, and the different spectrum is not reproduced within the diffusive shock acceleration theory. We argue that CRs are produced in the chemically enriched region, such as a superbubble, and the outward-decreasing abundance naturally leads to the hard spectrum of CR helium if CRs escape from the supernova remnant (SNR) shock in an energy-dependent way. We provide a simple analytical spectrum that also fits well the hardening (2) because of the decreasing Mach number in the hot superbubble with ~ 10^6 K. Our model predicts hard and concave spectra for heavier CR elements.

preprint2011arXiv

Gamma-Ray Burst without Baryonic and Magnetic Load?

We show that, contrary to common belief, internal shocks can arise in an accelerating radiation-dominated jet if it is confined even weakly to a converging opening angle because the acceleration declines. The radiation-dominated internal shock (RDIS) enables a very efficient yet highly nonthermal emission by Fermi-like photon acceleration, keeping the electron-positron ($e^{\pm}$) pair photosphere and inertia up to a high Lorentz factor >1000. In gamma-ray bursts (GRBs), a weak confinement would persist beyond the progenitor star or surrounding matter because of the fast cocoon accompanying the breakout jet. The simplest model predicts few high-energy cosmic rays and neutrinos, and a correlation between the early afterglow and the GeV-TeV prompt emission. The central engine allows a less fine-tuned baryon load than previously thought, even including pure-leptonic unmagnetized outflows.

preprint2011arXiv

Magnetar Asteroseismology with Long-Term Gravitational Waves

Magnetic flares and induced oscillations of magnetars (super-magnetized neutron stars) are promising sources of gravitational waves (GWs). We suggest that the GW emission, if any, would last longer than the observed X-ray quasi-periodic oscillations (X-QPOs), calling for the longer-term GW analyses from a day to months than the current searches. Like the pulsar timing, the oscillation frequency would also evolve with time because of the decay or reconfiguration of magnetic field, which is crucial for the GW detection. With the observed GW frequency and its time-derivatives, we can probe the interior magnetic field strength and its evolution to open a new GW asteroseismology with the next generation interferometers like advanced LIGO, advanced VIRGO, LCGT and ET.

preprint2011arXiv

Population III.1 and III.2 gamma-Ray Bursts: Constraints on the event rate for future radio and X-ray surveys

We calculate the theoretical event rate of gamma-ray bursts (GRBs) from the collapse of massive first-generation (Population III; Pop III) stars. The Pop III GRBs could be super-energetic with the isotropic energy up to $E_{\rm iso} \gtrsim 10^{55-57}$ ergs, providing a unique probe of the high-redshift Universe. We consider both the so-called Pop III.1 stars (primordial) and Pop III.2 stars (primordial but affected by radiation from other stars). We employ a semi-analytical approach that considers inhomogeneous hydrogen reionization and chemical evolution of the intergalactic medium. We show that Pop III.2 GRBs occur more than 100 times more frequently than Pop III.1 GRBs, and thus should be suitable targets for future GRB missions. Interestingly, our optimistic model predicts an event rate that is already constrained by the current radio transient searches. We expect $\sim10-10^4$ radio afterglows above $\sim$ 0.3 mJy on the sky with $\sim 1$ year variability and mostly without GRBs (orphans), which are detectable by ALMA, EVLA, LOFAR, and SKA, while we expect to observe maximum of $N < 20$ GRBs per year integrated over at $z>6$ for Pop III.2 and $N < 0.08$ per year integrated over at $z>10$ for Pop III.1 with EXIST, and $N < 0.2$ for Pop III.2 GRBs per year integrated over at $z > 6$ with \textit{Swift}.

preprint2011arXiv

TeV Electron Spectrum for Probing Cosmic-Ray Escape from a Supernova Remnant

One of the most essential but uncertain processes for producing cosmic-rays (CRs) and their spectra is how accelerated particles escape into the interstellar space. We propose that the CR electron spectra at >~TeV energy can provide a direct probe of the CR escape complementary to the CR nuclei and gamma-rays. We calculate the electron spectra from a young pulsar embedded in the supernova remnant (SNR), like Vela, taking into account the energy-dependent CR escape. SNRs would accelerate and hence confine particles with energy up to 10^{15.5}eV. Only energetic particles can escape first, while the lower energy particles are confined and released later. Then the observed electron spectrum should have a low energy cutoff whose position marks the age of the pulsar/SNR. The low energy cutoff is observable in the $\gtrsim {\rm TeV}$ energy window, where other contaminating sources are expected to be few due to the fast cooling of electrons. The spectrum looks similar to a dark matter annihilation line if the low energy cutoff is close to the high energy intrinsic or cooling break. The future experiments such as CALET and CTA are capable of directly detecting the CR escape features toward revealing the origin of CRs.

preprint2011arXiv

Three-Dimensional Simulations of MHD Turbulence Behind Relativistic Shock Waves and Their Implications for GRBs

Relativistic astrophysical phenomena such as gamma-ray bursts (GRBs) and active galactic nuclei often require long-lived strong magnetic field. Here, we report on three-dimensional special-relativistic magnetohydrodynamic (MHD) simulations to explore the amplification and decay of macroscopic turbulence dynamo excited by the so-called Richtmyer-Meshkov instability (RMI; a Rayleigh-Taylor type instability). This instability is an inevitable outcome of interactions between shock and ambient density fluctuations. We find that the magnetic energy grows exponentially in a few eddy turnover times, and then, following the decay of kinetic turbulence, decays with a temporal power-law exponent of -0.7. The magnetic-energy fraction can reach $epsilon_B \sim$ 0.1 but depends on the initial magnetic field strength. We find that the magnetic energy grows by at least two orders of magnitude compared to the magnetic energy immediately behind the shock. This minimum degree of the amplification does not depend on the amplitude of the initial density fluctuations, while the growth timescale and the maximum magnetic energy depend on the degree of inhomogeneity in the density. The transition from Kolmogorov cascade to MHD critical balance cascade occurs at $\sim$ 1/10th the initial inhomogeneity scale, which limits the maximum synchrotron polarization to less than 2%. New results include the avoidance of electron cooling with RMI turbulence, the turbulent photosphere model via RMI, and the shallow decay of the early afterglow from RMI. We also performed a simulation of freely decaying turbulence with relativistic velocity dispersion. We find that relativistic turbulence begins to decay much faster than one eddy-turnover time because of fast shock dissipation, which does not support the relativistic turbulence model by Narayan & Kumar.

preprint2010arXiv

A Gamma-Ray Burst/Pulsar for Cosmic-Ray Positrons with a Dark Matter-like Spectrum

We propose that a nearby gamma-ray burst (GRB) or GRB-like (old, single and short-lived) pulsar/supernova remnant/microquasar about 10^{5-6} years ago may be responsible for the excesses of cosmic-ray positrons and electrons recently observed by the PAMELA, ATIC/PPB-BETS, Fermi and HESS experiments. We can reproduce the smooth Fermi/HESS spectra as well as the spiky ATIC/PPB-BETS spectra. The spectra have a sharp cutoff that is similar to the dark matter predictions, sometimes together with a line (not similar), since higher energy cosmic-rays cool faster where the cutoff/line energy marks the source age. A GRB-like astrophysical source is expected to have a small but finite spread in the cutoff/line as well as anisotropy in the cosmic-ray and diffuse gamma-ray flux, providing a method for the Fermi and future CALET experiments to discriminate between dark matter and astrophysical origins.

preprint2010arXiv

Can Gamma-Ray Burst Jets Break Out the First Stars?

We show that a relativistic gamma-ray burst (GRB) jet can potentially pierce the envelope of very massive first generation star (Population III; Pop III) by using the stellar density profile to estimate both the jet luminosity (via accretion) and its penetrability. The jet breakout is possible even if the Pop III star has a supergiant hydrogen envelope without mass loss, thanks to the long-lived powerful accretion of the envelope itself. While the Pop III GRB is estimated to be energetic, E_{gamma,iso} ~ 10^{55} erg, the supergiant envelope hides the initial bright phase into the cocoon component, leading to a GRB with a long duration ~ 1000(1+z) sec and an ordinary isotropic luminosity ~ 10^{52} erg s^{-1} (~ 10^{-9} erg cm^{-2} s^{-1} at redshift z ~ 20). The neutrino-annihilation is not effective for Pop III GRBs because of a low central temperature, while the magnetic mechanism is viable. We also derive analytic estimates of the breakout conditions, which are applicable to various progenitor models. The GRB luminosity and duration are found to be very sensitive to the core and envelope mass, providing possible probes of the first luminous objects at the end of the high redshift dark ages.

preprint2010arXiv

Can we probe the Lorentz factor of gamma-ray bursts from GeV-TeV spectra integrated over internal shocks?

We revisit the high-energy spectral cutoff originating from the electron-positron pair creation in the prompt phase of gamma-ray bursts (GRBs) with numerical and analytical calculations. We show that the conventional exponential and/or broken power law cutoff should be drastically modified to a shallower broken power-law in practical observations that integrate emissions from different internal shocks. Since the steepening is tiny for observations, this "smearing" effect can generally reduce the previous estimates of the Lorentz factor of the GRB outflows. We apply our formulation to GRB 080916C, recently detected by the Large Area Telescope detector on the Fermi satellite, and find that the minimum Lorentz factor can be ~600 (or even smaller values), which is below but consistent with the previous result of ~900. Observing the steepening energy (so-called "pair-break energy") is crucial to diagnose the Lorentz factor and/or the emission site in the future observations, especially current and future Cherenkov telescopes such as MAGIC, VERITAS, and CTA.

preprint2010arXiv

Very High Lorentz Factor Fireballs and Gamma-Ray Burst Spectra

Collisionless entrainment of the surrounding matter imports the relativistic baryon component in the Gamma-Ray Burst (GRB) fireball frame. We show that half the fireball energy can be transferred from radiation to the comoving hot motions of baryons under the photosphere. The yet baryon-poor fireball can reexpand to a very high Lorentz factor (VHLF) Γ~ 10^3-10^6 by its own relativistic collisionless pressure beyond the photosphere (so-called collisionless bulk acceleration), leading to internal and external shocks. A simple synchrotron emission from the VHLF internal shocks produces (i) the extra power-law spectral component with variability observed in the Fermi GeV bursts, up to the TeV range for the future Cherenkov Telescope Array (CTA), (ii) the GeV onset delay with a weak luminosity dependence t_{delay} ~ L^{-1/5}, and (iii) the spectral break of GRB 090926 by the synchrotron cooling break or the maximum synchrotron cutoff limited by the dynamical time, not by the e+- creation cutoff. The relativistic baryon component could also heat the photospheric thermal photons into the main GRB Band spectrum via pp, pγ(Bethe-Heitler and photomeson), and Coulomb thermalization processes. In this hot photosphere-internal-external shock model, we can predict the anticorrelation of ~TeV neutrinos and GeV gamma-rays, which may be detectable using IceCube. The spectral peak and luminosity (Yonetoku) relation is also reproduced if the progenitor stars are nearly identical. We also discuss the steep/shallow decay of early X-ray afterglows and short GRBs.

preprint2010arXiv

White Dwarf Pulsars as Possible Cosmic Ray Electron-Positron Factories

We suggest that white dwarf (WD) pulsars can compete with neutron star (NS) pulsars for producing the excesses of cosmic ray electrons/positrons observed by the PAMELA, ATIC/PPB-BETS, Fermi and HESS experiments. A merger of two WDs leads to a rapidly spinning WD with a rotational energy comparable to the NS case. The birth rate is also similar, providing the right energy budget for the cosmic ray electrons/positrons. Applying the NS theory, we suggest that the WD pulsars can in principle produce electrons/positrons up to 10 TeV. In contrast to the NS model, the adiabatic and radiative energy losses of electrons/positrons are negligible since their injection continues after the expansion of the pulsar wind nebula, and hence it is enough that a fraction 1% of WDs are magnetized as observed. The long activity also increases the number of nearby sources, which reduces the Poisson fluctuation in the flux. The WD pulsars could dominate the quickly cooling electrons/positrons above TeV energy as a second spectral bump or even surpass the NS pulsars in the observing energy range 10 GeV - 1 TeV, providing a background for the dark matter signals and a nice target for the future AMS-02, CALET and CTA experiment.

preprint2009arXiv

Cosmic-Ray Electron Excess from Pulsars is Spiky or Smooth?: Continuous and Multiple Electron/Positron injections

We investigate the observed spectrum of cosmic-ray electrons and positrons from astrophysical sources, especially pulsars, and the physical processes for making the spectrum spiky or smooth via continuous and multiple electron/positron injections. We find that (1) the average electron spectrum predicted from nearby pulsars are consistent with PAMELA, Fermi and H.E.S.S. data. However, the ATIC/PPB-BETS peak around 500GeV is hard to produce by the sum of multiple pulsar contributions and requires a single (or a few) energetic pulsar(s). (2) A continuous injection produces a broad peak and a high energy tail above the peak, which can constrain the source duration ($\lesssim 10^5$yr with the current data). (3) The H.E.S.S. data in the TeV range suggest that young sources with age less than $\sim 6 \times 10^4$yr are less energetic than $\sim 10^{48}{\rm erg}$. (4) We also expect a large dispersion in the TeV spectrum due to the small number of sources, that may cause the high energy cutoff inferred by H.E.S.S. and potentially provide a smoking-gun for the astrophysical origin. These spectral diagnostics can be refined in the near future by the CALET experiments to discriminate different astrophysical and dark matter origins.

preprint2009arXiv

High-energy emission as a test of the prior emission model for gamma-ray burst afterglows

We study high-energy gamma-ray afterglow emission from gamma-ray bursts (GRBs) in the prior emission model, which is proposed to explain the plateau phase of the X-ray afterglow. This model predicts the high-energy gamma-ray emission when the prompt GRB photons from the main flow are up-scattered by relativistic electrons accelerated at the external shock due to the prior flow. The expected spectrum has the peak of 10-100 GeV at around the end time of the plateau phase for typical GRBs, and high-energy gamma rays from nearby and/or energetic GRBs can be detected by current and future Cherenkov telescopes such as MAGIC, VERITAS, CTA, and possibly Fermi. Multi-wavelength observations by ground-based optical telescopes as well as Fermi and/or Swift sattelites are important to constrain the model. Such external inverse-Compton emission may even lead to GeV-TeV gamma-ray signals with the delay time of 10-100 s, only if the plateau phase is short-lived.

preprint2009arXiv

Hypernova and Gamma-Ray Burst Remnants as TeV Unidentified Sources

We investigate hypernova (hyper-energetic supernova) and gamma-ray burst (GRB) remnants in our Galaxy as TeV gamma-ray sources, particularly in the role of potential TeV unidentified sources, which have no clear counterpart at other wavelengths. We show that the observed bright sources in the TeV sky could be dominated by GRB/hypernova remnants, even though they are fewer than supernova remnants (SNRs). If this is the case, TeV SNRs are more extended (and more numerous) than deduced from current observations. In keeping with their role as cosmic ray accelerators, we discuss hadronic gamma-ray emission from pi^0 decay, from beta decay followed by inverse Compton emission, and propose a third, novel process of TeV gamma-ray emission arising from the decay of accelerated radioactive isotopes such as 56Co entrained by relativistic or semi-relativistic jets in GRBs/hypernovae. We discuss the relevant observational signatures which could discriminate between these three mechanisms.

preprint2009arXiv

Is the PAMELA anomaly caused by the supernova explosions near the Earth?

We show that the anomaly of the positron fraction observed by the PAMELA experiment can be attributed to recent supernova explosion(s) in a dense gas cloud (DC) near the Earth. Protons are accelerated around the supernova remnant (SNR). Electrons and positrons are created through hadronic interactions inside the DC. Their spectrum is harder than that of the background because the SNR spends much time in a radiative phase. Our scenario predicts that the anti-proton flux dominates that of the background for >~100 GeV. We compare the results with observations (Fermi, HESS, PPB-BETS, and ATIC).

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

High-energy cosmic-ray nuclei from high- and low-luminosity gamma-ray bursts and implications for multi-messenger astronomy

Gamma-ray bursts (GRBs) are one of the candidates of ultra-high-energy (around 10^18.5 eV) cosmic-ray (UHECR) sources. We investigate high-energy cosmic-ray acceleration including heavy nuclei in GRBs by using Geant 4, and discuss its various implications, taking both of high-luminosity (HL) and low-luminosity (LL) GRBs into account. This is because LL GRBs may also make a significant contribution to the observed UHECR flux if they form a distinct population. We show that not only protons but also heavier nuclei can be accelerated up to ultra-high energies in the internal, (external) reverse and forward shock models. We also show that the condition for ultra-high-energy heavy nuclei such as iron to survive is almost the same as that for about TeV gamma-rays to escape from the source and for high-energy neutrinos not to be much produced. The multi-messenger astronomy by neutrino and GeV-TeV gamma-ray telescopes such as IceCube and KM3Net, GLAST and MAGIC will be important to see whether GRBs can be accelerators of ultra-high-energy heavy nuclei. We also demonstrate expected spectra of high-energy neutrinos and gamma rays, and discuss their detectabilities. In addition, we discuss implictaions of the GRB-UHECR hypothesis. We point out, since the number densities of HL-GRBs and LL-GRBs are quite different, its detemination by UHECR observations is also important.

preprint2008arXiv

Time-evolution of Peak Energy and Luminosity Relation within Pulses for GRB 061007: Probing Fireball Dynamics

We perform a time-resolved spectral analysis of bright, long Gamma-ray burst GRB 061007 using Suzaku/WAM and Swift/BAT. Thanks to the large effective area of the WAM, we can investigate the time evolution of the spectral peak energy, Et_peak and the luminosity Lt_iso with 1-sec time resolution, and we find that luminosity Lt_iso with 1-sec time resolution, and we find that the time-resolved pulses also satisfy the Epeak-Liso relation, which was found for the time-averaged spectra of other bursts, suggesting the same physical conditions in each pulse. Furthermore, the initial rising phase of each pulse could be an outlier of this relation with higher Et_peak value by about factor 2. This difference could suggest that the fireball radius expands by a factor of 2-4 and/or bulk Lorentz factor of the fireball is decelerated by a factor of 4 during the initial phase, providing a new probe of the fireball dynamics in real time.

preprint2007arXiv

A multi band study of the optically dark GRB 051028

Observations were made of the optical afterglow of GRB 051028 with the Lulin observatory's 1.0 m telescope and the WIDGET robotic telescope system. R band photometric data points were obtained on 2005 October 28 (UT), or 0.095-0.180 days after the burst. There is a possible plateau in the optical light curve around 0.1 days after the burst; the light curve resembles optically bright afterglows (e.g. GRB 041006, GRB 050319, GRB060605) in shape of the light curve but not in brightness. The brightness of the GRB 051028 afterglow is 3 magnitudes fainter than that of one of the dark events, GRB 020124. Optically dark GRBs have been attributed to dust extinction within the host galaxy or high redshift. However, the spectrum analysis of the X-rays implies that there is no significant absorption by the host galaxy. Furthermore, according to the model theoretical calculation of the Ly$α$ absorption to find the limit of GRB 051028's redshift, the expected $R$ band absorption is not high enough to explain the darkness of the afterglow. While the present results disfavor either the high-redshift hypothesis or the high extinction scenario for optically dark bursts, they are consistent with the possibility that the brightness of the optical afterglow, intrinsically dark.

preprint2007arXiv

Closure Relations for Electron-Positron Pair-Signatures in Gamma-Ray Bursts

We present recipes to diagnose the fireball of gamma-ray bursts (GRBs) by combining observations of electron-positron pair-signatures (the pair-annihilation line and the cutoff energy due to the pair-creation process). Our recipes are largely model-independent and extract information even from the non-detection of either pair-signature. We evaluate physical quantities such as the Lorentz factor, optical depth and pair-to-baryon ratio, only from the observable quantities. In particular, we can test whether the prompt emission of GRBs comes from the pair/baryonic photosphere or not. The future-coming Gamma-Ray Large Area Space Telescope (GLAST) satellite will provide us with good chances to use our recipes by detecting or non-detecting pair-signatures.

preprint2007arXiv

Low-Luminosity GRB 060218: A Collapsar Jet from a Neutron Star, Leaving a Magnetar as a Remnant?

The gamma-ray burst (GRB) 060218 has ~ 10^5 times lower luminosity than typical long GRBs, and is associated with a supernova (SN). The radio afterglow displays no jet break, so that this burst might arise from a mildly-relativistic spherical outflow produced by the SN shock sweeping the stellar surface. Since this model is energetically difficult, we propose that the radio afterglow is produced by a non-relativistic phase of an initially collimated outflow (jet). Our jet model is supported by the detection of optical linear polarization in the SN component. We also show analytically that the jet can penetrate a progenitor star. We analyzed the observational data of the prompt emission of this burst and obtained a smooth power-law light curve which might last longer than 10^6s. This behavior contrasts with the long intermittent activities with the X-ray flares of typical GRBs, implying that the central engine of this burst is different from those of typical GRBs. This argument is consistent with the analysis of the SN component of this burst, which suggests that the progenitor star was less massive and collapsed to a neutron star instead of a black hole. The collimation-corrected event rate of such low-luminosity GRBs is estimated to be ~ 10 times higher than that of typical long GRBs, and they might form a different GRB population: low-luminosity GRBs are produced by mildly-relativistic jets from neutron stars at the collapses of massive stars, while typical long GRBs by highly-relativistic jets from black holes. We suggest that the central engine of GRB 060218 is a pulsar (or a magnetar) with the initial rotation period P_0 \~ 10ms and the magnetic field B ~ 10^16G. A giant flare from the magnetar might be observed in future.

preprint2007arXiv

Probing the Efficiency of Electron-Proton Coupling in Relativistic Collisionless Shocks through the Radio Polarimetry of Gamma-Ray Burst Afterglows

The late-time optical/radio afterglows of $γ$-ray bursts (GRBs) are believed to be synchrotron emission of electrons accelerated in relativistic collisionless shocks propagating in the ambient medium of the sources. However, the fraction $f$ of electrons that are coupled to protons and accelerated remains unclear and a large number of thermal electrons that are not coupled to protons may be left behind. If $f<1$, the true explosion energies of GRBs are $f^{-1}$ times larger than those commonly estimated with $f=1$. Thus the value of $f$ gives an important constraint on the nature of the central engine of GRBs and the physics of collisionless shocks. Although early-time radio observations can probe the thermal electrons, they are difficult at present. We show that the Faraday rotation effects of the thermal electrons may suppress the linear polarization of the afterglow at frequencies higher than the absorption frequency in the late time, if the magnetic field is ordered at least in parts, and that $f$ can be constrained through the observation of the effects. We find that those effects may be detected with late-time, > 1 day, polarimetry with ALMA for a burst occurring within 1 Gpc (i.e., z ~ 0.2), if $f \sim 10^{-1}$.

preprint2007arXiv

Unstable GRB photospheres and electron-positron annihilation lines

We propose an emission mechanism of prompt gamma-ray bursts (GRBs) that can reproduce the observed non-thermal spectra with high radiative efficiencies, >50%. Internal dissipation below a photosphere can create a radiation-dominated thermal fireball. If electron-positron pairs outnumber protons, radiative acceleration of pairs drives the two-stream instabilities between pairs and protons, leading to the ``proton sedimentation'' in the accelerating pair frame. Pairs are continuously shock heated by proton clumps, scattering the thermal photons into the broken power-law shape, with the non-thermal energy that is comparable to the proton kinetic energy, consistent with observations. Pair photospheres become unstable around the radius of the progenitor star where strong thermalization occurs, if parameters satisfy the observed spectral (Yonetoku) relation. Pair annihilation lines are predicted above continua, which could be verified by GLAST.

preprint2006arXiv

High Energy Neutrinos and Cosmic-Rays from Low-Luminosity Gamma-Ray Bursts?

The recently discovered gamma-ray burst (GRB) GRB 060218/SN 2006aj is classified as an X-ray Flash with very long duration driven possibly by a neutron star. Since GRB 060218 is very near about 140 Mpc and very dim, one-year observation by Swift suggests that the rate of GRB 060218-like events might be very high so that such low luminosity GRBs (LL-GRBs) might form a different population from the cosmological high luminosity GRBs (HL-GRBs). We found that the high energy neutrino background from LL-GRBs could be comparable with that from HL-GRBs. If each neutrino event is detected by IceCube, later optical-infrared follow-up observations such as by Subaru and HST have possibilities to identify a Type Ibc supernova associated with LL-GRBs, even if gamma- and X-rays are not observed by Swift. This is in a sense a new window from neutrino astronomy, which might enable us to confirm the existence of LL-GRBs and to obtain information about their rate and origin. We also argue LL-GRBs as high energy gamma-ray and cosmic-ray sources.

preprint2006arXiv

Tail Emission of Prompt Gamma-Ray Burst Jets

Tail emission of the prompt gamma-ray burst (GRB) is discussed using a multiple emitting sub-shell (inhomogeneous jet, sub-jets or mini-jets) model, where the whole GRB jet consists of many emitting sub-shells. One may expect that such a jet with angular inhomogeneity should produce spiky tail emission. However, we found that the tail is not spiky but is decaying roughly monotonically. The global decay slope of the tail is not so much affected by the local angular inhomogeneity but affected by the global sub-shell energy distribution. The fact that steepening breaks appeared in some events prefers the structured jets. If the angular size of the emitting sub-shell is around 0.01-0.02 rad, some bumps or fluctuations appear in the tail emission observed frequently in long GRBs. If the parameter differences of sub-shell properties are large, the tail has frequent changes of the temporal slope observed in a few bursts. Therefore, the multiple emitting sub-shell model has the advantage of explaining the small-scale structure in the observed rapid decay phase.

preprint2005arXiv

Giant flare of SGR 1806-20 from a relativistic jet

Japanese magnetospheric explorer GEOTAIL recorded a detailed light curve during the initial 600 msec of a giant flare from SGR 1806-20 on December 27, 2004. We show that the observed light curve is well explained by an emission from relativistically expanding fireballs, like those of gamma-ray bursts (GRBs). Especially, the observed rapid fading after 500msec suggests that ejecta is collimated in a jet. We derive an upper limit on the jet opening half-angle of 0.2 radian that is as narrow as those of GRBs.

preprint2004arXiv

A Unified Model of Short and Long Gamma-Ray Bursts, X-Ray Rich Gamma-Ray Bursts, and X-Ray Flashes

We propose a possible unified model of short and long gamma-ray bursts (GRBs), X-ray rich GRBs, and X-ray flashes. The jet of a GRB is assumed to consist of multiple sub-jets or sub-shells (i.e., an inhomogeneous jet model). The multiplicity of the sub-jets along a line of sight n_s is an important parameter. If n_s is large (>> 1) the event looks like a long GRB, while if n_s=1, the event looks like a short GRB. Finally, when n_s=0, the event looks like an X-ray flash or an X-ray rich GRB. Furthermore, our model may also explain the bimodal distributions of T_{90} duration of BATSE-GRBs.

preprint2004arXiv

Gravitational Wave Memory of Gamma-Ray Burst Jets

Gamma-Ray Bursts (GRBs) are now considered as relativistic jets. We analyze the gravitational waves from the acceleration stage of the GRB jets. We show that (i) the point mass approximation is not appropriate if the opening half-angle of the jet is larger than the inverse of the Lorentz factor of the jet, (ii) the gravitational waveform has many step function like jumps, and (iii) the practical DECIGO and BBO may detect such an event if the GRBs occur in Local group of galaxy. We found that the light curve of GRBs and the gravitational waveform are anti-correlated so that the detection of the gravitational wave is indispensable to determine the structure of GRB jets.

preprint2003arXiv

Circular Polarization from Gamma-ray Burst Afterglows

We investigate the circular polarization (CP) from Gamma-Ray Burst (GRB) afterglows. We show that a tangled magnetic field cannot generate CP without an ordered magnetic field because there is always an oppositely directed field, so that no handedness exists. This implies the observation of CP could be a useful probe of an ordered field, which carries valuable information on the GRB central engine. By solving the transfer equation of polarized radiation, we find that the CP reaches 1% at radio frequencies and 0.01% at optical for the forward shock, and 10-1% at radio and 0.1-0.01% at optical for the reverse shock.

preprint2003arXiv

Delayed Flashes from Counter Jets of Gamma Ray Bursts

If X-ray flashes are due to the forward jet emissions from gamma ray bursts (GRBs) observed with large viewing angles, we show that a prompt emission from a counter jet should be observed as a delayed flash in the UV or optical band several hours to a day after the X-ray flash. Ultraviolet and Optical Telescope on Swift can observe the delayed flashes within ~13 Mpc, so that (double-sided) jets of GRBs may be directly confirmed. Since the event rate of delayed flashes detected by Swift may be as small as 6*10^{-5}events/yr, we require more sensitive detectors in future experiments.