Researcher profile

Fumio Takahara

Fumio Takahara contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
27works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2016arXiv

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

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

preprint2016arXiv

Evidence for a significant mixture of electron/positron pairs in FRII jets constrained by cocoon dynamics

We examine the plasma composition of relativistic jets in four FRII radio galaxies by analyzing the total cocoon pressure in terms of partial pressures of thermal and non-thermal electrons/positrons and protons. The total cocoon pressure is determined by cocoon dynamics via comparison of theoretical model with the observed cocoon shape. By inserting the observed number density of non-thermal electrons/positrons and the upper limit of thermal electron/positron number density into the equation of state, the number density of protons is constrained. We apply this method to four FRII radio galaxies (Cygnus A, 3C 219, 3C 223 and 3C 284), for which the total cocoon pressures have been already evaluated. We find that the positron-free plasma comprising of protons and electrons is ruled out, when we consider plausible particle distribution functions. In other words, the mixture of positrons is required for all four FRII radio galaxies; the number density ratio of electrons/positrons to protons is larger than two. Thus, we find that the plasma composition is independent of the jet power and the size of cocoons. We also investigate the additional contribution of thermal electrons/positrons and protons on the cocoon dynamics. When thermal electrons/positrons are absent, the cocoon is supported by the electron/ proton plasma pressure, while both electron/positron pressure supported and electron/proton plasma pressure supported cocoons are allowed if the number density of thermal electrons/positrons is about 10 times larger than that of non-thermal ones.

preprint2015arXiv

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

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

preprint2014arXiv

Electromotive Force in the Blandford-Znajek Process

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

preprint2014arXiv

General properties of the radiation spectra from relativistic electrons moving in a Langmuir turbulence

We examine the radiation spectra from relativistic electrons moving in a Langmuir turbulence expected to exist in high energy astrophysical objects by using numerical method. The spectral shape is characterized by the spatial scale λ, field strength σ, and frequency of the Langmuir waves, and in term of frequency they are represented by ω_0 = 2πc/λ, ω_st = eσ/mc, and ω_p, respectively. We normalize ω_st and ω_p by ω_0 as \a \equiv ω_st/ω_0 and \b \equivω_p/ω_0, and examine the spectral shape in the a-b plane. An earlier study based on Diffusive Radiation in Langmuir turbulence (DRL) theory by Fleishman and Toptygin showed that the typical frequency is γ^2ω_p and that the low frequency spectrum behaves as F_ω pronto ω^1 for b > 1 irrespective of a. Here, we adopt the first principle numerical approach to obtain the radiation spectra in more detail. We generate Langmuir turbulence by superposing Fourier modes,inject monoenergetic electrons, solve the equation of motion, and calculate the radiation spectra using Lienard-Wiechert potential. We find different features from the DRL theory for a>b>1. The peak frequency turns out to be γ^2ω_st which is higher than γ^2ω_p predicted in the DRL theory, and the spectral index of low frequency region is not 1 but 1/3. It is because the typical deflection angle of electrons is larger than the angle of the beaming cone sim 1/γ. We call the radiation for this case "Wiggler Radiation in Langmuir turbulence" (WRL).

preprint2013arXiv

Constraint on Pulsar Wind Properties from Induced Compton Scattering off Radio Pulses

Pulsar winds have longstanding problems in energy conversion and pair cascade processes which determine the magnetization $σ$, the pair multiplicity $κ$ and the bulk Lorentz factor $γ$ of the wind. We study induced Compton scattering by a relativistically moving cold plasma to constrain wind properties by imposing that radio pulses from the pulsar itself are not scattered by the wind as was first studied by Wilson & Rees. We find that relativistic effects cause a significant increase or decrease of the scattering coefficient depending on scattering geometry. Applying to the Crab, we consider uncertainties of an inclination angle of the wind velocity with respect to the radio beam $θ_{\rm pl}$ and the emission region size $r_{\rm e}$ which determines an opening angle of the radio beam. We obtain the lower limit $γ\gtrsim10^{1.7}r^{1/2}_{\rm e,3}θ^{-1}_{\rm pl}(1+σ)^{-1/4}$ ($r_{\rm e}=10^3r_{\rm e,3}$ cm) at the light cylinder $r_{\rm LC}$ for an inclined wind $θ_{\rm pl}>10^{-2.7}$. For an aligned wind $θ_{\rm pl}<10^{-2.7}$, we require $γ>10^{2.7}$ at $r_{\rm LC}$ and an additional constraint $γ>10^{3.4}r^{1/5}_{\rm e,3}(1+σ)^{-1/10}$ at the characteristic scattering radius $r_{\rm c}=10^{9.6}r^{2/5}_{\rm e,3}$ cm within which the `lack of time&#39; effect prevents scattering. Considering the lower limit $κ\gtrsim10^{6.6}$ suggested by recent studies of the Crab Nebula, for $r_{\rm e}=10^3$ cm, we obtain the most optimistic constraint $10^{1.7}\lesssimγ\lesssim10^{3.9}$ and $10^{6.6}\lesssimκ\lesssim10^{8.8}$ which are independent of $r$ when $θ_{\rm pl}\sim1$ and $1+σ\sim1$ at $r_{\rm LC}$.

preprint2013arXiv

Efficient Acceleration of Relativistic Magnetohydrodynamic Jets

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

preprint2013arXiv

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

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

preprint2012arXiv

A Leptonic Model of Steady High-Energy Gamma-Ray Emission from Sgr A$^*$

Recent observations of Sgr A$^*$ by Fermi and HESS have detected steady gamma-ray emission in the GeV and TeV bands. We present a new model to explain the GeV gamma-ray emission by inverse Compton scattering by nonthermal electrons supplied by the NIR/X-ray flares of Sgr A$^*$. The escaping electrons from the flare regions accumulate in a region with a size of $\sim 10^{18}$ cm and magnetic fields of $\lesssim 10^{-4}$ G. Those electrons produce gamma-rays by inverse Compton scattering off soft photons emitted by stars and dust around the central black hole. By fitting the GeV spectrum, we find constraints on the magnetic field and the energy density of optical-UV radiation in the central 1 pc region around the supermassive black hole. While the GeV spectrum is well fitted by our model, the TeV $γ$-rays, whose spectral index is different from that of the GeV emission, may be from different sources such as pulsar wind nebulae.

preprint2012arXiv

Baryon Loading of AGN Jets Mediated by Neutrons

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

preprint2012arXiv

Calorimetry of Active Galactic Nucleus jets: testing plasma composition in Cygnus A

We examine plasma composition of jets in active galactic nuclei through the comparison of the total pressure ($P$) with partial pressures of electrons and protons in a cocoon. The total pressure is estimated from the analysis of an expanding cocoon dynamics. We determine the average kinetic energy per particle for several representative cases of particle energy distribution such as one- and two-temperature thermal plasmas and non-thermal electrons by evaluating the dissipation of total kinetic energy of the jet into the internal energy of cocoon plasma. The number density of the total electrons/positrons ($n_{\pm}$) in the cocoon is constrained by using the particle supply from hot spots and the absence of thermal bremsstrahlung emission from radio lobes. By inserting $P$, $n_{\pm}$ and the particle energy of each population into the equation of state, the number density ($n_{p}$) and pressure ($P_{p}$) of protons in the cocoon can be constrained. Applying this method to Cygnus A, we find that (i) electron/positron ($e^{\pm}$) pairs always dominate in terms of number density, but that (ii) either an &#34;$e^{\pm}$-supported cocoon (i.e., $P_{\pm} >P_{p}$)&#34; or &#34;proton-supported one (i.e, $P_{\pm} <P_{p}$)&#34; is possible.

preprint2012arXiv

Inverse Compton Scattering Model for X-ray Emission of the Gamma-ray Binary LS 5039

We propose a model for the gamma-ray binary LS 5039 in which the X-ray emission is due to the inverse Compton (IC) process instead of the synchrotron radiation. Although the synchrotron model has been discussed in previous studies, it requires a strong magnetic field which leads to a severe suppression of the TeV gamma-ray flux in conflict with H.E.S.S. observations. In this paper, we calculate the IC emission by low energy electrons (γ_e \lesssim 10^3) in the Thomson regime. We find that IC emission of the low energy electrons can explain the X-ray flux and spectrum observed with Suzaku if the minimum Lorentz factor of injected electrons γ_min is around 10^3. In addition, we show that the Suzaku light curve is well reproduced if γ_min varies in proportion to the Fermi flux when the distribution function of injected electrons at higher energies is fixed. We conclude that the emission from LS 5039 is well explained by the model with the IC emission from electrons whose injection properties are dependent on the orbital phase. Since the X-ray flux is primarily determined by the total number of cooling electrons, this conclusion is rather robust, although some mismatches between the model and observations at the GeV band remain in the present formulation.

preprint2011arXiv

A Novel Emission Spectrum From A Relativistic Electron Moving In A Random Magnetic Field

We calculate numerically the radiation spectrum from relativistic electrons moving in small scale turbulent magnetic fields expected in high energy astrophysical sources. Such radiation spectrum is characterized by the strength parameter $a = λ_{B}} e|B|/mc^2$, where $λ_{B}}$ is the length scale of the turbulent field. When $a$ is much larger than the Lorentz factor of a radiating electron $γ$, synchrotron radiation is realized, while $a\ll 1$ corresponds to the so-called jitter radiation regime. Because for $1<a<γ$ we cannot use either approximations, we should have recourse to the Lienard-Wiechert potential to evaluate the radiation spectrum, which is performed in this paper. We generate random magnetic fields assuming Kolmogorov turbulence, inject monoenergetic electrons, solve the equation of motion, and calculate the radiation spectrum. We perform numerical calculations for several values of $a$ with $γ= 10$. We obtain various types of spectra ranging between jitter radiation and synchrotron radiation. For $a\sim 7$, the spectrum turns out to take a novel shape which has not been noticed up to now. It is like a synchrotron spectrum in the middle energy region, but in the low frequency region it is a broken power law and in the high frequency region an extra power law component appears beyond the synchrotron cutoff. We give a physical explanation of these features.

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.

preprint2011arXiv

Study of Four Young TeV Pulsar Wind Nebulae with a Spectral Evolution Model

We study four young Pulsar Wind Nebulae (PWNe) detected in TeV gamma-rays, G21.5-0.9, G54.1+0.3, Kes 75, and G0.9+0.1, using the spectral evolution model developed and applied to the Crab Nebula in our previous work. We model the evolution of magnetic field and particle distribution function inside a uniformly expanding PWN considering a time-dependent injection from the pulsar and radiative and adiabatic losses. Considering uncertainties in the interstellar radiation field (ISRF) and their distance, we study two cases for each PWN. Because TeV PWNe have a large TeV gamma-rays to X-rays flux ratio, the magnetic energy of the PWNe accounts for only a small fraction of the total energy injected (typically a few x 10^{-3}). The gamma-ray emission is dominated by inverse Compton scattering off the infrared photons of the ISRF. A broken power-law distribution function for the injected particles reproduces the observed spectrum well, except for G0.9+0.1. For G0.9+0.1, we do not need a low energy counterpart because adiabatic losses alone are enough to reproduce the radio observations. High energy power-law indices at injection are similar (2.5 -- 2.6), while low energy power-law indices range from 1.0 to 1.6. The lower limit of the particle injection rate indicates that the pair multiplicity is larger than 10^4. The corresponding upper limit of the bulk Lorentz factor of the pulsar winds is close to the break energy of the broken power-law injection, except for Kes 75.

preprint2010arXiv

A Model of the Spectral Evolution of Pulsar Wind Nebulae

We study the spectral evolution of PWNe taking into account the energy injected when they are young. We model the evolution of the magnetic field inside a uniformly expanding PWN. Considering time dependent injection from the pulsar and coolings by radiative and adiabatic losses, we solve the evolution of the particle distribution function. The model is calibrated by fitting the calculated spectrum to the observations of the Crab Nebula at an age of a thousand years. The spectral evolution of the Crab Nebula in our model shows that the flux ratio of TeV gamma-rays to X-rays increases with time, which implies that old PWNe are faint in X-rays, but not in TeV gamma-rays. The increase of this ratio is because the magnetic field decreases with time and is not because the X-ray emitting particles are cooled more rapidly than the TeV gamma-ray emitting particles. Our spectral evolution model matches the observed rate of the radio flux decrease of the Crab Nebula. This result implies that our magnetic field evolution model is close to the reality. Finally, from the viewpoint of the spectral evolution, only a small fraction of the injected energy from the Crab Pulsar needs to go to the magnetic field, which is consistent with previous studies.

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.

preprint2010arXiv

Modulation Mechanism of TeV, GeV, and X-ray Emission in LS5039

The emission mechanism of the gamma-ray binary LS5039 in energy bands of TeV, GeV, and X-ray is investigated. Observed light curves in LS5039 show that TeV and GeV fluxes anticorrelate and TeV and X-ray fluxes correlate. However, such correlated variations have not been explained yet reasonably at this stage. Assuming that relativistic electrons are injected constantly at the location of the compact object as a point source, and that they lose energy only by the inverse Compton (IC) process, we calculate gamma-ray spectra and light curves by the Monte Carlo method, including the full electromagnetic cascade process. Moreover, we calculated X-ray spectra and light curves by using the resultant electron distribution. As a result, we are able to reproduce qualitatively spectra and light curves observed by HESS, Fermi, and Suzaku for the inclination angle i = 30 dig and the index of injected electron distribution p = 2.5. We conclude that TeV-GeV anticorrelation is due to anisotropic IC scattering and anisotropic gamma-gamma absorption, and that TeV-X correlation is due to the dependence of IC cooling time on orbital phases. In addition, the constraint on the inclination angle implies that the compact object in LS5039 is a black hole.

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.

preprint2010arXiv

Synchrotron Blob Model of Infrared and X-ray Flares from Sagittarius A$^*$

Sagittarius A$^*$ in the Galactic center harbors a supermassive black hole and exhibits various active phenomena. Besides quiescent emission in radio and submillimeter radiation, flares in the near infrared (NIR) and X-ray bands are observed to occur frequently. We study a time-dependent model of the flares, assuming that the emission is from a blob ejected from the central object. Electrons obeying a power law with the exponential cutoff are assumed to be injected in the blob for a limited time interval. The flare data of 2007 April 4 were used to determine the values of model parameters. The spectral energy distribution of flare emission is explained by nonthermal synchrotron radiation in the NIR and X-ray bands. The model light curves suggest that electron acceleration is still underway during the rising phase of the flares. GeV gamma-rays are also emitted by synchrotron self-Compton scattering, although its luminosity is not strictly constrained by the current model. If the GeV emission is faint, the plasma blob is dominated by the magnetic energy density over the electron kinetic energy density. Observations in the GeV band will clarify the origin of the blob.

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.

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&#39;s mechanism.

preprint2009arXiv

Suzaku Observation of the Diffuse X-Ray Emission from the Open Cluster Westerlund 2: a Hypernova Remnant?

We present the analysis of Suzaku observations of the young open cluster Westerlund 2, which is filled with diffuse X-ray emission. We found that the emission consists of three thermal components or two thermal and one non-thermal components. The upper limit of the energy flux of the non-thermal component is smaller than that in the TeV band observed with H.E.S.S. This may indicate that active particle acceleration has stopped in this cluster, and that the accelerated electrons have already cooled. The gamma-ray emission observed with H.E.S.S. is likely to come from high-energy protons, which hardly cool in contrast with electrons. Metal abundances of the diffuse X-ray gas may indicate the explosion of a massive star in the past.

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 &#34;knee&#34; 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.

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.

preprint2003arXiv

Photon Emission in a Cascade from Relativistic Protons Initiated by Residual Thermal Photons in Gamma-Ray Bursts

Gamma-ray bursts are generally considered to be the result of internal shocks generated in an inhomogeneous relativistic outflow that arises from a fireball. In such shocks, the Fermi acceleration of protons is naturally expected to be at least as efficient as that of electrons. We investigate the consequences of proton acceleration in the residual thermal photon field of a fireball, especially those on the emission spectra of photons. In contrast to most other studies, we do not invoke a direct electron acceleration in shock waves. We show that the residual radiation field of the fireball can ignite the photopion production and subsequent cascade, and that the photons emitted in this process further enhance the photon-initiated cascade. We find that Fermi accelerated protons with $\gtrsim 10^{13}$ eV efficiently bring their energy into pion production and subsequent photon emission. The particles cascading from the pions emit photons over a wide energy range. The photons emitted from electron--positron pairs distribute continuously from the GeV range down to the X-ray range, while muons can emit gamma-rays by synchrotron radiation with a break at around 1--10 MeV. We also discuss several radiation processes which may possibly produce a break feature in the MeV range, as observed in addition to muon synchrotron radiation.

preprint1999arXiv

Global Relationships among Physical Properties of Galaxy Cluster Cores

Using X-ray data, we investigate interrelations between gas density rho_0, virial density rho_vir, core radius R, temperature T, entropy S_gas, and metal abundance Z in the core region of clusters of galaxies. First, we confirm that fundamental relations among rho_0, R, and T found by Fujita and Takahara are reproduced by another data catalogue. Second, we find that, when clusters have two components in their surface brightness distribution, the inner components also satisfy the same fundamental relations on the assumption that the average temperature of the inner component is the same as that of the outer component. These results strengthen our interpretation that clusters form a two parameter family in terms of mass and rho_vir; larger rho_vir corresponds to earlier formation epoch. We argue that the inner components represent distinct dark matter components which collapsed ahead of the outer components. Third, we also find a tight relation between S_gas and rho_vir both for the outer and inner components; S_gas is smaller for larger rho_vir but is larger than that produced through gravitational collapse alone for larger rho_vir. Although radiative cooling affects the thermal evolution, the tight relationship discovered suggests the existence of stable heating sources or stable energy transmission mechanisms. Finally, we find that the iron abundance at the centers of clusters is correlated with rho_vir for the inner components. This implies that iron produced by Type Ia supernovae has more accumulated for clusters formed earlier. We briefly discuss the implications of these findings.