Publ. Astron. Soc. Japan (2014) 00(0), 1–11 1 doi: 10.1093/pasj/xxx000 An application of the Ghosh & Lamb model to the accretion powered X-ray pulsar X Persei Fumiaki YATABE1,2 , Kazuo MAKISHIMA1 , Tatehiro MIHARA1, Motoki NAKAJIMA3 , Mutsumi SUGIZAKI4 , Shunji KITAMOTO2 , Yuki YOSHIDA2 and Toshihiro TAKAGI1 1High Energy Astrophysics Laboratory, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan 2Department of Physics, College of Science, Rikkyo University, 3-34-1 Nishi-Ikebukuro, Toshima, Tokyo 171-8501, Japan 3School of Dentistry at Matsudo, Nihon University, 2-870-1 Sakaecho-nishi, Matsudo, Chiba 101-8308, Japan 4Department of Physics, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8551, Japan ∗E-mail: [email protected] Received ; Accepted Abstract The accretion-induced pulse-period changes of the Be/X-ray binary pulsar X Persei were inves- tigated over a period of 1996 January to 2017 September. This study utilized the monitoring data acquired with the RXTE/ASM in 1.5–12 keV and MAXI/GSC in 2–20 keV. The source intensity changed by a factor of 5–6 over this period. The pulsar was spinning down for 1996- 2003, and has been spinning up since 2003, as already reported. The spin up/down rate and the 3–12 keV flux, determined every 250 d, showed a clear negative correlation, which can be successfully explained by the accretion torque model proposed by Ghosh & Lamb (1979). When the mass, radius and distance of the neutron star are allowed to vary over a range of 1.0–2.4 solar masses, 9.5–15 km, and 0.77–0.85 kpc, respectively, the magnetic field strength of B = (4 − 25) × 1013 G gave the best fits to the observation. In contrast, the observed results cannot be explained by the values of B ∼ 1012 G previously suggested for X Persei, as long as the mass, radius, and distance are required to take reasonable values. Assuming a distance of 0.81 ± 0.04 kpc as indicated by optical astrometry, the mass of the neutron star is estimated arXiv:1807.06252v2 [astro-ph.HE] 18 Jul 2018 as M =2.03 ± 0.17 solar masses. Key words: pulsars: individual (X Persei)stars: neutronX-rays: binaries 1 Introduction is the gravitational redshift at the neutron-star (hereafter NS) surface. No accreting pulsars with B ≥ 1013 G have yet been The magnetic fields of binary X-ray pulsars have been di- discovered. However, this could be due to a selection effect, rectly measured utilizing cyclotron resonance scattering fea- that the CRSF becomes difficult to detect in ≥ 100 keV (i.e, tures (hereafter CRSF) in their X-ray spectra (e.g., Makishima for B ≥ 1013 G) where the observational sensitivity becomes 2016). The magnetic field strength B measured in this way progressively lower. In order to search for such accreting NSs 12 is distributed in the range of B/(1 + zg) =(1 − 7) × 10 G with higher magnetic fields, we need to invoke other methods (Yamamoto et al. 2014 and references therein), where zg ≃ 0.24 c 2014. Astronomical Society of Japan. 2 Publications of the Astronomical Society of Japan, (2014), Vol. 00, No. 0 such as pulse timing analysis or modeling of continuum spectra Thus, the strong-field property of X Persei is suggested also by (Sasano 2015). its spectral properties. Using the MAXI Gas Slit Camera (GSC) data and some past In the present paper, we report the long-term continuous ob- measurements, Takagi et al. (2016) investigated long-term re- servations of the fluxes and P of X Persei by RXTE and MAXI. lations between the pulse period derivative P˙ and the fluxes Applying the GL79 model to the derived P˙ − L relation, we of the binary X-ray pulsar 4U 1626-67. They adopted the then estimate B and M of the NS in this binary. accretion torque theory proposed by Ghosh & Lamb (1979) (hereafter GL79), which describes P˙ as a function of the lu- minosity L, the pulse period P , the mass M, the radius R, 2 Observation and the surface magnetic field B of the NS. Then, employing 2.1 MAXI/GSC 12 B/(1+ zg) = 3.2 × 10 G measured with a CRSF (Orlandini Since August 2009, MAXI (Matsuoka et al. 2009) on board the et al. 1998), the GL79 model was confirmed to accurately ex- International Space Station (ISS) has been continuously scan- plain the observed relation between P˙ and L of 4U 1626-67, ning the X-ray sky every 92 min of the ISS orbital period. We including both the spin-up and spin-down phases. This means derived the data of X Persei from the MAXI/GSC (Mihara et al. that accurate measurements of P˙ over a sufficiently wide range 2011; Sugizaki et al. 2011), using On-Demand System 2 pro- of L of a binary X-ray pulsar would conversely allow us to esti- vided by the MAXI team. This data-analysis scheme enables mate its B. us to extract images, light curves, and energy spectra of a given X Persei (4U 0352+309) is a high-mass X-ray binary sys- source by specifying its coordinate and observation times. The tem, consisting of a Be-type star and a magnetized NS. Its pulse on-source events and background events, both in 2–20 keV were period of P ∼ 835 s, which is relatively longer than those of derived using standard regions. Every 92 min, we thus obtained other pulsars, was independently discovered by Ariel 5 and a background-subtracted data set (an energy spectrum), with an Copernicus (White et al. 1976). The distance to the source is integration time of ∼ 60 s corresponding to a single scan transit estimated as D = 0.81 ± 0.04 kpc by its optical parallax de- of the source. Since this is sufficiently shorter than compared to rived from the GAIA Astrometry 1. With X-ray observations, the pulse period of X Persei (∼ 835 s), the data obtained with Delgado-Mart´ı et al. (2001) determined the orbital period of a single-scan is just a snap shot at a certain pulse phase. These X Persei as P = 250.3 ± 0.6 d, and derived an eccentricity orb data sets have been combined into a 2–20 keV light curve from of e ∼0.11 which is considerably smaller than those of typical MJD 55054 to MJD 58049 and 12 energy spectra each summed Be/X-ray binaries, e ≥0.3. over the 250-d orbital period. The 250-d averaged spectra have X Persei was spinning up at a rate of ˙ −4 −1 P ∼−1.5 × 10 yr statistics high enough to yield accurate fluxes. until 1978, when it turned into a phase of spin down at a rate of P˙ ∼ 1.3 × 10−4 yr−1 (Delgado-Mart´ı et al. 2001). In 2002, the pulsar returned to the spin-up phase, and has since been 2.2 RXTE/ASM spinning up until today. This spin up/down alternation implies The All-Sky monitor (ASM; Levine et al. 1996) on board the that the source is close to a torque equilibrium, and hence the RXTE satellite (Bradt et al. 1993) continuously monitored the long P suggests a strong magnetic field, because we then expect 1 6 whole X-ray sky from MJD 50087 to MJD 55927. For this B ∼ L 2 P 7 (Makishima 2016) entire period, we retrieved the RXTE/ASM dwell-time light In the 2–20 keV range, the spectrum of X Persei can be fitted curves of X Persei in four energy ranges, namely, 1.5–12.0 by a powerlaw model with an exponential cutoff (Di Salvo et al. keV, 1.5–3.0 keV, 3.0–5.0 keV, and 5.0–12.0 keV3. These light 1998), which is typical of accreting X-ray pulsars. However, in curves were used to obtain P and the fluxes. the 20–100 keV band, the spectrum of X Persei is much flatter and harder than those of other X-ray pulsars, and lacks steep cutoff (Sasano 2015). Sasano (2015) empirically quantified this 2.3 RXTE/PCA shape of the spectrum in comparison with those of other pul- In order to cross check the validity of the fluxes of RXTE/ASM sars (mostly with B measured), and obtained an estimate of of and MAXI/GSC, RXTE/PCA data were also used. Since the 13 B ∼ 10 G. Although Coburn et al. (2001) noticed a spectral PCA data are generally too short and sparse for the accurate dip at ∼ 29 keV in an RXTE/PCA spectrum, and regarded it as a determination of P and calculation of P˙ , these data were not 12 CRSF corresponding to B/(1+ zg) = 2.5 × 10 G, the feature used for pulse period analysis. Through 16 years RXTE mis- is too shallow and broad for that interpretation, and can be ex- sion, total 156 pointing observations were performed on this plained away by a combination of two continuum components source by the Proportional Counter Array (PCA; Jahoda et al. without a local feature (Doroshenko et al. 2012; Sasano 2015). 2 https://maxi.riken.jp/mxondem/ 1 https://gea.esac.esa.int/archive/ 3 https://xte.mit.edu/ASM lc.html/ Publications of the Astronomical Society of Japan, (2014), Vol. 00, No. 0 3 0.04 Flux(Crab) 0.02 0 52000 54000 56000 58000 MJD Fig. 1. Long term (1996–2017) light curves of X Persei in 25-d bin, expressed in Crab units. Black filled circles and red open squares are data with the RXTE/ASM (1.5–12 keV) and the MAXI/GSC (2–20 keV), respectively.