arxiv: v1 [astro-ph.he] 18 Jan 2014

Transcription

1 PASJ: Publ. Astron. Soc. Japan,??, c 4. Astronomical Society of Japan. Firm Detection of a Cyclotron Resonance Feature with Suzaku in the X-ray Spectrum of GRO J87 during a Giant Outburst in arxiv:4.43v [astro-ph.he] 8 Jan 4 Takayuki YAMAMOTO, Tatehiro MIHARA, Mutsumi SUGIZAKI, Motoki NAKAJIMA, Kazuo MAKISHIMA,,3 and Makoto SASANO 3 MAXI team, RIKEN, - Hirosawa, Wako, Saitama 3-98 tyamamot@crab.riken.jp School of Dentistry at Matsudo, Nihon University, -87- Sakaecho-nishi, Matsudo, Chiba Department of Physics, The University of Tokyo, 7-3- Hongo, Bunkyo, Tokyo 3-33 (Received 3 November 9; accepted 4 January 8) Abstract 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 J87, achieved by the Suzaku Hard X-ray Detector during a giant outburst which was detected by the MAXI Gas Slit Camera in November. The Suzaku observation was carried out on November, outburst maximum when the X-ray flux reached.4 Crab in 4 kev, which corresponds to a luminosity of. 38 erg s in. kev at.8 kpc. The obtained broadband X-ray spectrum from. kev to 8 kev revealed a significant absorption feature, considered as the fundamental CRSF, at 76 kev. This unambiguously reconfirm the previously suggested 8 kev spectral feature in GRO J87. The implied surface magnetic field, 6.6 G, is the highest among binary X-ray pulsars from which CRSFs have ever been detected. Key words: stars: magnetic fields pulsars: individual (GRO J87) stars: neutron X-rays: binaries. Introduction Pulsars, which exhibit pulsating electromagnetic radiations in various wavelength, are strongly magnetized neutron stars. The rotation of the neutron star, combined with anisotropic radiation, causes the periodic pulsation. Although they are considered to be formed by supernova explosions of massive stars, the origin and time evolution of their magnetic fields are still open questions. X-ray binary pulsars (XBPs) are a group of X-ray binaries involving pulsating neutron stars. According to the type of the binary companion, they are classified into several subgroups including Super Giant XBPs and Be XBPs as major members (e.g. Reig ). Be XBPs produce recurrent outbursts synchronized with their binary orbital periods. The outbursts are considered to occur when the neutron star crosses a gaseous stellar disk of the Be star near the periastron passage. The outburst does not always appear every orbital cycle, and sometime arises in an irregular orbital phase, probably depending on the physical extent of the stellar disk. Surface magnetic fields of neutron stars in XBPs can be estimated from the cyclotron resonance scattering feature (CRSF), which has been observed as absorption features in their X- ray spectra. The CRSF is considered to appear at an energy of E a =.6(+z g ) B, where B is the magnetic field strength in G, and z g represents the gravitational redshift. Ginga/LAC observations in the 6 kev band detected the CRSFs from XBPs and showed that their surface magnetic fields are distributed in a very narrow range of (. 3.) G (Mihara et al. 998, Makishima et al. 999). Subsequently, ASCA, RXTE, BeppoSAX, INTEGRAL and Suzaku observations surveyed a wider energy band from. kev upto a few hundreds kev, and detected CRSFs from additional six XBPs (e.g. Coburn et al. ; Filippova et al. 7; Doroshenko et al. ; Yamamoto et al. ; Tsygankov et al. ; DeCesar et al. 3). However, the revised range of their surface magnetic fileds,(. 4.7) G, is still narrow. It is yet to be clarified whether this is intrinsic to XBPs, or a selection effect due to limited observations. GRO J87 is a Be XBP with a pulsation period of 93. s, discovered by the CGRO/BATSE in 993 (Stollberg et al. 993). Its optical counterpart was identified with a Be type star (Coe et al. 994) and the distance was estimated to be.8 kpc (Riquelme et al. ). Its X-ray outbursts have been monitored for about years by surveys with the CGRO/BATSE, RXTE/ASM, Swift/BAT, and MAXI/GSC. Since 3 January, the source has been in an active state exhibiting outbursts periodically (Kühnel et al. 3). From the recurrent outburst intervals and the pulsar period modulation, the binary orbital period was estimated as 47.8±.4 d (Coe et al. 7), which was recently refined to ±.4 d by the pulse arrival-time analysis (Kühnel et al. 3). Based on the CGRO/OSSE pointing observations performed in the 993 outburst and the BATSE earth-occultation data on that occasion, Shrader et al. (999) suggested a possible CRSF at around 88 kev in the X-ray spectra of GRO J87. In contrast, spectra of the 4 outburst obtained by the INTEGRAL/IBIS and JEM-X showed no feature in the 3 6 kev band (Coe et al. 7). Therefore, the possible absorption feature at 88 kev is considered to be the fundamental if it is real. Observations of the 7 November outburst by RXTE, Swift, and Suzaku were unable to confirm the suggestion, hampered by rather poor signal statistics (Naik et al. ; Kühnel et al. 3).

2 Yamamoto et al. [Vol., Flux [Crab] Dec Aug. Apr. Dec. Sep. Nov. 3 May Aug. 3 Jun. MJD Fig.. MAXI/GSC light curve of GRO J87 in the 4 kev band from 9 August to June. Vertical dash lines indicate the predicted epochs of periastron passage of the pulsar (Kühnel et al. 3). The inset zooms up the profile of the outburst in November, on which Suzaku data taken on November - are superposed with circles. In the present paper, we report the Suzaku observation performed at the peak of a giant outburst detected by MAXI in November, and the results of the spectral analysis for the CRSF. Unless otherwise specified, all errors hereafter refer to 9% confidence limits.. Outburst Activity Monitored by MAXI MAXI (Matsuoka et al. 9) Gas Slit Camear (GSC; Mihara et al. ) has been monitoring the X-ray flux of GRO J87 since the operation started on 9 August (MJD=8) (Sugizaki et al. ). Figure shows the obtained light curve until 3 June (MJD 64). By September (MJD 6), five outbursts were detected periodically at the same orbital phase close to the pulsar periastron passage. Their peak intensities in 4 kev are almost same at. Crab, which corresponds to a. kev luminosity of L X 37 erg s at.8 kpc assuming the same spectral shape as in Suzaku observation (section 3.4). Thus, these outbursts are categorized into the normal-type ones (Reig ). On November, the source exhibited unexpected brightening at an irregular orbital phase which is.3 cycle after the periastron (Nakajima et al. ). The 4 kev intensity reached.4 Crab at the maximum. Judging from the outburst phase and the peak luminosity, it is categorized into the giant-type outburst (Reig ). 3. Suzaku Observation of Giant Outburst and Data Analysis 3.. Observation and Data Reduction Triggered by the MAXI detection of the giant outburst from GRO J87, we requested a Suzaku ToO (Target of Opportunity) observation. It was performed on November, nearly coincident with the outburst maximum. Suzaku covers an energy band from. to kev with the X- ray Imaging Spectrometer (XIS: Koyama et al. 7) and the Hard X-ray Detector (HXD: Takahashi et al. 7, Kokubun et al. 7). The target was placed at the XIS nominal position on the focal-plane. The XIS was operated in the normal mode with /4-window and.3 s burst options, which affords a time resolution of s. The HXD was operated in the nominal mode. Table summarizes the Suzaku observations including exposure and count rate in each instrument. The data reduction and analysis were performed with the standard procedure using the Suzaku analysis software in HEASOFT version 6. and the CALDB files version 93, provided by the NASA/GSFC Suzaku GOF. All obtained data were first reprocessed by a Suzaku software tool, aepipeline to utilize the latest calibration. The net exposures after the standard event-screening process were 8. ks with the XIS and.4 ks with the HXD. Due to the.3 s burst option, the XIS exposure is about one third of that of the HXD. We started the XIS data analysis with the standard cleaned event files. On-source event data were collected from a circular region of 4 radius around the source position on the XIS CCD images, and background data from an annulus with the inner and outer radii of 3 and 4, respectively. The pileup effect on each image pixel was estimated by the Suzaku PileupTools. We excluded pixels on the image core in which the estimated pileup fraction is larger than % (Yamada et al. ). In the HXD data analysis, we created the background spectra with the standard procedure, using the archived background files provided by the Suzaku GOF. The obtained HXD- PIN background includes contribution from the Cosmic X- ray Background (CXB), while it is negligible in the HXD- GSO data (Fukazawa et al. 9). After subtracting the backgrounds, the source count rates became 4.6 ±. counts s in the PIN 6 kev band, and.9±. counts s in the GSO 6 kev band. 3.. Timing Analysis With the Suzaku analysis tool, aebarycen, we converted the photon arrival times of all events into those at the solarsystem barycenter and then searched the data for the coherent pulsation by epoch-folding analysis. The 93. s pulsation was detected significantly, both with the XIS and the HXD, and the best period was obtained as ±. s with the HXD-PIN data. Figure shows the folded pulse profiles in the XIS, HXD-PIN and HXD-GSO energy bands, where the phase φ = is set at the minimum in the HXD-PIN profile. We divided the HXD-GSO band into three, 7 kev, 7 8 kev and 8 kev, around the CRSF energy (section 3.3). The -7 kev profile is the same as that of HXD-PIN, while that in 7 8 kev is somewhat different. The pulsation is still significant in the highest energy band of 8- kev. The pulse profile in the XIS has two peaks at φ. and φ.6. The former tends to decrease towards higher energies. These double-peak profiles and their energy dependence are largely consistent with the results obtained in previous outbursts (Shrader et al. 999; Coe et al. 7; Naik et al. ; Kühnel et al. 3). However, details are rather different. Comparing the XIS-band profiles, the former peak ob- yamada/soft/ XISPileupDoc /XIS PileupDoc.html

3 No. ] Cyclotron Resonance Feature in GRO J87 3 Table. Log of Suzaku Observation of GRO J87 in the November Giant Outburst Date Obs Time XIS-FI (.8 kev) HXD-PIN ( 6 kev) HXD-GSO (6 kev) ( Start/End Exp. Rate Exp. Rate Exp. Rate Nov.) (UT) (ks) (counts s ) (ks) (counts s ) (ks) (counts s ) 4:44/: ± ±..38.9±. Observation ID = 976 Counts s.. Valley Rise Peak Fall XIS PIN GSO GSO kev kev 7 kev 7 8 kev. GSO.. 8 kev.. Phase Fig.. Folded pulse profiles by the XIS in kev, by HXD-PIN in kev, and by HXD-GSO in the 7 kev, 7 8 kev, and 8 kev bands. The dashed lines divide the pulse cycle into Valley, Rise, Peak, and Fall, referring to the -7 kev profile. tained here is apparently smaller than that in the 7 outburst. As illustrated in figure, we divided the pulse cycle into four phases, and named them Valley, Rise, Peak, and Fall, according to the profile in the GSO 7 kev band. They are used in the phase resolved spectral analysis in section Cyclotron Resonance Feature in Averaged Spectrum We examined a pulse-phase-averaged spectrum with the best photon statistics for the previously suggested CRSF signatures. All the spectral fitting attempts hereafter were carried out on Xspec version.7.. The cross normalization factor between the XIS and the HXD was fixed at :.6 according to the latest calibration information. We discard the energy bands of.7.9 kev around the silicon K edge and..4 kev around the gold M edges in the XIS data, where the calibration uncertainty is relatively larger. We did not use the XIS-BI data either in the spectral analysis, because it has larger calibration uncertainties than XIS-FI. Figure 3 (a) shows ratios of the spectra obtained with XIS-FI (.8 kev), HXD-PIN ( 6 kev) and HXD-GSO (6 kev) to those of the Crab nebula which has a simple powerlaw shape with a photon index of.. Figure 3 (b) shows the count-rate spectra without removing instrument responses. From the Crab ratios, the spectrum is found to be largely approximated by a smooth continuum with cutoffs below kev and above kev. In addition, iron-k emission lines at around 6. kev and an edge-like feature at around 7 8 kev are clearly seen. The latter looks like a typical CRSF observed in some XBP spectra, and its energy is close to those of the possible absorption features ( 88 kev) reported in past outbursts (Shrader et al. 999; Kühnel et al. 3). We fitted the spectrum above kev with typical XBP continuum models; cutoff power-law (CPL, cutoffpl in Xspec terminology), FDCO (Fermi-Dirac cutoff power-law; Tanaka 986), and NPEX (Negative and Positive power laws with EXponential cutoff: Mihara et al. 998) whose positive powerlaw index was fixed at.. However, as exemplified in figure 4 (a), none of these models alone were able to fit the data sufficiently, because of the feature at 7 8 kev. We thus applied a cyclotron absorption factor (CYAB; cyclabs in Xspec terminology, Mihara et al. 99, Makishima et al. 999) to the above continuum models. Since the width W of the CYAB factor cannot be constrained lower than the energy resolution, kev at 8 kev in HXD GSO, we set its lower limit at kev in the model fits. Then, all three continuum models became acceptable within 9% confidence limits, and the improvements of chi-squared ( ν ) for degree of freedom (ν) were estimated with the F-test to be significant above the 99% confidence limit. The case with NPEX * CYAB is shown in figure 4 (b). As listed in table, the best-fit CRSF energy, E a 7 8 kev, slightly depends on the continuum model. Since the CRSF energy, 8 kev, is rather high, we should examine the possibility that it is in reality the second harmonic. Actually, Vela X- has been sometimes reported to show a shallow absorption feature at kev (Makishima et al. 999), possibly interpreted as the fundamental CRSF, in addition to the more prominent feature at kev which is confirmed in many observations (Mihara et al. 998; Orlandini et al. 998; Makishima et al. 999; Kreykenbohm et al. 999; Kreykenbohm et al. ; Odaka et al. 3). We hence fitted the Suzaku spectra of GRO J87 by a pair of harmonic CYAB factors, with the fundamental resonance energy around 4 kev. Then, as shown in figure 4(c) and given in table, the best-fit ν slightly decreased to.8 from.4 of the initial single-cyab model, yielding E 37 kev. However, we consider this harmonic interpretation rather unlikely for the

4 4 Yamamoto et al. [Vol., Crab Model Ratio Counts s kev XIS FI phabs*(npex+bb+3gaus) phabs*(npex+bb+3gaus)*cyab PIN phabs*pcfabs*(npex+bb+3gaus)*cyab GSO Energy [kev] Fig. 3. Pulse-phase-averaged and background-subtracted broadband X-ray spectrum of GRO J87 by Suzaku XIS-FI (.8- kev), HXD-PIN ( 6 kev) and HXD-GSO (6 kev). (a) Ratio to the Crab Nebula spectrum. (b) Count-rate spectra and the folded best-fit model ofpcfabs phabs (NPEX+BB+3gaus) CYAB. (c) Residuals against phabs (NPEX + BB + 3gaus) model. (d) Residuals against phabs (NPEX+BB+3gaus) CYAB model. (e) Residuals againstpcfabs phabs (NPEX+BB+3gaus) CYAB model. following reasons. First, such a local feature at 4 kev is not visible in figure 4(b). Second, an F -test indicates that the fit improvement by introducing the second CRSF factor is less significant than 8%. Third, the derived ratio D /D =. of GRO J87 is even smaller than that of Vela X-, D /D =.7/.8=.9 (Makishima et al. 999). Finally, the obtained width W kev for the lower-energy feature is much wider than those of Vela X- (W. kev) and the higher-energy feature of GRO J87 (W =. kev). Therefore, we consider that the deep 7 8 kev feature GRO J87 is the fundamental resonance, although the alternate interpretation, that it is the second harmonic resonance, cannot be completely ruled out Broadband Spectral Model for Averaged Spectrum Now that the - kev HXD spectrum was successfully modeled and the CRSF was clearly detected, the next step is to search for broadband emission models that can explain the whole Suzaku spectrum from. kev to kev. We first (a) (b) (c) (d) (e) NPEX NPEX*CYAB NPEX*CYAB*CYAB Energy [kev] Fig. 4. Data-to-model residuals when the HXD data are fitted with (a) NPEX alone, (b) NPEX*CYAB, and (c) NPEX*CYAB*CYAB in kev band. tested the CPL, FDCO, and NPEX continuum models as used in the previous section, incorporating a CYAB factor at 7 8 kev and an interstellar absorption (phabs in Xspec terminology) whose hydrogen column density N H was set free. In any continuum model, however, the fit was far from acceptable, and the data-to-model residuals showed an excess in the soft X-ray band below 3 kev and iron K-lines at around 6. kev. This agrees with the results obtained in past outbursts (e.g. Naik et al. ; Kühnel et al. 3). We then added a blackbody (BB) model to account for the soft X-ray residuals, and three narrow gaussians (gaus) for K α lines from neutral iron (6.4 kev) and helium-like iron (6.7 kev), as well as K β line at 7. kev. Among the three continuum models, the NPEX-based composite model, expressed by phabs (NPEX + BB + 3gaus) CYAB, fit the data much better than the other two, but it is still unacceptable with ν =.8 for ν = 3 degree of freedom. The residuals, as shown in figure 3 (d), indicate that discrepancy remains at around 4 kev and around kev. To improve the fit, we tried to apply a partially covering absorption model (pcfabs in Xspec terminology). The fit became even better with ν =.4 forν=348 as shown in figure 3 (e). However, it is still outside the 9% confidence limit. This may be because the phase-averaged spectrum has complex features that can arise by averaging pulse-phase dependent spectra. Table 3 summarizes all the best-fit model parameters. 3.. Pulse-Phase Resolved Spectra As seen in figure, the folded pulse profiles from kev to kev are apparently energy dependent. This means that the energy spectrum depends on the pulse phase. We thus extracted four spectra, one from each of the four pulse phases defined in figure from the 7 kev GSO pulse profile. Figure shows the obtained four spectra in a form of their ratios to the phaseaveraged spectrum. While the spectrum in the Valley phase is softer than the average, that in the Peak phase is harder. No feature are apparent at the iron K-line band around 6 7 kev in these ratio plots. This indicates that the equivalent width of the iron lines does not depend on the pulse phases. We fitted each phase-resolved spectrum with the model which best described the phase averaged spectrum in section (a) (b) (c)

5 No. ] Cyclotron Resonance Feature in GRO J87 Table. Summary of model fits to Suzaku HXD kev spectrum with either CPL, FDCO, or NPEX continuum and zero, one, or two CYABs. Continuum CPL FDCO NPEX None CYAB CYAB None CYAB CYAB None CYAB CYAB α E cut (kev) E fold /kt (kev) A A D E a (kev) W (kev) D E a =E a W (kev) ν (ν) 6.(64).(6).4(9).9(63).9(6).(8).73(63).4(6).8(8) Units in photons s cm kev at kev. Table 3. Best-fit models for Suzaku.8 kev broadband spectrum Model Function phabs phabs phabs phabs pcfabs Component Parameter Cont(CPL) Cont(FDCO) Cont(NPEX) Cont(NPEX) CYAB CYAB CYAB CYAB phabs N H ( cm ) bbody kt BB (kev) I BB ( 3 ) gaus E Fe Kα (kev) I Fe K α ( 3 ) gaus E Fe 6.7 (kev) I Fe 6.7 ( 3 ) gaus3 E Fe Kβ (kev) I Fe K β ( 4 ) pcfabs N H ( cm ) f PCF NPEX α E cut (kev) 4.7. kt/e fold (kev) A ( ) A ( 4 ) CYAB D E a (kev) W (kev) ν (ν).3 (3) 3.9 (3).8 (3).4 (348) L. kev Units in photons s cm. Units in photons s cm kev at kev. Cont(CPL) is a CPL-based composite model including a soft BB and three iron-lines, expressed by CPL+BB+3gaus. Cont(FDCO) andcont(nepx) represent FDCO-based and NPEX-based composite models, respectivly. Units in 37 erg s.

6 6 Yamamoto et al. [Vol., Ratio.... Valley Rise Peak Fall Energy [kev] Fig.. Pulse-phase resolved spectra for Valley, Rise, Peak, and Fall phases in the form of the ratio to the model that best describes the pulse-phase averaged spectrum. GSO 7 kev [Counts s ] Depth E a [kev] Phase Fig. 6. Pulse-phase dependence of the CRSF depth D and energy E a. The vertical error bars ofd ande a represent 9% confidence limits of statitical uncertainties Here, we fixed the iron K β -line energy at 7. kev and the CYAB width at a typical value of. kev because they were poorly constrained by the data with lower statistics. The fits became acceptable except for the Peak phase. Table 4 summarizes the obtained best-fit parameters in each phase. Although these model parameters for the continuum are correlated in complex ways, those of the CRSF model are mostly free from them. As shown in figure 6, the derived CRSF parameters show some dependence on the pulse phase, but not more significantly than errors. Table. CRSF measurements in GRO J87. CRSF parameters Outburst Luminosity E a W D ( 37 erg s ) (kev) Nov (fix) Dec (fix) 993 Jul Calculated from the best spectral model in. kev. This work, Kühnel et al. (3), Shrader et al. (999), 4. Discussion 4.. Possible CRSF Energy Change We analyzed the broadband X-ray (.8 kev) spectrum of GRO J87 obtained by Suzaku, covering the peak of the giant outburst in detected by MAXI, and found a significant absorption signature at 7 8 kev (Yamamoto et al. 3). It can be interpreted as a fundamental CRSF, and reconfirms, with much higher significance, the previous suggestions (Shrader et al. 999, Kühnel et al. 3). Table compares the CRSF parameters and luminosity obtained in this work with those of the previous outbursts, and figure 7 gives its graphical plot. Thus, the CRSF energy might decrease towards higher luminosities, although the presently available information is very limited. The luminosity dependence of the CRSF energy has been observed in several XBPs. While some of them, 4U +63 (Mihara et al. 998; Mihara et al. 4; Nakajima et al. 6) and V 33+3 (Tsygankov et al. 6; Mowlavi et al. 6; Nakajima et al. ), showed negative correlations, others, Her X- (Gruber et al. ; Staubert et al. 7) and GX 34 (Yamamoto et al. ; Klochkov et al. ) showed positive. These are explained by variations of the cyclotronscattering photosphere; it increases by radiation pressure in the super-eddington luminosity regime (Mihara et al. 998), while it decreases due to dynamical pressure of the accretion in the sub-eddington luminosity (Staubert et al. 7). The luminosity of GRO J87 observed by Suzaku at the peak of the giant outburst,. 38 erg s, is close to the Eddington luminosity for the typical neutron-star mass of.4m. Therefore, the possible CRSF energy change suggests such a situation that the accretion mode changed from the sub- Eddington to the super-eddington regime at that time. 4.. Magnetic Fields in Binary Pulsars We now know 8 XBPs in which CRSFs are significantly detected and their parameters are determined well. The CRSF energy of 7 8 kev obtained here from GRO J87 is the highest among them. Therefore, the estimated surface magnetic field, 6.6 (+z g ) G, extends the highest end of their magnetic field distribution. Figure 8 shows the updated distribution of the XBP magnetic field strengths. It is still clustered in a very narrow range of (. 6.6) G, compared to the distribution of a larger number ( ) of single radio pulsars in the ATNF pulsar catalog (Manchester et al. ). Although the radio pulsars show considerably broader field distribution, this could be due to the much lower accuracy of their field determinations which assume spin down via mag-

7 No. ] Cyclotron Resonance Feature in GRO J87 7 Table 4. Best fit parameters of phase resolved spectra Pulse Phase Component Parameter Valley Rise Peak Fall phabs N H ( cm ) bbody kt BB (kev) I BB ( 3 ) gaus E Fe Kα (kev) I Fe Kα ( 3 ) gaus E Fe 6.7 (kev) I Fe 6.7 ( 3 ) gaus 3 I Fe Kβ ( 4 ) pcfabs N H ( cm ) f PCF NPEX α kt (kev) A ( ) A ( 4 ) CYAB D E a (kev) ν (ν).3 (7).99 (7).49 (7).7 (7) L. kev Spectal model function: phabs pcfabs (NPEX+BB+3gaus) CYAB Energy of iron K β is fixed to 7. kev. Width of CYAB is fixed to. kev. Units in photons s cm. Units in photons s cm kev at kev. Units in 37 erg s E a [kev] Luminosity [ 37 erg s ] Binary Pulsar Number Binary Pulsar Single Pulsar GRO J Single Pulsar Number Fig. 7. Luminosity (. kev) dependence of the CRSF energy from table. netic dipole radiation. In any case, the plots favor the scenario that the surface magnetic fields of neither XBPs nor radio pulsars would decay significantly within their lifetime of 8 yr (Itoh et al. 99; Makishima et al. 999). GRO J87 is known to have a large orbital eccentricity ofe=.68() (Coe et al. 7). In figure 9, we plot a relation between the surface magnetic fields and orbital eccentricities of XBPs whose CRSFs and binary orbital parameters are well determined. On this plot, GRO J87 locates at the upper right corner. Thus, the surface magnetic field and the orbital eccentricity of high mass X-ray binaries, including BeXBs, appear to have a positive correlation. This may suggest the evolutional relation between these parameters. Since the surface magmatic fields would not change as discussed above and the Magnetic Field Strength B [log G] Fig. 8. Distribution of magnetic field strengths on neutron stars in binary pulsars estimated from the CRSFs (red ordinate to the left) and in single pulsars from the period and the period derivative from ATNF pulsar catalog (Manchester et al. ) (black, to the right). It is updated from those in Mihara et al. (998) and Makishima et al. (999). orbital eccentricity would not change significantly within their lifetime, the correlation is considered to be formed when the XBPs are born. Further observational as well as theoretical studies are necessary. We thank the Suzaku operation team for arranging and carrying out the TOO observations. We are also grateful to all members of the MAXI and the ISS-operation teams. This research was partially supported by the Ministry of Education,

8 8 Yamamoto et al. [Vol., B [ G] 6 4 Her X 4U 66 BeXB HMXB LMXB X Per Cen X 3 38 Vela X Sw GRO J Orbital Eccentricity Fig. 9. The magnetic field strengths of XBPs measured with the CRSF technique, plotted against their orbital eccentricities. Meanings and references of labels to the data points are as follows : GRO J87 (This work; Coe et al. 7), 34 (GX 34 ; Yamamoto et al. ; Yamamoto 3; Yamamoto et al in prep.), 3 (A 3+6; Caballero et al. 3; Finger et al. 994), Her X (Vasco et al. ; Staubert et al. 9), 4U 66 (4U 66 67; Iwakiri et al. ; Chakrabarty et al. 997), 3 (GX 3 ; Suchy et al. ; Koh et al. 997), 946 (XTE J ; Maitra & Paul 3; Wilson et al. 3), X Per (Coburn et al. ; Delgado-Martí et al. ), Cen X 3 (Suchy et al. 8; Raichur & Paul ), 33 (X33+3; Nakajima et al. ; Raichur & Paul ), Vela X (Odaka et al. 3; Bildsten et al. 997), 38 (4U 38; Rodes-Roca et al. 9; Clark ), 97 (4U 97+9; Rivers et al. ; Baykal et al. 6), (4U +63; Nakajima et al. 6; Raichur & Paul ), Sw 66 (Swift J66.66; DeCesar et al. 3). Culture, Sports, Science and Technology (MEXT), Grant-in- Aid No and References Baykal, A., İnam, S. Ç., & Beklen, E. 6, MNRAS, 369, 76 Bildsten, L., et al. 997, ApJS, 3, 367 Caballero, I., et al. 3, ApJL, 764, L3 Chakrabarty, D., et al. 997, ApJ, 474, 44 Clark, G. W., ApJL, 4, L3 Coburn, W., Heindl, W. A., Gruber, D. E., Rothschild, R. E., Staubert, R., Wilms, J., Kreykenbohm, I.,, ApJ,, 738 Coburn, W., Heindl, W. A., Rothschild, R. E., Gruber, D. E., Kreykenbohm, I., Wilms, J., Kretschmar, P., & Staubert, R.,, ApJ, 8, 394 Coe, M. J., et al. 994, MNRAS, 7, L7 Coe, M. J., et al. 7, MNRAS, 378, 47 DeCesar, M. E., Boyd, P. T., Pottschmidt, K., Wilms, J., Suchy, S., Miller, M. C., 3, ApJ, 76, 6 Delgado-Martí, H., Levine, A. M., Pfahl, E., & Rappaport, S. A., ApJ, 46, 4 Doroshenko, V., Suchy, S., Santangelo, A., Staubert, R., Kreykenbohm, I., Rothschild, R. E., Pottschmidt, K., Wilms, J.,, A&A,, L Filippova, E. V., Tsygankov, S. S., Lutovinov, A. A., & Sunyaev, A. A. 7, ESA Special Publication, 6, 449 Finger, M. H., Wilson, R. B., & Hagedon, K. S. 994, IAU Circ., 93, Fukazawa, Y., et al. 9, PASJ, 6, 7 Gruber, D. E., Heindl, W. A., Rothschild, R. E., Coburn, W., Staubert, R., Kreykenbohm, I., & Wilms, J., ApJ, 6, 499 Itoh, N., Kotouda, T., & Hiraki, K. 99, ApJ, 4, 44 Iwakiri, W. B., et al., ApJ, 7, 3 Klochkov, D., et al., A&A, 4, L8 Koh, D. T., et al. 997, ApJ, 479, 933 Kokubun, M., et al. 7, PASJ, 9, 3 Koyama, K., et al. 7, PASJ, 9, 3 Kreykenbohm, I., Kretschmar, P., Wilms, J., Staubert, R., Kendziorra, E., Gruber, D. E., Heindl, W. A., Rothschild, R. E., 999, A&A, 34, 4 Kreykenbohm, I., Coburn, W., Wilms, J., Kretschmar, P., Staubert, R., Heindl, W. A., Rothschild, R. E.,, A&A, 39, 9 Kuehnel, M., et al., The Astronomer s Telegram, 477, Kühnel, M., et al. 3, A&A,, A9 Maitra, C., & Paul, B. 3, ApJ, 77, 96 Makishima, K., Mihara, T., Nagase, F. & Tanaka, Y., 999, ApJ,, 978 Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M., AJ, 9, 993 Matsuoka, M., et al. 9, PASJ, 6, 999 Mihara, T., Ph.D. thesis in University of Tokyo 99 Mihara, T., et al., PASJ, 63, 63 Mihara, T., Makishima, K., & Nagase, F., 998, Advances in Space Research,, 987 Mihara, T., Makishima, K., & Nagase, F., 4, ApJ, 6, 39 Mihara, T., Makishima, K., Ohashi, T., Sakao, T., & Tashiro, M. 99, Nature, 346, Mowlavi, N., et al. 6, A&A, 4, 87 Nakajima, M., et al., The Astronomer s Telegram, 46, Nakajima, M., Mihara, T., & Makishima, K., ApJ, 7, 7 Nakajima, M., Mihara, T., Makishima, K., & Niko, H. 6, ApJ, 646, Naik, S., Paul, B., Kachhara, C., Vadawale, S.,V., MNRAS, 43, 4 Odaka, H., Khangulyan, D., Tanaka, Y. T., et al. 3, ApJ, 767, 7 Orlandini, M., et al. 998, A&A, 33, Raichur, H., & Paul, B., MNRAS, 4, 3 Raichur, H., & Paul, B., MNRAS, 46, 663 Reig, P., Ap&SS, 33, Riquelme, M. S., Torrejón, J. M., & Negueruela, I., A&A, 39, A4 Rivers, E., et al., ApJ, 79, 79 Rodes-Roca, J. J., Torrejón, J. M., Kreykenbohm, I., Martinez, N. S., Camero-Arranz, A., Bernabeu, G., 9, A&A, 8, 39 Shrader, C. R., Sutaria, F. K., Singh, K. P., & Macomb, D. J. 999, ApJ,, 9 Staubert, R., Klochkov, D., & Wilms, J. 9, A&A,, 883 Staubert, R., Shakura, N. I., Postnov, K., Wilms, J., Rothschild, R. E., Coburn, W., Rodina, L., & Klochkov, D. 7, A&A, 46, L Stollberg, M. T., Finger, M. H., Wilson, R. B., et al. 993, IAU Circ., 836, Suchy, S., et al. 8, ApJ, 67, 487 Suchy, S., Fürst, F., Pottschmidt, K., Caballero, I., Kreykenbohm, I., Wilms, J., Markowitz, A., Rothschild, R. E.,, ApJ, 74, 4 Sugizaki, M., et al., PASJ, 63, 63 Takahashi, T., et al. 7, PASJ, 9, 3 Tanaka, Y. 986, IAU Colloq. 89: Radiation Hydrodynamics in Stars and Compact Objects,, 98 Truemper, J., Pietsch, W., Reppin, C., Voges, W., Staubert, R., Kendziorra, E., 978, ApJL, 9, L Tsygankov, S. S., Lutovinov, A. A., Churazov, E. M., & Sunyaev, R. A. 6, MNRAS, 37, 9

9 No. ] Cyclotron Resonance Feature in GRO J87 9 Tsygankov, S. S., Krivonos, R. A., & Lutovinov, A. A., MNRAS, 4, 47 Vasco, D., Klochkov, D., & Staubert, R., A&A, 3, A99 Verbunt, F. 993, ARA&A, 3, 93 White, N., Swank, J., & Holt, S. S., 983, ApJ, 7, 7 Wilson, C. A., Finger, M. H., Coe, M. J., & Negueruela, I. 3, ApJ, 84, 996 Yamada, S., Uchiyama, H., Dotani, T., et al., PASJ, 64, 3 Yamamoto, T., Ph.D. thesis in Nihon University 3 Yamamoto, T., et al., PASJ, 63, 7 Yamamoto, T., et al. 3, The Astron. Telegram, 479,