Neutrinos generated during solar flares remain elusive. However, after $50$ years of discussion and search, the potential knowledge unleashed by their discovery keeps the search crucial. Neutrinos associated with solar flares provide information on otherwise poorly known particle acceleration mechanisms during solar flare. For neutrino detectors, the separation between atmospheric neutrinos and solar flare neutrinos is technically encumbered by an energy band overlap. To improve differentiation from background neutrinos, we developed a method to determine the temporal search window for neutrino production during solar flares. Our method is based on data recorded by solar satellites, such as Geostationary Operational Environmental Satellite (GOES), Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI), and GEOTAIL. In this study, we selected 23 solar flares above the X5.0 class that occurred between 1996 and 2018. We analyzed the light curves of soft X-rays, hard X-rays, $\gamma$-rays, line $\gamma$-rays from neutron capture as well as the derivative of soft X-rays. The average search windows are determined as follows: $4,178$ s for soft X-ray, $700$ s for derivative of soft X-ray, $944$ s for hard X-ray ($100$-$800$ keV), $1,586$ s for line $\gamma$-ray from neutron captures, and $776$ s for hard X-ray (above $50$ keV). This method allows neutrino detectors to improve their sensitivity to solar flare neutrinos.

In this paper, we present the analysis and results of a direct measurement of the cosmic-ray proton spectrum with the CALET instrument onboard the International Space Station, including the detailed assessment of systematic uncertainties. The observation period used in this analysis is from October 13, 2015 to August 31, 2018 (1054 days). We have achieved the very wide energy range necessary to carry out measurements of the spectrum from 50 GeV to 10 TeV covering, for the first time in space, with a single instrument the whole energy interval previously investigated in most cases in separate subranges by magnetic spectrometers (BESS-TeV, PAMELA, and AMS-02) and calorimetric instruments (ATIC, CREAM, and NUCLEON). The observed spectrum is consistent with AMS-02 but extends to nearly an order of magnitude higher energy, showing a very smooth transition of the power-law spectral index from -2.81 +- 0.03 (50--500 GeV) neglecting solar modulation effects (or -2.87 +- 0.06 including solar modulation effects in the lower energy region) to -2.56 +- 0.04 (1--10 TeV), thereby confirming the existence of spectral hardening and providing evidence of a deviation from a single power law by more than 3 sigma.

The CALorimetric Electron Telescope (CALET) is a high-energy astroparticle physics space experiment installed on the International Space Station (ISS), developed and operated by Japan in collaboration with Italy and the United States. The CALET mission goals include the investigation of possible nearby sources of high-energy electrons, of the details of galactic particle acceleration and propagation, and of potential signatures of dark matter. CALET measures the cosmic-ray electron + positron flux up to 20 TeV, gamma-rays up to 10 TeV, and nuclei with Z=1 to 40 up to 1,000 TeV for the more abundant elements during a long-term observation aboard the ISS. Starting science operation in mid-October 2015, CALET performed continuous observation without major interruption with close to 20 million triggered events over 10 GeV per month. Based on the data taken during the first two-years, we present an overview of CALET observations: uses w/o major interruption 1) Electron + positron energy spectrum, 2) Nuclei analysis, 3) Gamma-ray observation including a characterization of on-orbit performance. Results of the electromagnetic counterpart search for LIGO/Virgo gravitational wave events are discussed as well.

There is growing evidence that a clear distinction between magnetars and radio pulsars may not exist, implying the population of neutron stars that exhibit both radio pulsations and bursting activities could be potentially large. In this situation, new insights into the burst mechanism could be gained by combining the temporal behavior of radio pulsations. We present a general model for radio suppression by relativistic $e^{\pm}$ plasma outflows at the onset of magnetar flares. A sudden ejection of magnetic energy into the magnetosphere would generate a fireball plasma, which is promptly driven to expand at relativistic speed. This would make the plasma cutoff frequency significantly higher than the rest frame radio frequency, resulting in the suppression of radio waves. We analytically show that any GHz radio emission arising from the magnetosphere is suppressed for $\sim100\ {\rm s}$, depending on the total fireball energy. On the other hand, thermal radiation is expected from the hot spot(s) on the stellar surface created by an inflow of dense plasma, which could be the origin of short bursts. Since our hypothesis predicts radio suppression in coincidence with short bursts, this could be an indirect method to constrain the occurrence rate of short bursts at the faint end that remain undetected by X-ray detectors. Furthermore, ultra-fast gamma-ray flashes from the fireball photosphere is also expected as a smoking gun, although the onboard detection is challenging due to its extremely short duration $\sim\mu$s. Finally, our model is applied to the radio pulsar with magnetar-like activities, PSR J1119-6127 in light of recent observations. Implications for fast radio bursts and the possibility of plasma lensing are also discussed.

We present results on searches for gamma-ray counterparts of the LIGO/Virgo gravitational-wave events using CALorimetric Electron Telescope ({\sl CALET}) observations. The main instrument of {\sl CALET}, CALorimeter (CAL), observes gamma-rays from $\sim1$ GeV up to 10 TeV with a field of view of nearly 2 sr. In addition, the {\sl CALET} gamma-ray burst monitor (CGBM) views $\sim$3 sr and $\sim2\pi$ sr of the sky in the 7 keV -- 1 MeV and the 40 keV -- 20 MeV bands, respectively, by using two different crystal scintillators. The {\sl CALET} observations on the International Space Station started in October 2015, and here we report analyses of events associated with the following gravitational wave events: GW151226, GW170104, GW170608, GW170814 and GW170817. Although only upper limits on gamma-ray emission are obtained, they correspond to a luminosity of $10^{49}\sim10^{53}$ erg s$^{-1}$ in the GeV energy band depending on the distance and the assumed time duration of each event, which is approximately the order of luminosity of typical short gamma-ray bursts. This implies there will be a favorable opportunity to detect high-energy gamma-ray emission in further observations if additional gravitational wave events with favorable geometry will occur within our field-of-view. We also show the sensitivity of {\sl CALET} for gamma-ray transient events which is the order of $10^{-7}$~erg\,cm$^{-2}$\,s$^{-1}$ for an observation of 100~s duration.

Extended results on the cosmic-ray electron + positron spectrum from 11 GeV to 4.8 TeV are presented based on observations with the Calorimetric Electron Telescope (CALET) on the International Space Station utilizing the data up to November 2017. The analysis uses the full detector acceptance at high energies, approximately doubling the statistics compared to the previous result. CALET is an all-calorimetric instrument with a total thickness of 30 $X_0$ at normal incidence and fine imaging capability, designed to achieve large proton rejection and excellent energy resolution well into the TeV energy region. The observed energy spectrum in the region below 1 TeV shows good agreement with Alpha Magnetic Spectrometer (AMS-02) data. In the energy region below $\sim$300 GeV, CALET's spectral index is found to be consistent with the AMS-02, Fermi Large Area Telescope (Fermi-LAT) and Dark Matter Particle Explorer (DAMPE), while from 300 to 600 GeV the spectrum is significantly softer than the spectra from the latter two experiments. The absolute flux of CALET is consistent with other experiments at around a few tens of GeV. However, it is lower than those of DAMPE and Fermi-LAT with the difference increasing up to several hundred GeV. The observed energy spectrum above $\sim$1 TeV suggests a flux suppression consistent within the errors with the results of DAMPE, while CALET does not observe any significant evidence for a narrow spectral feature in the energy region around 1.4 TeV. Our measured all-electron flux, including statistical errors and a detailed breakdown of the systematic errors, is tabulated in the Supplemental Material in order to allow more refined spectral analyses based on our data.

The CALorimetric Electron Telescope (CALET), launched for installation on the International Space Station (ISS) in August, 2015, has been accumulating scientific data since October, 2015. CALET is intended to perform long-duration observations of high-energy cosmic rays onboard the ISS. CALET directly measures the cosmic-ray electron spectrum in the energy range of 1 GeV to 20 TeV with a 2% energy resolution above 30 GeV. In addition, the instrument can measure the spectrum of gamma rays well into the TeV range, and the spectra of protons and nuclei up to a PeV. In order to operate the CALET onboard ISS, JAXA Ground Support Equipment (JAXA-GSE) and the Waseda CALET Operations Center (WCOC) have been established. Scientific operations using CALET are planned at WCOC, taking into account orbital variations of geomagnetic rigidity cutoff. Scheduled command sequences are used to control the CALET observation modes on orbit. Calibration data acquisition by, for example, recording pedestal and penetrating particle events, a low-energy electron trigger mode operating at high geomagnetic latitude, a low-energy gamma-ray trigger mode operating at low geomagnetic latitude, and an ultra heavy trigger mode, are scheduled around the ISS orbit while maintaining maximum exposure to high-energy electrons and other high-energy shower events by always having the high-energy trigger mode active. The WCOC also prepares and distributes CALET flight data to collaborators in Italy and the United States. As of August 31, 2017, the total observation time is 689 days with a live time fraction of the total time of approximately 84%. Nearly 450 million events are collected with a high-energy (E>10 GeV) trigger. By combining all operation modes with the excellent-quality on-orbit data collected thus far, it is expected that a five-year observation period will provide a wealth of new and interesting results.

In August 2015, the CALorimetric Electron Telescope (CALET), designed for long exposure observations of high energy cosmic rays, docked with the International Space Station (ISS) and shortly thereafter began tocollect data. CALET will measure the cosmic ray electron spectrum over the energy range of 1 GeV to 20 TeV with a very high resolution of 2% above 100 GeV, based on a dedicated instrument incorporating an exceptionally thick 30 radiation-length calorimeter with both total absorption and imaging (TASC and IMC) units. Each TASC readout channel must be carefully calibrated over the extremely wide dynamic range of CALET that spans six orders of magnitude in order to obtain a degree of calibration accuracy matching the resolution of energy measurements. These calibrations consist of calculating the conversion factors between ADC units and energy deposits, ensuring linearity over each gain range, and providing a seamless transition between neighboring gain ranges. This paper describes these calibration methods in detail, along with the resulting data and associated accuracies. The results presented in this paper show that a sufficient accuracy was achieved for the calibrations of each channel in order to obtain a suitable resolution over the entire dynamic range of the electron spectrum measurement.

First results of a cosmic-ray electron + positron spectrum, from 10 GeV to 3 TeV, is presented based upon observations with the CALET instrument on the ISS starting in October, 2015. Nearly a half million electron + positron events are included in the analysis. CALET is an all-calorimetric instrument with total vertical thickness of 30 $X_0$ and a fine imaging capability designed to achieve a large proton rejection and excellent energy resolution well into the TeV energy region. The observed energy spectrum over 30 GeV can be fit with a single power law with a spectral index of -3.152 $\pm$ 0.016 (stat.+ syst.). Possible structure observed above 100 GeV requires further investigation with increased statistics and refined data analysis.

To search for giant X-ray pulses correlated with the giant radio pulses (GRPs) from the Crab pulsar, we performed a simultaneous observation of the Crab pulsar with the X-ray satellite Hitomi in the 2 -- 300 keV band and the Kashima NICT radio observatory in the 1.4 -- 1.7 GHz band with a net exposure of about 2 ks on 25 March 2016, just before the loss of the Hitomi mission.The timing performance of the Hitomi instruments was confirmed to meet the timing requirement and about 1,000 and 100 GRPs were simultaneously observed at the main and inter-pulse phases, respectively, and we found no apparent correlation between the giant radio pulses and the X-ray emission in either the main or inter-pulse phases.All variations are within the 2 sigma fluctuations of the X-ray fluxes at the pulse peaks, and the 3 sigma upper limits of variations of main- or inter- pulse GRPs are 22\% or 80\% of the peak flux in a 0.20 phase width, respectively, in the 2 -- 300 keV band.The values become 25\% or 110\% for main or inter-pulse GRPs, respectively, when the phase width is restricted into the 0.03 phase.Among the upper limits from the Hitomi satellite, those in the 4.5-10 keV and the 70-300 keV are obtained for the first time, and those in other bands are consistent with previous reports.Numerically, the upper limits of main- and inter-pulse GRPs in the 0.20 phase width are about (2.4 and 9.3) $\times 10^{-11}$ erg cm$^{-2}$, respectively. No significant variability in pulse profiles implies that the GRPs originated from a local place within the magnetosphere and the number of photon-emitting particles temporally increases.However, the results do not statistically rule out variations correlated with the GRPs, because the possible X-ray enhancement may appear due to a $>0.02$\% brightening of the pulse-peak flux under such conditions.

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.

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.

We present upper limits in the hard X-ray and gamma-ray bands at the time of the LIGO gravitational-wave event GW 151226 derived from the CALorimetric Electron Telescope (CALET) observation. The main instrument of CALET, CALorimeter (CAL), observes gamma-rays from ~1 GeV up to 10 TeV with a field of view of ~2 sr. The CALET gamma-ray burst monitor (CGBM) views ~3 sr and ~2pi sr of the sky in the 7 keV - 1 MeV and the 40 keV - 20 MeV bands, respectively, by using two different scintillator-based instruments. The CGBM covered 32.5% and 49.1% of the GW 151226 sky localization probability in the 7 keV - 1 MeV and 40 keV - 20 MeV bands respectively. We place a 90% upper limit of 2 x 10^{-7} erg cm-2 s-1 in the 1 - 100 GeV band where CAL reaches 15% of the integrated LIGO probability (~1.1 sr). The CGBM 7 sigma upper limits are 1.0 x 10^{-6} erg cm-2 s-1 (7-500 keV) and 1.8 x 10^{-6} erg cm-2 s-1 (50-1000 keV) for one second exposure. Those upper limits correspond to the luminosity of 3-5 x 10^{49} erg s-1 which is significantly lower than typical short GRBs.

Recent $\gamma$-ray observations suggest that the particle acceleration occurs at the outer region of the pulsar magnetosphere. The magnetic field lines in the outer acceleration region (OAR) are connected to the neutron star surface (NSS). If copious electron--positron pairs are produced near the NSS, such pairs flow into the OAR and screen the electric field there. To activate the OAR, the electromagnetic cascade due to the electric field near the NSS should be suppressed. However, since a return current is expected along the field lines through the OAR, the outflow extracted from the NSS alone cannot screen the electric field just above the NSS. In this paper, we analytically and numerically study the electric-field screening at the NSS taking into account the effects of the back-flowing particles from the OAR. In certain limited cases, the electric field is screened without significant pair cascade if only ultrarelativistic particles ($\gamma\gg1$) flow back to the NSS. On the other hand, if electron--positron pairs with a significant number density and mildly relativistic temperature, expected to distribute in a wide region of the magnetosphere, flow back to the NSS, these particles adjust the current and charge densities, so that the electric field can be screened without pair cascade. We obtain the condition for the number density of particles to screen the electric field at the NSS. We also find that in ion-extracted case from the NSS, bunches of particles are ejected to the outer region quasi-periodically, which is a possible mechanism of observed radio emission.

We report on a very long baseline interferometry (VLBI) experiment on giant radio pulses (GPs) from the Crab pulsar in the radio 1.4 to 1.7 GHz range to demonstrate a VLBI technique for searching for fast radio bursts (FRBs). We carried out the experiment on 26 July 2014 using the Kashima 34 m and Usuda 64 m radio telescopes of the Japanese VLBI Network (JVN) with a baseline of about 200 km. During the approximately 1 h observation, we could detect 35 GPs by high-time-resolution VLBI. Moreover, we determined the dispersion measure (DM) to be 56.7585 +/- 0.0025 on the basis of the mean DM of the 35 GPs detected by VLBI. We confirmed that the sensitivity of a detection of GPs using our technique is superior to that of a single-dish mode detection using the same telescope.

The Square Kilometre Array will revolutionize pulsar studies with its wide field-of-view, wide-band observation and high sensitivity, increasing the number of observable pulsars by more than an order of magnitude. Pulsars are of interest not only for the study of neutron stars themselves but for their usage as tools for probing fundamental physics such as general relativity, gravitational waves and nuclear interaction. In this article, we summarize the activity and interests of SKA-Japan Pulsar Science Working Group, focusing on an investigation of modified gravity theory with the supermassive black hole in the Galactic Centre, gravitational-wave detection from cosmic strings and binary supermassive black holes, a study of the physical state of plasma close to pulsars using giant radio pulses and determination of magnetic field structure of Galaxy with pulsar pairs.

We study stochastic acceleration models for the Fermi bubbles. Turbulence is excited just behind the shock front via Kelvin--Helmholtz, Rayleigh--Taylor, or Richtmyer--Meshkov instabilities, and plasma particles are continuously accelerated by the interaction with the turbulence. The turbulence gradually decays as it goes away from the shock fronts. Adopting a phenomenological model for the stochastic acceleration, we explicitly solve the temporal evolution of the particle energy distribution in the turbulence. Our results show that the spatial distribution of high-energy particles is different from those for a steady solution. We also show that the contribution of electrons that escaped from the acceleration regions significantly softens the photon spectrum. The photon spectrum and surface brightness profile are reproduced by our models. If the escape efficiency is very high, the radio flux from the escaped low-energy electrons can be comparable to that of the WMAP haze. We also demonstrate hadronic models with the stochastic acceleration, but they are unlikely in the viewpoint of the energy budget.

The spectral shape of the prompt emissions of gamma-ray bursts (GRBs) is typically expressed by the Band function: smooth joining of two power-law functions for high-energy and low-energy regions. To reveal the origin of the Band function, we revisit the stochastic acceleration model, in which electrons are accelerated via scattering with turbulent waves in the GRB outflow. The balance between the acceleration and synchrotron cooling yields a narrow energy-distribution similar to the Maxwellian distribution. The synchrotron spectrum becomes consistent with the observed hard photon index for the low-energy region. On the other hand, the narrow electron energy distribution contradicts the power-law spectrum for the high-energy region. We consider an evolution of the electron energy distribution to solve this problem. The turbulence and magnetic field induced by a certain hydrodynamical instability gradually decay. According to this evolution, the typical synchrotron photon energy also decreases with time. The time-integrated spectrum forms the power-law shape for the high-energy region. We discuss the required evolutions of the turbulence and magnetic field to produce a typical Band function. Although the decay of the turbulence is highly uncertain, recent numerical simulations for decaying turbulence seem comparatively positive for the stochastic acceleration model. Another condition required to reconcile observations is a much shorter duration of the stochastic acceleration than the dynamical timescale.

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.

Fermi-LAT reveals two huge gamma-ray bubbles existing in the Galactic Center, called 'Fermi Bubbles'. The existence of two microwave bubbles at the same region are also reported by the observation by WMAP, dubbed 'WMAP haze'. In order to explain these components, It has been argued that the gamma-rays arise from Inverse-Compton scattering of relativistic electrons accelerated by plasma turbulence, and the microwaves are radiated by synchrotron radiation. But no previous research reproduces both the Fermi Bubbles and WMAP haze under typical magnetic fields in the galaxy. We assume that shocks present in the bubbles and the efficiency of the acceleration by plasma turbulence, 'stochastic acceleration', changes with the distance from the shock front. The distance from the shock front increases with time, accordingly the efficiency of the acceleration changes with time. We also consider the time development of the electrons escape from the turbulence by diffusive loss. Our model succeed to reproduce both the observed characteristics of the Fermi Bubbles and WMAP haze under typical magnetic fields.

In this paper we propose a new mechanism describing the initial spike of giant flares in the framework of the starquake model. We investigate the evolution of a plasma on a closed magnetic flux tube in the magnetosphere of a magnetar in the case of a sudden energy release and discuss the relationship with observations of giant flares. We perform one-dimensional numerical simulations of the relativistic magnetohydrodynamics in Schwarzschild geometry. We assume energy is injected at the footpoints of the loop by a hot star surface containing random perturbations of the transverse velocity. Alfv\'en waves are generated and propagate upward, accompanying very hot plasma that is also continuously heated by nonlinearly generated compressive waves. We find that the front edges of the fireball regions collide at the top of the tube with their symmetrically launched counterparts. This collision results in an energy release which can describe the light curve of initial spikes of giant flares.

We present the results of the search for a correlation between giant radio pulses (GRPs) at 1.4 GHz and hard X-rays at 15-75 keV from the Crab pulsar. We made simultaneous ground and satellite observations of the Crab pulsar over 12 hours in three occasions in April 2010, March and September 2011, and got a sample of 1.3*10^4 main-pulse phase GRPs. From these samples we have found statistically marginal enhancement (21.5%, 2.70 sigma) of hard X-ray flux within +/- 1.5 degree phase angle of the synchronous peak of main-pulse phase GRPs. This enhancement, if confirmed, implicates that GRPs may accompany plasma density increases in the pulsar magnetosphere.

A review is given for cyclotron resonant interactions in space plasmas. After giving a simple formulation for the test particle approach, illustrative examples for resonant interactions are given. It is shown that for obliquely propagating whistler waves, not only fundamental cyclotron resonance, but also other resonances, such as transit-time resonance, anomalous cyclotron resonance, higher-harmonic cyclotron resonance, and even subharmonic resonance can come into play. A few recent topics of cyclotron resonant interactions, such as electron injection in shocks, cyclotron resonant heating of solar wind heavy ions, and relativistic modifications, are also reviewed.

We study effects of particle re-acceleration (or heating) in the post-shock region via magnetohydrodynamic/plasma turbulence, in the context of a mixed hadronic-leptonic model for the prompt emission of gamma-ray bursts (GRBs), using both analytical and numerical methods. We show that stochastically accelerated (or heated) leptons, which are injected via pp and pg reactions and subsequent pair cascades, are plausibly able to reproduce the Band function spectra with alpha~1 and beta~2-3 in the ~MeV range. An additional hard component coming from the proton-induced cascade emission is simultaneously expected, which is compatible with observed extra power-law spectra far above the MeV range. We also discuss the specific implications of hadronic models for ongoing high-energy neutrino observations.

{\it Fermi} satellite has detected extra spectral components in GeV energy range in several GRBs. Those components have power-law shapes, which may contribute to also X-ray band. The limited photon statistics make it difficult to determine the origin of GeV photons, namely internal or external shocks. We try to explain the extra components with our numerical simulations based on internal dissipation picture. Our leptonic model may reproduce not only the GeV excess via SSC emission but also the low-energy excess via the late synchrotron emission from remnant electrons. The hadronic models also reproduce keV-GeV power-law components by synchrotron and SSC emissions from secondary electron-positron pairs. In most cases the hadronic models require a much larger energy of protons than gamma-rays. However, the keV-GeV flat spectra detected in GRB 090902B is well explained with a comparable energy in protons and gamma-rays Finally, we discuss both advantages and weaknesses for both the leptonic and hadronic models. To overcome difficulties in internal dissipation models, we propose introduction of continuous acceleration similar to the second-order Fermi acceleration.

We report on the first detection of ionospheric disturbances caused by short repeated gamma-ray bursts from the magnetar SGR J1550-5418. Very low frequency (VLF) radio wave data obtained in South America clearly show sudden amplitude and phase changes at the corresponding times of eight SGR bursts. Maximum amplitude and phase changes of the VLF signals appear to be correlated with the gamma-ray fluence. On the other hand, VLF recovery timescales do not show any significant correlation with the fluence, possibly suggesting that the bursts' spectra are not similar to each other. In summary, the Earth's ionosphere can be used as a very large gamma-ray detector and the VLF observations provide us with a new method to monitor high energy astrophysical phenomena without interruption such as Earth Occultation.

We present a theoretical model describing magnetar giant flares on the basis of solar flare/coronal mass ejection theory. In our model, a preflare activity plays a crucial role in driving evaporating flows and supplying baryonic matters into the magnetosphere. The loaded baryonic matter, that is called "prominence", is then gradually uplifted via crustal cracking with maintaining a quasi-force-free equilibrium of the magnetosphere. Finally the prominence is erupted by the magnetic pressure force due to the loss of equilibrium triggered by the explosive magnetic reconnection. The giant flare should be induced as a final outcome of the prominence eruption accompanied by large-scale field reconfigurations. An essential difference between solar and magnetar flares is the control process of their evolutionary dynamics. The flaring activity on magnetars is mainly controlled by the radiative process unlike the solar flare governed by the electron conduction. It is highly suggestive that our model is accountable for the physical properties of the extraordinary giant flare observed on 2004 December 27 from SGR1806-20, including the source of baryonic matters loaded in the expanding ejecta observed after the giant burst.

We propose a new mechanism for the prompt emission of gamma-ray burst. In our model electrons are continuously accelerated in the post shock region via plasma turbulence. Using the Monte Carlo technique, we mimic the second-order Fermi acceleration due to plasma turbulence and obtain photon spectra. Since the acceleration balances with the synchrotron cooling, the observed low-energy spectral index is naturally explained. The resultant spectra can be consistent with observed spectra at least below $\sim 1$ MeV. The model also predicts delayed GeV-TeV emission due to inverse Compton and broad pulse profile of optical emission in some cases. Although nontrivial assumptions are required to reproduce MeV-GeV power-law spectra, the model implies the possibility to explain various kinds of luminosity correlations.

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.

Earth's bow shock is known to produce non-thermal electrons which are generally observed as a `spike' in their flux profile. Here, in this paper, we present an analysis of electron and whistler wave properties for a quasi-perpendicular shock crossing that is supercritical, but subcritical to the so-called whistler critical Mach number, M$^w_{\rm crit}$, above which whistler waves cannot propagate upstream. We have found that the amplitudes of whistler waves increased exponentially as a function of time prior to the shock encounter, while the suprathermal ($>$ 2 keV) electron flux similarly increased with time, although with differing $e$-folding time scales. Comparison of the electron energy spectrum measured within the ramp with predictions from diffusive shock acceleration theory was poor, but the variation of pitch angle distribution showed scattering of non-thermal electrons in the upstream region. While not finding a specific mechanism to account for the electron diffusion, we suggest that the whistlers seen probably account for the differences observed between this `gradual' event and the `spike' events seen at shocks with no upstream whistlers.

We present the unsaturated peak profile of SGR 1900+14 giant flare on 1998 August 27. This was obtained by particle counters of the Low Energy Particle instrument onboard the GEOTAIL spacecraft. The observed peak profile revealed four characteristic structures: initial steep rise, intermediate rise to the peak, exponential decay and small hump in the decay phase. From this light curve, we found that the isotropic peak luminosity was $2.3\times10^{46}$ erg s$^{-1}$ and the total energy was $4.3 \times 10^{44}$ erg s$^{-1}$ ($E\gtrsim$ 50 keV), assuming that the distance to SGR 1900+14 is 15 kpc and that the spectrum is optically thin thermal bremsstrahlung with $kT =$ 240 keV. These are consistent with the previously reported lower limits derived from Ulysses and Konus-Wind observations. A comparative study of the initial spikes of SGR 1900+14 giant flare in 1998 and SGR 1806-20 in 2004 is also presented. The timescale of the initial steep rise shows the magnetospheric origin, while the timescale of the intermediate rise to the peak indicates that it originates from the crustal fracturing. Finally, we argue that the four structures and their corresponding timescales provide a clue to identify extragalactic SGR giant flares among short GRBs.

On December 27, 2004, plasma particle detectors on the GEOTAIL spacecraft detected an extremely strong signal of hard X-ray photons from the giant flare of SGR1806-20, a magnetar candidate. While practically all gamma-ray detectors on any satellites were saturated during the first ~500 ms interval after the onset, one of the particle detectors on GEOTAIL was not saturated and provided unique measurements of the hard X-ray intensity and the profile for the first 600 ms interval with 5.48 ms time resolution. After ~50 ms from the initial rapid onset, the peak photon flux (integrated above ~50 keV) reached the order of 10^7 photons sec^{-1} cm^{-2}. Assuming a blackbody spectrum with kT=175 keV, we estimate the peak energy flux to be 21 erg sec^{-1} cm^{-2} and the fluence (for 0-600 ms) to be 2.4 erg cm^{-2}. The implied energy release comparable to the magnetic energy stored in a magnetar (~10^{47} erg) suggests an extremely efficient energy release mechanism.

The outer shells of young supernova remnants (SNRs) are the most plausible acceleration sites of high-energy electrons with the diffusive shock acceleration (DSA) mechanism. We studied spatial and spectral properties close to the shock fronts in four historical SNRs (Cas A, Kepler's remnant, Tycho's remnant, and RCW 86) with excellent spatial resolution of {\it Chandra}. In all of the SNRs, hard X-ray emissions were found on the rims of the SNRs, which concentrate in very narrow regions (so-called "filaments"); apparent scale widths on the upstream side are below or in the order of the point spread function of {\it Chandra}, while 0.5--40 arcsec (0.01--0.4 pc) on the downstream side with most reliable distances. The spectra of these filaments can be fitted with both thermal and nonthermal (power-law and {\tt SRCUT}) models. The former requires unrealistic high temperature ($\ga$2 keV) and low abundances ($\la$1 solar) for emission from young SNRs and may be thus unlikely. The latter reproduces the spectra with best-fit photon indices of 2.1--3.8, or roll-off frequencies of (0.1--28)$\times 10^{17}$ Hz, which reminds us of the synchrotron emission from electrons accelerated via DSA. We consider various physical parameters as functions of the SNR age, including the previous results on SN 1006 \citep{bamba2003b}; the filament width on the downstream side increases with the SNR age, and the spectrum becomes softer keeping a nonthermal feature. It was also found that a function, that is the roll-off frequency divided by the square of the scale width on the downstream side, shows negative correlation with the age, which might provide us some information on the DSA theory.

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