Source author record

Ryo Yamazaki

Ryo Yamazaki 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

66works
13topics
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

66 published item(s)

preprint2024arXiv

GRB 080710: A narrow, structured jet showing a late, achromatic peak in the optical and infrared afterglow?

We present a possible theoretical interpretation of the observed afterglow emission of long gamma-ray burst GRB 080710. While its prompt GRB emission properties are normal, the afterglow light curves in the optical and infrared bands are exceptional in two respects. One is that the observed light curves of different wavelengths have maximum at the same time, and that the achromatic peak time, $2.2\times10^3$ s after the burst trigger, is about an order of magnitude later than typical events. The other is that the observed flux before the peak increases more slowly than theoretically expected so far. Assuming that the angular distribution of the outflow energy is top-hat or Gaussian-shaped, we calculate the observed light curves of the synchrotron emission from the relativistic jets and explore the model parameters that explain the observed data. It is found that a narrowly collimated Gaussian-shaped jet with large isotropic-equivalent energy is the most plausible model for reproducing the observed afterglow behavior. Namely, an off-axis afterglow scenario to the achromatic peak is unlikely. The inferred values of the opening angle and the isotropic equivalent energy of the jet are possibly similar to those of GRB 221009A, but the jet of GRB 080710 has a much smaller efficiency of the prompt gamma-ray emission. Our results indicate a greater diversity of the GRB jet properties than previously thought.

preprint2022arXiv

Gamma-ray Diagnostics of r-process Nucleosynthesis in the Remnants of Galactic Binary Neutron-Star Mergers

We perform a full nuclear-network numerical calculation of the $r$-process nuclei in binary neutron-star mergers (NSMs), with the aim of estimating $γ$-ray emissions from the remnants of Galactic NSMs up to $10^6$ years old. The nucleosynthesis calculation of 4,070 nuclei is adopted to provide the elemental composition ratios of nuclei with an electron fraction $Y_{\rm e}$ between 0.10 and 0.45 . The decay processes of 3,237 unstable nuclei are simulated to extract the $γ$-ray spectra. As a result, the NSMs have different spectral color in $γ$-ray band from various other astronomical objects at less than $10^5$ years old. In addition, we propose a new line-diagnostic method for $Y_{\rm e}$ that uses the line ratios of either $^{137{\rm m}}$Ba/$^{85}$K or $^{243}$Am/$^{60{\rm m}}$Co, which become larger than unity for young and old $r$-process sites, respectively, with a low $Y_{\rm e}$ environment. From an estimation of the distance limit for $γ$-ray observations as a function of the age, the high sensitivity in the sub-MeV band, at approximately $10^{-9}$ photons s$^{-1}$ cm$^{-2}$ or $10^{-15}$ erg s$^{-1}$ cm$^{-2}$, is required to cover all the NSM remnants in our Galaxy if we assume that the population of NSMs by \citet{2019ApJ...880...23W}. A $γ$-ray survey with sensitivities of $10^{-8}$--$10^{-7}$ photons s$^{-1}$ cm$^{-2}$ or $10^{-14}$--$10^{-13}$ erg s$^{-1}$ cm$^{-2}$ in the 70--4000 keV band is expected to find emissions from at least one NSM remnant under the assumption of NSM rate of 30 Myr$^{-1}$. The feasibility of $γ$-ray missions to observe Galactic NSMs are also studied.

preprint2022arXiv

High-power laser experiment forming a supercritical collisionless shock in a magnetized uniform plasma at rest

We present a new experimental method to generate quasi-perpendicular supercritical magnetized collisionless shocks. In our experiment, ambient nitrogen (N) plasma is at rest and well-magnetized, and it has uniform mass density. The plasma is pushed by laser-driven ablation aluminum (Al) plasma. Streaked optical pyrometry and spatially resolved laser collective Thomson scattering clarify structures of plasma density and temperatures, which are compared with one-dimensional particle-in-cell simulations. It is indicated that just after the laser irradiation, the Al plasma is magnetized by a self-generated Biermann battery field, and the plasma slaps the incident N plasma. The compressed external field in the N plasma reflects N ions, leading to counter-streaming magnetized N flows. Namely we identify the edge of the reflected N ions. Such interacting plasmas form a magnetized collisionless shock.

preprint2022arXiv

Synchrotron Self-Compton Emission in the Two-Component Jet Model for Gamma-Ray Bursts

Gamma-ray bursts (GRBs) are intense bursts of high-energy photons (prompt emissions) caused by relativistic jets. After the emissions, multi-wavelength afterglows, from radio to very-high-energy (VHE) gamma-ray, last for more than a few days. In the past three years, the VHE gamma-ray photons from four GRBs (GRBs 180720B, 190114C, 190829A and 201216C) were detected by ground-based Imaging Atmospheric Cherenkov Telescopes, such as the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) telescopes and the High Energy Stereoscopic System (H.E.S.S.). One of them, GRB 190829A, had some peculiar features of showing achromatic peaks in X-ray and optical bands at $1.4\times10^3$ s and being classified as low-luminosity GRBs. Previously, we proposed a two-component jet model, which has `narrow jet' with a small initial opening half-angle $θ_0=0.015$ rad and large bulk Lorentz factor $Γ_0=350$, and `wide jet' with $θ_0=0.1$ rad and $Γ_0=20$. The narrow jet explained the early X-ray and optical emissions and apparently small isotropic gamma-ray energy and peak energy in the off-axis viewing case. Furthermore, the late X-ray and radio (1.3 and 15.5 GHz) afterglows were emitted from the wide jet. Here, we calculate the VHE gamma-ray flux by the synchrotron self-Compton (SSC) emission. The multi-wavelength afterglows of GRB 190829A including the VHE gamma-ray emission are well explained by our two-component jet model. The afterglow emissions from our two-component jet are also consistent with the observational results of GRBs 180720B, 190114C and 201216C, when the jets are viewed on-axis. Furthermore, we discuss the detectability of off-axis orphan afterglows by the Cherenkov Telescope Array (CTA).

preprint2022arXiv

The Gamma-Ray Emission from the Supernova Remnant RX J1713.7-3946 Interacting with Two-phase Medium

We study the origin of gamma rays from the supernova remnant (SNR) RX J1713.7-3946. Using an analytical model, we calculate the distribution of cosmic rays (CRs) around the SNRs. Motivated by the results of previous studies, we assume that the SNR is interacting with two-phase interstellar medium (ISM), where dense clumps are surrounded by tenuous interclump medium. We also assume that only higher-energy protons (~> TeV) can penetrate the dense clumps. We find that pi^0-decay gamma rays produced by protons reproduce the observed gamma-ray spectrum peaked at ~TeV. On the other hand, it has recently been indicated that the observed ISM column density (N_p), the X-ray surface brightness (I_X), and the gamma-ray surface brightness (I_g) at grid points across the SNR form a plane in the three-dimensional (3D) space of (N_p, I_X, I_g). We find that the planar configuration is naturally reproduced if the ISM or the CR electron-to-proton ratio is not spherically uniform. We show that the shift of the observed data in the 3D space could be used to identify which of the quantities, the ISM density, the CR electron-to-proton ratio, or the magnetic field, varies in the azimuthal direction of the SNR.

preprint2022arXiv

X-ray Line Diagnostics of Ion Temperature at Cosmic-Ray Accelerating Collisionless Shocks

A novel collisionless shock jump condition is suggested by modeling the entropy production at the shock transition region. We also calculate downstream developments of the atomic ionization balance and the ion temperature relaxation in supernova remnants (SNRs). The injection process and subsequent acceleration of cosmic-rays (CRs) in the SNR shocks are closely related to the formation process of the collisionless shocks. The formation of the shock is caused by wave-particle interactions. Since the wave-particle interactions result in energy exchanges between electromagnetic fields and charged particles, the randomization of particles associated with the shock transition may occur with the rate given by the scalar product of the electric field and current. We find that order-of-magnitude estimates of the randomization with reasonable strength of the electromagnetic fields in the SNR constrain the amount of the CR nuclei and ion temperatures. The constrained amount of the CR nuclei can be sufficient to explain the Galactic CRs. The ion temperature becomes significantly lower than in the case of no CRs. To distinguish the case without CRs, we perform synthetic observations of atomic line emissions from the downstream region of the SNR RCW~86. Future observations by {\it XRISM} and {\it Athena} can distinguish whether the SNR shock accelerates the CRs or not from the ion temperatures.

preprint2021arXiv

Observational constraints on the maximum energies of accelerated particles in supernova remnants: low maximum energies and a large variety

Supernova remnants (SNRs) are thought to be the most promising sources of Galactic cosmic rays. One of the principal questions is whether they are accelerating particles up to the maximum energy of Galactic cosmic rays ($\sim$ PeV). In this work, a systematic study of gamma-ray emitting SNRs is conducted as an advanced study of Suzuki et al. 2021. Our purpose is to newly measure the evolution of maximum particle energies with increased statistics and better age estimates. We model their gamma-ray spectra to constrain the particle-acceleration parameters. Two candidates of the maximum energy of freshly accelerated particles, the gamma-ray cutoff and break energies, are found to be well below PeV. We also test a spectral model that includes both the freshly accelerated and escaping particles to estimate the maximum energies more reliably, but no tighter constraints are obtained with current statistics. The average time dependences of the cutoff energy ($\propto t^{-0.81 \pm 0.24}$) and break energy ($\propto t^{-0.77 \pm 0.23}$) cannot be explained with the simplest acceleration condition (Bohm limit), and requires shock-ISM (interstellar medium) interaction. The average maximum energy during lifetime is found to be $\lesssim 20$ TeV $(t_{\rm M}/1~{\rm kyr})^{-0.8}$ with $t_{\rm M}$ being the age at the maximum, which reaches PeV if $t_{\rm M} \lesssim 10$ yr. The maximum energies during lifetime are suggested to have a variety of 1.1--1.8 dex from object to object. Although we cannot isolate the cause of this variety, this work provides an important clue to understand the microphysics of particle acceleration in SNRs.

preprint2020arXiv

Fast Particle Acceleration at Perpendicular Shocks with Uniform Upstream Magnetic Field and Strong Downstream Turbulence

The shock waves produce relativistic particles via the diffusive shock acceleration (DSA) mechanism. Among various circumstances, fast acceleration has been expected for perpendicular shocks. We investigate the acceleration time and the energy spectrum of particles accelerated at a perpendicular shock. In our model, the upstream perpendicular magnetic field has no fluctuation, and the downstream region is highly turbulent. Then, the particle motion is the gyration in the upstream region and Bohm-like diffusion downstream. Under this situation, we derive an analytical form of the acceleration time. Using test particle simulations, the validity of our formula is verified. In addition, the energy spectrum of particles is the same as those predicted by standard DSA. Therefore, presently proposed mechanism simultaneously realizes the rapid acceleration and the canonical spectrum, $dN/dp \propto p^{-2}$, even if there is no upstream magnetic amplification.

preprint2020arXiv

Less noticeable shallow decay phase in early X-ray afterglows of GeV/TeV-detected gamma-ray bursts

The nature of the shallow decay phase in the X-ray afterglow of the gamma-ray burst (GRB) is not yet clarified. We analyze the data of early X-ray afterglows of 26 GRBs triggered by Burst Alert Telescope onboard Neil Gehrels Swift Observatory and subsequently detected by Fermi Large Area Telescope (LAT) and/or Imaging Atmospheric Cherenkov Telescopes. It is found that 9 events (including 2 out of 3 very-high-energy gamma-ray events) have no shallow decay phase and that their X-ray afterglow light curves are well described by single power-law model except for the jet break at later epoch. The rest are fitted by double power-law model and have a break in the early epoch (around ks), however, 8 events (including a very-high-energy gamma-ray event) have the pre-break decay index larger than 0.7. We also analyze the data of well-sampled X-ray afterglows of GRBs without LAT detection, and compare their decay properties with those of high-energy and very-high-energy gamma-ray events. It is found that for the GeV/TeV bursts, the fraction of events whose X-ray afterglows are described by single power-law is significantly larger than those for non GeV/TeV GRBs. Even if the GeV/TeV GRBs have shallow decay phase, their decay slope tends to be steeper than non GeV/TeV bursts, that is, they have less noticeable shallow decay phase in the early X-ray afterglow. A possible interpretation along with the energy injection model is briefly discussed.

preprint2020arXiv

Study on the escape timescale of high-energy particles from supernova remnants through thermal X-ray properties

In this decade, GeV/TeV gamma-ray observations of several supernova remnants (SNRs) have implied that accelerated particles are escaping from their acceleration sites. However, when and how they escape from the SNR vicinities are yet to be understood. Recent studies have suggested that the particle escape might develop with thermal plasma ages of the SNRs. In this paper, we present a systematic study on time evolution of particle escape using thermal X-ray properties and gamma-ray spectra. We used 38 SNRs which associate with GeV/TeV gamma-ray emissions. We conducted spectral fittings on the gamma-ray spectra using exponential cutoff power law and broken power law models to estimate the exponential cutoff or the break energies, both of which are indicators of particle escape. The plots of the gamma-ray cutoff/break energies over the plasma ages show similar tendencies to those predicted by simple theories of the particle escape under conditions in which a shock is interacting with thin interstellar medium or clouds. The particle escape timescale is estimated as $\sim$100 kyr from decreasing trends of the total energy of the confined protons with the plasma age. The large dispersions of the cutoff/break energies of the data may suggest an intrinsic variety of particle escape environments. This might be the cause of the complicated Galactic cosmic-ray spectral shape measured on Earth.

preprint2020arXiv

Toward experimental observations of induced Compton scattering by high-power laser facilities

Induced Compton scattering (ICS) is a nonlinear interaction between intense electromagnetic radiation and a rarefied plasma. Although the magnetosphere of pulsars is a potential cite at which ICS occurs in nature, the ICS signatures have not been discovered so far. One of the reasons for non-detection of the ICS signatures is that we still do not attain the concrete understanding of such nonlinear plasma interactions because of their nonlinear nature and of the lack of experimental confirmations. Here, we propose a possible approach to understand ICS experimentally in laboratories, especially, with the use of the up-to-date short-pulse lasers. We find that the scattered light of ICS has characteristic signatures in the spectrum. The signatures will be observed in some current laser facilities. The characteristic spectrum is quantitatively predictable and we can diagnose the properties of the scattering plasma from the signatures.

preprint2019arXiv

Neutrino Emissions from Tidal Disruption Remnants

We study high-energy neutrino emissions from tidal disruption remnants (TDRs) around supermassive black holes. The neutrinos are produced by the decay of charged pions originating in ultrarelativistic protons that are accelerated there. In the standard theory of tidal disruption events (TDEs), there are four distinct phases from debris circularization of stellar debris to super- and sub-Eddington to radiatively inefficient accretion flows (RIAFs). In addition, we consider the magnetically arrested disk (MAD) state in both the super-Eddington accretion and RIAF phases. We find that there are three promising cases to produce neutrino emissions: the super-Eddington accretion phase of the MAD state and the RIAF phases of both the non-MAD and MAD states. In the super-Eddington MAD state, the enhanced magnetic field makes it possible to accelerate the protons to $E_{p,max}~0.35 PeV (M_bh/10^{7.7}M_\odot)^{41/48}$ with the other given appropriate parameters. The neutrino energy is then $E_{ν,pk}~67 TeV (M_bh/10^{7.7}M_\odot)^{41/48}$ at the peak of the energy spectrum. For $M_bh\gtrsim10^{7.7} M_\odot$, the neutrino light curve is proportional to $t^{-65/24}$, while it follows the standard $t^{-5/3}$ decay rate for $M_bh<10^{7.7} M_\odot$. In both cases, the large luminosity and characteristic light curves diagnose the super-Eddington MAD state in TDEs. In the RIAF phase of the non-MAD state, we find $E_{p, max}~0.45 PeV (M_bh/10^7M_\odot)^{5/3}$ and $E_{ν,pk}~0.35 PeV (M_bh/10^7M_\odot)^{5/3}$, and its light curve is proportional to $t^{-10/3}$. This indicates that one can identify whether the existing RIAFs are the TDE origin or not. TDRs are potentially a population of hidden neutrino sources invisible in gamma rays.

preprint2019arXiv

Suzaku and Fermi view of the supernova remnant 3C 396

3C 396 is a composite supernova remnant (SNR), consisting of a central pulsar wind nebula (PWN) and a bright shell in the west, which is known to be interacting with molecular clouds (MCs). We present a study of X-ray emission from the shell and the PWN of the SNR 3C 396 using archival Suzaku data. The spectrum of the SNR shell is clearly thermal, without a signature of a non-thermal component. The abundances of Al and Ca from the shell are slightly enhanced, which indicates the presence of metal-enriched supernova ejecta. The PWN spectra are well described by a power-law model with a photon index of $\sim$1.97 and a thermal component with an electron temperature of $\sim$0.93 keV. The analysis of about 11-years of Fermi data revealed an 18 sigma-detection of gamma-ray emission from the location overlapping with the position of 3C 396 / 4FGL J1903.8+0531. The spectrum of 3C 396 / 4FGL J1903.8+0531 is best-fitted with a log-parabola function with parameters of $α$ = 2.66 and $β$ = 0.16 in the energy range of 0.2$-$300 GeV. The luminosity of 3C 396 / 4FGL J1903.8+0531 was found to be $>$10$^{35}$ erg s$^{-1}$ at 6.2 kpc, which rules out the inverse Compton emission model. Possible scenarios of gamma-ray emission are hadronic emission and bremsstrahlung processes, due to the fact that the SNR is expanding into dense MCs in the western and northern regions of the SNR.

preprint2019arXiv

Weibel-Mediated Shocks Propagating into Inhomogeneous Electron-Positron Plasmas

The external forward shock emitting the gamma-ray burst (GRB) afterglow is collisionless, and it is mediated by the Weibel instability which generates the magnetic field. The GRB afterglow shows that the magnetic field in the large downstream region is much stronger than the shock-compressed pre-shock field. However, particle-in-cell (PIC) simulations of relativistic shocks propagating into homogeneous media show that the Weibel generated field decays near the shock front. Some GRB observations and theoretical studies suggest that the preshock medium is inhomogeneous. We perform the PIC simulation of a relativistic shock propagating into inhomogeneous plasma. It is found that the post-shock magnetic field decays slowly compared with the homogeneous case. Sound waves and entropy waves are also generated by the shock-wave interaction, and temperature anisotropy is produced by the sound wave in the downstream region. The free energy of the temperature anisotropy is large enough to explain the observed field strength. Our results show that the upstream density fluctuation has a significant effect in the downstream region of collisionless shocks even if the wavelength of the upstream inhomogeneity is much larger than the kinetic scale.

preprint2018arXiv

The Transition from Young to Middle-aged Supernova Remnants: Thermal and Nonthermal Aspects of SNR N132D

Supernova remnants (SNRs) are the primary candidate of Galactic cosmic-ray accelerators. It is still an open issue when and how young SNRs, which typically exhibit strong synchrotron X-rays and GeV and TeV gamma-rays, undergo the state transition to middle-aged SNRs dominated by thermal X-rays and GeV gamma-rays. The SNR N132D in the Large Magellanic Cloud is an ideal target to study such a transition, exhibiting bright X-rays and gamma-rays, and with the expected age of ~2500 yrs. In this paper we present results of NuSTAR and Suzaku spectroscopy. We reveal that N132D has a nearly equilibrium plasma with a temperature of > 5 keV or a recombining plasma with a lower temperature (~1.5 keV) and a recombining timescale (net) of 8.8 (7.0--10.0)e12 cm^-3s. Together with the center filled morphology observed in the iron K line image, our results suggest that N132D is now at transition stage from a young SNR to middle-aged. We have constrained the tight upper-limit of nonthermal X-rays. Bright gamma-rays compared to faint nonthermal X-rays suggest that the gamma-rays are hadronic in origin. The spectral energy distribution from radio to gamma-rays shows a proton cut-off energy of ~30 TeV. These facts confirm that N132D is in the transition from young to middle-aged SNR. The large thermal energy of > 10^51 erg and accelerated proton energy of ~ 10^50 erg suggest the supernova explosion might have been very energetic.

preprint2016arXiv

Electromagnetic Afterglows Associated with Gamma-Ray Emission Coincident with Binary Black Hole Merger Event GW150914

The Fermi Gamma-ray Burst Monitor reported the possible detection of the gamma-ray counterpart of a binary black hole merger event, GW150914. We show that the gamma-ray emission is caused by a relativistic outflow with Lorentz factor larger than 10. Subsequently, debris outflow pushes the ambient gas to form a shock, which is responsible for the afterglow synchrotron emission. We find that the 1.4 GHz radio flux peaks at $\sim10^5$ sec after the burst trigger. If the ambient matter is dense enough with density larger than $\sim10^{-2}$ cm$^{-3}$, then the peak radio flux is $\sim0.1$ mJy, which is detectable with radio telescopes such as the Very Large Array. The optical afterglow peaks earlier than the radio, and if the ambient matter density is larger than $\sim0.1$ cm$^{-3}$, the optical flux is detectable with large telescopes such as the Subaru Hyper Suprime-Cam. To reveal the currently unknown mechanisms of the outflow and its gamma-ray emission associated with the binary black hole merger event, follow-up electromagnetic observations of afterglows are important. Detection of the afterglow will localize the sky position of the gravitational wave and the gamma-ray emissions, and it will support the physical association between them.

preprint2016arXiv

Investigating the X-ray and Gamma-ray Properties of the Galactic Supernova Remnants Kes 69, 3C 396, 3C 400.2

Kes 69, 3C 396, and 3C 400.2 are mixed-morphology (MM) Galactic supernova remnants (SNRs), where Kes 69 and 3C 396 are interacting with molecular clouds (MCs). Previous X-ray studies showed that the emission from these SNRs is thermal. It has been suggested that MM SNRs interacting with MCs are potential candidates for recombining plasma (RP) in X-rays and hadronic gamma-ray emission. Recently, Chandra observations revealed signs of RP in 3C 400.2. Our preliminary analyses show that the X-ray emission of NW and SE region of 3C 400.2 arises from recombining plasma. We detected GeV gamma-ray emission from Kes 69 and 3C 396 above 5$σ$.

preprint2016arXiv

Searching for the Time Variation in Supernova Remnant RX J1713.7-3946

Supernova Remnant RX J1713.7-3946 emits synchrotron X-rays and very high energy $γ$-rays. Recently, thermal X-ray line emission is detected from ejecta plasma. CO and HI observations indicate that a highly inhomogeneous medium surrounding the SNR. It is interacting with dense molecular clouds in the northwest and the southwest of the remnant. The origin of the $γ$-ray emission from RX J1713.7-3946 is still uncertain. Detection of rapid variability in X-ray emission from RX J1713.7-3946 indicates the magnetic field $B$ $\sim$ mG. In this work, we investigate the time variation in X-ray flux, luminosity and photon index of RX J1713.7-3946. For this investigation, we study the northwest part of the remnant using Suzaku data in 2006 and 2010. We present preliminary results based on our analysis and interpretations about these X-ray time variability.

preprint2016arXiv

Spatially-Resolved Spectroscopy of a Balmer-Dominated Shock in the Cygnus Loop: An Extremely Thin Cosmic-Ray Precursor?

We present high-resolution long-slit spectroscopy of a Balmer-dominated shock in the northeastern limb of the Cygnus Loop with the Subaru high dispersion spectrograph. By setting the slit angle along the shock normal, we investigate variations of the flux and profile of the H-alpha line from preshock to postshock regions with a spatial resolution of about 4 times 10^{15} cm. The H-alpha line profile can be represented by a narrow (28.9+/-0.7 km/s) Gaussian in a diffuse region ahead of the shock, i.e., a photoionization precursor, and narrow (33.1+/-0.2 km/s) plus broad (130-230 km/s) Gaussians at the shock itself. We find that the width of the narrow component abruptly increases up to 33.1+/-0.2 km/s, or 38.8+/-0.4 km/s if we eliminate projected emission originating from the photoionization precursor, in an unresolved thin layer (< 4 times 10^{15} cm at a distance of 540 pc) at the shock. We show that the sudden broadening can be best explained by heating via damping of Alfven waves in a thin cosmic-ray precursor, although other possibilities are not fully ruled out. The thickness of the cosmic-ray precursor in the Cygnus Loop (a soft gamma-ray emitter) is an order of magnitude thinner than that in Tycho's Knot g (a hard gamma-ray emitter), which may be caused by different energy distribution of accelerated particles between the two sources. In this context, systematic studies might reveal a positive correlation between the thickness of the cosmic-ray precursor and the hardness of the cosmic-ray energy distribution.

preprint2016arXiv

Suzaku observations of the hard X-ray spectrum of Vela Jr

We report the results of Suzaku observations of the young supernova remnant, Vela Jr.\ (RX J0852.0$-$4622), which is known to emit synchrotron X-rays, as well as TeV gamma-rays. Utilizing 39 Suzaku mapping observation data from Vela Jr., a significant hard X-ray emission is detected with the hard X-ray detector (HXD) from the north-west TeV-emitting region. The X-ray spectrum is well reproduced by a single power-law model with the photon index of 3.15$^{+1.18}_{-1.14}$ in the 12--22 keV band. Compiling this with the soft X-ray spectrum simultaneously observed with the X-ray imaging spectrometer (XIS) onboard Suzaku, we find that the wide-band X-ray spectrum in the 2--22 keV band is reproduced with a single power-law or concave broken power-law model, which are statistically consistent with each other. Whichever the model of a single or broken power-law is appropriate, clearly the spectrum has no rolloff structure. Applying this result to the method introduced in \citet{yama2014}, we find that one-zone synchrotron model with electron spectrum having a power-law plus exponential cutoff may not be applicable to Vela Jr.

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

Electron acceleration with improved Stochastic Differential Equation method: cutoff shape of electron distribution in test-particle limit

We develop a method of stochastic differential equation to simulate electron acceleration at astrophysical shocks. Our method is based on Itô's stochastic differential equations coupled with a particle splitting, employing a skew Brownian motion where an asymmetric shock crossing probability is considered. Using this code, we perform simulations of electron acceleration at stationary plane parallel shock with various parameter sets, and studied how the cutoff shape, which is characterized by cutoff shape parameter $a$, changes with the momentum dependence of the diffusion coefficient $β$. In the age-limited cases, we reproduce previous results of other authors, $a\approx2β$. In the cooling-limited cases, the analytical expectation $a\approxβ+1$ is roughly reproduced although we recognize deviations to some extent. In the case of escape-limited acceleration, numerical result fits analytical stationary solution well, but deviates from the previous asymptotic analytical formula $a\approxβ$.

preprint2015arXiv

Extremely Soft X-ray Flash as the indicator of off-axis orphan GRB afterglow

We verified the off-axis jet model of X-ray flashes (XRFs) and examined a discovery of off-axis orphan gamma-ray burst (GRBs) afterglows. The XRF sample was selected on the basis of the following three factors: (1) a constraint on the lower peak energy of the prompt spectrum $E^{src}_{obs}$, (2) redshift measurements, and (3) multi-color observations of an earlier (or brightening) phase. XRF020903 was the only sample selected basis of these criteria. A complete optical multi-color afterglow light curve of XRF020903 obtained from archived data and photometric results in literature showed an achromatic brightening around 0.7 days. An off-axis jet model with a large observing angle (0.21 rad, which is twice the jet opening half-angle, $θ_{jet}$) can naturally describe the achromatic brightening and the prompt X-ray spectral properties. This result indicates the existence of off-axis orphan GRB afterglow light curves. Events with a larger viewing angle ($>\sim2θ_{jet}$) could be discovered using an 8-m class telescope with wide field imagers such as Subaru Hyper-Suprime-Cam and the Large Synoptic Survey Telescope.

preprint2015arXiv

On cosmic-ray production efficiency at supernova remnant shocks propagating into realistic diffuse interstellar medium

Using three-dimensional magnetohydrodynamics simulations, we show that the efficiency of cosmic-ray (CR) production at supernova remnants (SNRs) is over-predicted if it could be estimated based on proper motion measurements of H$α$ filaments in combination with shock-jump conditions. Density fluctuations of upstream medium make shock waves rippled and oblique almost everywhere. The kinetic energy of the shock wave is transferred into that of downstream turbulence as well as thermal energy which is related to the shock velocity component normal to the shock surface. Our synthetic observation shows that the CR acceleration efficiency as estimated from a lower downstream plasma temperature, is overestimated by 10-40%, because rippled shock does not immediately dissipate all upstream kinetic energy.

preprint2015arXiv

Optical inverse-Compton emission from clusters of galaxies

Shocks around clusters of galaxies accelerate electrons which upscatter the Cosmic Microwave Background photons to higher-energies. We use an analytical model to calculate this inverse Compton (IC) emission, taking into account the effects of additional energy losses via synchrotron and Coulomb scattering. We find that the surface brightness of the optical IC emission increases with redshift and halo mass. The IC emission surface brightness, 32--34~mag~arcsec$^{-2}$, for massive clusters is potentially detectable by the newly developed Dragonfly Telephoto Array.

preprint2015arXiv

Turbulent Cosmic Ray Reacceleration at Radio Relics and Halos in Clusters of Galaxies

Radio relics are synchrotron emission found on the periphery of galaxy clusters. From the position and the morphology, it is often believed that the relics are generated by cosmic-ray (CR) electrons accelerated at shocks through a diffusive shock acceleration (DSA) mechanism. However, some radio relics have harder spectra than the prediction of the standard DSA model. One example is observed in the cluster 1RXS J0603.3+4214, which is often called the ``Toothbrush Cluster''. Interestingly, the position of the relic is shifted from that of a possible shock. In this study, we show that these discrepancies in the spectrum and the position can be solved if turbulent (re)acceleration is very effective behind the shock. This means that for some relics turbulent reacceleration may be the main mechanism to produce high-energy electrons, contrary to the common belief that it is the DSA. Moreover, we show that for efficient reacceleration, the effective mean free path of the electrons has to be much smaller than their Coulomb mean free path. We also study the merging cluster 1E 0657-56 or the ``Bullet Cluster'', in which a radio relic has not been found at the position of the prominent shock ahead of the bullet. We indicate that a possible relic at the shock is obscured by the observed large radio halo that is generated by strong turbulence behind the shock. We propose a simple explanation if the morphological differences of radio emission among the Toothbrush, the Bullet, and the Sausage (CIZA J2242.8+5301) Clusters.

preprint2014arXiv

A Hadronic-Leptonic Model for the Fermi Bubbles: Cosmic-Rays in the Galactic Halo and Radio Emission

We investigate non-thermal emission from the Fermi bubbles on a hadronic model. Cosmic-ray (CR) protons are accelerated at the forward shock of the bubbles. They interact with the background gas in the Galactic halo and create $π^0$-decay gamma-rays and secondary electrons through proton-proton interaction. We follow the evolution of the CR protons and electrons by calculating their distribution functions. We find that the spectrum and the intensity profile of $π^0$-decay gamma-rays are consistent with observations. We predict that the shock front is located far ahead of the gamma-ray boundary of the Fermi bubbles. This naturally explains the fact that a clear temperature jump of thermal gas was not discovered at the gamma-ray boundary in recent Suzaku observations. We also consider re-acceleration of the background CRs in the Galactic halo at the shock front. We find that it can significantly affect the gamma-rays from the Fermi bubbles, unless the density of the background CRs is $\lesssim 10$\% of that in the Galactic disk. We indicate that secondary electrons alone cannot produce the observed radio emission from the Fermi bubbles. However, the radio emission from the outermost region of the bubbles can be explained, if electrons are directly accelerated at the shock front with an efficiency of $\sim 0.1$\% of that of protons.

preprint2014arXiv

Detailed Investigation of the Gamma-Ray Emission in the Vicinity of SNR W28 with FERMI-LAT

We present a detailed investigation of the $γ$-ray emission in the vicinity of the supernova remnant (SNR) W28 (G6.4$-$0.1) observed by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope. We detected significant $γ$-ray emission spatially coincident with TeV sources HESS J1800$-$240A, B, and C, located outside the radio boundary of the SNR. Their spectra in the 2-100 GeV band are consistent with the extrapolation of the power-law spectra of the TeV sources. We also identified a new source of GeV emission, dubbed Source W, which lies outside the boundary of TeV sources and coincides with radio emission from the western part of W28. All of the GeV $γ$-ray sources overlap with molecular clouds in the velocity range from 0 to 20 km s$^{-1}$. Under the assumption that the $γ$-ray emission towards HESS J1800-240A, B, and C comes from $π^0$ decay due to the interaction between the molecular clouds and cosmic rays (CRs) escaping from W28, they can be naturally explained by a single model in which the CR diffusion coefficient is smaller than the theoretical expectation in the interstellar space. The total energy of the CRs escaping from W28 is constrained through the same modeling to be larger than $\sim$ 2 $\times$ 10$^{49}$ erg. The emission from Source W can also be explained with the same CR escape scenario.

preprint2014arXiv

Dynamics and microinstabilities at perpendicular collisionless shock: A comparison of large-scale two-dimensional full particle simulations with different ion to electron mass ratio

Large-scale two-dimensional (2D) full particle-in-cell simulations are carried out for studying the relationship between the dynamics of a perpendicular shock and microinstabilities generated at the shock foot. The structure and dynamics of collisionless shocks are generally determined by Alfven Mach number and plasma beta, while microinstabilities at the shock foot are controlled by the ratio of the upstream bulk velocity to the electron thermal velocity and the ratio of the plasma-to-cyclotron frequency. With a fixed Alfven Mach number and plasma beta, the ratio of the upstream bulk velocity to the electron thermal velocity is given as a function of the ion-to-electron mass ratio. The present 2D full PIC simulations with a relatively low Alfven Mach number (M_A ~ 6) show that the modified two-stream instability is dominant with higher ion-to-electron mass ratios. It is also confirmed that waves propagating downstream are more enhanced at the shock foot near the shock ramp as the mass ratio becomes higher. The result suggests that these waves play a role in the modification of the dynamics of collisionless shocks through the interaction with shock front ripples.

preprint2014arXiv

Synchrotron X-ray diagnostics of cutoff shape of nonthermal electron spectrum at young supernova remnants

Synchrotron X-rays can be a useful tool to investigate electron acceleration at young supernova remnants (SNRs). At present, since the magnetic field configuration around the shocks of SNRs is uncertain, it is not clear whether electron acceleration is limited by SNR age, synchrotron cooling, or even escape from the acceleration region. We study whether the acceleration mechanism can be constrained by the cutoff shape of the electron spectrum around the maximum energy. We derive analytical formulae of the cutoff shape in each case where the maximum electron energy is determined by SNR age, synchrotron cooling and escape from the shock. They are related to the energy dependence of the electron diffusion coefficient. Next, we discuss whether information on the cutoff shape can be provided by observations in the near future which will simply give the photon indices and the flux ratios in the soft and hard X-ray bands. We find that if the power-law index of the electron spectrum is independently determined by other observations, then we can constrain the cutoff shape by comparing theoretical predictions of the photon indices and/or the flux ratios with observed data which will be measured by NuSTAR and/or ASTRO-H. Such study is helpful in understanding the acceleration mechanism. In particular, it will supply another independent constraint on the magnetic field strength around the shocks of SNRs.

preprint2014arXiv

X-ray Spectroscopy of the Mixed Morphology Supernova Remnant W28 with XMM-Newton

We report on spatially resolved X-ray spectroscopy of the north-eastern part of the mixed morphology supernova remnant (SNR) W28 with {\it XMM-Newton}. The observed field of view includes a prominent and twisted shell emission forming the edge of this SNR as well as part of the center-filled X-ray emission brightening toward the south-west edge of the field of view. The shell region spectra are in general represented by an optically thin thermal plasma emission in collisional ionization equilibrium with a temperature of $\sim$0.3~keV and a density of $\sim$10~cm$^{-3}$, which is much higher than the density obtained for inner parts. In contrast, we detected no significant X-ray flux from one of the TeV $γ$-ray peaks with an upper-limit flux of 2.1$\times$10$^{-14}$ erg cm$^{-2}$ s$^{-1}$ in the 2--10~keV band. The large flux ratio of TeV to X-ray, larger than 16, and the spatial coincidence of the molecular cloud and the TeV $γ$-ray emission site indicate that the TeV $γ$-ray of W28 is $π^{0}$-decay emission originating from collisions between accelerated protons and molecular cloud protons. Comparing the spectrum in the TeV band and the X-ray upper limit, we obtained a weak upper limit on the magnetic field strength $B\lesssim$ 1500~$μ$G.

preprint2013arXiv

A critical shock Mach number for particle acceleration in the absence of pre-existing cosmic rays: $M=\sqrt 5$

It is shown that, under some generic assumptions, shocks cannot accelerate particles unless the overall shock Mach number exceeds a critical value M > sqrt(5). The reason is that for M <= sqrt(5) the work done to compress the flow in a particle precursor requires more enthalpy flux than the system can sustain. This lower limit applies to situations without significant magnetic field pressure. In case that the magnetic field pressure dominates the pressure in the unshocked medium, i.e. for low plasma beta, the resistivity of the magnetic field makes it even more difficult to fulfil the energetic requirements for the formation of shock with an accelerated particle precursor and associated compression of the upstream plasma. We illustrate the effects of magnetic fields for the extreme situation of a purely perpendicular magnetic field configuration with plasma beta = 0, which gives a minimum Mach number of M = 5/2. The situation becomes more complex, if we incorporate the effects of pre-existing cosmic rays, indicating that the additional degree of freedom allows for less strict Mach number limits on acceleration. We discuss the implications of this result for low Mach number shock acceleration as found in solar system shocks, and shocks in clusters of galaxies.

preprint2013arXiv

Entropy at the Outskirts of Galaxy Clusters as Implications for Cosmological Cosmic-Ray Acceleration

Recently, gas entropy at the outskirts of galaxy clusters attracts much attention. We propose that the entropy profiles could be used to study cosmic-ray (CR) acceleration around the clusters. If the CRs are effectively accelerated at the formation of clusters, the kinetic energy of infalling gas is consumed by the acceleration and the gas entropy should decrease. As a result, the entropy profiles become flat at the outskirts. If the acceleration is not efficient, the entropy should continue to increase outwards. By comparing model predictions with X-ray observations with Suzaku, which show flat entropy profiles, we find that the CRs have carried ~<7% of the kinetic energy of the gas away from the clusters. Moreover, the CR pressure at the outskirts can be ~<40% of the total pressure. On the other hand, if the entropy profiles are not flat at the outskirts as indicated by combined Plank and ROSAT observations, the carried energy and the CR pressure should be much smaller than the above estimations.

preprint2013arXiv

Origin of Radially Aligned Magnetic Fields in Young Supernova Remnants

It has been suggested by radio observations of polarized synchrotron emissions that downstream magnetic field in some young supernova remnants are oriented radially. We study magnetic field distribution of turbulent supernova remnant driven by the Richtmyer-Meshkov instability -- in other words, the effect of rippled shock -- by using three-dimensional magnetohydrodynamics simulations. We find that the induced turbulence has radially biased anisotropic velocity dispersion that leads to a selective amplification of the radial component of the magnetic field. The Richtmyer-Meshkov instability is induced by the interaction between the shock and upstream density fluctuations. Future high-resolution polarization observations can distinguish the following candidates responsible for the upstream density fluctuations: (i) inhomogeneity caused by the cascade of large-scale turbulence in the ISM so-called the big-power-law-in-the-sky, (ii) structures generated by the Drury instability in the cosmic-ray modified shock, and (iii) fluctuations induced by the non-linear feedback of the cosmic-ray streaming instability.

preprint2013arXiv

The Fermi Bubbles as a Scaled-up Version of Supernova Remnants

In this study, we treat the Fermi bubbles as a scaled-up version of supernova remnants (SNRs). The bubbles are created through activities of the super-massive black hole (SMBH) or starbursts at the Galactic center (GC). Cosmic-rays (CRs) are accelerated at the forward shocks of the bubbles like SNRs, which means that we cannot decide whether the bubbles were created by the SMBH or starbursts from the radiation from the CRs. We follow the evolution of CR distribution by solving a diffusion-advection equation, considering the reduction of the diffusion coefficient by CR streaming. In this model, gamma-rays are created through hadronic interaction between CR protons and the gas in the Galactic halo. In the GeV band, we can well reproduce the observed flat distribution of gamma-ray surface brightness, because some amount of gas is left behind the shock. The edge of the bubbles is fairly sharp owing to the high gas density behind the shock and the reduction of the diffusion coefficient there. The latter also contributes the hard gamma-ray spectrum of the bubbles. We find that the CR acceleration at the shock has started when the bubbles were small, and the time-scale of the energy injection at the GC was much smaller than the age of the bubbles. We predict that if CRs are accelerated to the TeV regime, the apparent bubble size should be larger in the TeV band, which could be used to discriminate our hadronic model from other leptonic models. We also present neutrino fluxes.

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

Microinstabilities at perpendicular collisionless shocks: A comparison of full particle simulations with different ion to electron mass ratio

A full particle simulation study is carried out for studying microinstabilities generated at the shock front of perpendicular collisionless shocks. The structure and dynamics of shock waves are determined by Alfven Mach number and plasma beta, while microinstabilities are controlled by the ratio of the upstream bulk velocity to the electron thermal velocity and the plasma-to-cyclotron frequency. Thus, growth rates of microinstabilities are changed by the ion-to-electron mass ratio, even with the same Mach number and plasma beta. The present two-dimensional simulations show that the electron cyclotron drift instability is dominant for a lower mass ratio, and electrostatic electron cyclotron harmonic waves are excited. For a higher mass ratio, the modified two-stream instability is dominant and oblique electromagnetic whistler waves are excited, which can affect the structure and dynamics of collisionless shocks by modifying shock magnetic fields.

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

Search for Diffuse X-rays from the Bow Shock Region of Runaway Star BD+43$^\circ$3654 with Suzaku

The bow shocks of runaway stars with strong stellar winds of over 2000 km s$^{-1}$ can serve as particle acceleration sites. The conversion from stellar wind luminosity into particle acceleration power has an efficiency of the same order of magnitude as those in supernova remnants, based on the radio emission from the bow shock region of runaway star BD+43$^\circ$3654 \citep{Benaglia10}.If this object exhibits typical characteristics, then runaway star systems can contribute a non-negligible fraction of Galactic cosmic-ray electrons. To constrain the maximum energy of accelerated particles from measurements of possible non-thermal emissions in the X-ray band, Suzaku observed BD+43$^\circ$3654 in April 2011 with an exposure of 99 ks. Because the onboard instruments have a stable and low background level, Suzaku detected a possible enhancement over the background of $7.6\pm 3.4$ cnt arcmin$^{-2}$ at the bow shock region, where the error represents the 3 sigma statistics only. However, the excess is not significant within the systematic errors of non-X-ray and cosmic-ray backgrounds of the X-ray Imaging Spectrometer, which are $\pm 6.0$ and $\pm 34$ cnt arcmin$^{-2}$, respectively, and the 3-sigma upper limit in the X-ray luminosity from the shock region, which is $1.1 \times 10^{32}$ erg s$^{-1}$ per 41.2 arcmin$^2$ in the 0.5 to 10 keV band. This result leads to three conclusions: (1) a shock-heating process is inefficient on this system; (2) the maximum energy of electrons does not exceed $\sim$ 10 TeV, corresponding to a Lorentz factor of less than $10^7$; and (3) the magnetic field in the shock acceleration site might not be as turbulent as those in pulsar wind nebulae and supernova remnants.

preprint2012arXiv

Suzaku Observations of the Non-thermal Supernova Remnant HESS J1731-347

A detailed analysis of the nonthermal X-ray emission from the North-Western and Southern parts of the supernova remnant (SNR) HESS J1731$ - $347 with {\it Suzaku} is presented. The shell portions covered by the observations emit hard and line-less X-rays. The spectrum can be reproduced by a simple absorbed power-law model with a photon index $Γ$ of 1.8-2.7 and an absorption column density $N_{\rm H}$ of (1.0-2.1)$\times 10^{22}$ cm$^{-2}$. These quantities change significantly from region to region; the North-Western part of the SNR has the hardest and most absorbed spectrum. The Western part of the X-ray shell has a smaller curvature than North-Western and Southern shell segments. A comparison of the X-ray morphology to the Very High Energy (VHE) gamma-ray and radio images was performed. The efficiency of electron acceleration and emission mechanism in each portion of the shell are discussed. Thermal X-ray emission from the SNR was searched for but could not be detected at a significant level.

preprint2011arXiv

Cross-scale coupling at a perpendicular collisionless shock

A full particle simulation study is carried out on a perpendicular collisionless shock with a relatively low Alfven Mach number (M_A=5). In the present study, we have performed a two-dimensional (2D) electromagnetic full particle simulation with a "shock-rest-frame model". The simulation domain is taken to be larger than the ion inertial length in order to include full kinetics of both electrons and ions. The present simulation result has confirmed the transition of shock structures from the cyclic self-reformation to the quasi-stationary shock front. During the transition, electrons and ions are thermalized in the direction parallel to the shock magnetic field. Ions are thermalized by low-frequency electromagnetic waves (or rippled structures) excited by strong ion temperature anisotropy at the shock foot, while electrons are thermalized by high-frequency electromagnetic waves (or whistler mode waves) excited by electron temperature anisotropy at the shock overshoot. Ion acoustic waves are also excited at the shock overshoot where the electron parallel temperature becomes higher than the ion parallel temperature. We expect that ion acoustic waves are responsible for parallel diffusion of both electrons and ions, and that a cross-scale coupling between an ion-scale mesoscopic instability and an electron-scale microscopic instability is important for structures and dynamics of a collisionless perpendicular shock.

preprint2011arXiv

Evolution of Synchrotron X-rays in Supernova Remnants

A systematic study of the synchrotron X-ray emission from supernova remnants (SNRs) has been conducted. We selected a total of 12 SNRs whose synchrotron X-ray spectral parameters are available in the literature with reasonable accuracy, and studied how their luminosities change as a function of radius. It is found that the synchrotron X-ray luminosity tends to drop especially when the SNRs become larger than ~5 pc, despite large scatter. This may be explained by the change of spectral shape caused by the decrease of the synchrotron roll-off energy. A simple evolutionary model of the X-ray luminosity is proposed and is found to reproduce the observed data approximately, with reasonable model parameters. According to the model, the total energy of accelerated electrons is estimated to be 10^(47-48) ergs, which is well below the supernova explosion energy. The maximum energies of accelerated electrons and protons are also discussed.

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

Relativistic Electron Shock Drift Acceleration in Low Mach Number Galaxy Cluster Shocks

An extreme case of electron shock drift acceleration in low Mach number collisionless shocks is investigated as a plausible mechanism of initial acceleration of relativistic electrons in large-scale shocks in galaxy clusters where upstream plasma temperature is of the order of 10 keV and a degree of magnetization is not too small. One-dimensional electromagnetic full particle simulations reveal that, even though a shock is rather moderate, a part of thermal incoming electrons are accelerated and reflected through relativistic shock drift acceleration and form a local nonthermal population just upstream of the shock. The accelerated electrons can self-generate local coherent waves and further be back-scattered toward the shock by those waves. This may be a scenario for the first stage of the electron shock acceleration occurring at the large-scale shocks in galaxy clusters such as CIZA J2242.8+5301 which has well defined radio relics.

preprint2011arXiv

Toward Understanding the Cosmic-ray Acceleration at Young Supernova Remnants Interacting with Interstellar Clouds: Possible Applications to RX J1713.7-3946

Using three-dimensional magnetohydrodynamics simulations, we investigate general properties of a blast wave shock interacting with interstellar clouds. The shock-cloud interaction generates a turbulent shell through the vorticity generations. In the turbulent shell, the magnetic field is amplified as a result of turbulent dynamo action. In the case of a young supernova remnant, the corresponding strength of the magnetic field is approximately 1 mG. The propagation speed of the shock wave is significantly stalled in the clouds because of the high density, while the shock maintains a high velocity in the diffuse surrounding. In addition, when the shock wave hits the clouds, reflection shock waves are generated that propagate back into the shocked shell. From these simulation results, many observational characteristics of a young SNR RX J1713.7-3946 that is suggested to be interacting with molecular clouds, can be explained as follows: The reflection shocks can accelerate particles in the turbulent downstream region where the magnetic field strength reaches 1mG, which causes short-time variability of synchrotron X-rays. Since the shock velocity is stalled locally in the clouds, the temperature in the shocked cloud is suppressed below 1 keV. Thus, thermal X-ray line emissions would be faint even if the SNR is interacting with molecular clouds. We also find that the photon index of the $π^0$-decay gamma rays generated by cosmic-ray protons can be 1.5, because the penetration depth of high-energy particles into the clumpy clouds depends on their energy, which is consistent with the recent gamma-ray observation by Fermi space telescope.

preprint2010arXiv

Gamma rays from molecular clouds illuminated by cosmic rays escaping from interacting supernova remnants

Recently, gamma-ray telescopes AGILE and Fermi observed several middle-aged supernova remnants (SNRs) interacting with molecular clouds. It is likely that their gamma rays arise from the decay of neutral pions produced by the inelastic collision between cosmic rays (CRs) and nucleons, which suggests that SNRs make the bulk of Galactic CRs. In this letter, we provide the analytical solution of the distribution of CRs that have escaped from a finite-size region, which naturally explains observed broken power-law spectra of the middle-aged SNRs. In addition, the typical value of the break energy of the gamma-ray spectrum, 1-10 GeV, is naturally explained from the fact that the stellar wind dynamics shows the separation between the molecular clouds and the explosion center of about 10 pc. We find that a runaway-CR spectrum of the four middle-aged SNRs (W51C, W28, W44 and IC 443) interacting with molecular clouds could be the same, even though it leads to different gamma-ray spectra. This result is consistent with that of recent studies of the Galactic CR propagation, and supports that SNRs are indeed the sources of Galactic CRs.

preprint2010arXiv

The radiative X-ray and gamma-ray efficiencies of rotation powered pulsars

We present a statistical analysis of the X-ray luminosity of rotation powered pulsars and their surrounding nebulae using the sample of Kargaltsev & Pavlov (2008) and we complement this with an analysis of the gamma-ray-emission of Fermi detected pulsars. We report a strong trend in the efficiency with which spin-down power is converted to X-ray and gamma-ray emission with characteristic age: young pulsars and their surrounding nebulae are efficient X-ray emitters, whereas in contrast old pulsars are efficient gamma-ray emitters. We divided the X-ray sample in a young (Tau < 1.7x10^4 yr) and old sample and used linear regression to search for correlations between the logarithm of the X-ray and gamma-ray luminosities and the logarithms of the periods and period derivatives. The X-ray emission from young pulsars and their nebulae are both consistent with L_X ~ Pdot^3/P^6. For old pulsars and their nebulae the X-ray luminosity is consistent with a more or less constant efficiency eta = L_X/Edot = ~ 8x10^-5. For the gamma-ray luminosity we confirm that L_gamma ~ Edot^(1/2). We discuss these findings in the context of pair production inside pulsar magnetospheres and the striped wind model. We suggest that the striped wind model may explain the similarity between the X-ray properties of the pulsar wind nebulae and the pulsars themselves, which according to the striped wind model may both find their origin outside the light cylinder, in the pulsar wind zone.

preprint2010arXiv

Two-Step Acceleration Model of Cosmic Rays at Middle-Aged SNR

Recent gamma-ray observations of middle-aged supernova remnants revealed a mysterious broken power-law spectrum. Using three-dimensional magnetohydrodynamics simulations, we show that the interaction between a supernova blast wave and interstellar clouds formed by thermal instability generates multiple reflected shocks. The typical Mach numbers of the reflected shocks are shown to be M ~ 2 depending on the density contrast between the diffuse intercloud gas and clouds. These secondary shocks can further energize cosmic-ray particles originally accelerated at the blast-wave shock. This "two-step" acceleration scenario reproduces the observed gamma-ray spectrum and predicts the high-energy spectral index ranging approximately from 3 to 4.

preprint2009arXiv

Escape-limited Model of Cosmic-ray Acceleration Revisited

The spectrum of cosmic rays (CRs) is affected by their escape from an acceleration site. This may have been observed not only in the gamma-ray spectrum of young supernova remnants (SNRs) such as RX J1713.7-3946, but also in the spectrum of CRs showering on the Earth. The escape-limited model of cosmic-ray acceleration is studied in general. We discuss the spectrum of CRs running away from the acceleration site. The model may also constrain the spectral index at the acceleration site and the ansatz with respect to the unknown injection process into the particle acceleration. We apply our model to CR acceleration in SNRs and in active galactic nuclei (AGN), which are plausible candidates of Galactic and extragalactic CRs, respectively. In particular, for young SNRs, we take account of the shock evolution with cooling of escaping CRs in the Sedov phase. The spectrum of escaping CRs generally depends on the physical quantities at the acceleration site, such as the spectral index, the evolution of the maximum energy of CRs and the evolution of the number of CRs. It is found that the spectrum of run-away particles can be both softer and harder than that of the acceleration site. The model could explain spectral indices of both Galactic and extragalactic CRs produced by SNRs and AGNs, respectively, suggesting the unified picture of CR acceleration.

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

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.

preprint2009arXiv

The Nature of a Cosmic-Ray Accelerator, CTB37 B, Observed with Suzaku and Chandra

We report on Suzaku and Chandra observations of the young supernova remnant CTB37B, from which TeV gamma-rays were detected by the H.E.S.S. Cherenkov telescope. The 80 ks Suzaku observation provided us with a clear image of diffuse emission and high-quality spectra. The spectra revealed that the diffuse emission is comprised of thermal and non-thermal components. The thermal component can be represented by an NEI model with a temperature, a pre-shock electron density and an age of 0.9(0.7-1.1) keV, 0.4(0.3-0.5) cm^{-3} and 650(350-3150) yr, respectively. This suggests that the explosion of CTB37B occurred in a low-density space. A non-thermal power-law component was found from the southern region of CTB37B. Its photon index of ~1.5 and a high roll-off energy (>15 keV) indicate efficient cosmic-ray acceleration. A comparison of this X-ray spectrum with the TeV gamma-ray spectrum leads us to conclude that the TeV gamma-ray emission seems to be powered by either multi-zone Inverse Compton scattering or the decay of neutral pions. The point source resolved by Chandra near the shell is probably associated with CTB37B, because of the common hydrogen column density with the diffuse thermal emission. Spectral and temporal characteristics suggest that this source is a new anomalous X-ray pulsar.

preprint2009arXiv

Turbulence and Magnetic Field Amplification in Supernova Remnants: Interactions Between A Strong Shock Wave and Multi-Phase Interstellar Medium

We examine MHD simulations of the propagation of a strong shock wave through the interstellar two-phase medium composed of small-scale cloudlets and diffuse warm neutral medium in two-dimensional geometry. The pre-shock two-phase medium is provided as a natural consequence of the thermal instability that is expected to be ubiquitous in the interstellar medium. We show that the shock-compressed shell becomes turbulent owing to the preshock density inhomogeneity and magnetic field amplification takes place in the shell. The maximum field strength is determined by the condition that plasma beta ~ 1, which gives the field strength on the order of 1 mG in the case of shock velocity ~ 1,000 km/s. The strongly magnetized region shows filamentary and knot-like structures in two-dimensional simulations. The spatial scale of the regions with magnetic field of 1 mG in our simulation is roughly 0.05 pc which is comparable to the spatial scale of the X-ray hot spots recently discovered in supernova remnants where the magnetic field strength is indicated to be amplified up to the order of 1 mG. This result may also suggest that the turbulent region with locally strong magnetic field is expected to be spread out in the region with frequent supernova explosions, such as in the Galactic center and starburst galaxies.

preprint2008arXiv

Electron acceleration at a low-Mach-number perpendicular collisionless shock

A full particle simulation study is carried out on the electron acceleration at a collisionless, relatively low Alfven Mach number (M_A=5), perpendicular shock. Recent self-consistent hybrid shock simulations have demonstrated that the shock front of perpendicular shocks has a dynamic rippled character along the shock surface of low-Mach-number perpendicular shocks. In this paper, the effect of the rippling of perpendicular shocks on the electron acceleration is examined by means of large-scale (ion-scale) two-dimensional full particle simulations. It has been shown that a large-amplitude electric field is excited at the shock front in association with the ion-scale rippling, and that reflected ions are accelerated upstream at a localized region where the shock-normal electric field of the rippled structure is polarized upstream. The current-driven instability caused by the highly-accelerated reflected ions has a high growth rate to large-amplitude electrostatic waves. Energetic electrons are then generated by the large-amplitude electrostatic waves via electron surfing acceleration at the leading edge of the shock transition region. The present result suggests that the electron surfing acceleration is also a common feature at low-Mach-number perpendicular collisionless shocks.

preprint2008arXiv

Gamma-ray spectrum of RX J1713.7-3946 in the Fermi era and future detection of neutrinos

The recently launched satellite, Fermi Gamma-ray Space Telescope, is expected to find out if cosmic-ray (CR) protons are generated from supernova remnants (SNRs), especially RX J1713.7-3946, by observing the GeV-to-TeV gamma-rays. The GeV emission is thought to be bright if the TeV emission is hadronic, i.e., of proton origin, while dim if leptonic. We reexamine the above view using a simple theoretical model of nonlinear acceleration of particles to calculate the gamma-ray spectrum of Galactic young SNRs. If the nonlinear effects of CR acceleration are considered, it may be impossible to distinguish the evidence of proton acceleration from leptonic in the gamma-ray spectrum of Galactic young SNRs like RX J1713.7-3946. On the other hand, future km^3-class neutrino observations will likely find a clear evidence of the proton acceleration there.

preprint2008arXiv

X-Ray Observation of Very High Energy Gamma-Ray Source, HESS J1745-303, with Suzaku

Suzaku observations of a TeV unidentified (unID) source, HESS J1745-303, are presented. A possible excess of neutral iron line emission is discovered, and is likely associated with the main part of HESS J1745-303, named "region A". It may be an X-ray reflection nebula where the X-rays from previous Galactic Center (GC) activity are reflected by a molecular cloud. This result further strengthens the assumption that the molecular cloud which is spatially coincident with region A of HESS J1745-303 is located in the GC region. The TeV emission from molecular clouds is reminiscent of the diffuse TeV gamma-rays from the GC giant molecular clouds, and it could have the same emission mechanism. With deep exposure mapping observations by Suzaku, a tight upper-limit on the 2-10 keV continuum diffuse emission from region A is obtained, as 2.1x10^-13ergs s^-1cm^-2. The flux ratio between 1-10 TeV and 2-10 keV is larger than 4. Possible scenarios to reproduce wide-band spectra from keV to TeV are examined. Thermal X-rays from nearby two old supernova remnants, G359.0-0.9 and G359.1-0.5, are detected, and their emission properties are well determined in the present study with deep exposure.

preprint2006arXiv

Full particle simulation of a perpendicular collisionless shock: A shock-rest-frame model

The full kinetic dynamics of a perpendicular collisionless shock is studied by means of a one-dimensional electromagnetic full particle simulation. The present simulation domain is taken in the shock rest frame in contrast to the previous full particle simulations of shocks. Preliminary results show that the downstream state falls into a unique cyclic reformation state for a given set of upstream parameters through the self-consistent kinetic processes.

preprint2006arXiv

Possibility of Ultra High-Energy Cosmic Rays from the Giant Flare in SGR 1806-20

On 2004 December 27, a giant flare from the soft gamma repeater 1806$-$20 was observed. The radiation mechanism of the initial peak of the flare would be controversial. In this letter we point out that very high-energy cosmic rays would be produced in the case that the flare was caused by internal shocks, as is usually considered for gamma-ray bursts. The highest energy of cosmic rays can reach $10^{19}$ eV, if the Lorentz factor of the shocks is sufficiently high. Future observations of cosmic rays will inform us about the mechanism of the giant flare.

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.

preprint2006arXiv

The X-ray synchrotron emission of RCW 86 and the implications for its age

We report here X-ray imaging spectroscopy observations of the northeastern shell of the supernova remnant RCW 86 with Chandra and XMM-Newton. Along this part of the shell the dominant X-ray radiation mechanism changes from thermal to synchrotron emission. We argue that both the presence of X-ray synchrotron radiation and the width of the synchrotron emitting region suggest a locally higher shock velocity of V_s = 2700 km/s and a magnetic field of B = 24+/-5 microGauss. Moreover, we also show that a simple power law cosmic ray electron spectrum with an exponential cut-off cannot explain the broad band synchrotron emission. Instead a concave electron spectrum is needed, as predicted by non-linear shock acceleration models. Finally, we show that the derived shock velocity strengthens the case that RCW 86 is the remnant of SN 185.

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

Constraints On the Diffusive Shock Acceleration From the Nonthermal X-ray Thin Shells In SN1006 NE Rim

Characteristic scale lengths of nonthermal X-rays from the SN1006 NE rim, which are observed by Chandra, are interpreted in the context of the diffusive shock acceleration on the assumption that the observed spatial profile of nonthermal X-rays corresponds to that of accelerated electrons with energies of a few tens of TeV. To explain the observed scale lengths, we construct two simple models with a test particle approximation, where the maximum energy of accelerated electrons is determined by the age of SN1006 (age-limited model) or the energy loss (energy loss-limited model), and constrain the magnetic field configuration and the diffusion coefficients of accelerated electrons. When the magnetic field is nearly parallel to the shock normal, the magnetic field should be in the range of 20-85 micro Gauss and highly turbulent both in upstream and downstream, which means that the mean free path of accelerated electrons is on the order of their gyro-radius (Bohm limit). This situation can be realized both in the age-limited and energy loss-limited model. On the other hand, when the magnetic field is nearly perpendicular to the shock normal, which can exist only in the age-limited case, the magnetic field is several micro Gauss in the upstream and 14-20 micro Gauss in the downstream, and the upstream magnetic field is less turbulent than the downstream.

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.