Search

Pulsar wind nebulae (PWNe) are clouds of the magnetized relativistic electron/positron plasma supplied from the central pulsar. However, the number of radio-emitting particles inside a PWN is larger than the expectation from the study of pulsar magnetospheres and then their origin is still unclear. A stochastic acceleration of externally injected particles by a turbulence inside the PWN is proposed by our previous studies. In this paper, the previous stochastic acceleration model of the PWN broadband spectra is improved by taking into account the time evolution of the turbulent energy and then the total energy balance inside a PWN is maintained. The turbulent energy supplied from the central pulsar is wasted by the backreaction from the stochastic particle acceleration and the adiabatic cooling according the PWN expansion. The model is applied to the Crab Nebula and reproduce the current broadband emission spectrum, especially the flat radio spectrum although time evolution of the turbulent energy (diffusion coefficient) is a bit complicated compared with our previous studies, where we assumed an exponential behavior of the diffusion coefficient., Comment: Accepted for publication in Prog. Theor. Exp. Phys., 16 pages, 4 figures and 1 table

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

The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very-high-energy (VHE, > 100GeV) photons which have yet to be detected in coincidence with a gravitational wave signal. The Cherenkov Telescope Array (CTA) is a next-generation VHE observatory which aims to be indispensable in this search, with an unparalleled sensitivity and ability to slew anywhere on the sky within a few tens of seconds. New observing modes and follow-up strategies are being developed for CTA to rapidly cover localization areas of gravitational wave events that are typically larger than the CTA field of view. This work will evaluate and provide estimations on the expected number of of gravitational wave events that will be observable with CTA, considering both on- and off-axis emission. In addition, we will present and discuss the prospects of potential follow-up strategies with CTA., Comment: Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 (arXiv:2309.08219)

Over three thousand pulsars have been discovered, but none have been confirmed to be younger than a few hundred years. Observing a pulsar after a supernova explosion will help us understand the properties of newborn ones, including their capability to produce gamma-ray bursts and fast radio bursts. Here, the possible youngest pulsar wind nebula (PWN) at the center of the SN 1986J remnant is studied. We demonstrate that the 5 GHz flux of 'PWN 1986J', increasing with time, is consistent with a stochastic acceleration model of PWNe developed to explain the flat radio spectrum of the Crab Nebula. We obtain an acceleration time-scale of electrons/positrons and a decay time-scale of the turbulence responsible for the stochastic acceleration as about 10 and 70 years, respectively. Our findings suggest that efficient stochastic acceleration and rising radio/submm light curves are characteristic signatures of the youngest PWNe. Follow-up ${\it ALMA}$ observations of decades-old supernovae within a few tens of Mpc, including SN 1986J, are encouraged to reveal the origin of the flat radio spectrum of PWNe., Comment: Accepted for publication in MNRAS, 8 pages, 2 figures and 1 table

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 $\theta_0=0.015$ rad and large bulk Lorentz factor $\Gamma_0=350$, and `wide jet' with $\theta_0=0.1$ rad and $\Gamma_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)., Comment: 16 pages, 9 figures, 3 tables, JHEAp, in press

Pulsar wind nebulae (PWNe) are one of the most energetic galactic sources with bright emissions from radio waves to very high-energy gamma-rays. We perform wideband X-ray spectroscopy of four energetic PWNe, N157B, PSR J1813-1749, PSR J1400-6325, and G21.5-0.9, with the Suzaku, Chandra, NuSTAR, and Hitomi observatories. A significant spectral break or cutoff feature is found in the hard X-ray band for all the samples, except for N157B. The break energies in the broken power-law fitting are in the range of 4--14 keV, whereas the cutoff energies in the cutoff power-law fitting are at 22 keV or higher. The break or cutoff energy does not show a significant correlation with either the spin-down energy or characteristic age of the hosting pulsars. A possible correlation is found between the photon index change in the broken power-law fitting and the X-ray emitting efficiency of the pulsars, although its significance is not high enough to be conclusive. We discuss what determines the break parameters based on simple models., Comment: 21 pages, 22 figures, PASJ, accepted

Recently, the Large High Altitude Air Shower Observatory (LHAASO) reported the discovery of 12 ultrahigh-energy (UHE; $\mathrm{\varepsilon} \ge 100$ TeV) gamma-ray sources located in the Galactic plane. A few of these UHE gamma-ray emitting regions are in spatial coincidence with pulsar wind nebulae (PWNe). We consider a sample of five sources; two of them are LHAASO sources (LHAASO J1908+0621 and LHAASO J2226+6057) and the remaining three are GeV-TeV gamma-ray emitters. In addition, X-rays, radio observations, or upper limits are also available for these objects. We study multiwavelength radiation from these sources by considering a PWN origin, where the emission is powered by the spin-down luminosity of the associated pulsars. In this leptonic emission model, the electron population is calculated at different times under the radiative (synchrotron and inverse-Compton) and adiabatic cooling. We also include the onset of the reverberation phase for the PWN, by assuming radially symmetric expansion. However, in this work, we find that multiwavelength emission can be interpreted before the onset of this phase. The maximum energy of the electrons based on the spectral fit is found to be above 0.1 PeV and close to 1 PeV. For LHAASO J2226+6057, using its observations in radio to UHE gamma-rays, we find that UHE gamma-rays can be interpreted using electrons with maximum energy of 1 PeV. We estimate the upper limits on the minimum Lorentz factor of the electrons and it also infers the minimum value of the pair-multiplicity of charged pairs., Comment: 13 pages, 11 figures, Revised version, Accepted for publication in MNRAS, Comments welcome

We present a systematic study of particle transport by diffusion in young pulsar wind nebulae (PWNe). We selected nine bright sources that are well resolved with the Chandra X-ray Observatory. We analyzed archival data to obtain their radial profiles of photon index ($\Gamma$) and surface brightness ($\Sigma$) in a consistent way. These profiles were then fit with a pure diffusion model that was tested on Crab, 3C 58 and G21.5$-$0.9 before. In addition to the spectral softening due to the diffusion, we calculated the synchrotron power and built up the theoretical surface brightness profile. For each source, we performed separate fits to the $\Gamma$ and the $\Sigma$ profiles. We found these two profiles of most PWNe are similar, except for Crab and Vela. Both profiles can be well described by our model, suggesting that diffusion dominates the particle transport in most sampled PWNe. The discrepancy of parameters between the $\Gamma$ profile and $\Sigma$ profiles is relatively large for 3C 58 and G54.1+0.3. This difference could be attributed to the elongated shape, reflecting boundary, and the non-uniform magnetic field. Finally, we found no significant correlations between diffusion parameters and physical parameters of PWN and pulsar., Comment: 15 pages, 6 figures, accepted for publication in ApJ

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., Comment: 18 pages, 14 figures, 3 tables, accepted for publication in PASJ

Giant radio pulses (GRPs) are sporadic bursts emitted by some pulsars, lasting a few microseconds. GRPs are hundreds to thousands of times brighter than regular pulses from these sources. The only GRP-associated emission outside radio wavelengths is from the Crab Pulsar, where optical emission is enhanced by a few percent during GRPs. We observed the Crab Pulsar simultaneously at X-ray and radio wavelengths, finding enhancement of the X-ray emission by $3.8\pm0.7\%$ (a 5.4$\sigma$ detection) coinciding with GRPs. This implies that the total emitted energy from GRPs is tens to hundreds of times higher than previously known. We discuss the implications for the pulsar emission mechanism and extragalactic fast radio bursts., Comment: 63 pages, 23 figures, 7 Tables. Published in Science

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

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

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 $\alpha$ = 2.66 and $\beta$ = 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., Comment: 9 pages, 5 figures; accepted for publication in MNRAS

A pulsar wind is a relativistic outflow dominated by Poynting energy at its base. Based on the standard ideal magnetohydrodynamic (MHD) model of pulsar wind nebulae (PWNe) with the ordered magnetic field, the observed slow expansion $v_{\rm PWN} \ll c$ requires the wind to be dominated by kinetic energy at the upstream of its termination shock, which conflicts with the pulsar wind theory ($\sigma$-problem). In this paper, we extend the standard model of PWNe by phenomenologically taking into account conversion of the ordered to turbulent magnetic field and dissipation of the turbulent magnetic field. Disordering of the magnetic structure is inferred from the recent three-dimensional relativistic ideal MHD simulations, while magnetic dissipation is a non-ideal MHD effect requiring a finite resistivity. We apply this model to the Crab Nebula and find that the conversion effect is important for the flow deceleration, while the dissipation effect is not. Even for Poynting-dominated pulsar wind, we obtain the Crab Nebula's $v_{\rm PWN}$ by adopting a finite conversion time-scale of $\sim 0.3$ yr. Magnetic dissipation primarily affects the synchrotron radiation properties. Any values of the pulsar wind magnetization $\sigma_{\rm w}$ are allowed within the present model of the PWN dynamics alone, and even a small termination shock radius of $\ll 0.1$ pc reproduces the observed dynamical features of the Crab Nebula. In order to establish a high-$\sigma_{\rm w}$ model of PWNe, it is important to extend the present model by taking into account the broadband spectrum and its spacial profiles., Comment: Accepted for publication in MNRAS, 13 pages, 4 figures, 1 table

Low-mass population III (PopIII) stars of $\lesssim 0.8 M_{\odot}$ could survive up until the present. Non-detection of low-mass PopIII stars in our Galaxy has already put a stringent constraint on the initial mass function (IMF) of PopIII stars, suggesting that PopIII stars have a top-heavy IMF. On the other hand, some claims that the lack of such stars stems from metal enrichment of their surface by accretion of heavy elements from interstellar medium (ISM). We investigate effects of the stellar wind on the metal accretion onto low-mass PopIII stars because accretion of the local ISM onto the Sun is prevented by the solar wind even for neutrals. The stellar wind and radiation of low-mass PopIII stars are modeled based on knowledge of nearby low-mass stellar systems including our Sun. We find that low-mass PopIII stars traveling across the Galaxy forms the stellar magnetosphere in most of their life. Once the magnetosphere is formed, most of neutral interstellar particles are photoionized before reaching to the stellar surface and are blown away by the wind. Especially, the accretion abundance of iron will be reduced by a factor of $< 10^{-12}$ compared with Bondi-Hoyle-Lyttleton accretion. The metal accretion can enhance iron abundance [Fe/H] only up to $\sim -14$. This demonstrates that low-mass PopIII stars remain pristine and will be found as metal free stars and that further searches for them are valuable to constrain the IMF of PopIII stars., Comment: Accepted for publication in ApJ, 12 pages, 6 figures, 2 tables

The broadband emission of Pulsar Wind Nebulae (PWNe) is well described by non-thermal emissions from accelerated electrons and positrons. However, the standard shock acceleration model of PWNe does not account for the hard spectrum in radio wavelengths. The origin of the radio-emitting particles is also important to determine the pair production efficiency in the pulsar magnetosphere. Here, we propose a possible resolution for the particle energy distribution in PWNe; the radio-emitting particles are not accelerated at the pulsar wind termination shock but are stochastically accelerated by turbulence inside PWNe. We upgrade our past one-zone spectral evolution model including the energy diffusion, i.e., the stochastic acceleration, and apply to the Crab Nebula. A fairly simple form of the energy diffusion coefficient is assumed for this demonstrative study. For a particle injection to the stochastic acceleration process, we consider the continuous injection from the supernova ejecta or the impulsive injection associated with supernova explosion. The observed broadband spectrum and the decay of the radio flux are reproduced by tuning the amount of the particle injected to the stochastic acceleration process. The acceleration time-scale and the duration of the acceleration are required to be a few decades and a few hundred years, respectively. Our results imply that some unveiled mechanisms, such as back reaction to the turbulence, are required to make the energies of stochastically and shock accelerated particles comparable., Comment: Accepted for publication in ApJ, 11 pages, 4 figures

The observed radial profiles of the X-ray emission from Pulsar Wind Nebulae (PWNe) have been claimed to conflict with the standard one-dimensional (1-D) steady model. However, the 1-D model has not been tested to reproduce both the volume-integrated spectrum and the radial profile of the surface brightness, simultaneously. We revisit the 1-D steady model and apply it to PWNe 3C 58 and G21.5-0.9. We find that the parameters of the pulsar wind, the radius of the termination shock $r_{\rm s}$ and magnetization $\sigma$, greatly affect both the photon spectrum and radial profile of the emission. We have shown that the parameters constrained by the entire spectrum lead to a smaller X-ray nebula than observed one. We have also tested the case that reproduces only the observations in X and gamma-rays, ignoring the radio and optical components. In this case, there are parameter sets that reproduce both the spectrum and emission profile, but the advection time to the edge of the nebula becomes much smaller than the age. Our detailed discussion clarifies that the standard 1-D steady model has severe difficulty to reproduce both the volume-integrated spectrum and the surface brightness simultaneously. This implies that the model should be improved by taking into account extra physical processes such as spatial diffusion of particles. Additionally, we calculate the surface brightness profile of the radio, optical and TeV gamma-rays. The future observations in these wavelengths are also important to probe the spatial distributions of the relativistic plasma and the magnetic field of PWNe., Comment: 15 pages, 10 figures, accepted for publication in ApJ

Synchrotron radiation is widely considered as the origin of the pulsed non-thermal emissions from rotation-powered pulsars in optical and X-ray bands. In this paper, we study the synchrotron radiation emitted by the created electron and positron pairs in the pulsar magnetosphere to constrain on the energy conversion efficiency from the Poynting flux to the particle energy flux. We model two pair creation processes, two-photon collision which efficiently works in young $\gamma$-ray pulsars ($\lesssim10^6$ yr), and magnetic pair creation which is the dominant process to supply pairs in old pulsars ($\gtrsim10^6$ yr). Using the analytical model, we derive the maximum synchrotron luminosity as a function of the energy conversion efficiency. From the comparison with observations, we find that the energy conversion efficiency to the accelerated particles should be an order of unity in the magnetosphere, even though we make a number of the optimistic assumptions to enlarge the synchrotron luminosity. In order to explain the luminosity of the non-thermal X-ray/optical emission from pulsars with low spin-down luminosity $L_{\rm sd}\lesssim10^{34}$ erg s$^{-1}$, non-dipole magnetic field components should be dominant at the emission region. For the $\gamma$-ray pulsars with $L_{\rm sd}\lesssim10^{35}$ erg s$^{-1}$, observed $\gamma$-ray to X-ray and optical flux ratios are much higher than the flux ratio between curvature and the synchrotron radiations. We discuss some possibilities such as the coexistence of multiple accelerators in the magnetosphere as suggested from the recent numerical simulation results. The obtained maximum luminosity would be useful to select observational targets in X-ray and optical bands., Comment: 18 pages, 6 figures. Accepted for publication in Astrophysical Journal

We present the results of the simultaneous observation of the giant radio pulses (GRPs) from the Crab pulsar at 0.3, 1.6, 2.2, 6.7, and 8.4 GHz with four telescopes in Japan. We obtain 3194 and 272 GRPs occurring at the main pulse and the interpulse phases, respectively. A few GRPs detected at both 0.3 and 8.4 GHz are the most wide-band samples ever reported. In the frequency range from 0.3 to 2.2 GHz, we find that about 70\% or more of the GRP spectra are consistent with single power laws and the spectral indices of them are distributed from $-4$ to $-1$. We also find that a significant number of GRPs have such a hard spectral index (approximately $-1$) that the fluence at 0.3 GHz is below the detection limit ("dim-hard" GRPs). Stacking light curves of such dim-hard GRPs at 0.3 GHz, we detect consistent enhancement compared to the off-GRP light curve. Our samples show apparent correlations between the fluences and the spectral hardness, which indicates that more energetic GRPs tend to show softer spectra. Our comprehensive studies on the GRP spectra are useful materials to verify the GRP model of fast radio bursts in future observations., Comment: 19 pages, 25 figures, accepted for publication in ApJ

Angular momentum loss by the plasma wind is considered as a universal feature of isolated neutron stars including magnetars. The wind nebulae powered by magnetars allow us to compare the wind properties and the spin-evolution of magnetars with those of rotation-powered pulsars (RPPs). In this paper, we construct a broadband emission model of magnetar wind nebulae (MWNe). The model is similar to past studies of young pulsar wind nebulae (PWNe) around RPPs, but is modified for the application to MWNe that have far less observational information than the young PWNe. We apply the model to the MWN around the youngest ($\sim$ 1kyr) magnetar 1E 1547.0-5408 that has the largest spin-down power $L_{\rm spin}$ among all the magnetars. However, the MWN is faint because of low $L_{\rm spin}$ of 1E 1547.0-5408 compared with the young RPPs. Since most of parameters are not well constrained only by an X-ray flux upper limit of the MWN, we adopt the model parameters from young PWN Kes 75 around PSR J1846-0258 that is a peculiar RPP showing magnetar-like behaviors. The model predicts $\gamma$-ray flux that will be detected in a future TeV $\gamma$-ray observation by {\it CTA}. The MWN spectrum does not allow us to test hypothesis that 1E 1547.0-5408 had milliseconds period at its birth because the particles injected early phase of evolution are suffered from severe adiabatic and synchrotron losses. Further both observational and theoretical studies of the wind nebulae around magnetars are required to constrain the wind and spin-down properties of magnetars., Comment: Accepted for publication in ApJ, 27 pages, 3 figures

Induced Compton scattering (ICS) is an interaction between intense electro-magnetic radiations and plasmas, where ICS transfers the energy from photons to plasmas. Although ICS is important for laser plasma interactions in laboratory experiments and for radio emission from pulsars propagating in pulsar wind plasmas, the detail of photon cooling process has not been understood. The problem is that, when ICS dominates, evolution of photon spectra is described as a nonlinear convection equation, which makes photon spectra to be multi-valued. Here, we propose a new approach to treat evolution of photon spectra affected by ICS. Starting from the higher-order Kompaneets equation, we find a new equation that resolves the unphysical behavior of photon spectra. In addition, we find the steady-state analytic solution, which is linearly stable. We also successfully simulate the evolution of photon spectra without artificial viscosity. We find that photons rapidly lose their energy by ICS with continuously forming solitary structures in frequency-space. The solitary structures have the logarithmically same width characterized by an electron temperature. The energy transfer from photons to plasma is more effective for broader spectrum of photons such as expected in astrophysical situations., Comment: Accepted for publication in PTEP, 15 pages, 3 figures

Electrons/positrons produced in a pulsar magnetosphere emit synchrotron radiation, which is widely believed as the origin of the non-thermal X-ray emission detected from pulsars. Particles are produced by curvature photons emitted from accelerated particles in the magnetosphere. These curvature photons are detected as pulsed $\gamma$-ray emissions from pulsars with age $\lesssim10^6$ yr. Using $\gamma$-ray observations and analytical model, we impose severe constraints on the synchrotron radiation as a mechanism of the non-thermal X-ray emission. In most middle-aged pulsars ($\sim10^5-10^6$ yr) which photon-photon pair production is less efficient in their magnetosphere, we find that the synchrotron radiation model is difficult to explain the observed non-thermal X-ray emission., Comment: 2014 Fermi Symposium proceedings - eConf C14102.1

We study the synchrotron radiation as the observed non-thermal X-ray emission from old pulsars ($\gtrsim1-10$Myr) to investigate the particle acceleration in their magnetospheres. We assume that the power-law component of the observed X-ray spectra is caused by the synchrotron radiation from electrons and positrons in the magnetosphere. We consider two pair production mechanisms of X-ray emitting particles, the magnetic and the photon-photon pair productions. High-energy photons, which ignite the pair production, are emitted via the curvature radiation of the accelerated particles. We use the analytical description for the radiative transfer and estimate the luminosity of the synchrotron radiation. We find that for pulsars with the spin-down luminosity $L_{\rm sd}\lesssim10^{33}$ erg s$^{-1}$, the locations of the particle acceleration and the non-thermal X-ray emission are within $\lesssim10^7$cm from the centre of the neutron star, where the magnetic pair production occurs. For pulsars with the spin-down luminosity $L_{\rm sd}\lesssim10^{31}$ erg s$^{-1}$ such as J0108-1431, the synchrotron radiation is difficult to explain the observed non-thermal component even if we consider the existence of the strong and small-scale surface magnetic field structures., Comment: 25 pages, 7 figures, 2 tables, accepted for publication in MNRAS

Pulsar winds have longstanding problems in energy conversion and pair cascade processes which determine the magnetization $\sigma$, the pair multiplicity $\kappa$ and the bulk Lorentz factor $\gamma$ 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 $\theta_{\rm pl}$ and the emission region size $r_{\rm e}$ which determines an opening angle of the radio beam. We obtain the lower limit $\gamma\gtrsim10^{1.7}r^{1/2}_{\rm e,3}\theta^{-1}_{\rm pl}(1+\sigma)^{-1/4}$ ($r_{\rm e}=10^3r_{\rm e,3}$ cm) at the light cylinder $r_{\rm LC}$ for an inclined wind $\theta_{\rm pl}>10^{-2.7}$. For an aligned wind $\theta_{\rm pl}<10^{-2.7}$, we require $\gamma>10^{2.7}$ at $r_{\rm LC}$ and an additional constraint $\gamma>10^{3.4}r^{1/5}_{\rm e,3}(1+\sigma)^{-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' effect prevents scattering. Considering the lower limit $\kappa\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\gamma\lesssim10^{3.9}$ and $10^{6.6}\lesssim\kappa\lesssim10^{8.8}$ which are independent of $r$ when $\theta_{\rm pl}\sim1$ and $1+\sigma\sim1$ at $r_{\rm LC}$., Comment: 26 pages, 11 figures, 4 tables, accepted for publication in PTEP

Among dozens young pulsar wind nebulae, some have been detected in TeV \gamma-rays (TeV PWNe), while others have not (non-TeV PWNe). The TeV emission detectability is not correlated either with the spin-down power or with the characteristic age of their central pulsars, and it is an open problem what determines the detectability. To study this problem, we investigate spectral evolution of five young non-TeV PWNe, 3C58, G310.6-1.6, G292.0+1.8, G11.2-0.3 and SNR B0540-69.3. We use a spectral evolution model which has been developed to be applied to young TeV PWNe in our previous works. TeV \gamma-ray flux upper limits of non-TeV PWNe give upper or lower limits on parameters, such as the age of the PWN and the fraction of the spin-down power going to the magnetic energy injection (the fraction parameter). Combined with other independent observational and theoretical studies, we can guess a plausible value of the parameters for each object. For 3C58, we prefer the parameters with an age of 2.5 kyr old and the fraction parameter of 3.0x10^{-3}, although the spectral modeling alone does not rule out a shorter age and a higher fraction parameter. The fraction parameter of 3.0x10^{-3} is also consistent for other non-TeV PWNe and then the value is regarded as common to young PWNe including TeV PWNe. Moreover, we find that the intrinsic properties of the central pulsars are similar, 10^{48-50}erg for the initial rotational energy and 10^{42-44}erg for the magnetic energy (2x10^{12} - 3x10^{13}G for the dipole magnetic field strength at their surfaces). The TeV detectability is correlated with the total injected energy and the energy density of the interstellar radiation field around PWNe. Except for G292.0+1.8, a broken power-law injection of the particles well reproduces the broadband emission from non-TeV PWNe., Comment: Accepted for publication in MNRAS

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., Comment: Accepted for publication in ApJ; 56pages, 15figures

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., Comment: Accepted for publication in ApJ; 37 pages, 7 figures.

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., Comment: Accepted for publication in ApJ Letters