The Suzaku data of the highly variable magnetar 1E 1547.0$-$5408, obtained during the 2009 January activity, were reanalyzed. The 2.07 s pulsation of the 15--40 keV emission detected with the HXD was found to be phase modulated, with a period of $36.0^{+4.5}_{-2.5}$ ks and an amplitude of $0.52 \pm 0.14$ s. The modulation waveform is suggested to be more square-wave like rather than sinusoidal. While the effect was confirmed with the 10--14 keV XIS data, the modulation amplitude decreased towards lower energies, becoming consistent with 0 below 4 keV. After the case of 4U 0142+61, this makes the 2nd example of this kind of behavior detected from magnetars. The effect can be interpreted as a manifestation of free precession of this magnetar, which is suggested to be oblately deformed under the presence of strong toroidal field of $\sim 10^{16}$ G.

Active shielding is an effective technique to reduce background signals in hard X-ray detectors and to enable observing darker sources with high sensitivity in space. Usually the main detector is covered with some shield detectors made of scintillator crystals such as BGO (Bi$_4$Ge$_3$O$_{12}$), and the background signals are filtered out using anti-coincidence among them. Japanese X-ray observing satellites "Suzaku" and "ASTRO-H" employed this technique in their hard X-ray instruments observing at > 10 keV. In the next generation X-ray satellites, such as the NGHXT proposal, a single hybrid detector is expected to cover both soft (1-10 keV) and hard (> 10 keV) X-rays for effectiveness. However, present active shielding is not optimized for the soft X-ray band, 1-10 keV. For example, Bi and Ge, which are contained in BGO, have their fluorescence emission lines around 10 keV. These lines appear in the background spectra obtained by ASTRO-H Hard X-ray Imager, which are non-negligible in its observation energy band of 5-80 keV. We are now optimizing the design of active shields for both soft and hard X-rays at the same time. As a first step, we utilized a BGO crystal as a default material, and measured the L lines of Bi and K lines of Ge from it using the X-ray SOIPIX, "XRPIX".

As the endpoints of massive star evolution, neutron stars are enigmatic celestial objects characterized by extremely dense and exotic nuclear matter, magnetospheres with positrons (antimatter), rapid rotation and ultra-strong magnetic fields. Such an extreme environment has provided an accessible astrophysical laboratory to study physics under conditions unattainable on Earth and to tackle a range of fundamental questions related to: the aftermath of stellar evolution and the powerful explosions of massive stars, the equation of state and physics of some of the most exotic and magnetic stars in the Universe, the workings of the most powerful particle accelerators in our Galaxy and beyond, and the sources of gravitational waves that are yet to be detected. Recent observations revealed a great diversity of neutron stars, including ultra-strongly magnetized pulsars, referred to as "magnetars", and unusual types of accreting X-ray pulsars. In this white paper, we highlight the prospects of the upcoming X-ray mission, ASTRO-H, in studying these highly magnetized neutron stars.

The dipping low-mass X-ray binary 4U 1915-05 was observed by Suzaku on 2007 November 8 for a net exposure of 39 ksec. It was detected by the XIS with a 0.8-10 keV signal rate of 9.84+\-0.01 cts/s per camera, and HXD-PIN with a 12-45 keV signal rate of 0.29+/-0.01 cts/s. After removing the periodic dips and an X-ray burst, the 0.8 - 45 keV continuum was successfully described by an optically thick disk emission with an inner-disk temperature ~ 0.7 keV and a neutron-star blackbody emission with a temperature ~ 1.3 keV, on condition that the blackbody component, or possibly the disk emission too, is significantly Comptonized. This successful modeling is consistent with 4U 1915-05 being in a high-soft state in this observation, and implies that its broadband spectrum can be interpreted in the same scheme as for many non-dipping Low-mass X-ray binaries in the soft state. Its bolometric luminosity (~ 0.02 times the Eddington limit) is relatively low for the soft state, but within a tolerance, if considering the distance and inclination uncertainties. As a high-inclination binary, this source exhibited stronger Comptonization effect, with a larger Comptonizing y-parameter, compared to low and medium inclination binaries. This suggests that the Comptonizing coronae of these objects in the soft state is in an oblate (rather than spherical) shape, extending along the accretion disk plane, because the y-parameter would not depend on the inclination if the corona were spherical.

Magnetars are a special type of neutron stars, considered to have extreme dipole magnetic fields reaching ~1e+11 T. The magnetar 4U 0142+61, one of prototypes of this class, was studied in broadband X-rays (0.5-70 keV) with the Suzaku observatory. In hard X-rays (15-40 keV), its 8.69 sec pulsations suffered slow phase modulations by +/-0.7 sec, with a period of ~15 hours. When this effect is interpreted as free precession of the neutron star, the object is inferred to deviate from spherical symmetry by ~1.6e-4 in its moments of inertia. This deformation, when ascribed to magnetic pressure, suggests a strong toroidal magnetic field, ~1e+12 T, residing inside the object. This provides one of the first observational approaches towards toroidal magnetic fields of magnetars.

The symbiotic X-ray binary 4U 1954+319 is a rare system hosting a peculiar neutron star (NS) and an M-type optical companion. Its ~5.4h NS spin period is the longest among all known accretion-powered pulsars and exhibited large (~7%) fluctuations over 8 years. A spin trend transition was detected with Swift/BAT around an X-ray brightening in 2012. The source was in quiescent and bright states before and after this outburst based on 60 ks Suzaku observations in 2011 and 2012. The observed continuum is well described by a Comptonized model with the addition of a narrow 6.4 keV Fe Kalpha line during the outburst. Spectral similarities to slowly rotating pulsars in high-mass X-ray binaries, its high pulsed fraction (~60-80%), and the location in the Corbet diagram favor high B-field (>~1e+12 G) over a weak field as in low-mass X-ray binaries. The observed low X-ray luminosity (1e+33-1e+35 erg/s), probable wide orbit, and a slow stellar wind of this SyXB make quasi-spherical accretion in the subsonic settling regime a plausible model. Assuming a ~1e+13 G NS, this scheme can explain the ~5.4 h equilibrium rotation without employing the magnetar-like field (~1e+16 G) required in the disk accretion case. The time-scales of multiple irregular flares (~50 s) can also be attributed to the free-fall time from the Alfven shell for a ~1e+13 G field. A physical interpretation of SyXBs beyond the canonical binary classifications is discussed.

We report on the firm detection of a cyclotron resonance scattering feature (CRSF) in the X-ray spectrum of the Be X-ray binary pulsar, GRO J1008-57, achieved by the Suzaku Hard X-ray Detector during a giant outburst which was detected by the MAXI Gas Slit Camera in 2012 November. The Suzaku observation was carried out on 2012 November 20, outburst maximum when the X-ray flux reached $\sim 0.45$ Crab in 4-10 keV, which corresponds to a luminosity of $1.1 \times 10^{38}$ erg s$^{-1}$ in 0.5--100 keV at 5.8 kpc. The obtained broadband X-ray spectrum from 0.5 keV to 118 keV revealed a significant absorption feature, considered as the fundamental CRSF, at $\sim 76$ keV. This unambiguously reconfirm the previously suggested $\sim$ 80 keV spectral feature in GRO J1008$-$57. The implied surface magnetic field, $6.6\times 10^{12}$ G, is the highest among binary X-ray pulsars from which CRSFs have ever been detected.

The dipping X-ray source 4U 1822$-$37 was observed by Suzaku on 2006 Octrober 20 for a net exposure of 37 ks. The source was detected with the XIS at a 1-10 keV flux of 5.5$\times10^{-10}$ erg cm$^{-2}$ s$^{-1}$, and with the HXD (HXD-PIN) at a 10-50 keV flux of 8.9$\times10^{-10}$ erg cm$^{-2}$ s$^{-1}$. With HXD-PIN, the pulsation was detected at a barycentric period of 0.592437 s, and its change rate was reconfirmed as $-2.43\times$10$^{-12}$ s s$^{-1}$. The 1-50 keV spectra of 4U 1822-37 were found to be very similar to those of Her X-1 in the slopes, cutoff and iron lines. Three iron lines (Fe K$\alpha$, Fe XXV, and Fe XXVI) were detected, on top of a 1-50 keV continuum that is described by an NPEX model plus a soft blackbody. In addition, a cyclotron resonance scattering feature was detected significantly ($>99\%$ confidence), at an energy of 33$\pm$2 keV with a depth of 0.4$^{+0.6}_{-0.3}$. Therefore, the neutron star in this source is concluded to have a strong magnetic field of 2.8$\times10^{12}$ G. Further assuming that the source has a relatively high intrinsic luminosity of several times 10$^{37}$ erg s$^{-1}$, its spectral and timing properties are consistently explained.

Improvements of in-orbit calibration of GSO scintillators in the Hard X-ray Detector on board Suzaku are reported. To resolve an apparent change of the energy scale of GSO which appeared across the launch for unknown reasons, consistent and thorough re-analyses of both pre-launch and in-orbit data have been performed. With laboratory experiments using spare hardware, the pulse height offset, corresponding to zero energy input, was found to change by ~0.5 of the full analog voltage scale, depending on the power supply. Furthermore, by carefully calculating all the light outputs of secondaries from activation lines used in the in-orbit gain determination, their energy deposits in GSO were found to be effectively lower, by several percent, than their nominal energies. Taking both these effects into account, the in-orbit data agrees with the on-ground measurements within ~5%, without employing the artificial correction introduced in the previous work (Kokubun et al. 2007). With this knowledge, we updated the data processing, the response, and the auxiliary files of GSO, and reproduced the HXD-PIN and HXD-GSO spectra of the Crab Nebula over 12-300 keV by a broken powerlaw with a break energy of ~110 keV.

Using data from more than ten-years of observations with the Akeno Giant Air Shower Array (AGASA), we published a result that the energy spectrum of ultra-high energy cosmic rays extends beyond the cutoff energy predicted by Greisen, and Zatsepin and Kuzmin. In this paper, we reevaluate the energy determination method used for AGASA events with respect to the lateral distribution of shower particles, their attenuation with zenith angle, shower front structure, delayed particles observed far from the core and other factors. The currently assigned energies of AGASA events have an accuracy of $\pm$25% in event-reconstruction resolution and $\pm$18% in systematic errors around 10$^{20}$eV. This systematic uncertainty is independent of primary energy above 10$^{19}$eV. Based on the energy spectrum from 10$^{14.5}$eV to a few times 10$^{20}$eV determined at Akeno, there are surely events above 10$^{20}$eV and the energy spectrum extends up to a few times 10$^{20}$eV without a GZK-cutoff.

We have developed an atmospheric monitoring system for the Telescope Array experiment at Akeno Observatory. It consists of a Nd:YAG laser with an alt-azimuth shooting system and a small light receiver. This system is installed inside an air conditioned weather-proof dome. All parts, including the dome, laser, shooter, receiver, and optical devices are fully controlled by a personal computer utilizing the Linux operating system. It is now operated as a back-scattering LIDAR System. For the Telescope Array experiment, to estimate energy reliably and to obtain the correct shower development profile, the light transmittance in the atmosphere needs to be calibrated with high accuracy. Based on observational results using this monitoring system, we consider this LIDAR to be a very powerful technique for Telescope Array experiments. The details of this system and its atmospheric monitoring technique will be discussed.

Anisotropy in the arrival directions of cosmic rays around 10^{18}eV is studied using data from the Akeno 20 km^2 array and the Akeno Giant Air Shower Array (AGASA), using a total of about 216,000 showers observed over 15 years above 10^{17}eV. In the first harmonic analysis, we have found significant anisotropy of $\sim$ 4 % around 10^{18}eV, corresponding to a chance probability of $\sim 10^{-5}$ after taking the number of independent trials into account. With two dimensional analysis in right ascension and declination, this anisotropy is interpreted as an excess of showers near the directions of the Galactic Center and the Cygnus region. This is a clear evidence for the existence of the galactic cosmic ray up to the energy of 10^{18}eV. Primary particle which contribute this anisotropy may be proton or neutron.

The Telescope Array prototype detectors were installed at Akeno Observatory and at the Utah Fly's Eye site. Using these detectors, we have observed the Crab Nebula and AGN's since the end of 1995. The successful detections of TeV gamma rays from Crab Nebula and Mkn501 are reported.