Source author record

Kenji Toma

Kenji Toma 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

34works
12topics
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

34 published item(s)

preprint2022arXiv

High-energy Gamma-rays from Magnetically Arrested Disks in Nearby Radio Galaxies

The origins of the GeV gamma-rays from nearby radio galaxies are unknown. Hadronic emission from magnetically arrested disks (MADs) around central black holes (BHs) is proposed as a possible scenario. Particles are accelerated in the MAD by magnetic reconnection and stochastic turbulence acceleration. We pick out the fifteen brightest radio galaxies in the GeV band from the Fermi 4LAC-DR2 catalog and apply the MAD model. We find that we can explain the data in the GeV bands by the MAD model if the accretion rate is lower than 0.1% of the Eddington rate. For a higher accretion rate, GeV gamma-rays are absorbed by two-photon interaction due to copious low-energy photons. If we assume another proposed prescription of the electron heating rate by magnetic reconnection, the MAD model fails to reproduce the GeV data for the majority of our sample. This indicates that the electron heating rate is crucial. We also apply the MAD model to Sgr A* and find that GeV gamma-rays observed at the Galactic center do not come from the MAD of Sgr A*. We estimate the cosmic ray intensity from Sgr A*, but it is too low to explain the high-energy cosmic ray intensity on Earth.

preprint2022arXiv

Interaction of a Relativistic Magnetized Collisionless Shock with a Dense Clump

The interactions between a relativistic magnetized collisionless shock and dense clumps have been expected to play a crucial role on the magnetic field amplification and cosmic-ray acceleration. We investigate this process by two-dimensional Particle-In-Cell (PIC) simulations for the first time, where the clump size is much larger than the gyroradius of downstream particles. We also perform relativistic magnetohydrodynamic (MHD) simulations for the same condition to see the kinetic effects. We find that particles escape from the shocked clump along magnetic field lines in the PIC simulations, so that the vorticity is lower than that in the MHD simulations. Moreover, in both the PIC and MHD simulations, the shocked clump quickly decelerates because of relativistic effects. Owing to the escape and the deceleration, the shocked clump cannot amplify the downstream magnetic field in relativistic collisionless shocks. This large-scale PIC simulation opens a new window to understand large-scale behaviors in collisionless plasma systems.

preprint2022arXiv

The Photon Ring in M87*

We report measurements of the gravitationally lensed secondary image -- the first in an infinite series of so-called "photon rings" -- around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical expectations, providing clear evidence that such measurements probe spacetime and a striking confirmation of the models underlying the first set of EHT results. The residual diffuse emission evolves on timescales comparable to one week. We are able to detect with high significance a southwestern extension consistent with that expected from the base of a jet that is rapidly rotating in the clockwise direction. This result adds further support to the identification of the jet in M87* with a black hole spin-driven outflow, launched via the Blandford-Znajek process. We present three revised estimates for the mass of M87* based on identifying the modeled thin ring component with the bright ringlike features seen in simulated images, one of which is only weakly sensitive to the astrophysics of the emission region. All three estimates agree with each other and previously reported values. Our strongest mass constraint combines information from both the ring and the diffuse emission region, which together imply a mass-to-distance ratio of $4.20^{+0.12}_{-0.06}~μ{\rm as}$ and a corresponding black hole mass of $(7.13\pm0.39)\times10^9M_\odot$, where the error on the latter is now dominated by the systematic uncertainty arising from the uncertain distance to M87*.

preprint2021arXiv

Hadronic High-energy Emission from Magnetically Arrested Disks in Radio Galaxies

We propose a novel interpretation that gamma-rays from nearby radio galaxies are hadronic emission from magnetically arrested disks (MADs) around central black holes (BHs). The magnetic energy in MADs is higher than the thermal energy of the accreting plasma, where the magnetic reconnection or turbulence may efficiently accelerate non-thermal protons. They emit gamma-rays via hadronic processes, which can account for the observed gamma-rays for M87 and NGC 315. Non-thermal electrons are also accelerated with protons and produce MeV gamma-rays, which is useful to test our model by proposed MeV satellites. The hadronic emission from the MADs may significantly contribute to the GeV gamma-ray background and produce the multi-PeV neutrino background detectable by IceCube-Gen2. In addition, gamma-rays from MADs provide electron-positron pairs through two-photon pair production at the BH magnetosphere. These pairs can screen the vacuum gap, which affects high-energy emission and jet-launching mechanisms in radio galaxies.

preprint2021arXiv

The response of black hole spark gaps to external changes: A production mechanism of rapid TeV flares?

We study the response of a starved Kerr black hole magnetosphere to abrupt changes in the intensity of disk emission and in the global magnetospheric current, by means of 1D general relativistic particle-in-cell simulations. Such changes likely arise from the intermittency of the accretion process. We find that in cases where the pair production opacity contributed by the soft disk photons is modest, as in, e.g., M87, such changes can give rise to delayed, strong TeV flares, dominated by curvature emission of particles accelerated in the gap. The flare rise time, and the delay between the external variation and the onset of the flare emitted from the outer gap boundary, are of the order of the light crossing time of the gap. The rapid, large amplitude TeV flares observed in M87 and, perhaps, other AGNs may be produced by such a mechanism.

preprint2020arXiv

A Thermodynamic Approach towards the Question "What is Cellular Life?"

The question "What is life?" has been asked and studied by the researchers of various fields. Nevertheless, no global theory which unified various aspects of life has been proposed so far. Considering that the physical principle for the theory of birth should be the one known for the unanimated world, and that the life processes are irreversibly selective, we showed by a deductive inference that the maximum entropy production principle plays an essential role for the birth and the evolution of life in a fertile environment. In order to explain the survival strategy of life in a barren period of environment, we also proposed that life had simultaneously developed a reversible on and off switching mechanism of the chemical reactions by the dynamics of equilibrium thermodynamics. Thus, the birth and evolution of life have been achieved by the cooperation between the driving force due to the non-equilibrium thermodynamics and the protective force due to the equilibrium thermodynamics in the alternating environmental conditions.

preprint2020arXiv

Efficient acceleration of cylindrical jets: effects of radiative cooling and tangled magnetic field

Diverging supersonic flows are accelerating, as in the case of a de Laval nozzle, and the same concept has been applied for acceleration of magnetohydrodynamic flows in the universe. Here, we study the dynamics of "non-diverging" cylindrical supersonic flows and show that they can be accelerated by effects of radiative cooling and the tangled magnetic field. In addition to radiative cooling of the jet materials (cooling effect), conversion of the ordered magnetic field into the turbulent one (conversion effect) and dissipation of the turbulent magnetic field (dissipation effect) are formulated according to our study on pulsar wind nebulae. Although each of the cooling and conversion effects is an ineffective acceleration process, the terminal velocity of magnetized cylindrical jets attains about half of the maximum possible value when the cooling, conversion and dissipation effects work simultaneously. The radiation efficiency is also about half of the total luminosity of the jet in the case of maximal acceleration. The concept for flow acceleration by the non-ideal MHD effects may be useful for studying relativistic jets in active galactic nuclei, in which the region near the jet axis is expected to be cylindrical and kink unstable.

preprint2016arXiv

Causal production of the electromagnetic energy flux and role of the negative energies in Blandford-Znajek process

Blandford-Znajek process, the steady electromagnetic energy extraction from a rotating black hole (BH), is widely believed to work for driving relativistic jets in active galactic nuclei, gamma-ray bursts and Galactic microquasars, although it is still under debate how the Poynting flux is causally produced and how the rotational energy of the BH is reduced. We generically discuss the Kerr BH magnetosphere filled with a collisionless plasma screening the electric field along the magnetic field, extending the arguments of Komissarov and our previous paper, and propose a new picture for resolving the issues. For the magnetic field lines threading the equatorial plane in the ergosphere, we find that the inflow of particles with negative energy as measured in the coordinate basis is generated near that plane as a feedback from the Poynting flux production, which appears to be a similar process to the mechanical Penrose process. For the field lines threading the event horizon, we first show that the concept of the steady inflow of negative electromagnetic energy is not physically essential, partly because the sign of the electromagnetic energy density depends on the coordinates. Then we build an analytical toy model of a time-dependent process both in the Boyer-Lindquist and Kerr-Schild coordinate systems in which the force-free plasma injected continuously is filling a vacuum, and suggest that the structure of the steady outward Poynting flux is causally constructed by the displacement current and the cross-field current at the in-going boundary between the plasma and the vacuum. In the steady state, the Poynting flux is maintained without any electromagnetic source.

preprint2016arXiv

Gamma-ray bursts and Population III stars

Gamma-ray bursts (GRBs) are ideal probes of the epoch of the first stars and galaxies. We review the recent theoretical understanding of the formation and evolution of the first (so-called Population III) stars, in light of their viability of providing GRB progenitors. We proceed to discuss possible unique observational signatures of such bursts, based on the current formation scenario of long GRBs. These include signatures related to the prompt emission mechanism, as well as to the afterglow radiation, where the surrounding intergalactic medium might imprint a telltale absorption spectrum. We emphasize important remaining uncertainties in our emerging theoretical framework.

preprint2016arXiv

Stochastic Particle Acceleration in Turbulence Generated by the Magnetorotational Instability

We investigate stochastic particle acceleration in accretion flows. It is believed that the magnetorotational instability (MRI) generates turbulence inside accretion flows and that cosmic rays (CRs) are accelerated by the turbulence. We calculate equations of motion for CRs in the turbulent fields generated by MRI with the shearing box approximation without back reaction to the field. The results show that the CRs randomly gain or lose their energies through the interaction with the turbulent fields. The CRs diffuse in the configuration space anisotropically: The diffusion coefficient in direction of the unperturbed flow is about twenty times higher than the Bohm coefficient, while those in the other directions are only a few times higher than the Bohm. The momentum distribution is isotropic, and its evolution can be described by the diffusion equation in momentum space where the diffusion coefficient is a power-law function of the CR momentum. We show that the shear acceleration efficiently works for energetic particles. We also cautiously note that in the shearing box approximation, particles that cross the simulation box many times along the radial direction suffer unphysical runaway acceleration by the Lorentz transformation, which needs to be taken with special care.

preprint2015arXiv

Neutrino and Cosmic-Ray Emission and Cumulative Background from Radiatively Inefficient Accretion Flows in Low-Luminosity Active Galactic Nuclei

We study high-energy neutrino and cosmic-ray (CR) emission from the cores of low-luminosity active galactic nuclei (LLAGN). In LLAGN, the thermalization of particles is expected to be incomplete in radiatively inefficient accretion flows (RIAFs), allowing the existence of non-thermal particles. In this work, assuming stochastic particle acceleration due to turbulence in RIAFs, we solve the Fokker-Planck equation and calculate spectra of escaping neutrinos and CRs. The RIAF in LLAGN can emit CR protons with $\gtrsim10$ PeV energies and TeV-PeV neutrinos generated via $pp$ and/or $pγ$ reactions. We find that, if $\sim1$\% of the accretion luminosity is carried away by non-thermal ions, the diffuse neutrino intensity from the cores of LLAGN may be as high as $E_ν^2Φ_ν\sim3\times{10}^{-8} {\rm GeV} {\rm cm}^{-2} {\rm s}^{-1}{\rm sr}^{-1}$, which can be compatible with the observed IceCube data. This result does not contradict either of the diffuse gamma-ray background observed by {\it Fermi} or observed diffuse cosmic-ray flux. Our model suggests that, although very-high-energy gamma rays may not escape, radio-quiet AGN with RIAFs can emit GeV gamma-rays, which could be used for testing the model. We also calculate the neutron luminosity from RIAFs of LLAGN, and discuss a strong constraint on the model of jet mass loading mediated by neutrons from the diffuse neutrino observation.

preprint2015arXiv

Synchrotron Self-Compton Emission by Relativistic Electrons under Stochastic Acceleration: Application to Mrk 421 and Mrk 501

We examine the applicability of the stochastic electron acceleration to two high synchrotron peaked blazars, Mrk 421 and Mrk 501, assuming synchrotron self-Compton emission of gamma-rays. Our model considers an emitting region moving at relativistic speed, where non-thermal electrons are accelerated and attain a steady-state energy spectrum together with the photons they emit. The kinetic equations of the electrons and photons are solved numerically, given a stationary wave number spectrum of the magnetohydrodynamic (MHD) disturbances, which are responsible for the electron acceleration and escape. Our simple formulation appears to reproduce the two well-sampled, long-term averaged photon spectra. In order to fit the model to the emission component from the radio to the X-ray bands, we need both a steeper wave spectral index than the Kolmogorov spectrum and efficient particle escape. Although the model provides a natural explanation for the high-energy cutoff of the electron energy distribution, the derived physical parameters raise a problem with an energy budget if the MHD waves with the Alfv{é}n velocity are assumed to be the acceleration agent.

preprint2015arXiv

X-Ray Polarimetry with the Polarization Spectroscopic Telescope Array (PolSTAR)

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

preprint2014arXiv

Electromotive Force in the Blandford-Znajek Process

One of the mechanisms widely considered for driving relativistic jets in active galactic nuclei, galactic microquasars, and gamma-ray bursts is the electromagnetic extraction of the rotational energy of a central black hole, i.e., the Blandford-Znajek process, although the origin of the electromotive force in this process is still under debate. We study this process as the steady unipolar induction in the Kerr black hole magnetosphere filled with a collisionless plasma screening the electric field (the D field) along the magnetic field (the B field), i.e., D dot B = 0. We extend the formulations and arguments made by Komissarov, and generally show that the origin of the electromotive force is ascribed to the ergosphere. It is explicitly shown that open magnetic field lines penetrating the ergosphere have a region where the D field is stronger than the B field in the ergosphere, and it keeps driving the poloidal currents and generating the electromotive force and the outward Poynting flux. The range of the possible value of the so-called angular velocity of the magnetic field line Omega_F is deduced for the field lines threading the equatorial plane in the ergosphere. We briefly discuss the relation between our conclusion and the ideal magnetohydrodynamic condition.

preprint2013arXiv

Efficient Acceleration of Relativistic Magnetohydrodynamic Jets

Relativistic jets in active galactic nuclei, galactic microquasars, and gamma-ray bursts are widely considered to be magnetohydrodynamically driven by black hole accretion systems, although conversion mechanism from Poynting into particle kinetic energy flux is still open. Recent detailed numerical and analytical studies of global structures of steady, axisymmetric magnetohydrodynamic (MHD) flows with specific boundary conditions have not reproduced as rapid an energy conversion as required by observations. In order to find more suitable boundary conditions, we focus on the flow along a poloidal magnetic field line just inside the external boundary, without treating transfield force balance in detail. We find some examples of the poloidal field structure and corresponding external pressure profile for an efficient and rapid energy conversion as required by observations, and that the rapid acceleration requires a rapid decrease of the external pressure above the accretion disk. We also clarify the differences between the fast magnetosonic point of the MHD flow and the sonic point of de Laval nozzle.

preprint2013arXiv

Time-Dependent Models for Blazar Emissions with the Second-Order Fermi Acceleration

The second-order Fermi acceleration (Fermi-II) driven by turbulence may be responsible for the electron acceleration in blazar jets. We test this model with time-dependent simulations. The hard electron spectrum predicted by the Fermi-II process agrees with the hard photon spectrum of 1ES 1101-232. For other blazars that show softer spectra, the Fermi-II model requires radial evolution of the electron injection rate and/or diffusion coefficient in the outflow. Such evolutions can yield a curved electron spectrum, which can reproduce the synchrotron spectrum of Mrk 421 from the radio to the X-ray regime.The photon spectrum in the GeV energy range of Mrk 421 is hard to fit with a synchrotron self-Compton model. However, if we introduce an external radio photon field with a luminosity of $4.9 \times 10^{38}~\mbox{erg}~\mbox{s}^{-1}$, GeV photons are successfully produced via inverse Compton scattering. The temporal variability of the diffusion coefficient or injection rate causes flare emission. The observed synchronicity of X-ray and TeV flares implies a decrease of the magnetic field in the flaring source region.

preprint2012arXiv

Baryon Loading of AGN Jets Mediated by Neutrons

Plasmas of geometrically thick, black hole (BH) accretion flows in active galactic nuclei (AGNs) are generally collisionless for protons, and involve magnetic field turbulence. Under such conditions a fraction of protons can be accelerated stochastically and create relativistic neutrons via nuclear collisions. These neutrons can freely escape from the accretion flow and decay into protons in dilute polar region above the rotating BH to form relativistic jets. We calculate geometric efficiencies of the neutron energy and mass injections into the polar region, and show that this process can deposit luminosity as high as L_j ~ 2e-3 dot{M} c^2 and mass loading dot{M}_j ~ 6e-4 dot{M} for the case of the BH mass M ~ 1e8 M_sun, where dot{M} is mass accretion rate. The terminal Lorentz factors of the jets are Gamma ~ 3, and they may explain the AGN jets having low luminosities. For higher luminosity jets, which can be produced by additional energy inputs such as Poynting flux, the neutron decay still can be a dominant mass loading process, leading to e.g., Gamma ~ 50 for L_{j,tot} ~ 3e-2 dot{M}c^2.

preprint2012arXiv

Magnetic Structures in Gamma-Ray Burst Jets Probed by Gamma-Ray Polarization

We report polarization measurements in two prompt emissions of gamma-ray bursts, GRB 110301A and GRB 110721A, observed with the Gamma-ray burst polarimeter (GAP) aboard IKAROS solar sail mission. We detected linear polarization signals from each burst with polarization degree of $Π= 70 \pm 22$% with statistical significance of $3.7 σ$ for GRB 110301A, and $Π= 84^{+16}_{-28}$% with $3.3 σ$ confidence level for GRB 110721A. We did not detect any significant change of polarization angle. These two events had shorter durations and dimmer brightness compared with GRB 100826A, which showed a significant change of polarization angle, as reported in Yonetoku et al. (2011). Synchrotron emission model can be consistent with all the data of the three GRBs, while photospheric quasi-thermal emission model is not favorable. We suggest that magnetic field structures in the emission region are globally-ordered fields advected from the central engine.

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

Strict Limit on CPT Violation from Polarization of Gamma-Ray Bursts

We report the strictest observational verification of CPT invariance in the photon sector, as a result of gamma-ray polarization measurement of distant gamma-ray bursts (GRBs), which are brightest stellar-size explosions in the universe. We detected the gamma-ray polarization of three GRBs with high significance, and the source distances may be constrained by a well-known luminosity indicator for GRBs. For the Lorentz- and CPT-violating dispersion relation E_{\pm}^2=p^2 \pm 2ξp^3/M_{Pl}, where \pm denotes different circular polarization states of the photon, the parameter ξis constrained as |ξ|<O(10^{-15}). Barring precise cancellation between quantum gravity effects and dark energy effects, the stringent limit on the CPT-violating effect leads to the expectation that quantum gravity presumably respects the CPT invariance.

preprint2011arXiv

A Photosphere-Internal Shock Model of Gamma-Ray Bursts: Case Studies of Fermi/LAT Bursts

Radially inhomogeneous gamma-ray burst (GRB) jets release variable photospheric emission and can have internal shocks occurring above the photosphere. We generically formulate a photospheric emission model of GRBs including Compton up-scattered photospheric (UP) emission off the electrons (and positrons) in the internal shocks, and find that the photospheric emission may correspond to the traditional (Band) component at <~1 MeV and the UP emission to the high-energy emission observed by Fermi/LAT for some GRBs at >~ 10 MeV. The two components can be separate in the spectrum in some cases or can mimic a smooth broad Band spectrum in other cases. We apply our formulation to the well-studied long and short LAT GRBs, GRB 080916C, GRB 090902B, and GRB 090510, and typically find reasonable parameters for fitting the time-binned spectra, although fine tuning of several parameters is required. The observed delays of the high-energy emission with respect to the MeV emission which are large compared to the variability times are unlikely to be due to simple kinematic effects of a non-evolving jet. These delays may be attributed to the temporal evolution of the physical parameters of the jet, and thus the delay timescales could provide a potential tool for investigating the structures of GRB jets themselves and their progenitors. The difference of the delay timescales of long and short GRBs inferred from the Fermi data might be due to the differences in the progenitors of long and short GRBs. Some other properties and consequences of this model are discussed, including temporal correlations among the prompt optical, the soft X-ray, and the distinct high-energy component as well as the Band component.

preprint2011arXiv

Detection of Gamma-Ray Polarization in Prompt Emission of GRB 100826A

We report the polarization measurement in prompt $γ$-ray emission of GRB 100826A with the Gamma-Ray Burst Polarimeter (GAP) aboard the small solar power sail demonstrator IKAROS. We detected the firm change of polarization angle (PA) during the prompt emission with 99.9% ($3.5 σ$) confidence level, and the average polarization degree ($Π$) of $27 \pm 11$% with 99.4% ($2.9 σ$) confidence level. Here the quoted errors are given at 1 $σ$ confidence level for two parameters of interest. The systematic errors have been carefully included in this analysis, unlike any previous reports. Such a high $Π$ can be obtained in several emission models of gamma-ray bursts (GRBs), including synchrotron and photospheric models. However, it is difficult to explain the observed significant change of PA within the framework of axisymmetric jet as considered in many theoretical works. The non-axisymmetric (e.g., patchy) structures of the magnetic fields and/or brightness inside the relativistic jet are therefore required within the observable angular scale of $\sim Γ^{-1}$. Our observation strongly indicates that the polarization measurement is a powerful tool to constrain the GRB production mechanism, and more theoretical works are needed to discuss the data in more details.

preprint2011arXiv

High Energy Neutrino Emission from the Earliest Gamma-Ray Bursts

We discuss the high energy neutrino emission from gamma-ray bursts resulting from the earliest generation (`population III') stars forming in the Universe, whose core collapses into a black hole. These gamma-ray bursts are expected to produce a highly relativistic, magnetically dominated jet, where protons can be accelerated to ultra-high energies. These interact with the photons produced by the jet, leading to ultra-high energy photo-meson neutrinos as well as secondary leptons and photons. The photon luminosity and the shock properties, and thus the neutrino spectrum, depend on the mass of the black holes as well as on the density of the surrounding external gas. We calculate the individual source neutrino spectral fluxes and the expected diffuse neutrino flux for various source parameters and evolution scenarios. Both the individual and diffuse signals appear detectable in the 1-300 PeV range with current and planned neutrino detectors such as IceCube and ARIANNA, provided the black hole mass is in excess of 30-100 solar masses. This provides a possible test for the debated mass of the progenitor stellar objects, as well as a probe for the early cosmological environment and the formation rate of the earliest structures.

preprint2011arXiv

On the High Energy Emission of the Short GRB 090510

Long-lived high-energy (>100MeV) emission, a common feature of most Fermi-LAT detected gamma-ray burst, is detected up to \sim 10^2 s in the short GRB 090510. We study the origin of this long-lived high-energy emission, using broad-band observations including X-ray and optical data. We confirm that the late > 100 MeV, X-ray and optical emission can be naturally explained via synchrotron emission from an adiabatic forward shock propagating into a homogeneous ambient medium with low number density. The Klein-Nishina effects are found to be significant, and effects due to jet spreading and magnetic field amplification in the shock appear to be required. Under the constraints from the low-energy observations, the adiabatic forward shock synchrotron emission is consistent with the later-time (t>2s) high-energy emission, but falls below the early-time (t < 2s) high energy emission. Thus we argue that an extra high energy component is needed at early times. A standard reverse shock origin is found to be inconsistent with this extra component. Therefore, we attribute the early part of the high-energy emission (t< 2s) to the prompt component, and the long-lived high energy emission (t>2s) to the adiabatic forward shock synchrotron afterglow radiation. This avoids the requirement for an extremely high initial Lorentz factor.

preprint2011arXiv

On the Implications of Late Internal Dissipation for Shallow-Decay Afterglow Emission and Associated High-Energy Gamma-Ray Signals

The origin of the shallow-decay emission during early X-ray afterglows has been an open issue since the launch of the Swift satellite. One of the appealing models is the late internal dissipation model, where X-ray emission during the shallow-decay phase is attributed to internal dissipation, analogous to the prompt gamma-ray emission. We discuss possible scenarios of the late prompt emission, such as late internal shocks, magnetic reconnection, and photospheric dissipation. We also consider the consequences of late dissipation and a two-component (early and late) jet model for the high-energy (GeV-TeV) emission. We study not only synchrotron self-Compton (SSC) emission from the early and late jets but also external inverse-Compton (EIC) emission, which is naturally predicted in the late dissipation model. For the latter, we perform numerical calculations taking into account the equal-arrival-time surface of EIC photons and show that the EIC component typically has a peak at ~1-100 GeV which may dominate over the SSC components. We demonstrate that very high energy gamma rays from both these components are detectable for nearby and/or energetic gamma-ray bursts, with current and future Cherenkov detectors such as MAGIC, VERITAS, CTA and HAWC, and possibly Fermi. Although the expected event rate would not be large, detections should be useful as a test of the model. Multi-wavelength observations using both the ground-based telescopes and the Swift and/or Fermi satellites are also important to constrain the models.

preprint2011arXiv

Population III GRB Afterglows: Constraints on Stellar Masses and External Medium Densities

Population III stars are theoretically expected to be prominent around redshifts z ~ 20, consisting of mainly very massive stars with M_* >~ 10 M_sun$, but there is no direct observational evidence for these objects. They may produce collapsar gamma-ray bursts (GRBs), with jets driven by magnetohydrodynamic processes, whose total isotropic-equivalent energy could be as high as E_iso >~ 10^{57} erg over a cosmological-rest-frame duration of t_d >~ 10^4 s, depending on the progenitor mass. Here we calculate the afterglow spectra of such Pop. III GRBs based on the standard external shock model, and show that they will be detectable with the Swift BAT/XRT and Fermi LAT instruments. We find that in some cases a spectral break due to electron-positron pair creation will be observable in the LAT energy range, which can put constraints on the ambient density of the pre-collapse Pop. III star. Thus, high redshift GRB afterglow observations could be unique and powerful probes of the properties of Pop. III stars and their environments. We examine the trigger threshold of the BAT instrument in detail, focusing on the image trigger system, and show that the prompt emission of Pop. III GRBs could also be detected by BAT. Finally we briefly show that the late-time radio afterglows of Pop. III GRBs for typical parameters, despite the large distances, can be very bright: ~ 140 mJy at 1 GHz, which may lead to a constraint on the Pop. III GRB rate from the current radio survey data, and ~2.4 mJy at 70 MHz, which implies that Pop. III GRB radio afterglows could be interesting background source candidates for 21 cm absorption line detections.

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

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.

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

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.

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.