The anomalous X-ray pulsar 4U 0142+61 was observed with Suzaku on 2007 August 15 for a net exposure of -100 ks, and was detected in a 0.4 to ~70 keV energy band. The intrinsic pulse period was determined as 8.68878 \pm 0.00005 s, in agreement with an extrapolation from previous measurements. The broadband Suzaku spectra enabled a first simultaneous and accurate measurement of the soft and hard components of this object by a single satellite. The former can be reproduced by two blackbodies, or slightly better by a resonant cyclotron scattering model. The hard component can be approximated by a power-law of photon index \Gamma h ~0.9 when the soft component is represented by the resonant cyclotron scattering model, and its high-energy cutoff is constrained as >180 keV. Assuming an isotropic emission at a distance of 3.6 kpc, the unabsorbed 1-10 keV and 10-70 keV luminosities of the soft and hard components are calculated as 2.8e+35 erg s^{-1} and 6.8e+34 erg s^{-1}, respectively. Their sum becomes ~10^3 times as large as the estimated spin-down luminosity. On a time scale of 30 ks, the hard component exhibited evidence of variations either in its normalization or pulse shape.

We report the results of a Suzaku observation of the anomalous X-ray pulsar (AXP) 1E 1841-045 at a center of the supernova remnant Kes 73. We confirmed that the energy-dependent spectral models obtained by the previous separate observations were also satisfied over a wide energy range from 0.4 to ~70 keV, simultaneously. Here, the models below ~10 keV were a combination of blackbody (BB) and power-law (PL) functions or of two BBs wit h different temperatures at 0.6 - 7.0 keV (Morii et al. 2003), and that above ~20 keV was a PL function (Kuiper Hermsen Mendez 2004). The combination BB + PL + PL was found to best represent the phase-averaged spectrum. Phase-resolved spectroscopy indicated the existence of two emission regions, one with a thermal and the other with a non-thermal nature. The combination BB + BB + PL was also found to represent the phase-averaged spectrum well. However, we found that this model is physically unacceptable due to an excessively large area of the emission region of the blackbody. Nonetheless, we found that the temperatures and radii of the two blackbody components showed moderate correlations in the phase-resolved spectra. The fact that the same correlations have been observed between the phase-averaged spectra of various magnetars (Nakagawa et al. 2009) suggests that a self-similar function can approximate the intrinsic energy spectra of magnetars below ~10 keV.

Spectral studies of quiescent emission and bursts of magnetar candidates using XMM-Newton, Chandra and Swift data are presented. Spectra of both the quiescent emission and the bursts for most magnetar candidates are reproduced by a photoelectrically absorbed two blackbody function (2BB). There is a strong correlation between lower and higher temperatures of 2BB (kT_LT and kT_HT) for the magnetar candidates of which the spectra are well reproduced by 2BB. In addition, a square of radius for kT_T (R_LT^2) is well correlated with a square of radius for kT_HT (R_HT^2). A ratio kT_LT/kT_HT ~ 0.4 is nearly constant irrespective of objects and/or emission types (i.e., the quiescent emission and the bursts). This would imply a common emission mechanism among the magnetar candidates. The relation between the quiescent emission and the bursts might be analogous to a relation between microflares and solar flares of the sun. Three AXPs (4U 0142+614, 1RXS J170849.0-400910 and 1E 2259+586) seem to have an excess above ~7 keV which well agrees with a non-thermal hard component discovered by INTEGRAL.

We have searched for pulsation of the anomalous X-ray pulsar (AXP) 4U 0142+61 in the K' band ($\lambda_{\rm eff} = 2.11$ $\mu$m) using the fast-readout mode of IRCS at the Subaru 8.2-m telescope. We found no significant signal at the pulse frequency expected by the precise ephemeris obtained by the X-ray monitoring observation with RXTE. Nonetheless, we obtained a best upper limit of 17% (90% C.L.) for the root-mean-square pulse fraction in the K' band. Combined with i' band pulsation (Dhillon et al. 2005), the slope of the pulsed component ($F_\nu \propto \nu^\alpha$) was constrained to $\alpha > -0.87$ (90% C.L.) for an interstellar extinction of $A_{V} = 3.5$.

Anomalous X-ray pulsars (AXPs) are thought to be magnetars, which are neutron stars with ultra strong magnetic field of $10^{14}$-- $10^{15}$ G. Their energy spectra below $\sim$10 keV are modeled well by two components consisting of a blackbody (BB) ($\sim$0.4 keV) and rather steep power-law (POW) function (photon index $\sim$2-4). Kuiper et al.(2004) discovered hard X-ray component above $\sim$10 keV from some AXPs. Here, we present the Suzaku observation of the AXP 1E 1841-045 at the center of supernova remnant Kes 73. By this observation, we could analyze the spectrum from 0.4 to 50 keV at the same time. Then, we could test whether the spectral model above was valid or not in this wide energy range. We found that there were residual in the spectral fits when fit by the model of BB + POW. Fits were improved by adding another BB or POW component. But the meaning of each component became ambiguous in the phase-resolved spectroscopy. Alternatively we found that NPEX model fit well for both phase-averaged spectrum and phase-resolved spectra. In this case, the photon indices were constant during all phase, and spectral variation seemed to be very clear. This fact suggests somewhat fundamental meaning for the emission from magnetars.

We study thermal structure and evolution of magnetars as cooling neutron stars with a phenomenological heat source in a spherical internal layer. We explore the location of this layer as well as the heating rate that could explain high observable thermal luminosities of magnetars and would be consistent with the energy budget of neutron stars. We conclude that the heat source should be located in an outer magnetar's crust, at densities rho < 5e11 g/cm^3, and should have the heat intensity of the order of 1e20 erg/s/cm^3. Otherwise the heat energy is mainly emitted by neutrinos and cannot warm up the surface.

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.

Japanese magnetospheric explorer GEOTAIL recorded a detailed light curve during the initial 600 msec of a giant flare from SGR 1806-20 on December 27, 2004. We show that the observed light curve is well explained by an emission from relativistically expanding fireballs, like those of gamma-ray bursts (GRBs). Especially, the observed rapid fading after 500msec suggests that ejecta is collimated in a jet. We derive an upper limit on the jet opening half-angle of 0.2 radian that is as narrow as those of GRBs.

We report a glitch-like pulse frequency deviation from the simple spin-down law in an anomalous X-ray pulsar (AXP) 4U 0142+61 detected by ASCA observations. We also found a significant pulse profile change after the putative glitch. The glitch parameters resemble those found in another AXP 1RXS J170849.0$-$400910, in the Vela pulsar, and in other radio pulsars. This suggests that the radio pulsars and AXPs have the same internal structure and glitch mechanism. It must be noted, however, that the pulse frequency anomaly can also be explained by a gradual change of the spin-down rate ($\dot{P}$) without invoking a glitch.

PSR 1259-63 is a radio pulsar orbiting a Be star in a highly eccentric orbit. Soft and hard X-rays are observed from this binary system. We apply the shock powered emission model to this system. The collision of the pulsar and Be star winds forms a shock, which accelerates electrons and positrons to the relativistic energies. We derive the energy distribution of relativistic electrons and positrons as a function of the distance from the shock in the pulsar nebula. We calculate the X-rays and $\gamma$-rays emitted from the relativistic electrons and positrons in the nebula at various orbital phases, taking into account the Klein-Nishina effect fully. The shock powered emission model can explain the observed X-ray properties approximately. We obtain from the comparison with observations that a fraction of $\sim 0.1$ of the pulsar spin-down luminosity should be transformed into the relativistic electrons and positrons. We find that the magnetization parameter of the pulsar wind, the ratio of the Poynting flux to the kinetic energy flux, is $\sim 0.1$ immediately upstream of the termination shock of the pulsar wind, and may decrease with distance from the pulsar. We predict the flux of 10 MeV - 100 GeV $\gamma$-rays which may be nearly equal to the detection threshold in the future projects.

The rotating neutron superfluid in the inner crust of a neutron star is threaded by quantized vortex lines. The pinning force from lattice nuclei and the Magnus force from neutron superfluid act onto a vortex line that has a finite tension. The configuration of a vortex line in equilibrium comprises a number of pinned straight lengths separated by unpinned kinks when the rotation axis is slightly inclined to the major axis of a crystal lattice and the Magnus force is not so strong. Energy of ~ 4 MeV is required to form a kink at densities of ~ 3.4 x 10^13 g cm^-3. Magnus force makes the vortex line parabolic on global scale. There exist two modes, rotational and helical (Kelvin), for oscillations allowed on a vortex line. The vortex oscillations in the pinned straight segments are possible only above a minimum frequency. We find no unstable mode that grows with time. Hence, the vortex configurations with kinks may be stable and yield the strongly pinned state. The essence of our results may not be altered even when we consider the vortex behaviors in a polycrystalline structure. Our studies suggest that pinning of a vortex line in a polycrystalline structure would be still strong enough to explain the large glitches.

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.

The rotating neutron superfluid in the inner crust of a neutron star is threaded by quantized vortex lines. The pinning force from lattice nuclei and the Magnus force from neutron superfluid act onto a vortex line that has a finite tension. The configuration of a vortex line in equilibrium comprises a number of pinned straight lengths separated by unpinned kinks when the rotation axis is slightly inclined to the major axis of a crystal lattice and the Magnus force is not so strong. Energy of $\sim 4 {MeV}$ is required to form a kink at densities of $\sim 3.4 \times 10^{13}{g cm}^{-3}$. Magnus force makes the vortex line parabolic on global scale. There exist two modes, rotational and helical (Kelvin), for oscillations allowed on a vortex line. The vortex oscillations in the pinned straight segments are possible only above a minimum frequency. We find no unstable mode that grows with time. Hence, the vortex configurations with kinks may be stable and yield the strongly pinned state. The essence of our results may not be altered even when we consider the vortex behaviors in a polycrystalline structure. Our studies suggest that pinning of a vortex line in a polycrystalline structure would be still strong enough to explain the large glitches.

We study the structure and emission of the SS 433 jets in the X-ray emitting region and in the inner and hotter portion inside the X-ray emitting region. In order to consider the jet structure from the inner to outer regions we develop the hybrid model combining the conical beam and the model beam whose cross section grows with the distance more slowly. We find that the jet beams in the inner and hotter portion are of two-temperature and emit a large amount of high energy gamma photons. Our analyses suggest the thick absorbing envelope to exist in the SS 433 system. Based on our results, we discuss the possible acceleration mechanism for the SS 433 jets.