BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Yutaka Ohira

Yutaka Ohira appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-ph.COastro-ph.GAhep-phphysics.plasm-phastro-ph

Catalog footprint

What is connected

56works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

56 published item(s)

preprint2026arXiv

Cosmic-ray Acceleration in Core-Collapse Supernova Remnants with the Wind Termination Shock

We investigate the attainable maximum energy of particles accelerated in the core-collapse supernova remnant (SNR) shock propagating in the free wind region with the Parker-spiral magnetic field, current sheet, and the wind termination shock (WTS) by using test particle simulations. This work focuses on Wolf-Rayet stars as progenitors. The magnetic field amplification in the free wind region (shock upstream region) is not considered in this work. Test particle simulations show that particles escaped from the core-collapse SNR reach and move along the WTS, and eventually return to the SNR shock from the poles or equator of the WTS. The particle attainable energy can be boosted by this cyclic motion between the SNR shock and WTS and can be larger than the particle energy that is limited by escape from the SNR shock. The particle energy limited by the cyclic motion between the SNR shock and WTS is about $10-100~{\rm TeV}$. Thus, the core-collapse SNR without upstream magnetic field amplification can be the origin of the break around $10~{\rm TeV}$ of the energy spectrum of observed cosmic ray protons and helium.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Biermann Battery Powered by Resistive Heating Induced by Cosmic Ray Streaming

It is recently proposed that cosmic rays generate a seed magnetic field in the early universe. In this paper, we propose another generation mechanism of magnetic fields by cosmic rays, which is the Biermann battery driven by resistive heating induced by the streaming of cosmic rays. This mechanism is dominant in small-scale, low-temperature, and strongly-ionized regions, compared with other previously proposed mechanisms. Because cosmic rays are expected to be accelerated after the death of the first stars, this mechanism can work during structure formation in the early universe. We show that it makes the seed magnetic field with sufficient strength for the subsequent dynamo to amplify it to the micro Gauss level in the current galaxies.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Magnetic field generation by charge exchange in a supernova remnant in the early universe

We present new generation mechanisms of magnetic fields in supernova remnant shocks propagating to partially ionized plasmas in the early universe. Upstream plasmas are dissipated at the collisionless shock, but hydrogen atoms are not dissipated because they do not interact with electromagnetic fields. After the hydrogen atoms are ionized in the shock downstream region, they become cold proton beams that induce the electron return current. The injection of the beam protons can be interpreted as an external force acting on the downstream proton plasma. We show that the effective external force and the electron return current can generate magnetic fields without any seed magnetic fields. The magnetic field strength is estimated to be $B\sim 10^{-14}-10^{-11}~{\rm G}$, where the characteristic lengthscale is the mean free path of charge exchange, $\sim 10^{15}~{\rm cm}$. Since protons are marginally magnetized by the generated magnetic field in the downstream region, the magnetic field could be amplified to larger values and stretched to larger scales by turbulent dynamo and expansion.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Self-discharge by streaming cosmic rays

A new nonthermal phenomenon caused by streaming cosmic rays (CRs) in the universe is proposed. The streaming CRs drive the return current of thermal electrons to compensate for the CR current. Then, electric fields are induced by the resistivity of the return current. It is shown that the resistive electric fields can accelerate secondary electrons generated by the streaming CRs. This is the self-discharge by streaming CRs. In this work, the self-discharge condition and the condition for runaway acceleration of secondary electrons are presented. The self-discharge makes high-energy secondary electrons, resulting in enhancements of ionization and nonthermal emission including K$α$ emission line of neutral irons. After the self-discharge, the return current of thermal electrons is replaced by the electric current of secondary electrons. Since some magnetic field generations and amplifications are driven by the return current of thermal electrons, the self-discharge can significantly influence them.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetic field generation by an inhomogeneous return current

A new generation mechanism of the magnetic field in an inhomogeneous collisionless plasma with a beam component is proposed. We show that even though the current and charge neutralities are initially satisfied, the current neutrality is eventually violated if there is an inhomogeneity, so that the magnetic field is generated. By conducting ab initio two-dimensional particle-in-cell simulations, we demonstrate that the magnetic field is generated as expected. The new generation mechanism of the magnetic field can play an important role in the current universe because cosmic rays can be regarded as the beam component in the astrophysical plasma. We propose that the first cosmic rays generate the magnetic field with a large scale at the redshift of z ~ 20.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Particle acceleration in a shock wave propagating to an inhomogeneous medium

We investigate effects of upstream density fluctuations on the diffusive shock acceleration by Monte Carlo simulations. The simulations show that particles are reaccelerated in the shock downstream region by a sound wave generated at the shock front. The time scale of turbulent acceleration by the sound wave is estimated. We propose a new back reaction of particles accelerated around the shock front. The accelerated particles generate the upstream density fluctuations by the Drury instability, which are converted to the downstream sound waves by the shock. The downstream sound waves modify the momentum spectrum of particles accelerated around the shock front. This new back reaction affects emission from the accelerated particles, which gives a new constraint on the acceleration efficiency of the diffusive shock acceleration.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cosmic-Ray Hardenings in the Light of AMS-02

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Injection to rapid diffusive shock acceleration at perpendicular shocks in partially ionized plasmas

We present a three-dimensional hybrid simulation of a collisionless perpendicular shock in a partially ionized plasma for the first time. In this simulation, the shock velocity and the upstream ionization fraction are Vsh ~ 1333 km/s and fi ~ 0.5, that are typical values for isolated young supernova remnants in the interstellar medium. We confirm previous two-dimensional simulation results that downstream hydrogen atoms leak into the upstream region, they are accelerated by the pickup process in the upstream region, and large magnetic field fluctuations are generated both in the upstream and downstream regions. In addition, we find that the magnetic field fluctuations have three-dimensional structures and the leaking hydrogen atoms are injected to the diffusive shock acceleration at the perpendicular shock after the pickup process. The observed diffusive shock acceleration can be interpreted as the shock drift acceleration with scattering. Particles are accelerated to v ~ 100 Vsh ~ 0.3c within ~ 100 gyroperiods in this simulation. The acceleration time scale is faster than that of the diffusive shock acceleration in parallel shocks. Our simulation results suggest that supernova remnants can accelerate cosmic rays to 10^{15.5} eV (the knee) during the Sedov phase.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Weibel instability driven by spatially anisotropic density structures

Observations of afterglows of gamma-ray bursts suggest (GRBs) that post-shock magnetic fields are strongly amplified to about 100 times the shock-compressed value. The Weibel instability appears to play an important role in generating of the magnetic field. However, recent simulations of collisionless shocks in homogeneous plasmas show that the magnetic field generated by the Weibel instability rapidly decays. There must be some density fluctuations in interstellar and circumstellar media. The density fluctuations are anisotropically compressed in the downstream region of relativistic shocks. In this paper, we study the Weibel instability in electron--positron plasmas with the spatially anisotropic density distributions by means of two-dimensional particle-in-cell simulations. We find that large magnetic fields are maintained for a longer time by the Weibel instability driven by the spatially anisotropic density structure. Particles anisotropically escape from the high density region, so that the temperature anisotropy is generated and the Weibel instability becomes unstable. Our simulation results suggest that the Weibel instability driven by an anisotropic density structure can generate sufficiently large magnetic fields and they can cover sufficiently large regions to explain the afterglow emission of GRBs.

0 citations0 reviewsNo code linkedOpen access
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β$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Evidence for Thermal X-Ray Line Emission from the Synchrotron-Dominated Supernova Remnant RX J1713.7-3946

We report the first detection of thermal X-ray line emission from the supernova remnant (SNR) RX J1713.7-3946, the prototype of the small class of synchrotron dominated SNRs. A softness-ratio map generated using XMM-Newton data shows that faint interior regions are softer than bright shell regions. Using Suzaku and deep XMM-Newton observations, we have extracted X-ray spectra from the softest area, finding clear line features at 1 keV and 1.35 keV. These lines can be best explained as Ne Ly-alpha and Mg He-alpha from a thermal emission component. Since the abundance ratios of metals to Fe are much higher than solar values in the thermal component, we attribute the thermal emission to reverse-shocked SN ejecta. The measured Mg/Ne, Si/Ne, and Fe/Ne ratios of 2.0-2.6, 1.5-2.0, and <0.05 solar suggest that the progenitor star of RX J1713.7-3946 was a relatively low-mass star (<~20 M_sun), consistent with a previous inference based on the effect of stellar winds of the progenitor star on the surrounding medium. Since the mean blastwave speed of ~6000 km/s (the radius of 9.6 pc divided by the age of 1600 yr) is relatively fast compared with other core-collapse SNRs, we propose that RX J1713.7-3946 is a result of a Type Ib/c supernova whose progenitor was a member of an interacting binary. While our analysis provides strong evidence for X-ray line emission, our interpretation of its nature as thermal emission from SN ejecta requires further confirmation especially through future precision spectroscopic measurements using ASTRO-H.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Magnetic field amplification by collisionless shocks in partially ionized plasmas

In this paper, we study shock structures of collisionless shocks in partially ionized plasmas by means of two-dimensional hybrid simulations, where the shock is a perpendicular shock with shock velocity Vsh ~ 40 Va ~ 1333 km/s and the upstream ionization fraction is 0.5. We find that large density fluctuations and large magnetic fields fluctuations are generated both in the upstream and downstream regions. In addition, we find that the velocity distribution of downstream hydrogen atoms has three components. Observed shock structures suggest that diffusive shock acceleration can operate at perpendicular shocks propagating into partially ionized plasmas in real three-dimensional systems.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Acoustic instability in the neutral precursor region of collisionless shocks propagating into partially ionized plasmas

Recent studies about collisionless shocks in partially ionized plasmas showed that some of neutral particles leak into the shock upstream region from the downstream region. In this paper, we perform a linear analysis and show that acoustic waves are unstable in the neutral precursor region. The acoustic instability amplifies fluctuations of magnetic field and density in the upstream region. The fluctuations are indispensable for the diffusive shock acceleration and could be important for the downstream turbulence.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamics of X-Ray-Emitting Ejecta in the Oxygen-Rich Supernova Remnant Puppis A Revealed by the XMM-Newton RGS

Using the unprecedented spectral resolution of the reflection grating spectrometer (RGS) onboard XMM-Newton, we reveal dynamics of X-ray--emitting ejecta in the oxygen-rich supernova remnant Puppis A. The RGS spectrum shows prominent K-shell lines, including O VII He_alpha forbidden and resonance, O VIII Ly_alpha, O VIII Ly_beta, and Ne IX He_alpha resonance, from an ejecta knot positionally coincident with an optical oxygen-rich filament (the so-called Omega filament in the northeast of the remnant). We find that the line centroids are blueshifted by 1480+/-140+/-60 km/s (the first and second term errors are measurement and calibration uncertainties, respectively), which is fully consistent with that of the optical Omega filament. Line broadening at 654 eV (corresponding to O VIII Ly_alpha) is obtained to be less than ~0.9 eV, indicating an oxygen temperature of less than ~30 keV. Analysis of XMM-Newton MOS spectra shows an electron temperature of ~0.8 keV and an ionization timescale of ~2e10 cm^-3 s. We show that the oxygen and electron temperatures as well as the ionization timescale can be reconciled if the ejecta knot was heated by a collisionless shock whose velocity is ~600-1200 km/s and was subsequently equilibrated due to Coulomb interactions. The RGS spectrum also shows relatively weak K-shell lines of another ejecta feature located near the northeastern edge of the remnant, from which we measure redward Doppler velocities of 650+/-70+/-60 km/s.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulations of Collisionless Perpendicular Shocks in Partially Ionized Plasmas

Perpendicular collisionless shocks propagating into partially ionized plasmas are investigated by two-dimensional hybrid particle simulations. It is shown that some neutral particles leak into the upstream region from the downstream region, the leaking neutral particles become pickup ions in the upstream region and modify the shock structure, the pickup ions are preferentially accelerated, and plasma instabilities are excited by the pickup ions in the upstream and downstream regions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stability Analysis for Cosmic-Ray Heating of Cool Cores in Galaxy Clusters

We study the heating of the cool cores in galaxy clusters by cosmic-rays (CRs) accelerated by the central active galactic nuclei (AGNs). We especially focus on the stability of the heating. The CRs stream with Alfvén waves in the intracluster medium (ICM) and heat the ICM. First, assuming that the heating and radiative cooling is balanced, we search steady state solutions for the ICM and CR profiles of clusters by solving a boundary value problem. The boundary conditions are set so that the solutions are consistent with observations of clusters. We find steady state solutions if the magnetic fields are strong enough and the association between the magnetic fields and the ICM is relatively weak. Then, we analyse the stability of the solutions via a Lagrangian perturbation analysis and find that the solutions are globally stable. We confirm the results by numerical simulations. Using the steady state solutions as the initial conditions, we follow the evolution of the profiles for 100 Gyr. We find that the profiles do not evolve on time scales much larger than cluster lifetimes. These results, as well as consistency with observations of radio mini-halos, suggest that the CR heating is a promising mechanism to solve the so-called "cooling flow problem".

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Turbulent Shear Acceleration

We consider particle acceleration by large-scale incompressible turbulence with a lengthscale larger than the particle mean free path. We derive an ensemble-averaged transport equation of energetic charged particles from an extended transport equation which contains the shear acceleration. The ensemble-averaged transport equation describes particle acceleration by incompressible turbulence (turbulent shear acceleration). We find that for Kolmogorov turbulence, the turbulent shear acceleration becomes important in small scale. Moreover, by Monte Carlo simulations, we confirm that the ensemble-averaged transport equation describes the turbulent shear acceleration.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effects of Leakage Neutral Particles on Shocks

In this paper we investigate effects of neutral particles on shocks propagating into the partially ionized medium. We find that for 120 km/s < u_{sh} < 3000 km/s (u_{sh} is the shock velocity), about ten percent of upstream neutral particles leak into the upstream region from the downstream region. Moreover, we investigate how the leakage neutral particles affect the upstream structure of the shock and particle accelerations. Using four fluid approximations (upstream ions, upstream neutral particles, leakage neutral particles and pickup ions), we provide analytical solutions of the precursor structure due to leakage neutral particles. It is shown that the upstream flow is decelerated in the precursor region and the shock compression ratio becomes smaller than without leakage neutral particles, but the total compression ratio does not change. Even if leakage of neutral particles is small (a few percents of total upstream particles), this smaller compression ratio of the shock can explain steep gamma-ray spectra from young supernova remnants. Furthermore, leakage neutral particles could amplify the magnetic field and heat the upstream region.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Positron annihilation as a cosmic-ray probe

Recently, the gamma-ray telescopes AGILE and Fermi observed several middle-aged supernova remnants (SNRs) interacting with molecular clouds. A plausible emission mechanism of the gamma rays is the decay of neutral pions produced by cosmic ray (CR) nuclei (hadronic processes). However, observations do not rule out contributions from bremsstrahlung emission due to CR electrons. TeV gamma-ray telescopes also observed many SNRs and discovered many unidentified sources. It is still unclear whether the TeV gamma-ray emission is produced via leptonic processes or hadronic processes. In this Letter, we propose that annihilation emission of secondary positrons produced by CR nuclei is a diagnostic tool of the hadronic processes. We investigate MeV emissions from secondary positrons and electrons produced by CR protons in molecular clouds. The annihilation emission of the secondary positrons from SNRs can be robustly estimated from the observed gamma-ray flux. The expected flux of the annihilation line from SNRs observed by AGILE and Fermi is sufficient for the future Advanced Compton Telescope to detect. Moreover, synchrotron emission from secondary positrons and electrons and bremsstrahlung emission from CR protons can be also observed by the future X-ray telescope NuSTAR and ASTRO-H.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radio Mini-Halo Emission from Cosmic Rays in Galaxy Clusters and Heating of the Cool Cores

It has been proposed that the cool cores of galaxy clusters are stably heated by cosmic rays (CRs). If this is the case, radio mini-halos, which are often found in the central regions of cool core clusters, may be attributed to the synchrotron emission from the CRs. Based on this idea, we investigate the radial profiles of the mini-halos. First, using numerical simulations, we confirm that it is appropriate to assume that radiative cooling of the intracluster medium (ICM) is balanced with the heating by CR streaming. In these simulations, we assume that the streaming velocity of the CRs is the sound velocity of the ICM, and indicate that the heating is even more stable than the case where the streaming velocity is the Alfven velocity. Then, actually assuming the balance between cooling and heating, we estimate the radial profiles of CR pressure in six clusters only from X-ray observations. Since the CR protons interact with the ICM protons, we can predict the radial profiles of the resultant synchrotron radiation. We compare the predictions with the observed radial profiles of the mini-halos in the six clusters and find that they are consistent if the momentum spectra of the CRs are steep. These results may indicate that the cores are actually being heated by the CRs. We also predict broad-band spectra of the six clusters, and show that the non-thermal fluxes from the clusters are small in hard X-ray and gamma-ray bands.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Alfven Wave Amplification and Self-Containment of Cosmic-Rays Escaping from a Supernova Remnant

We study the escape of cosmic-ray (CR) protons accelerated at a supernova remnant (SNR) by numerically solving a diffusion-convection equation from the vicinity of the shock front to the region far away from the front. We consider the amplifications of Alfven waves generated by the escaping CR particles and their effects on CR escape into interstellar medium (ISM). We find that the amplification of the waves significantly delays the escape of the particles even far away from the shock front (on a scale of the SNR). This means that the energy spectrum of CR particles measured through gamma-ray observations at molecular clouds around SNRs is seriously affected by the particle scattering by the waves.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Cosmic Ray Helium Hardening

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gamma-Ray Burst without Baryonic and Magnetic Load?

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Non-thermal Emissions from Cool Cores Heated by Cosmic-Rays in Galaxy Clusters

We study non-thermal emissions from cool cores in galaxy clusters. We adopted a recent model, in which cosmic-rays (CRs) prevail in the cores and stably heat them through CR streaming. The non-thermal emissions come from the interaction between CR protons and intracluster medium (ICM). Comparison between the theoretical predictions and radio observations shows that the overall CR spectra must be steep, and most of the CRs in the cores are low-energy CRs. Assuming that the CRs are injected through AGN activities, we study the nature of the shocks that are responsible for the CR acceleration. The steep CR spectra are likely to reflect the fact that the shocks travel in hot ICM with fairly small Much numbers. We also study the dependence on the CR streaming velocity. The results indicate that synchrotron emissions from secondary electrons should be observed as radio mini-halos in the cores. In particular, low-frequency observations (e.g. LOFAR) are promising. On the other hand, the steepness of the spectra makes it difficult to detect non-thermal X-ray and gamma-ray emissions from the cores. The low-energy CRs may be heating optical filaments observed in the cores.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stable Heating of Cluster Cooling Flows by Cosmic-Ray Streaming

We study heating of cool cores in galaxy clusters by cosmic-ray (CR) streaming using numerical simulations. In this model, CRs are injected by the central active galactic nucleus (AGN) and move outward with Alfven waves. The waves are excited by the streaming itself and become non-linear. If magnetic fields are large enough, CRs can prevail in and heat the entire core because of a large Alfven velocity. We find that the CR streaming can stably heat both high and low temperature clusters for a long time without the assistance of thermal conduction, and it can prevent the development of massive cooling flows. If there is even minor contribution of thermal conduction, the heating can be more stabilized. We discuss the reason of the stability and indicate that the CR pressure is insensitive to the change of intracluster medium (ICM), and that the density dependence of the heating term is similar to that of the radiative cooling.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Effects of Neutral Particles on Modified Shocks at Supernova Remnants

H-alpha emission from supernova remnants (SNRs) implies the existence of neutral hydrogens in the ambient medium. In the precursor of an SNR shock modified by cosmic rays (CRs), upstream plasmas are pushed by the CR pressure, but neutral particles are not, so that the relative velocity appears and some neutral particles become pickup ions by the charge exchange process in the precursor. We investigate how the pickup ions generated in the precursor affect the shock structure and the particle acceleration. If the CR pressure is larger than 20% of the shock ram pressure, the compression of the subshock becomes smaller than that without pickup ions because of the pressure of the pickup ions. Moreover, even if the shock is modified by CRs, the total compression ratio can be smaller than 4. In addition, the pickup ions play an important role for the injection into the particle acceleration. If the shock is a quasi-perpendicular shock and if the the multiply reflected ion acceleration occurs, the CR spectrum can be harder than that of the test particle diffusive shock acceleration below GeV.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Slow Diffusion of Cosmic-Rays around a Supernova Remnant

We study the escape of cosmic-ray protons accelerated at a supernova remnant (SNR). We are interested in their propagation in interstellar medium (ISM) after they leave the shock neighborhood where they are accelerated, but when they are still near the SNR with their energy density higher than that in the average ISM. Using Monte-Carlo simulations, we found that the cosmic-rays with energies of <~TeV excite Alfven waves around the SNR on a scale of the SNR itself if the ISM is highly ionized. Thus, even if the cosmic-rays can leave the shock, scattering by the waves prevents them from moving further away from the SNR. The cosmic-rays form a slowly expanding cosmic-ray bubble, and they spend a long time around the SNR. This means that the cosmic-rays cannot actually escape from the SNR until a fairly late stage of the SNR evolution. This is consistent with some results of Fermi and H.E.S.S. observations.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Molecular Clouds as a Probe of Cosmic-Ray Acceleration in a Supernova Remnant

We study cosmic-ray acceleration in a supernova remnant (SNR) and the escape from it. We model nonthermal particle and photon spectra for the hidden SNR in the open cluster Westerlund 2, and the old-age mixed-morphology SNR W 28. We assume that the SNR shock propagates in a low-density cavity, which is created and heated through the activities of the progenitor stars and/or previous supernova explosions. We indicate that the diffusion coefficient for cosmic-rays around the SNRs is less than ~1% of that away from them. We compare our predictions with the gamma-ray spectra of molecular clouds illuminated by the cosmic-rays (Fermi and H.E.S.S.). We found that the spectral indices of the particles are ~2.3. This may be because the particles were accelerated at the end of the Sedov phase, and because energy dependent escape and propagation of particles did not much affect the spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Plasma instabilities as a result of charge exchange in the downstream region of supernova remnant shocks

H-alpha emission from supernova remnants (SNRs) implies the existence of neutral hydrogen in the circumstellar medium. Some of the neutral particles penetrating the shock are ionized by the charge exchange process and make a cold ion beam in the shock downstream region. We perform linear analyses of collisionless plasma instabilities between the cold beam and the hot downstream plasma. We find that, under typical SNR conditions, either the resonant instability or the Weibel instability is the most unstable. This mechanism may amplify the magnetic field to more than 100 micro G and changes the shock structure. As a result, the radio spectrum and the large magnetic field can be explained, apart from the widely discussed Bell's mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Two-Dimensional particle-in-cell simulations of the nonresonant, cosmic-ray driven instability in SNR shocks

In supernova remnants, the nonlinear amplification of magnetic fields upstream of collisionless shocks is essential for the acceleration of cosmic rays to the energy of the "knee" at 10^{15.5}eV. A nonresonant instability driven by the cosmic ray current is thought to be responsible for this effect. We perform two-dimensional, particle-in-cell simulations of this instability. We observe an initial growth of circularly polarized non-propagating magnetic waves as predicted in linear theory. It is demonstrated that in some cases the magnetic energy density in the growing waves, can grow to at least 10 times its initial value. We find no evidence of competing modes, nor of significant modification by thermal effects. At late times we observe saturation of the instability in the simulation, but the mechanism responsible is an artefact of the periodic boundary conditions and has no counterpart in the supernova-shock scenario.

0 citations0 reviewsNo code linkedOpen access