Low energy ground-based cosmic ray air shower experiments generally have energy threshold in the range of a few tens to a few hundreds of TeV. The shower observables are measured indirectly with an array of detectors. The atmospheric absorption of low energy secondaries limits their detection frequencies at the Earth's surface. However, due to selection effects, a tiny fraction of low energy showers, which are produced in the lower atmosphere can reach the observational level. But, due to less information of shower observables, the reconstruction of these showers are arduous. Hence, it is believed that direct measurements by experiments aboard on satellites and balloon flights are more reliable at low energies. Despite having very small efficiency ($\sim$0.1%) at low energies, the large acceptance ($\sim$5 m$^2$sr) of GRAPES-3 experiment allows observing primary cosmic rays down below to $\sim$1 TeV and opens up the possibility to measure primary energy spectrum spanning from a few TeV to beyond cosmic ray knee (up to 10$^{16}$ eV), covering five orders of magnitude. The GRAPES-3 energy threshold for primary protons through Monte Carlo simulations are calculated, which gives reasonably good agreement with data. Furthermore, the total efficiencies and acceptance are also calculated for protons primaries. The ability of GRAPES-3 experiment to cover such a broader energy range may provide a unique handle to bridge the energy spectrum between direct measurements at low energies and indirect measurements at ultra-high energies.

The angular resolution of an extensive air shower (EAS) array plays a critical role in determining its sensitivity for the detection of point $\gamma$-ray sources in the multi-TeV energy range. The GRAPES-3 an EAS array located at Ooty in India (11.4$^{\circ}$N, 76.7$^{\circ}$E, 2200 m altitude) is designed to study $\gamma$-rays in the TeV-PeV energy range. It comprises of a dense array of 400 plastic scintillators deployed over an area of 25000 m$^2$ and a large area (560 m$^2$) muon telescope. A new statistical method allowed real time determination of the propagation delay of each detector in the GRAPES-3 array. The shape of shower front is known to be curved and here the details of a new method developed for accurate measurement of the shower front curvature is presented. These two developments have led to a sizable improvement in the angular resolution of GRAPES-3 array. It is shown that the curvature depends on the size and age of an EAS. By employing two different techniques, namely, the odd-even and the left-right methods, independent estimates of the angular resolution are obtained. The odd-even method estimates the best achievable resolution of the array. For obtaining the angular resolution, the left-right method is used after implementing the size and age dependent curvature corrections. A comparison of the angular resolution as a function of EAS energy by these two methods shows them be virtually indistinguishable. The angular resolution of GRAPES-3 array is 47$^{\prime}$ for energies E$>$5 TeV and improves to 17$^{\prime}$ at E$>$100 TeV and finally approaching 10$^{\prime}$ at E$>$500 TeV.

Context.Recent studies suggest that many giant exoplanets are highly enriched with heavy elements compared to their host star andcontain several tens of Earth masses or more of heavy elements. Such enrichment is considered to have been brought by accretionof planetesimals in late formation stages. Previous dynamical simulations, however, show that planets are unable to collect so muchheavy elements throughin situplanetesimal accretion. Aims.We investigate whether a giant planet migrating inward can capture planetesimals efficiently to significantly increase its metal-licity. Methods.We performed orbital integrations of a migrating giant planet and planetesimals in a protoplanetary gas disc to infer theplanetesimal mass that is accreted by the planet. Results.We find that the two shepherding processes of mean motion resonances trapping and aerodynamic gas drag inhibit plan-etesimal capture of a migrating planet. However, the amplified libration allows the highly-excited planetesimals in the resonances toescape from the resonance trap and be accreted by the planet. Consequently, we show that a migrating giant planet captures planetes-imals with total mass of several tens of Earth masses, if the planet forms at a few tens of AU in a relatively massive disc. We alsofind that planetesimal capture occurs efficiently in a limited range of semi-major axis, and that the total captured planetesimal massincreases with increasing migration distances. Our results have important implications for understanding the relation between giantplanet metallicity and mass, as we suggest that it reflects the formation location of the planet, or more precisely, the location whererunaway gas accretion occurred. We also suggest the observed metal-rich close-in Jupiters migrated to their present locations fromafar, where they formed.

Observations of the molecular gas in galaxies are vital to understanding the evolution and star-forming histories of galaxies. However, galaxies with molecular gas maps of their whole discs having sufficient resolution to distinguish galactic structures are severely lacking. Millimeter wavelength studies at a high angular resolution across multiple lines and transitions are particularly needed, severely limiting our ability to infer the universal properties of molecular gas in galaxies. Hence, we conducted a legacy project with the 45 m telescope of the Nobeyama Radio Observatory, called the CO Multi-line Imaging of Nearby Galaxies (COMING), which simultaneously observed 147 galaxies with high far-infrared flux in $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=1-0$ lines. The total molecular gas mass was derived using the standard CO-to-H$_2$ conversion factor and found to be positively correlated with the total stellar mass derived from the WISE $3.4 \mu$m band data. The fraction of the total molecular gas mass to the total stellar mass in galaxies does not depend on their Hubble types nor the existence of a galactic bar, although when galaxies in individual morphological types are investigated separately, the fraction seems to decrease with the total stellar mass in early-type galaxies and vice versa in late-type galaxies. No differences in the distribution of the total molecular gas mass, stellar mass, and the total molecular gas to stellar mass ratio was observed between barred and non-barred galaxies, which is likely the result of our sample selection criteria, in that we prioritized observing FIR bright (and thus molecular gas-rich) galaxies.

The SEDA-FIB is a detector designed to measure solar neutrons. This solar neutron detector was operated onboard the ISS on July 16, 2009 and March 31, 2018. Eighteen large solar flares were later observed by the GOES satellite in solar active region 12673 that appeared on September 4 and lasted until September 10, 2017, with intensity higher than > M2. In nine of those solar flares, the SEDA-FIB detected clear signals of solar neutrons, along with five minor excesses. Among these events, we focus on two associated with the flares of X2.2 (SOL2017-09-06) and X8.2 (SOL2017-09-10) that share a common feature: a process of accelerating electrons into high energies as clearly recorded by the FERMI-GBM detector. These events may provide us with useful information to elucidate the ion acceleration process. The X8.2 event was a limb flare that proved adequate for fixing the parameters needed to explain the process of particle acceleration into high energies. According to our analysis, the electron acceleration process may possibly be explained by the shock acceleration model. However, we found that it would be difficult to explain the simultaneous acceleration of ions with electrons, unless the ions were preheated prior to their rapid acceleration.

Studies of internal structure of gas giant planets suggest that their envelopes are enriched with heavier elements than hydrogen and helium relative to their central stars. Such enrichment likely occurred by solid accretion during late formation stages of gas giant planets in which gas accretion dominates protoplanetary growth. Some previous studies performed orbital integration of planetesimals around a growing protoplanet with the assumption of an uniform circumstellar disc to investigate how efficiently the protoplanet captures planetesimals. However, not only planetesimals but also disc gas are gravitationally perturbed by the protoplanet in its late formation stages, resulting in gap opening in the circumstellar disc. In this study, we investigate the effects of such gap formation on the capture of planetesimals by performing dynamical simulations of planetesimals around a growing proto-gas giant planet. Gap formation reduces the surface density of disc gas, makes a steep pressure gradient and limits the growth rate of the protoplanet. We find that the first effect enhances the capture of planetesimals, while the others reduce it. Consequently the amount of planetesimals captured during the gas accretion is estimated to be at most ~3 Mearth. We conclude that the in-situ capture of planetesimals needs the initial solid surface density more than five times higher than that of the minimum mass solar nebula for explaining the inferred large amount of heavy element in Jupiter. For highly dense warm Jupiters, we would need additional processes enhancing the capture and/or supply of planetesimals.

We present the X-ray survey results of high-magnetic-field radio pulsars (high-B PSRs) with Swift/XRT. X-ray observations of the rotation-powered pulsars with the dipole magnetic field ${\it B_{\rm d}}$ near the quantum critical field ${\it B_{\rm q}}=4.4\times10^{13}$ G is of great importance for understanding the transition between the rotation-powered pulsars and the magnetars, because there are a few objects that have magnetar-like properties. Out of the 27 high-B PSRs that are in the ATNF pulsar catalogue but have not been reported or have no effective upper-limits in the X-ray bands, we analyze the Swift/XRT data for 21 objects, where 6 objects are newly observed and 15 objects are taken from the archival data. As a result, we have new $3\sigma$ upper-limits for all the 21 objects. Since the upper-limits are tight, we conclude that we do not find any magnetar-like high-B PSRs such as PSR J1819$-$1458. The probability of the high X-ray efficiency in the high-B PSRs is obtained to be $11\%-29\%$. Combining the previous observations, we discuss which parameter causes magnetar-like properties. It may be suggested that the magnetar-like properties appear only when ${\it B_{\rm d}}\gtrsim10^{13.5}$ G for the radio pulsar population. This is true even if the radio-quiet high-B RPP are included.

We observe the magnetar CXOU J171405.7-381031 with XMM-Newton and obtain the most reliable X-ray spectral parameters for this magnetar. After removing the flux from the surrounding supernova remnant CTB~37B, the radiation of CXOU J171405.7-381031 is best described by a two-component model, consisting of a blackbody and power law. We obtain a blackbody temperature of 0.58^{+0.03}_{-0.03} keV, photon index of 2.15^{+0.62}_{-0.68}, and unabsorbed 2-10 keV flux of 2.33^{+0.02}_{-0.02} x 10^{-12} erg cm^{-2} s^{-1}. These new parameters enable us to compare CXOU 171405.7-381031 with other magnetars, and it is found that the luminosity, temperature and the photon index of CXOU J171405.7-381031 are aligned with the known trend among the magnetar population with a slightly higher temperature, which could be caused by its young age. All the magnetars with a spin-down age of less than 1~kyr show time variation or bursts except for CXOU J171405.7-381031. We explore the time variability for ten observations in between 2006 and 2015, but there is no variation larger than \sim 10%.

We present the results of $^{12}$CO($J$=1-0) and $^{13}$CO($J$=1-0) simultaneous mappings toward the nearby barred spiral galaxy NGC 4303 as a part of the CO Multi-line Imaging of Nearby Galaxies (COMING) project. Barred spiral galaxies often show lower star-formation efficiency (SFE) in their bar region compared to the spiral arms. In this paper, we examine the relation between the SFEs and the volume densities of molecular gas $n(\rm{H}_2)$ in the eight different regions within the galactic disk with CO data combined with archival far-ultraviolet and 24 $\mu$m data. We confirmed that SFE in the bar region is lower by 39% than that in the spiral arms. Moreover, velocity-alignment stacking analysis was performed for the spectra in the individual regions. The integrated intensity ratios of $^{12}$CO to $^{13}$CO ($R_{12/13}$) range from 10 to 17 as the results of stacking. Fixing a kinetic temperature of molecular gas, $n(\rm{H}_2)$ was derived from $R_{12/13}$ via non-local thermodynamic equilibrium (non-LTE) analysis. The density $n(\rm{H}_2)$ in the bar is lower by 31-37% than that in the arms and there is a rather tight positive correlation between SFEs and $n(\rm{H}_2)$, with a correlation coefficient of $\sim 0.8$. Furthermore, we found a dependence of $n(\rm{H}_2)$ on the velocity dispersion of inter-molecular clouds ($\Delta V/ \sin i$). Specifically, $n(\rm{H}_2)$ increases as $\Delta V/ \sin i$ increases when $\Delta V/ \sin i < 100$ km s$^{-1}$. On the other hand, $n(\rm{H}_2)$ decreases as $\Delta V/ \sin i$ increases when $\Delta V/ \sin i > 100$ km s$^{-1}$. These relations indicate that the variations of SFE could be caused by the volume densities of molecular gas, and the volume densities could be governed by the dynamical influence such as cloud-cloud collisions, shear and enhanced inner-cloud turbulence.

In the solar flare observed on June 3, 2012, high energy gamma-rays and neutrons were observed. The event includes a remarkable feature of a high neutron/gamma-ratio in the secondary particles. We have examined whether this high n/$\gamma$-ratio can be explained by simulation. As a result of simulations using the GEANT4 program, the high n/$\gamma$-ratio may be reproduced for the case that helium and other heavy ions were dominantly accelerated in the flare.

We examined radial variations in molecular-gas based star formation efficiency (SFE), which is defined as star formation rate per unit molecular gas mass, for 80 galaxies selected from the CO Multi-line Imaging of Nearby Galaxies project (Sorai et al. 2019). The radial variations in SFE for individual galaxies are typically a factor of 2 -- 3, which suggests that SFE is nearly constant along galactocentric radius. We found the averaged SFE in 80 galaxies of $(1.69 \pm 1.1) \times 10^{-9}$ yr$^{-1}$, which is consistent with Leroy et al. 2008 if we consider the contribution of helium to the molecular gas mass evaluation and the difference in the assumed initial mass function between two studies. We compared SFE among different morphological (i.e., SA, SAB, and SB) types, and found that SFE within the inner radii ($r/r_{25} < 0.3$, where $r_{25}$ is $B$-band isophotal radius at 25 mag arcsec$^{-2}$) of SB galaxies is slightly higher than that of SA and SAB galaxies. This trend can be partly explained by the dependence of SFE on global stellar mass, which probably relates to the CO-to-H$_2$ conversion factor through the metallicity. For two representative SB galaxies in our sample, NGC 3367 and NGC 7479, the ellipse of $r/r_{25}$ = 0.3 seems to cover not only the central region but also the inner part of the disk, mainly the bar. These two galaxies show higher SFE in the bar than in spiral arms. However, we found an opposite trend in NGC 4303; SFE is lower in the bar than in spiral arms, which is consistent with earlier studies (e.g., Momose et al. 2010). These results suggest diversity of star formation activities in the bar.

The $^{12}$CO $(J=1\rightarrow0)$ velocity fields of a sample of 20 nearby spiral galaxies, selected from the CO Multi-line Imaging of Nearby Galaxies (COMING) legacy project of Nobeyama Radio Observatory, have been analyzed by Fourier decomposition to determine their basic kinematic properties, such as circular and noncircular velocities. On average, the investigated barred (SAB and SB) galaxies exhibit a ratio of noncircular to circular velocities of molecular gas larger by a factor of 1.5-2 than non-barred (SA) spiral galaxies at radii within the bar semimajor axis $a_\mathrm{b}$ at 1 kpc resolution, with a maximum at a radius of $R/a_\mathrm{b}\sim0.3$. Residual velocity field images, created by subtracting model velocity fields from the data, reveal that this trend is caused by kpc-scale streaming motions of molecular gas in the bar region. Applying a new method based on radial velocity reversal, we estimated the corotation radius $R_\mathrm{CR}$ and bar pattern speed $\Omega_\mathrm{b}$ in seven SAB and SB systems. The ratio of the corotation to bar radius is found to be in a range of $\mathcal{R}\equiv R_\mathrm{CR}/a_\mathrm{b}\sim0.8\mathrm{-}1.6$, suggesting that intermediate (SBb-SBc), luminous barred spiral galaxies host fast and slow rotator bars. Tentative negative correlations are found for $\Omega_\mathrm{b}$ vs. $a_\mathrm{b}$ and $\Omega_\mathrm{b}$ vs. total stellar mass $M_\ast$, indicating that bars in massive disks are larger and rotate slower, possibly a consequence of angular momentum transfer. The kinematic properties of SAB and SB galaxies, derived from Fourier decomposition, are compared with recent numerical simulations that incorporate various rotation curve models and galaxy interactions.

The GRAPES-3 muon telescope in Ooty, India had claimed detection of a 2 hour (h) high-energy ($\sim$20 GeV) burst of galactic cosmic-rays (GCRs) through a $>$50$\sigma$ surge in GeV muons, was caused by reconnection of the interplanetary magnetic field (IMF) in the magnetosphere that led to transient weakening of Earth's magnetic shield. This burst had occurred during a G4-class geomagnetic storm (storm) with a delay of $\frac{1}{2}$h relative to the coronal mass ejection (CME) of 22 June 2015 (Mohanty et al., 2016). However, recently a group interpreted the occurrence of the same burst in a subset of 31 neutron monitors (NMs) to have been the result of an anisotropy in interplanetary space (Evenson et al., 2017) in contrast to the claim in (Mohanty et al., 2016). A new analysis of the GRAPES-3 data with a fine 10.6$^{\circ}$ angular segmentation shows the speculation of interplanetary anisotropy to be incorrect, and offers a possible explanation of the NM observations. The observed 28 minutes (min) delay of the burst relative to the CME can be explained by the movement of the reconnection front from the bow shock to the surface of Earth at an average speed of 35 km/s, much lower than the CME speed of 700 km/s. This measurement may provide a more accurate estimate of the start of the storm.

Recent observations found that electrons are accelerated to $\sim$10 GeV and emit synchrotron hard X-rays in two magnetic white dwarfs (WDs), also known as cataclysmic variables (CVs). In nova outbursts of WDs, multi-GeV gamma-rays were detected inferring that protons are accelerated to 100 GeV or higher. In recent optical surveys, the WD density is found to be higher near the Sun than in the Galactic disk by a factor $\sim$2.5. The cosmic rays (CR) produced by local CVs and novae will accumulate in the local bubble for $10^6$ - $10^7$ yrs. On these findings, we search for CRs from historic CVs and novae in the observed CR spectra. We model the CR spectra at the heliopause as sums of Galactic and local components based on observational data as much as possible. The initial Galactic CR electron and proton spectra are deduced from the gamma-ray emissivity, the local electron spectrum from the hard X-ray spectra at the CVs, and the local proton spectrum inferred by gamma-ray spectrum at novae. These spectral shapes are then expressed in a simple set of polynomial functions of CR energy and regressively fitted until the high-energy ($>$100 GeV) CR spectra near Earth and the Voyager-1 spectra at the heliopause are reproduced. We then extend the modeling to nuclear CR spectra and find that one spectral shape fits all local nuclear CRs and the apparent hardening of the nuclear CR spectra is caused by the roll-down of local nuclear spectra around 100 - 200 GeV. All local CR spectra populate in a limited energy band below 100 - 200 GeV and enhance gamma-ray emissivity below $\sim$10 GeV. Such an enhancement is observed in the inner Galaxy, suggesting the CR fluxes from CVs and novae are substantially higher there.

Recent localization of the repeating Fast Radio Burst (FRB) 121102 revealed the distance of its host galaxy and luminosities of the bursts. We investigated constraints on the young neutron star (NS) model, that (a) the FRB intrinsic luminosity is supported by the spin-down energy, and (b) the FRB duration is shorter than the NS rotation period. In the case of a circular cone emission geometry, conditions (a) and (b) determine the NS parameters within very small ranges, compared with that from only condition (a) discussed in previous works. Anisotropy of the pulsed emission does not affect the area of the allowed parameter region by virtue of condition (b). The determined parameters are consistent with those independently limited by the properties of the possible persistent radio counterpart and the circumburst environments such as surrounding materials. Since the NS in the allowed parameter region is older than the spin-down timescale, the hypothetical GRP-like model expects a rapid radio flux decay of $\lesssim1$ Jy within a few years as the spin-down luminosity decreases. The continuous monitoring will give a hint of discrimination of the models. If no flux evolution will be seen, we need to consider an alternative model, e.g., the magnetically powered flare.

At the 33rd ICRC, we reported the possible detection of solar gamma rays by a ground level detector and later re-examined this event. On March 7, 2011, the solar neutron telescope (SNT) located at Mt. Sierra Negra, Mexico (4,600 m) observed enhancements of the counting rate from 19:49 to 20:02 UT and from 20:50 to 21:01 UT. The statistical significance was 9.7sigma and 8.5sigma, respectively. This paper discusses the possibility of using this mountain detector to detect solar gamma rays. In association with this event, the solar neutron detector SEDA-FIB onboard the International Space Station has also detected solar neutrons with a statistical significance of 7.5sigma. The FERMI-LAT detector also observed high-energy gamma rays from this flare with a statistical significance of 6.7sigma. We thus attempted to make a unified model to explain this data. In this paper, we report on another candidate for solar gamma rays detected on September 25th, 2011 by the SNT located in Tibet (4,300 m) from 04:37 to 04:47 UT with a statistical significance of 8.0sigma (by the Li-Ma method).

Studies were made of the 1-70 keV persistent spectra of fifteen magnetars as a complete sample observed with Suzaku from 2006 to 2013. Combined with early NuSTAR observations of four hard X-ray emitters, nine objects showed a hard power-law emission dominating at $\gtrsim$10 keV with the 15--60 keV flux of $\sim$1-$11\times 10^{-11}$ ergs s$^{-1}$ cm$^{-2}$. The hard X-ray luminosity $L_{\rm h}$, relative to that of a soft-thermal surface radiation $L_{\rm s}$, tends to become higher toward younger and strongly magnetized objects. Updated from the previous study, their hardness ratio, defined as $\xi=L_{\rm h}/L_{\rm s}$, is correlated with the measured spin-down rate $\dot{P}$ as $\xi=0.62 \times (\dot{P}/10^{-11}\,{\rm s}\,{\rm s}^{-1})^{0.72}$, corresponding with positive and negative correlations of the dipole field strength $B_{\rm d}$ ($\xi \propto B_{\rm d}^{1.41}$) and the characteristic age $\tau_{\rm c}$ ($\xi \propto \tau_{\rm c}^{-0.68}$), respectively. Among our sample, five transients were observed during X-ray outbursts, and the results are compared with their long-term 1-10 keV flux decays monitored with Swift/XRT and RXTE/PCA. Fading curves of three bright outbursts are approximated by an empirical formula used in the seismology, showing a $\sim$10-40 d plateau phase. Transients show the maximum luminosities of $L_{\rm s}$$\sim$$10^{35}$ erg s$^{-1}$, which is comparable to those of the persistently bright ones, and fade back to $\lesssim$$10^{32}$ erg s$^{-1}$. Spectral properties are discussed in a framework of the magnetar hypothesis.

We reported a photon detection at the arrival time of the gravitational wave on December 26th, 2015. According to the LIGO-Virgo collaboration, the gravitational wave was produced by the merging process of the two black holes. The merged time was determined at 03:38:53.6 UT on December 26th. At 03:38:54.05 GPS time, one of the detectors of SEDA-FIB on board the International Space Station (ISS) detected a photon, arriving from the direction of Corona Australis near the Galactic Bulge (GB). The energy of the photon was about (35+-7) MeV. A 3.0 sigma level detection of the photon is discussed.

Previous works have suggested a correlation between the X-ray luminosity Lx and the rotational luminosity Lrot of the radio pulsars.However, none of the obtained regression lines are statistically acceptable due to large scatters. We construct a statistical model which has an intrinsic Lx-Lrot relation and reproduces the observed Lx distribution about it by using a Monte Carlo simulator, which takes into account the effects obscuring the intrinsic relation,i.e., the anisotropy of radiation, additional heating, uncertainty in distance and detection limit of the instruments. From the ATNF pulsar catalog we collect 57 `ordinary radio pulsars' with significant detection and 42 with upper limits.The sample does not include the high-magnetic field pulsars (>10^{13} G), which are separately analyzed. We obtain a statistically acceptable relation Lx (0.5 - 10 keV)= 10^{31.69} (Lrot / L_0)^{c_1} with c_1 = 1.03 \pm 0.27 and L_0 =10^{35.38}. The distribution about the obtained Lx-Lrot relation is reproduced well by the simulator. Pulsars with abnormally high Lx fall into two types: one is the soft gamma-ray pulsars, and the other is thermally bright pulsars in comparison with the standard cooling curve. On the other hand, pulsars showing low Lx are found to have dim pulsar wind nebulae. We argue that there is an unknown mechanism that governs both the magnetospheric emission and the pulsar wind nebulae, and it might involve the production rate of electron-positron pairs. The high-field pulsars form a distinctive population other than the ordinary pulsars due to their excess luminosities.

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.

An M6.5-class flare was observed at N12E56 of the solar surface at 16:06 UT on July 8, 2014. In association with this flare, solar neutron detectors located on two high mountains, Mt. Sierra Negra and Chacaltaya and at the space station observed enhancements in the neutral channel. The authors analysed these data and a possible scenario of enhancements produced by high-energy protons and neutrons is proposed, using the data from continuous observation of a solar surface by the ultraviolet telescope onboard the Solar Dynamical Observatory (SDO).

The solar neutron detector SEDA-FIB onboard the International Space Station (ISS) has detected several events from the solar direction associated with three large solar flares observed on March 5th (X1.1), 7th (X5.4), and 9th (M6.3) of 2012. In this study, we present the time profiles of those neutrons and discuss the physics that may be related to a possible acceleration scenario for ions over the solar surface. We compare our data with the dynamical pictures of the flares obtained by the ultra-violet telescope of the space-based Solar Dynamics Observatory.

We analyzed a highly impulsive solar flare observed on June 3, 2012. In association with this flare, emissions of hard X-rays, high-energy gamma rays, and neutrons were detected by the detectors onboard the FERMI, RHESSI satellites and the International Space Station. We compared those results with the pictures taken by the UV telescope onboard the Solar Dynamics Observatory satellite and found the crossing structure of two magnetic ropes at two positions on the solar surface almost at the same time. High-energy gamma rays were detected by the Fermi Large Area Telescope satellite, implying that the impulsive flare was one of a major source of proton acceleration processes on the solar surface. At the beginning of research, impulsive solar flares were considered to be the main source of particle acceleration processes; our current observations have confirmed this hypothesis.

We develop a time-dependent multi-group multidimensional relativistic radiative transfer code, which is required to numerically investigate radiation from relativistic fluids involved in, e.g., gamma-ray bursts and active galactic nuclei. The code is based on the spherical harmonic discrete ordinate method (SHDOM) that evaluates a source function including anisotropic scattering in spherical harmonics and implicitly solves the static radiative transfer equation with a ray tracing in discrete ordinates. We implement treatments of time dependence, multi-frequency bins, Lorentz transformation, and elastic Thomson and inelastic Compton scattering to the publicly available SHDOM code. Our code adopts a mixed frame approach; the source function is evaluated in the comoving frame whereas the radiative transfer equation is solved in the laboratory frame. This implementation is validated with various test problems and comparisons with results of a relativistic Monte Carlo code. These validations confirm that the code correctly calculates intensity and its evolution in the computational domain. The code enables us to obtain an Eddington tensor that relates first and third moments of intensity (energy density and radiation pressure) and is frequently used as a closure relation in radiation hydrodynamics calculations.

We investigate the fraction of metal nuclei in the relativistic jets of gamma-ray bursts associated with core-collapse supernovae. We simulate the fallback in jet-induced explosions with two-dimensional relativistic hydrodynamics calculations and the jet acceleration with steady, radial, relativistic magnetohydrodynamics calculations, and derive detail nuclear composition of the jet by postprocessing calculation. We found that if the temperature at the jet launch site is above $4.7\times 10^9$K, quasi-statistical equilibrium (QSE) is established and heavy nuclei are dissociated to light particles such as $^4$He during the acceleration of the jets. The criterion for the survival of metal nuclei is written in terms of the isotropic jet luminosity as $L_{\rm j}^{\rm iso} \lesssim 3.9\times 10^{50}(R_{\rm i}/10^7{\rm cm})^2 (1+\sigma_{\rm i})~{\rm erg~s^{-1}}$, where $R_{\rm i}$ and $\sigma_{\rm i}$ are the initial radius of the jets and the initial magnetization parameter, respectively. If the jet is initially dominated by radiation field (i.e., $\sigma_{\rm i} \ll 1$) and the isotropic luminosity is relatively high ($L_{\rm j}^{\rm iso}\gtrsim 4\times 10^{52}~{\rm erg s^{-1}}$), the metal nuclei cannot survive in the jet. On the other hand, if the jet is mainly accelerated by magnetic field (i.e., $\sigma_{\rm i} \gg1$), metal nuclei initially contained in the jet can survive without serious dissociation even for the case of high luminosity jet. If the jet contains metal nuclei, the dominant nuclei are $^{28}$Si, $^{16}$O, and $^{32}$S and the mean mass number can be $\langle$A$\rangle\sim25$.

The 2.1-s anomalous X-ray pulsar 1E 1547.0-5408 exhibited an X-ray outburst on 2009 January 22, emitting a large number of short bursts. The wide-band all- sky monitor (WAM) on-board Suzaku detected at least 254 bursts in the 160keV-6.2MeV band over the period of January 22 00:57-17:02 UT from the direction of 1E 1547.0-5408. One of these bursts, which occurred at 06:45:13, produced the brightest fluence in the 0.5-6.2MeV range, with an averaged 0.16-6.2MeV flux and extrapolated 25 keV-2 MeV fluence of about 3x10-6 erg cm-2 s-1 and about 3x10-4 erg cm-2, respectively. After pile-up corrections, the time-resolved WAM spectra of this burst were well-fitted in the 0.16-6.2MeV range by two-component models; specifically, a blackbody plus an optically thin thermal bremsstrahlung or a combination of a blackbody and a power-law component with an exponential cutoff. These results are compared with previous works reporting the persistent emission and weaker short bursts followed by the same outburst.

Multi-wavelength observations of the black-widow binary system 2FGL J2339.6-0532 are reported. The Fermi gamma-ray source 2FGL J2339.6-0532 was recently categorized as a black widow in which a recycled millisecond pulsar (MSP) is evaporating up the companion star with its powerful pulsar wind. Our optical observations show clear sinusoidal light curves due to the asymmetric temperature distribution of the companion star. Assuming a simple geometry, we constrained the range of the inclination angle of the binary system to 52$^{\circ}$ < i < 59$^{\circ}$, which enables us to discuss the interaction between the pulsar wind and the companion in detail. The X-ray spectrum consists of two components: a soft, steady component that seems to originate from the surface of the MSP, and a hard variable component from the wind-termination shock near the companion star. The measured X-ray luminosity is comparable to the bolometric luminosity of the companion, meaning that the heating efficiency is less than 0.5. In the companion orbit, 10$^{11}$ cm from the pulsar, the pulsar wind is already in particle dominant-stage, with a magnetization parameter of $\sigma$ < 0.1. In addition, we precisely investigated the time variations of the X-ray periodograms and detected a weakening of orbital modulation. The observed phenomenon may be related to an unstable pulsar-wind activity or a weak mass accretion, both of which can result in the temporal extinction of radio-pulse.

We investigate the random walk process in relativistic flow. In the relativistic flow, photon propagation is concentrated in the directions of the flow velocity due to relativistic beaming effect. We show that, in the pure scattering case, the number of scatterings is proportional to the size parameter $\xi\equiv L/l_0$ if the flow velocity $\beta\equiv v/c$ satisfies $\beta/\Gamma\gg \xi^{-1}$, while it is proportional to $\xi^2$ if $\beta/\Gamma\ll \xi^{-1}$ where $L$ and $l_0$ are the size of the system in the observer frame and the mean free path in the comoving frame, respectively. We also examine the photon propagation in the scattering and absorptive medium. We find that, if the optical depth for absorption $\tau_{\rm a}$ is considerably smaller than the optical depth for scattering $\tau_{\rm s}$ ($\tau_{\rm a}/\tau_{\rm s} \ll 1$) and the flow velocity satisfies $\beta\gg \sqrt{2\tau_{\rm a}/\tau_{\rm s}}$, the effective optical depth is approximated by $\tau_*\simeq\tau_{\rm a}(1+\beta)/\beta$. Furthermore, we perform Monte Carlo simulations of radiative transfer and compare the results with the analytic expression for the number of scattering. The analytic expression is consistent with the results of the numerical simulations. The expression derived in this Letter can be used to estimate the photon production site in relativistic phenomena, e.g., gamma-ray burst and active galactic nuclei.

The relationship between solar activity and the global climate is not only an academically interesting issue, but also an important problem for human beings. Lean and Rind have analyzed a considerable amount of climate data from around the world from 1889 to 2006. According to their analysis, the global effect was estimated to be 0.17 plus-minus 0.01K between the solar maximum and minimum. However, they noticed that the effect strongly appeared in the zones between 70N and 30N, and between 25S and 50S. At its peak latitude (near 40), the effect was estimated to be 0.5 K. Therefore, we analyzed a tree that survived at the Murooji temple in Nara prefecture (34d32mN, 136d2mE) Japan, for 391 years. Quite surprisingly, Fourier analysis of the annual growth rate identified two cycles with periodicities of 12 and 25 years during the Maunder minimum. We have continued the analysis, using different samples from the Nagusa shrine in Hyogo Prefecture, Itayanagi City, Aomori Prefecture, Japan, and from Yaku Island in Kyusyu, Japan. An evidence of solar activity was found in all of the samples. In particular, we have discovered a correlation between Swiss glacier fluctuation and the growth rate of the Yaku tree ring.

A new type of solar neutron detector (SEDA-FIB) was launched on board the Space Shuttle Endeavor on July 16 2009, and began collecting data at the International Space Station (ISS) on August 25 2009. This paper summarizes four years of observations with the solar neutron detector SEDA-FIB (Space Environment Data Acquisition using the FIBer detector). The solar neutron detector FIB can determine both the energy and arrival direction of solar neutrons. In this paper, we first present the angular distribution of neutron induced protons obtained in Monte Carlo simulations. The results are compared with the experimental results. Then we provide the angular distribution of background neutrons during one full orbit of the ISS (90 minutes). Next, the angular distribution of neutrons during the flare onset time from 20:02 to 20:10 UT on March 7 2011 is presented. It is compared with the distribution when a solar flare is not occurring. Observed solar neutrons possibly originated from the M-class solar flares that occurred on March 7 (M3.7), June 7 (M2.5), September 24 (M3.0) (weak signal) and November 3 (X1.9) of 2011 and January 23 of 2012 (M8.7). This marked the first time that neutrons have been observed from M-class solar flares. A possible interpretation of the neutron production process will be also provided.

On March 7, 2011 from 19:48:00 to 20:03:00 UT, the solar neutron telescope located at Mt. Sierra Negra, Mexico (4,600m) observed a 8.8sigma enhancement. In this paper, we would like to try to explain this enhancement by a hypothesis that a few GeV gamma-rays arrived at the top of the mountain produced by the Sun. We postulate that protons were accelerated at the shock front. They precipitate at the solar surface and produced those gamma-rays. If hypothesis is confirmed, this enhancement is the first sample of GeV gamma-rays observed by a ground level detector.

A new type of solar neutron detector (FIB) was launched onboard the Space Shuttle Endeavour on July 16, 2009, and it began collecting data at the International Space Station (ISS) on August 25, 2009. This paper summarizes the three years of observations obtained by the solar neutron detector FIB until the end of July 2012. The solar neutron detector FIB can determine both the energy and arrival direction of neutrons. We measured the energy spectra of background neutrons over the SAA region and elsewhere, and found the typical trigger rates to be 20 counts/sec and 0.22 counts/sec, respectively. It is possible to identify solar neutrons to within a level of 0.028 counts/sec, provided that directional information is applied. Solar neutrons were observed in association with the M-class solar flares that occurred on March 7 (M3.7) and June 7 (M2.5) of 2011. This marked the first time that neutrons were observed in M-class solar flares. A possible interpretaion of the prodcution process is provided.

We report on an 105 ks Suzaku observation of the supernova remnant CTA 1 (G119.5+10.2). The Suzaku soft X-ray observation was carried out with both timing mode and imaging mode. A ~ 10' extended feature, which is interpreted as a bow-shock component of the pulsar wind nebula (PWN), is revealed in this deep observation for the first time. The nebular spectrum can be modelled by a power-law with a photon index of ~ 1.8 which suggests a slow synchrotron cooling scenario. The photon index is approximately constant across this extended feature. We compare and discuss our observations of this complex nebula with previous X-ray investigations. We do not obtain any significant pulsation from the central pulsar in the soft (0.2-12 keV) and hard (10-60 keV) X-ray data. The non-detection is mainly due to the loss of the precise imaging ability to accurately determine the source contribution. The spectra of XIS and HXD can be directly connected without a significant spectral break according to our analysis. Future observations of NuSTAR and Astro-H would be able to resolve the contamination and provide an accurate hard X-ray measurement of CTA 1.

To explain the pulsed emission of the rotation powered pulsars from radio to gamma-ray, the polar cap models, the slot gap models, and the outer gap models are proposed. The recent observations suggest that these models are likely to co-exist in the same magnetosphere. If so, their mutual relation is known to be troublesome (Harding 2009) due to the boundary conditions and the direction of the current which are properly assumed in each acceleration models. We performed a particle simulation for the global magnetospheric structure. Based on the simulation, we present a new picture of the global structure of the pulsar magnetosphere. It is found that a new dead zone is formed along the current neutral line which separates the oppositely directed current. We shall call this the current- neutral zone. We suggest that the polar cap accelerators and the slot gaps locate above the current-neutral zone, and the outer gap exist between the current neutral zone and the traditional dead zone. We also give an estimate of the super-rotation region.

The main issue of the pulsar magnetosphere is how the rotation power is converted into both particle beams which causes pulsed emissions, and a highly relativistic wind of electron-positron plasmas which forms surrounding nebulae shining in X-rays and TeV gamma-rays. As a sequel of the first paper (Wada & Shibata 2007), we carried out a three dimensional particle simulation for the axisymmetric global magnetosphere. We present the results of additional calculations, which are higher resolution model and higher pair creation rate cases, and a detailed analysis for the solution. We confined to demonstrate the cases of low pair creation rate, i.e., the magnetic field is fixed dipole. The radiation drag of the plasma is taken in a form with the curvature radius along the dipole magnetic field. The electrostatic interactions are calculated by a programmable special purpose computer, GRAPE-DR (Makino et al 2007). Once pair creation is onset in the outer gaps, the both signed particles begin to drift across the closed magnetic field due to radiation drag, and they create outflow. Eventually, the steady magnetosphere has outer gaps, both signed outflow of plasma and a region in which the electric field is dominant extending from the equator. In the steady state, the magnetic field made by magnetospheric current is comparable to the dipole magnetic field outside of several light radii from the star. In much more pair creation rate model, the effect of modification of the magnetic field will bring about modification of the outflow of the plasma, requiring further study with higher pair creation rate model in a subsequent paper.

Broad-band (0.8-70 keV) spectra of the persistent X-ray emission from 9 magnetars were obtained with Suzaku, including 3 objects in apparent outburst. The soft X-ray component was detected from all of them, with a typical blackbody temperature of kT ~ 0.5 keV, while the hard-tail component, dominating above ~10 keV, was detected at ~1 mCrab intensity from 7 of them. Therefore, the spectrum composed of a soft emission and a hard-tail component may be considered to be a common property of magnetars, both in their active and quiescent states. Wide-band spectral analyses revealed that the hard-tail component has a 1-60 keV flux, Fh, comparable to or even higher than that carried by the 1-60 keV soft component, Fs. The hardness ratio of these objects, defined as xi=Fh/Fs, was found to be tightly anti-correlated with their characteristic age tau as xi=(3.3+/-0.3)x(tau/1 kyr)^(-0.67+/-0.04) with a correlation coefficient of -0.989, over the range from xi~10 to xi~0.1. Magnetars in outburst states were found to lie on the same correlation as relatively quiescent ones. This hardness ratio is also positively correlated with their surface magnetic fields with a correlation coefficient of 0.873. In addition, the hard-tail component becomes harder towards sources with older characteristic ages, with the photon index changing from ~1.7 to ~0.4.

Recent observations in X-ray and Gamma-ray suggest that the emission region of the pulsar magnetosphere can be multifold. In particular, the open-close boundary of the magnetic field, so-called the Y-point, can be the new candidate place where magnetic field energy converts into the plasma heat and/or flow energy. Here, we present a new Particle-in-Cell code, which can be applied to the Y-point of the pulsar magnetosphere in axisymmetric geometry. The electromagnetic solver is used in the two-dimensional grid points with the cylindrical coordinate (R, z), while the particle solver operates in the three-dimensional Cartesian coordinate (x, y, z), where the Buneman-Boris method is used. The particle motion is treated in special relativity. The inner boundary conditions are set up to generate rotation of the magnetosphere by use of the force-free semi-analytic solution given by Uzdensky (2003, ApJ, 598, 446). The code has been verified by dispersion relations of all the wave modes in electron-positron plasmas. The initial test run is also presented to demonstrate the Y-shaped structure at the top of the dead zone on the light cylinder. We suggest that the structure is variable with quasi-periodicity associated with magnetic reconnection and that plasma will be accelerated and/or heated. In a time-averaged point of view, break up of the ideal-MHD condition takes place in the vicinity of the Y-point.

The results from a systematic study of eleven pulsar wind nebulae with a torus structure observed with the Chandra X-ray observatory are presented. A significant observational correlation is found between the radius of the tori, r, and the spin-down luminosity of the pulsars, Edot. A logarithmic linear fit between the two parameters yields log r = (0.57 +- 0.22) log Edot -22.3 +- 8.0 with a correlation coefficient of 0.82, where the units of r and Edot are pc and ergs s^-1, respectively. The value obtained for the Edot dependency of r is consistent with a square root law, which is theoretically expected. This is the first observational evidence of this dependency, and provides a useful tool to estimate the spin-down energies of pulsars without direct detections of pulsation. Applications of this dependency to some other samples are also shown.

A Chandra study of pulsar wind nebula around the young energetic pulsar PSR B1509-58 is presented. The high resolution X-ray image with total exposure time of 190 ks reveals a ring like feature 10'' apart from the pulsar. This feature is analogous to the inner ring seen in the Crab nebula and thus may correspond to a wind termination shock. The shock radius enables us to constrain the wind magnetization, sigma>= 0.01. The obtained sigma is one order of magnitude larger than that of the Crab nebula. In the pulsar vicinity, the southern jet appears to extend beyond the wind termination shock, in contrast to the narrow jet of the Crab. The revealed morphology of the broad jet is coincident with the recently proposed theoretical model in which a magnetic hoop stress diverts and squeezes the post-shock equatorial flow towards the poloidal direction generating a jet.

A two-dimensional electrodynamic model is used to study particle acceleration and non-thermal emission mechanisms in the pulsar magnetospheres. We solve distribution of the accelerating electric field with the emission process and the pair-creation process in meridional plane, which includes the rotational axis and the magnetic axis. By solving the evolutions of the Lorentz factor, and of the pitch angle, we calculate spectrum in optical through $\gamma$-ray bands with the curvature radiation, synchrotron radiation, and inverse-Compton process not only for outgoing particles, but also for ingoing particles, which were ignored in previous studies. We apply the theory to the Vela pulsar. We find that the curvature radiation from the outgoing particles is the major emission process above 10 MeV bands. In soft $\gamma$-ray to hard X-ray bands, the synchrotron radiation from the ingoing primary particles in the gap dominates in the spectrum. Below hard X-ray bands, the synchrotron emissions from both outgoing and ingoing particles contribute to the calculated spectrum. The calculated spectrum is consistent with the observed phase-averaged spectrum of the Vela pulsar. We show that the observed five-peak pulse profile in the X-ray bands of the Vela pulsar is reproduced by the inward and outward emissions, and the observed double-peak pulse profile in $\gamma$-ray bands is explained by the outward emissions.

In this paper, we present an expanding disc model to derive polarization properties of the Crab nebula. The distribution function of the plasma and the energy density of the magnetic field are prescribed as function of the distance from the pulsar by using the model by Kennel and Coroniti (1984) with $\sigma = 0.003$, where $\sigma$ is the ratio of Poynting flux to the kinetic energy flux in the bulk motion just before the termination shock. Unlike previous models, we introduce disordered magnetic field, which is parameterized by the fractional energy density of the disordered component. Flow dynamics is not solved. The mean field is toroidal. Averaged polarization degree over the disc is obtained as a function of inclination angle and fractional energy density of the disordered magnetic field. It is found for the Crab that the disordered component has about 60 percent of the magnetic field energy. This value is also supported by the facts that the disc appears not `lip-shape' but as `rings' in the intensity map as was observed, and that the highest polarization degree of $\sim 40$ percent is reproduced for rings, being consistent with the observation. We suggest that because the disordered field contributes rather pressure than tension, the pinch force may be over-estimated in previous relativistic magnetohydrodynamic simulations. Disruption of the toroidal magnetic field with alternating direction, which is proposed by Lyubarski (2003), may actually takes place. The relativistic flow speed, which is indicated by the front-back contrast, can be detected in asymmetry in distributions of the position angle and depolarization.

Soon after the discovery of radio pulsars in 1967, the pulsars are identified as strongly magnetic (typically $10^{12}$G) rapidly rotating ($\sim 10^{2}-0.1$ Hz) neutron stars. However, the mechanism of particle acceleration in the pulsar magnetosphere has been a longstanding problem. The central problem is why the rotation power manifests itself in both gamma-ray beams and a highly relativistic wind of electron-positron plasmas, which excites surrounding nebulae observed in X-ray. Here we show with a three dimensional particle simulation for the global axisymmetric magnetosphere that a steady outflow of electron-positron pairs is formed with associated pair sources, which are the gamma-ray emitting regions within the light cylinder. The magnetic field is assumed to be dipole, and to be consistent, pair creation rate is taken to be small, so that the model might be applicable to old pulsars such as Geminga. The pair sources are charge-deficient regions around the null surface, and we identify them as the outer gap. The wind mechanism is the electromagnetic induction which brings about fast azimuthal motion and eventually trans-field drift by radiation drag in close vicinity of the light cylinder and beyond. The wind causes loss of particles from the system. This maintains charge deficiency in the outer gap and pair creation. The model is thus in a steady state, balancing loss and supply of particles. Our simulation implies how the wind coexists with the gamma-ray emitting regions in the pulsar magnetosphere.

We report here on the first pointed X-ray observation of the supernova remnant (SNR) G330.2+1.0. The X-ray morphology is characterized by an extended shell. Its X-ray spectrum is well represented by a single power-law function with a photon index of $\gamma\simeq 2.8$ and interstellar absorption of $n_{\rm H}\simeq2.6\times 10^{22}$[cm$^{-2}$]. We interpret this emission as synchrotron radiation from accelerated electrons at the SNR shock, as seen in SN 1006. The surface brightness of the X-ray emission is anti-correlated with the radio emission, and the power-law spectrum is dominated at the western shell where the radio emission is weak. The co-existence of two distinct (radio bright/X-ray faint and radio faint/X-ray bright) shells in a single supernova remnant challenges our understanding of the particle acceleration and radiation mechanisms in different interstellar environments. The object may be a good target for searching TeV gamma-rays and molecular gas surrounding the blast shock. We also report on the nature of a bright point-like source (AX J1601-5143) to the south of the SNR.

A two-dimensional electrodynamical model is used to study particle acceleration in the outer magnetosphere of a pulsar. The charge depletion from the Goldreich-Julian charge density causes a large electric field along the magnetic field lines. The charge particles are accelerated by the electric field and emit $\gamma$-rays via the curvature process. Some of the emitted $\gamma$-rays may collide with $X$-ray photons to make new pairs, which are accelerated again on the different field lines in the gap and proceed similar processes. We simulate the pair creation cascade in the meridional plane using the pair creation mean-free path, in which the $X$-ray photon number density is proportional to inverse square of radial distance. With the space charge density determined by the pair creation simulation, we solve the electric structure of the outer gap in the meridional plane and calculate the curvature spectrum. Because the two-dimensional model can link both gap width along the magnetic field line and trans-field thickness with the spectral cut-off energy and flux, we can diagnose both the current through the gap and inclination angle between the rotational and magnetic axes. We apply the theory to the Vela pulsar. By comparing the results with the $EGRET$ data, we rule out any cases that have a large particle injection at the outer boundary. We also suggest the inclination angle of $\alpha_{inc}\geq65^{\circ}$. The present model predicts the outer gap starting from near the conventional null charge surface for the Vela pulsar.

We carried out a deep subarcsecond BRI imaging of the two middle-aged pulsars to establish their properties in the optical range. Both pulsars are detected at >10 sigma level. Geminga is for the first time reliably detected in the I band with a magnitude of 25.10+/-0.14. We also reanalyze archival ESO/NTT and HST broadband data and find that some published fluxes for Geminga were estimated inaccurately. The resulting dereddened broadband spectra of both pulsars are remarkably similar to each other and show significant flux increases towards the far-UV and near-IR, and a wide flux excess in V-I bands. This suggests a multicomponent structure of the optical emission. The nonthermal power law component of the pulsar magnetospheric origin dominates in the most part of the optical range. For PSR B0656+14 it is compatible with a low energy extension of the power law tail seen in hard X-rays. For Geminga the respective extension overshoots by a factor of 100 the nonthermal optical flux, which has a less steep spectral slope than in X-rays. This implies a spectral break at a photon energy of about 1 keV. The flux increases towards the far-UV are compatible with contributions of the Rayleigh-Jeans parts of the blackbody components from whole surfaces of the neutron stars dominating in soft X-rays. The V-I excess, which is most significant for PSR B0656+14, suggests a third spectral component of still unidentified origin. Faint, a few arcseconds in size nebulae extended perpendicular to the proper motion directions of the pulsars, are seen around both objects in our deepest I band images. They can be optical counterparts of the bow-shock head of Geminga and of the tentative pulsar wind nebula of PSR B0656+14 observed in X-rays.

During the period when the Sun was intensely active on October-November 2003, two remarkable solar neutron events were observed by the ground-based neutron monitors. On October 28, 2003, in association with an X17.2 large flare, solar neutrons were detected with high statistical significance (6.4 sigma) by the neutron monitor at Tsumeb, Namibia. On November 4, 2003, in association with an X28 class flare, relativistic solar neutrons were observed by the neutron monitors at Haleakala in Hawaii and Mexico City, and by the solar neutron telescope at Mauna Kea in Hawaii simultaneously. Clear excesses were observed at the same time by these detectors, with the significance calculated as 7.5 sigma for Haleakala, and 5.2 sigma for Mexico City. The detector onboard the INTEGRAL satellite observed a high flux of hard X-rays and gamma-rays at the same time in these events. By using the time profiles of the gamma-ray lines, we can explain the time profile of the neutron monitor. It appears that neutrons were produced at the same time as the gamma-ray emission.

We investigate the electrodynamics of an outer gap in the meridional plane of the aligned-rotator. The charge depletion from the Goldreich-Julian charge density causes a large electric field along the magnetic field line. The electrons or the positrons are accelerated by the field-aligned electric field and radiate the $\gamma$-rays tangentially to the local magnetic field line. Some of such $\gamma$-rays collide with $X$-rays to materialize as the electron-positron pairs on different field lines from the field line on which they were emitted. As a result, the electric field structure is expected to change across the field lines. Including these trans-field effects, we solve the formation of the electric field self-consistently with the curvature radiation and the pair creation processes. The $\gamma$-ray emission and the pair creation are treated by use of Monte Carlo technique. We demonstrate that the distribution of the electric field along the field lines is affected by both the gap geometry and the external currents coming into the gap through the boundaries. In the electrodynamical model, it has been known that the solution disappears if the current density carried by the electron-positron pairs produced in the gap exceeds a critical value. We show that the critical current density is significantly increased when the trans-field structure is taken into account. We also find that the location of the inner boundary of the gap shifts toward the stellar surface from the conventional null surface as the current density increases. The reason for the shift is derived from the stability condition of the inner boundary. We also argue that the ideal-MHD condition holds outside of the gap only when the low energy particles coexist with the high energy particles migrating from the gap.

We present spectral analysis of the Crab Nebula obtained with the {\it Chandra} X-ray observatory. The X-ray spectrum is characterized by a power-law whose index varies across the nebula. The variation can be discussed in terms of the particle injection from the pulsar in two different directions: the equatorial plane containing the torus and the symmetry axis along the jet. In the equatorial plane, spectra within the torus are the hardest, with photon index $\alpha \approx 1.9$, and are almost independent of the surface brightness. At the periphery of the torus, the spectrum gradually softens in the outer, lower surface brightness regions, up to $\alpha \approx 3.0$. This indicates that synchrotron losses become significant to X-ray emitting particles at the outer boundary of the torus. We discuss the nature of the torus, incorporating information from observations at other wavelengths. Spectral variations are also seen within the southern jet. The core of the jet is the hardest with $\alpha \approx$ 2.0, and the outer sheath surrounding the core becomes softer with $\alpha$ up to 2.5 at the outermost part. Based on the similarity between the spectra of the jet core and the torus, we suggest that the electron spectra of the particles injected from the pulsar are also similar in these two different directions. The brightness ratio between the near and far sides of the torus can be explained by Doppler boosting and relativistic aberration; however, the observed ratio cannot be derived from the standard weakly magnetized pulsar wind model. We also found a site where an optical filament comprised of supernova ejecta is absorbing the soft X-ray emission ($<$ 2 keV).

We investigate the $\gamma$-ray emission from a outer-gap accelerator, which is located in the outer-magnetosphere of a pulsar. The charge depletion from the Goldreich-Julian density causes a large electric field along the magnetic field lines. The electric field accelerates electrons and positrons to high energies. We solve the formation of the electric field self-consistently with the radiation and the pair-creation processes in one-dimensional geometry along the magnetic field lines, and calculate curvature spectrum to compare with the observations. We find that because the particles escape from the accelerating region with large Lorentz factors as $10^{7.5}$, curvature radiation from the outside of the gap, in which there is no electric field, still has a important contribution to the spectrum. Including the outside emission, the model spectra are improved between 100 MeV and 1 GeV for the Vela-type pulsars. Since the field line curvature close to the light cylinder affects the spectrum in several hundred MeV, there is a possibility to diagnose the field structure by using the spectrum in these bands. The difference found in the $\gamma$-ray spectra of the twin pulsars, PSR B0833-45 (Vela) and B1706-44, is deduced from the dimensionless electric current running through the gap. The extension to the three-dimensional modeling with self-consistent electrodynamics is necessary to understand the difference in the observed flux between the twin pulsars.

We present year-scale morphological variations of the Crab Nebula revealed by Chandra X-ray observatory. Observations have been performed about every 1.7 years over the 3 years from its launch. The variations are clearly recognized at two sites: the torus and the southern jet. The torus, which had been steadily expanding until 1.7 years ago, now appears to have shrunk in the latest observation. Additionally, the circular structures seen to the northeast of the torus have decayed into several arcs. On the other hand, the southern jet shows the growth of its overall kinked-structure. We discuss the nature of these variations in terms of the pulsar wind nebula (PWN) mechanism.

An axisymmetric model for the Crab Nebula is constructed to examine the flow dynamics in the nebula. The model is based on that of Kennel and Coroniti (1984), although we assume that the kinetic-energy-dominant wind is confined in an equatorial region. The evolution of the distribution function of the electron-positron plasma flowing out in the nebula is calculated. Given viewing angles, we reproduce an image of the nebula and compare it with Chandra observation. The reproduced image is not a ring-like but rather 'lip-shaped'. It is found that the assumption of toroidal field does not reproduce the Chandra image. We must assume that there is disordered magnetic field with an amplitude as large as the mean toroidal field. In addition, the brightness contrast between the front and back sides of the ring cannot be reproduced if we assume that the magnetization parameter $\sigma$ is as small as $\sim 10^{-3}$. The brightness profile along the semi-major axis of the torus is also examined. The non-dissipative, ideal-MHD approximation in the nebula appears to break down. We speculate that if the magnetic energy is released by some process that produce turbulent field in the nebula flow and causes heating and acceleration, e.g. by magnetic reconnection, then the present difficulties may be resolved (i.e. we can reproduce a ring image, and a higher brightness contrast). Thus, the magnetization parameter $\sigma$ can be larger than previously expected.

We investigate a stationary pair production cascade in the outer magnetosphere of a spinning neutron star. The charge depletion due to global flows of charged particles, causes a large electric field along the magnetic field lines. Migratory electrons and/or positrons are accelerated by this field to radiate gamma-rays via curvature and inverse-Compton processes. Some of such gamma-rays collide with the X-rays to materialize as pairs in the gap. The replenished charges partially screen the electric field, which is self-consistently solved together with the energy distribution of particles and gamma-rays at each point along the field lines. By solving the set of Maxwell and Boltzmann equations, we demonstrate that an external injection of charged particles at nearly Goldreich-Julian rate does not quench the gap but shifts its position and that the particle energy distribution cannot be described by a power-law. The injected particles are accelerated in the gap and escape from it with large Lorentz factors. We show that such escaping particles migrating outside of the gap contribute significantly to the gamma-ray luminosity for young pulsars and that the soft gamma-ray spectrum between 100 MeV and 3 GeV observed for the Vela pulsar can be explained by this component. We also discuss that the luminosity of the gamma-rays emitted by the escaping particles is naturally proportional to the square root of the spin-down luminosity.

We develop a model for gamma-ray emission from the outer magnetosphere of pulsars (the outer-gap model). The charge depletion causes a large electric field which accelerates electrons and positrons. We solve the electric field with radiation and pair creation processes self-consistently, and calculate curvature spectrum and Inverse-Compton (IC) spectrum. We apply this theory to PSR B0833-45 (Vela) and B1706-44 for which their surface magnetic fields, observed thermal X-rays are similar to each other. We find that each observed cut-off energies of the gamma-rays are well explained. By inclusion of emission outside the gap, the spectrum is in better agreement with the observations than the spectrum arising only from the inside of the gap. The expected TeV fluxes are much smaller than that observed by CANGAROO group in the direction of B1706-44.

We construct an axisymmetric model for the Crab Nebula. The flow dyamics is based on the model by Kennel and Coroniti (1984), but we assume that the kinetic-energy-dominant wind is confined in the equatorial region. We reproduce spacially resolved spectra which agree well with the Chandra results. However, the reproduced image is not a ring but lip-shaped. In addition, brightness contrast between fore and back sides of the ring cannot be reproduced if we assumes that the sigma-parameter is as small as 10**-3. We suggest that the nebula field is highly disordered due to for instance magnetic reconnection. The estimate of sigma can be larger than previously expected.

We reconstruct a 3+1 formalism of general relativistic electromagnetism, and derive the equations of motion of charged particles in the pulsar magnetosphere, taking account of the inclination between the rotation axis and the magnetic axis. Apart from the previous works where space charge is evaluated by assuming the flow velocity being the light speed, we analyze particle motion in the polar cap, finding that gravity changes significantly its dynamics and the condition for acceleration.

We study self-consistency for the pulsar polar cap model. In general, there will remain an unscreened electric field at the pair production front (PPF). We have derived the condition of electric field screening by pairs in the presence of returning particles. A previous belief that pairs with a density higher than the Goldreich-Julian density immediately screen out the electric field beyond PPF is unjustified. Pairs have little contribution on screening the electric field, and the geometrical screening is the only known way of getting a finite potential drop.

We study on the self-consistency of the pulsar polar cap model, i.e., the problem of whether the field-aligned electric field is screened by electron-positron pairs that are injected beyond the pair production front. We solve the one-dimensional Poisson equation along a magnetic field line, both analytically and numerically, for a given current density incorporating effects of returning positrons, and we obtain the conditions for the electric-field screening. The formula which we obtained gives the screening distance and the return flux for given primary current density, field geometry and pair creation rate at the pair production front. If the geometrical screening is not possible, for instance, on field lines with a super-Goldreich-Julian current, then the electric field at the pair production front is constrained to be fairly small in comparison with values expected typically by the conventional polar cap models. This is because (1) positive space charge by pair polarization is limited to a small value, and (2) returning of positrons leave pair electrons behind. A previous belief that pair creation with a pair density higher than the Goldreich-Julian density immediately screens out the electric field is unjustified at least for for a super Goldreich-Julian current density. We suggest some possibilities to resolve this difficulty.

This paper deals with the Crab Nebula problem to suggest that particle acceleration takes place not only at the inner shock but also over a larger region in the nebula with disordered magnetic field. Kennel and Cornoniti (1984) constructed a spherically symmetric model of the Crab Nebula and concluded that the pulsar wind which excites the nebular is kinetic-energy dominant (KED) because the nebula flow induced by KED wind is favorable to explain the nebula spectrum and expansion speed. This is true even with new Chandra observation, which provides newly the spatially resolved spectra. We have shown below with 3D modelling and the Chandra image that pure toroidal magnetic field and KED wind are incompatible with the Chandra observation.

We developed a magneto-turbulent model for the cosmic ray (CR) electrons seen in the radio halo clusters of galaxies. Steady state momentum distribution functions of the CR electrons are calculated for given spectra of the turbulent Alfven waves. The radio spectrum produced by the obtained CR electron distribution is compared to the observed radio spectrum of the Coma radio halo. We find that the observed radio spectrum of the Coma cluster is well reproduced when the spectral index of the turbulent Alfven waves is 2.8. The obtained energy spectrum of the turbulent Alfven waves is steeper than that expected from the turbulence theory, suggesting back reaction of the particle acceleration. The fitting procedure constraints the amplitude of the turbulent Alfven waves as well as the spectral index. Then we estimate the dissipation of the turbulent Alfven waves, which is found to be less than the cooling rate by X-ray radiation. We suggest that the turbulence which is sufficient for particle acceleration is developed in the clusters containing the radio halo. It is most likely that cluster mergers create the turbulence and seed relativistic electrons.

We report the B band optical observations of an old (17.5 Myr) radiopulsar PSR B0950+08 obtained with the Suprime-Cam at the Subaru telescope. We detected a faint object, B=27.07(16). Within our astrometrical accuracy it coincides with the radio position of the pulsar and with the object detected earlier by Pavlov et al. (1996) in UV with the HST/FOC/F130LP. The positional coincidence and spectral properties of the object suggest that it is the optical counterpart of PSR B0950+08. Its flux in the B band is two times higher than one would expect from the suggested earlier Rayleigh-Jeans interpretation of the only available HST observations in the adjacent F130LP band. Based on the B and F130LP photometry of the suggested counterpart and on the available X-ray data we argue in favour of nonthermal origin of the broad-band optical spectrum of PSR B0950+08, as it is observed for the optical emission of the younger, middle-aged pulsars PSR B0656+14 and Geminga. At the same time, the optical efficiency of PSR B0950+08, estimated from its spin-down power and the detected optical flux, is by several orders of magnitude higher than for these pulsars, and comparable with that for the much younger and more energetic Crab pulsar. We cannot exclude the presence of a compact, about 1'', faint pulsar nebula around PSR B0950+08, elongated perpendicular to the vector of its proper motion, unless it is not a projection of a faint extended object on the pulsar position.

We solve the one-dimensional Poisson equation along a magnetic field line, both analytically and numerically, for a given current density incorporating effects of returning positrons. We find that the number of returning positrons per one primary electrons should be smaller than unity, and the returning of positrons occurs only in a very short braking distance scale. As a result, for realistic polar cap parameters, the accelerating electric field will not be screened out; thus, the model fails to be self-consistent. A previous belief that pair creation with a pair density higher than the Goldreich-Julian density immediately screens out the electric field is unjustified. We suggest some possibilities to resolve this difficulty.

We investigate a stationary pair production cascade in the outer magnetosphere of a spinning neutron star. The charge depletion due to global flows of charged particles, causes a large electric field along the magnetic field lines. Migratory electrons and/or positrons are accelerated by this field to radiate curvature gamma-rays, some of which collide with the X-rays to materialize as pairs in the gap. The replenished charges partially screen the electric field, which is self-consistently solved, together with the distribution functions of particles and gamma-rays. By solving the Vlasov equations describing this pair production cascade, we demonstrate the existence of a stationary gap in the outer magnetosphere of PSR B1055-52 for a wide range of current densities flowing in the accelerator: From sub to super Goldreich-Julian values. However, we find that the expected GeV spectrum becomes very soft if the current density exceeds the Goldreich-Julian value. We also demonstrate that the GeV spectrum softens with decreasing magnetic inclination and with increasing distance to this pulsar. We thus conclude that a sub-Goldreich-Julian current, a large magnetic inclination, and a small distance (500 pc, say) are plausible to account for EGRET observations. Furthermore, it is found that the TeV flux due to inverse Compton scatterings of infrared photons whose spectrum is inferred from the Rayleigh-Jeans side of the soft blackbody spectrum is much less than the observational upper limit.

We investigate a stationary pair production cascade in the outer magnetosphere of a spinning neutron star. The charge depletion due to global flows of charged particles, causes a large electric field along the magnetic field lines. Migratory electrons and/or positrons are accelerated by this field to radiate curvature gamma-rays, some of which collide with the X-rays to materialize as pairs in the gap. The replenished charges partially screen the electric field, which is self-consistently solved, together with the distribution functions of particles and gamma-rays. By solving the Vlasov equations describing this pair production cascade, we demonstrate the existence of a stationary gap in the outer magnetosphere of PSR B1055-52 for a wide range of current densities flowing in the accelerator: From sub to super Goldreich-Julian values. However, we find that the expected GeV spectrum becomes very soft if the current density exceeds the Goldreich-Julian value. We also demonstrate that the GeV spectrum softens with decreasing magnetic inclination and with increasing distance to this pulsar. We thus conclude that a sub-Goldreich-Julian current, a large magnetic inclination, and a small distance (500 pc, say) are plausible to account for EGRET observations. Furthermore, it is found that the TeV flux due to inverse Compton scatterings of infrared photons whose spectrum is inferred from the Rayleigh-Jeans side of the soft blackbody spectrum is much less than the observational upper limit.

ASCA's long look at the 288 millisecond radio pulsar, PSR J0631+1036, reveals coherent X-ray pulsation from this source for the first time. The source was first detected in the serendipitous Einstein observation and later identified as a radio pulsar. Possible pulsation in the gamma-ray band has been detected from the CGRO EGRET data (Zepka, et al. 1996). The X-ray spectrum in the ASCA band is characterized by a hard power-law type emission with a photon index of about 2.3, when fitted with a single power-law function modified with absorption. An additional blackbody component of about 0.14 keV increases the quality of the spectral fit. The observed X-ray flux is 2.1e-13 ergs/s/cm2 in the 1-10 keV band. We find that many characteristics of PSR J0631+1036 are similar to those of middle-aged gamma-ray pulsars such as PSR B1055-52, PSR B0633+17 (Geminga), and PSR B0656+14.

We investigate a stationary pair production cascade in the outer magnetosphere of a spinning neutron star. The charge depletion due to global flows of charged particles, causes a large electric field along the magnetic field lines. Migratory electrons and/or positrons are accelerated by this field to radiate curvature gamma-rays, some of which collide with the X-rays to materialize as pairs in the gap. The replenished charges partially screen the electric field, which is self-consistently solved together with the distribution functions of particles and gamma-rays. If no current is injected at neither of the boundaries of the accelerator, the gap is located around the conventional null surface, where the local Goldreich-Julian charge density vanishes. However, we first find that the gap position shifts outwards (or inwards) when particles are injected at the inner (or outer) boundary. Applying the theory to the Crab pulsar, we demonstrate that the pulsed TeV flux does not exceed the observational upper limit for moderate infrared photon density and that the gap should be located near to or outside of the conventional null surface so that the observed spectrum of pulsed GeV fluxes may be emitted via a curvature process. Some implications of the existence of a solution for a super Goldreich-Julian current are discussed.

We investigate a stationary particle acceleration zone in the outer magnetosphere of a spinning neutron star. The charge depletion due to a global current causes a large electric field along the magnetic field lines. Migratory electrons and/or positrons are accelerated by this field to radiate curvature gamma-rays, some of which collide with the X-rays to materialize as pairs in the gap. As a result of this pair-production cascade, the replenished charges partially screen the electric field, which is self-consistently solved together with the distribution of particles and gamma-rays. If no current is injected at neither of the boundaries of the accelerator, the gap is located around the so-called null surface, where the local Goldreich-Julian charge density vanishes. However, we find that the gap position shifts outwards (or inwards) when particles are injected at the inner (or outer) boundary. We apply the theory to the nine pulsars of which X-ray fields are known from observations. We show that the gap should be located near to or outside of the null surface for the Vela pulsar and PSR B1951+32, so that their expected GeV spectrum may be consistent with observations. We then demonstrate that the intrinsically large TeV flux from the outer gap of PSR B0540-69 is absorbed by the magnetospheric infrared photons to be undetectable. We also point out that the electrodynamic structure and the resultant GeV emission properties of millisecond pulsars are similar to young pulsars.

We investigate a stationary pair production cascade in the outer magnetosphere of a spinning neutron star. The charge depletion due to a global current, causes a large electric field along the magnetic field lines. Migratory electrons and/or positrons are accelerated by this field to radiate curvature gamma-rays, some of which collide with the X-rays to materialize as pairs in the gap. The replenished charges partially screen the electric field, which is self-consistently solved together with the distribution functions of particles and gamma-rays. If no current is injected at neither of the boundaries of the accelerator, the gap is located around the so-called null surface, where the local Goldreich-Julian charge density vanishes. However, we first find that the gap position shifts outwards (or inwards) when particles are injected at the inner (or outer) boundary. Applying the theory to the Crab pulsar, we demonstrate that the pulsed TeV flux does not exceed the observational upper limit for moderate infrared photon density and that the gap should be located near to or outside of the null surface so that the observed spectrum of pulsed GeV fluxes may be emitted via curvature process.

We present the results of a BeppoSAX observation of PSR J0218+4232 which, for the first time, provides detailed information on the pulsar's temporal and spectral emission properties in the broad band 1-10 keV. We detected X-ray pulses with a pulsed fraction of 73 +/- 12 %. The pulse profile is characterized by two peaks phase separated by 0.47 +/-0.05. The pulsed spectrum is best described by a power-law of photon index 0.61 +/- 0.32 with an unabsorbed (2-10 keV) X-ray flux of 4.1E-13 erg/cm^2/s implying a luminosity of Lx = 1.3E32 Theta (d/5.7 kpc})^2 erg/s and an X-ray efficiency of 4.8E-4 Theta (d/5.7 kpc})^2 where Theta is the solid angle spanned by the emission beam.

First I emphasize the magnetospheric loop current system induced by the wind activity. Next, I address the formation of the field-aligned electric field. An interesting aspect is the reaction of the inner magnetosphere to an induced current. In the traditional outer gap, the current running through the gap is all carried by the particles created in the gap. However, if external current-carrying particles are added, the gap moves inward or outward; it can even appear near the star. It is suggested that the gaps appear at various altitudes, depending on the current distribution on field lines, the field line geometry, and the source of the current-carrying particles. For definiteness of sign of charge and current direction, $\vec{\Omega} \cdot \vec{\mu} >0$ is assumed throughout.

Several weeks after the explosion of supernova (SN) 1987A, the UV flash of the SN illuminated a ring-like structure in the circumstellar material at about 0.65 ly from the SN. The interaction between the stellar winds from the SN progenitor is considered to be the candidate for the formation of the circumstellar structure. In the case that the stellar winds are spherically symmetric, the interaction should result in a shell-like structure. However, in this paper we show that the magnetic field in the winds causes an anisotropy which leads to the formation of a ring-like structure. When the fast wind of the blue supergiant phase of the progenitor sweeps up the surrounding slow wind of the red-supergiant phase, the magnetic field as well as the wind material are piled up in the interaction region. Since the magnetic energy increases in proportion to the square of the amplitude, the magnetic field exhibits its effect prominently at the interaction region; due to the magnetic pressure force the material at lower latitudes is compressed into a ring-like structure. It is suggested that this magnetic process can also explain the newly observed pair of rings of the SN 1987A nebula.