Arxiv:1012.2817V1

arXiv:1012.2817v1 [astro-ph.HE] 13 Dec 2010 o,A 6 G,Canada 2G7, T6G AB ton, ekly A97075,UA(-al [email protected] (e-mail: USA 94720-7450, CA Berkeley, slowl very 2005; from pulsations Zand long-period in’t or 2006) (SFXTs, al. th Transients et of Smith X-ray flaring X-ray Fast amplitude Supergiant high the including properties, 04 one l 07 ht ta.20) oemembers Some 2008). of Chaty & al. group (Filliatre et Chaty the wind object 2007; of compact stellar al. the et strong to Moon a due 2004; in is this embedded In that being HMXBs. appears obscured it called cases, commonly many are these and 2006), xesa eetn ore nte2–0kVbn,i only it 1 band, keV to 20–50 sources the While the in sources localizes 2010). detecting sources al. at et “IGR” Bird excels studied 2007; poorly al. previously et (Bodaghee or new of number (INTE- Laboratory GRAL) Astrophysics Gamma-Ray International A522 USA 52242, IA 18CACR/nvriéPrs7DnsDdrt,CASacla CEA Me des Diderot), L’Orme Denis 709, France Cedex, Gif-sur-Yvette 7 191 Bât. FR-91 Paris d’Astrophysique, DSM/IRFU/Service CEA/CNRS/Université 7158 MB Bdge ta.20)hv ag eeso intrin- of of levels dozens large the with have absorption of 2007) of sic Many proper- al. number the et (Bodaghee as large systems. HMXBs well these the as of is (HMXBs) ties Binaries surprise X-ray biggest Mass High the possibly 2008; fied, al. et Tomsick na- 2006; the 2009). al. Chaty et of & (Walter Tomsick optica determination sources Rodriguez, a secure the for obtain of allowing ture to counterparts, satellites IR X-ray or other with vations 2t tet e ok Y10760,USA 10027-6601, NY York, New Street, 120th 3 2 5 4 1 h adXryiaigo h aatcpaeb the by plane Galactic the of imaging X-ray hard The smr n oeo hs ore aebe identi- been have sources these of more and more As rpittpstuigL using typeset Preprint T eateto hsc,Uiest fAbra om28CEB, 238 Room Alberta, of University Physics, of Department PERI THE IN APPEAR O eateto hsc n srnm,Uiest fIw,Io Iowa, of University Astronomy, and Physics of Department oubaAtohsc aoaoy oubaUiest,5 University, Columbia Laboratory, Astrophysics Columbia pc cecsLbrtr,7GusWy nvriyo Calif of University Way, Gauss 7 Laboratory, Sciences Space I srpyiu ntuetto oéiain(UMR Modélisation Instrumentation Astrophysique - AIM OFRAINO G 03361 SAB -A IAYWT VER WITH BINARY X-RAY BE A AS J01363+6610 IGR OF CONFIRMATION aelt Wnlre l 03 a noee large a uncovered has 2003) al. et (Winkler satellite ujc headings: Subject J01363+6610. IGR of properties physical uigthe during h pe ii ntelmnst a vnlwr( lower even was luminosity the on limit upper the ntefil,mksteascainwt h esa eylikel 9 very star of Be luminosity the X-ray with an association the indicates makes field, the in 10 ossetwt adpwrlwwt htnidxof index photon a with power-law hard a with consistent h esa,w n aibesuc XM J013549.5+661243 (XMMU source variable a find we star, Be the rscm ergs n20 uy3 o 8k.AB trwspeiul ugse as suggested previously was star Be A although ks. and 28 for source, IGR 31 Binary X-ray July Mass 2009 High on candidate the containing field The 22 cm INTEGRAL − − 2 2 A A T J Chandra 9%cndneerr) hs eut,aogwt u optic our with along results, These errors). confidence (90% E s OHN N STROPHYSICAL tl mltajv 7/15/03 v. emulateapj style X − 1. 1 H 021 e,uasre)a hspsto with position this at unabsorbed) keV, (0.2–12 ∼ INTRODUCTION .T A. 10 ′ tr:nurn—Xry:sas—sas msinln,B b — Be emission-line, stars: — stars X-rays: — neutron stars: niiul(G J01363+6610) (IGR individual –5 23 bevto.W ics osbeipiain ftevr l very the of implications possible discuss We observation. OMSICK MB xii te extreme other exhibit HMXBs ′ –10 Chandra eurn olwu obser- follow-up requiring , J 24 OURNAL 1 cm C , uiga20 bevto,ddntdtc nXrysuc at source X-ray an detect not did observation, 2007 a during , J − RAIG EROME 2 eg,Wle tal. et Walter (e.g., . 1 H × EINKE oapa nteAtohsclJournal Astrophysical the in appear To R ley.edu) 10 INTEGRAL INTEGRAL ODRIGUEZ 31 2 0West 50 rss ergs aCity, wa J , Edmon- risiers, LUMINOSITY ornia, ULES ABSTRACT y, y e − 5 l , 1 < hsi oe hntpclfraB -a iay and binary, X-ray Be a for typical than lower is This . H AND 1 ALPERN . 4 Γ 04.Tepa u bevdfo h orews1 mcrab 17 was et source (Grebenev the from transient ( observed is flux peak source The the 2004). that detected not indicating was later, but 19 during April discovered 2004 on was observations source The (Corbet binary. star X-ray neutron the of period spin the the 1986). between correlation and the observed period the and bital is systems class cases, prop this remaining interesting of the 42 the of erties One in for in- 2009). Ziolkowski star unknown & pulsations is neutron (Belczynski the X-ray type a Of object systems, of compact presence binary periastron. the X-ray outbursts dicate approaches Be X-ray object known periodic tran- 64 compact to the cases, of lead the class most nature orbits when In binary eccentric binary the as demonstrates disk. system X-ray emis- emission circumstellar with X-ray a Be star sient from early-type from flux an the lines by optical sion dominated The in is systems 1982). HMXBs these Heuvel den van more & (Rappaport of tion 2007). al. et (Patel stars neutron rotating ca ( mcrab ad2wesltrwas later weeks 2 3- band The 2004). al. et (Grebenev n h orehsntbe eetdi other in detected been not has source the and evtosee huh23M of Ms 2.3 though even servations × r ubrti 04(rvnse l 07 ide l 2007, al. et Bird discov- 2007; the al. et during (Krivonos detected 2010). 2004 only in was outburst indi source ery catalogs and the 2006) that al. t cating et of 2010) Hartog studies al. (den detailed et region both (Bird Cassiopeia include source non-detection this of of Reports position the at accumulated a pcrsoy n thsbe ugse ob h likely the be to suggested been has it and spectroscopy, cal ag ro ein esa a on ihnteerror the within found was star H Be narrow-band a using region, circle error large a A 1 = ∼ RASH G 03361 a entnaieycasfida Be a as classified tentatively been has J01363+6610 IGR INTEGRAL lhuhthe Although 10 2 ,adte2kcdsac siaefrteB star Be the for estimate distance kpc 2 the and y, . . 4 6 31 XMM-Newton 3 P , ± × 03361 a bevdby observed was J01363+6610 rss ergs B ∼ 0 HILIP 10 h osbecutrato the of counterpart possible the . ODAGHEE 9 n oundniyof density column a and 3 iha vrg -a u f2 of flux X-ray average an with ) . − 1 10 − × livsiaino te -a sources X-ray other of investigation al K 1 rscm ergs bevtoshv loldt h addi- the to led also have observations INTEGRAL AARET 10 suigtesm pcrlmodel) spectral same the assuming − 1 h pcrmo hssuc is source this of spectrum The . 11 wqisetlmnst o the for luminosity quiescent ow rscm ergs 4 − O UECN X-RAY QUIESCENT LOW Y S , 2 < akhl hsc stars: — physics hole lack s YLVAIN − 1mrb(rbnve l 2004), al. et (Grebenev mcrab 11 1 oiinucranyo 3 of uncertainty position α nte1–5kVbn n 9 and band keV 17–45 the in ) σ − 2 mgn n olwu opti- follow-up and imaging pe ii nte1–5keV 17–45 the in limit upper s C INTEGRAL − HATY 1 N nte81 e band keV 8–15 the in ) h oiinof position the H XMM-Newton (1 = INTEGRAL 5 , . × 5 INTEGRAL − + iehv been have time 0 0 10 . . 5 7 ) − INTEGRAL × 13 ∼ ′ . weeks 2 leaves 7 ob- or- al. he - - . 2 Tomsick et al. counterpart (Reig et al. 2004, 2005). However, a sensitive X- (2008), and the procedures used for data reduction are de- ray observation taken with the Chandra X-ray Observatory in scribed in that work. 2007 failed to detect the Be star (Tomsick et al. 2008). At the 2 kpc distance estimated for the Be star (Reig et al. 2005), the 2.3. MDM Optical Observations non-detection implies an upper limit on the X-ray luminosity of <2×1031 ergs cm−2 s−1 (Tomsick et al. 2008), which is Optical spectroscopic observations were made with the lower than quiescent luminosities for other Be X-ray binaries Boller & Chivens CCD spectrograph on the 2.4 m Hiltner (Campana et al. 2002). This luminosity approaches the level Telescope of the MDM Observatory on 2009 August 23 UT. that has been seen during quiescent periods from transient Conditions were clear, but with high humidity. The 150 g/mm ′′ Low Mass X-ray Binaries (LMXBs) due to thermal emission grating, blazed at 4700 Å, was used with a 1 wide slit, which from the neutron star surface (Brown, Bildsten & Rutledge gave a wavelength coverage of 3700−7360 Å at 7.6 Å reso- 1998; Campanaet al. 1998). Although neutron stars in lution. Spectra of six stars coinciding with X-ray sources in HMXBs might not be heated to the high levels seen for the field of IGR J01363+6610 were obtained, each with an LMXBs, it is interesting that these observations probe this lu- exposure time of 120 s. Of these, three are inside the refined minosity regime. INTEGRAL error circle and are presented here. Spectral re- Currently, there are significant uncertainties about the na- duction was performedusing standard IRAF procedures. Flux ture of IGR J01363+6610. While we know that it is an X-ray calibration used Oke & Gunn (1983) standard stars, although transient, and the hard X-ray emission makes it likely that it is the spectrophotometry is not expected to be accurate because a binary (although its orbital period is unknown), the Chandra of the narrow slit. non-detection makes it unclear whether it is really a Be X-ray binary. Finding a Be star in the relatively large INTEGRAL 3. RESULTS error circle is not convincing because Be stars are more com- 3.1. X-ray Sources in the INTEGRAL Error Circle monly found as single stars rather than being part of a binary system (Porter & Rivinius 2003). Furthermore, we cannot be For the XMM-Newton analysis, we used the SAS tool certain that the compact object in the system is a neutron star edetect_chain to search for sources in the ∼30′ diame- since pulsations have not been detected. ter field-of-view (FOV) of the MOS1 and MOS2 instruments, In this paper, we report on a second sensitive X-ray ob- and included photons in the 0.1–10 keV energy band. All servation of the IGR J01363+6610 field with XMM-Newton seven of the MOS2 CCDs were operational, but, for MOS1, along with optical spectroscopy of the Be star as well as other CCD#6 was not active because it was switched off after an X-ray sources with optical counterparts in the field. With anomaly that occurred in 2005. While this means that some XMM-Newton, we confirm that the Be star is an X-ray source. of the MOS2 FOV is not covered by MOS1, the central CCD, We also re-analyze the Chandra observation and discuss the which includes the entire INTEGRAL error circle, is covered 2007 results in the contextof the new informationfrom XMM- by both MOS units. Newton. We found a total of 21 sources, including 14 detected by both MOS units and seven detected by MOS2 in the part 2. OBSERVATIONS of the FOV that MOS1 did not cover. We determined the ′′ 2.1. XMM-Newton number of counts for each source using 25 radius apertures and subtracted the background using circular apertures of the We observed the IGR J01363+6610 field with XMM- same size in 16 source-free regions of the detector. The Newton on 2009 July 31 from 14.4 h to 22.3 h UT. The obser- mean number of background counts per aperture is 47.1±1.8 vation (ObsID 0603850101) occurred during XMM-Newton for MOS1 and 45.9 ± 1.7 for MOS2 (in 28,180 s of expo- revolution 1766. The EPIC pn, MOS1, and MOS2 instru- sure time). The sources range in brightness from 12 ± 8 to ments (Strüder et al. 2001; Turner et al. 2001) were all op- 156 ± 13 counts (both of these sources are in the region cov- erated in Full Frame mode with a medium filter. We used ered only by MOS2). Five of the sources are in the 3′.7 IN- the XMM-Newton Science Analysis Software (SAS-10.0.0) TEGRAL 90% confidence error circle for IGR J01363+6610 package to process the raw data files. We used the SAS tool given in Bird et al. (2010). These five sources include the sec- cifbuild to obtain the necessary “current calibration files” ond brightest source in the field (132 ± 9 counts) as well as for the observation. Then, we reprocessed the pn, MOS1, and sources with 74 ± 8, 31 ± 6, 25 ± 6, and 19 ± 6 counts, where MOS2 data using epproc and emproc, yielding photon these numbers of counts are averages of the two MOS detec- event lists for the three instruments. Although proton flares tors. sometimes cause portions of XMM-Newton observations to The XMM-Newton source names and positions are given in have very high backgrounds,we did not find evidence for pro- Table 1. In addition to MOS1 and MOS2, we also determined ton flares, and we were able to use the full exposure time. the position measured by the pn instrument for each source, and we report the average position measured by the three in- 2.2. Chandra X-ray Observatory struments (except for XMMU J013632.4+660924, which fell We also used an observation of the IGR J01363+6610 field between two pn CCD chips). In averaging the positions, we made by the Chandra X-ray Observatory. The observation weighted each measurement by its statistical uncertainty. In (ObsID 7533), which had an exposure time of 4,976 s, was Table 1, we also report the overall 90% confidence position made on 2007 June 8. We used the Advanced CCD Imaging uncertainties, and, in each case, the error is dominated by the ′′ Spectrometer (ACIS; Garmire et al. 2003), and the aimpoint systematic pointing uncertainty of 3 .46 was placed on the ACIS-S3 chip. The 90% confidence IN- 6 TEGRAL error circle for IGR J01363+6610 (Bird et al. 2010) In the document entitled “EPIC status of calibration and data analysis” (XMM-SOC-CAL-TN-0018), Guainazzi et al. report an rms value for the is contained on the ACIS-S3 and ACIS-S4 chips. Results systematic pointing uncertainty of 2′′.0, and we have multiplied this by 1.7 to of the observation were previously reported in Tomsick et al. obtain the 90% confidence value. TheBeX-raybinaryIGRJ01363+6610inquiescence 3

TABLE 1 XMM-Newton SOURCESINTHE INTEGRAL ERROR CIRCLE

Name R.A. (J2000)a Decl. (J2000)a Position MOS Uncertaintyb countsc XMMU J013549.5+661243 01h35m49s.53 +66◦12′43′′.1 ±3′′.4 132 ± 9 XMMU J013606.5+661304 01h36m06s.54 +66◦13′04′′.5 ±3′′.9 25 ± 6 XMMU J013620.8+660851 01h36m20s.80 +66◦08′51′′.0 ±3′′.7 31 ± 6 XMMU J013632.4+660924d 01h36m32s.48 +66◦09′24′′.0 ±4′′.1 19 ± 6 XMMU J013644.2+661302 01h36m44s.26 +66◦13′02′′.2 ±3′′.5 74 ± 8 aThe position is the weighted average of the positions measured by the three XMM-Newton instruments. bThe 90% conﬁdence uncertainty in the XMM-Newton position. This includes a systematic contribution of 3′′.4 due to the absolute pointing uncertainty and a statistical contribution. We have added the two contributions in quadrature. cThe average of the 0.1–10 keV count rates measured by MOS1 and MOS2. dThis source is detected by MOS1 and MOS2, but not by the pn instrument. It falls on a gap between pn CCD chips.

Using the 2007 observation of the field by Chandra, we N previously reported three sources in the INTEGRAL error circle (Tomsick et al. 2008). However, in the most recent IN- TEGRAL source catalog (Bird et al. 2010), the best estimate of the IGR J01363+6610 position has shifted by 1′.8 rela- tive to the value used in Tomsick et al. (2008), and a fourth 7 2 Chandra source, CXOU J013644.5+661301,is now also con- 1 sistent with the INTEGRAL position for IGR J01363+6610, 3 and its Chandra position is R.A. (J2000) = 01h36m44s 54, . E Decl. (J2000) = +66◦13′01′′.6 with a 90% confidence uncertainty of 0′′.64. Of the four Chandra sources, two are coincident with XMM-Newton sources. Thus, merging the 2007 Chandra source list with the list of 2009 XMM-Newton leaves a list of seven X-ray sources in the INTEGRAL error circle, 6 5 and these sources are listed in Table 2. 4 3.2. Optical and IR Counterparts and Optical Spectroscopy For the seven X-ray sources listed in order of R.A. in Ta- ble 2, we searched for optical and infrared (IR) counterparts in the United States Naval Observatory (USNO-B1.0) and 2 Mi- FIG. 1.— Red optical image from the Digitized Sky Survey. The large circle cron All-Sky Survey (2MASS) catalogs. Four of the sources is the 3′.7 INTEGRAL error circle from Bird et al. (2010). The seven X-ray (designated as #1, #2, #5, and #7) have optical and IR coun- sources detected by XMM-Newton and Chandra are labeled. The arrows point terparts, and the names of these counterparts and their magni- to the optical counterparts for four of the X-ray sources, and their optical and tudes are given in Table 2. For the other three X-ray sources IR magnitudes are given in Table 2. The small circles mark the locations of (#3, #4, and #6), the nearest USNO and 2MASS sources are the X-ray sources without known optical counterparts. >7′′ away, which indicates that these X-ray sources are not associated with any sources in the USNO or 2MASS catalogs. Figure 1 shows a red (close to R-band) optical image from the X-ray source #2 is coincident with USNO-B1.0 1562- Digitized Sky Survey with the revised INTEGRAL error circle 0030364, which is also known as the bright (V = 11.5) op- and the seven X-ray sources labeled. tical source TYC 4043-860-1. There is some confusion about X-ray source #1 is the Be star that was previously iden- whether this is an emission line star based on a catalog of stars tified using optical imaging and spectroscopy taken in 2004 with Hα in emission (González & González 1956). However, (Reig et al. 2004, 2005). The optical spectrum of this star Hα images and optical spectra taken, respectively, 2 months (USNO-B1.0 1562-0030282) is shown in Figure 2 and has and 5 months after the X-ray outburst from IGR J01363+6610 the blue continuum and the Hα and Hβ emission lines in- do not show evidence for Hα emission and the spectrum dicative of a Be star. The equivalent width (EW) of Hα is shows Hα in absorption(Reig et al. 2004). The XMM-Newton –54 ± 3 Å, which is consistent with the value measured by observation provides the first evidence that this is also an X- Reig et al. (2004), suggesting that the Be star’s circumstellar ray source, but the lack of an Hα emission line in 2004 makes disk is stable. Interstellar absorption features in the spectrum it very unlikely to be the correct IGR J01363+6610 counter- can be used to estimate the extinction E(B −V) following the part. correlations in Herbig (1975). Most commonly, the 4430 Å We also obtained optical spectra for the only other two diffuse interstellar band is used for this purpose, but it falls X-ray sources in the INTEGRAL error circle with opti- in a poorly exposed region of our spectrum. Instead we use cal/IR counterparts. The spectra for sources #5 and #7 the 5780 Å feature. With an EW of 1.0 Å, it corresponds to are shown in Figure 3, and we identify them as G-type E(B−V) in the range 1.5−2.0, which is consistent with the es- and early M-type stars, respectively. Like Be X-ray bina- timate of Reig et al. (2005) based on the spectral classification ries, LMXBs usually show Hα in emission during outbursts and photometry, E(B −V)=1.6. (Fender et al. 2009) and during quiescence (Orosz et al. 2002; 4 Tomsick et al.

TABLE 2 OPTICAL/INFRARED IDENTIFICATIONS

Catalog/Sourcea Separationb Magnitudes (1) XMMU J013549.5+661243 USNO-B1.0 1562-0030282 2′′.0 ± 3′′.4 B = 14.9 ± 0.3 R = 12.4 ± 0.3 I = 11.1 ± 0.3 ′′ ′′ 2MASS J01354986+6612433 2 .0 ± 3 .4 J = 10.04 ± 0.02 H = 9.57 ± 0.03 Ks = 9.12 ± 0.02 (2) XMMU J013606.5+661304 USNO-B1.0 1562-0030364 2′′.5 ± 3′′.9 B = 12.1 ± 0.3 R = 10.7 ± 0.3 I = 10.1 ± 0.3 ′′ ′′ 2MASS J01360684+6613021 3 .0 ± 3 .9 J = 10.42 ± 0.02 H = 10.10 ± 0.03 Ks = 10.00 ± 0.02 (3) CXOU J013609.9+661157 The closest USNO-B1.0 source is 9′′.1 ± 0′′.6 away. The closest 2MASS source is 9′′.0 ± 0′′.6 away. (4) XMMU J013620.8+660851 The closest USNO-B1.0 source is 7′′.2 ± 3′′.7 away. The closest 2MASS source is 8′′.7 ± 3′′.7 away. (5) CXOU J013621.2+660928 USNO-B1.0 1561-0031113 0′′.52 ± 0′′.64 B = 17.2 ± 0.3 R = 13.6 ± 0.3 I = 12.0 ± 0.3 ′′ ′′ 2MASS J01362115+6609286 0 .38 ± 0 .64 J = 10.42 ± 0.03 H = 9.54 ± 0.03 Ks = 9.26 ± 0.02 (6) XMMU J013632.4+660924/CXOU J013632.8+660924 The closest USNO-B1.0 source is 9′′.1 ± 0′′.6 away. The closest 2MASS source is 7′′.5 ± 0′′.6 away. (7) XMMU J013644.2+661302/CXOU J013644.5+661301 USNO-B1.0 1562-0030589 0′′.73 ± 0′′.64 B = 17.4 ± 0.3 R = 14.4 ± 0.3 I = 13.1 ± 0.3 ′′ ′′ 2MASS J01364458+6613014 0 .33 ± 0 .64 J = 12.77 ± 0.02 H = 12.01 ± 0.03 Ks = 11.83 ± 0.03 aThe catalogs are the 2 Micron All-Sky Survey (2MASS) and the United States Naval Observatory (USNO-B1.0). bThe uncertainty on the separation is the addition (in quadrature) of the X-ray source position error and the 2MASS or USNO-B1.0 error.

FIG. 2.— Optical spectrum of the Be star XMMU J013549.5+661243 taken on 2009 August 23, which is a few weeks after the XMM-Newton observation. The optical spectrum was taken with the MDM 2.4 meter telescope and shows strong Hα and Hβ lines, indicating the presence of a circumstellar disk. The lower spectrum is the upper spectrum divided by 6, allowing for all of the Hα line to be visible. This source was previously suggested as the most likely counterpart to IGR J01363+6610, and the detection of the Be star with XMM-Newton provides further conﬁrmation of this association.

FIG. 3.— Optical spectra of two of the X-ray sources in the INTEGRAL error circle for IGR J01363+6610. The stars are listed as source #5 (top Charles & Coe 2006). Also, symbiotic systems often show panel) and source #7 (bottom panel) in Table 2. We identify source #5 as a strong Hα in emission (Chakrabarty & Roche 1997). How- highly reddened G-type star and source #7 as an early M-type star. ever, there are counter-examples for both LMXBs and symbi- otics. In any case, as sources #5 and #7 do not have emission lines, there is no reason to consider that they might be the we focus exclusively on the X-ray properties of the Be star IGR J01363+6610 counterpart. XMMU J013549.5+661243 in the remainder of this paper. Thus, these observations show that the Be star that was previously considered to be the likely coun- 3.3. XMMU J013549.5+661243 terpart of IGR J01363+6610 is coincident with 3.3.1. Energy Spectrum XMMU J013549.5+661243, which was the brightest X- We used the SAS tool xmmselect to produce MOS1, ray source in the INTEGRAL error circle during the 2009 MOS2, and pn energy spectra for XMMU J013549.5+661243 XMM-Newton observation. The fact that we have now and considered the recommendations given in a recent EPIC conclusively shown that the Be star is a transient X-ray (MOS and pn) calibration document7 for event ﬁltering and source and that none of the other fainter X-ray sources energy ranges. For MOS, we used event ﬁltering with the in the INTEGRAL error circle have optical properties that would be expected of a hard X-ray source strengthens the 7 The document entitled “EPIC status of calibration and data analysis” by association between the Be star and IGR J01363+6610, and Guainazzi et al. (XMM-SOC-CAL-TN-0018) is based on results obtained TheBeX-raybinaryIGRJ01363+6610inquiescence 5

(a) XMM/MOS1, 2009, 0.1-10 keV (b) Chandra/ACIS, 2007, 0.3-10 keV

60 arcsec 60 arcsec

FIG. 4.— (a) X-ray image from the 2009 July 31 XMM-Newton observation showing the detection of XMMU J013549.5+661243, which is the very likely counterpart to IGR J01363+6610. The image uses 0.1–10 keV photons from the MOS1 instrument and was rebinned so that each pixel is 4′′ wide. The inner circle is the source region used for making the spectrum and the light curve. The annulus (outer two circles) is the region used to estimate the background. (b) X-ray image from a 2007 Chandra observation. The image uses 0.3–10 keV photons from the ACIS instrument and was rebinned so that the pixels are 1′′. The circle shows the 90% conﬁdence XMM-Newton error circle, highlighting the fact that IGR J01363+6610 was not detected. In both images, North is up and East is to the left.

expression “#XMMEA_EM && PATTERN<=12” and included events within an aperture with a 25′′ radius. A back- 10-3 ground spectrum was extracted from an annulus centered on the source, and the source and background regions are shown -4 in Figure 4a. For the pn instrument, we used event ﬁlter- 10 )

ing with the expression “FLAG=0 && PATTERN<=4” and -1 s extracted spectra from source and background regions. For -2 -5 all three instruments, we used the SAS tools rmfgen and 10 arfgen to make response matrices. For the MOS detectors, (keV cm

we used the 0.1–10 keV bandpass, and for the pn detector, we E 10-6

used the 0.2–12 keV bandpass. We rebinned each spectrum, E F and we ﬁtted the spectra using χ2 statistics.

We used the XSPEC version 12 software for spectral fit- -7 ting, and we tried several models. The results of using an ab- 10 sorbed power-law are shown in Table 3 and Figure 5. For absorption, we used the photoelectric absorption cross sections 0.1 1 10 from Balucinska-Church & McCammon (1992) and elemen- Energy (keV) tal abundances from Wilms, Allen & McCray (2000), which correspond to the estimated abundances for the interstellar FIG. 5.— XMM-Newton energy spectrum for XMMU J013549.5+661243 medium. This model gives a good fit (χ2/ν = 23.1/26), fitted with an absorbed power-law model. The model parameters are given in Table 3, indicating a relatively hard spectrum (Γ = 1.4 ± 0.3) but not a highly and requires a relatively hard power-law photon index of absorbed spectrum. The pn data are shown in black, and the MOS data are Γ = 1.4 ± 0.3. The power-law fit gives a 0.2–12 keV un- shown in grey with diamonds plotted on MOS1 points and squares plotted on +0.3 × −13 −2 −1 MOS2 points. absorbed flux of (1.9−0.2) 10 ergs cm s and a col- +0.7 × 22 −2 umn density of NH = (1.5−0.5) 10 cm , which is some- what higher than the value through the Galaxy along the line 21 −2 of sight, NH = 5.2 × 10 cm (Kalberla et al. 2005). Al- though this could indicate a small amount of absorption lo- which is consistent with the column density determined from cal to the source, it is much lower than the values near 1023 the X-ray spectrum. cm−2 seen for the sources usually considered to be obscured We also fitted the energy spectrum with blackbodyand ther- HMXBs. Another argumentagainst local absorption of the X- mal Bremsstrahlung models, and the parameters are given in ray source is that the value of E(B −V)=1.5–2.0 corresponds Table 3. A fit with a 1.3 keV blackbody is similar in qual- 22 −2 2 to an NH in the range (1.0–1.4)×10 cm (Ryter 1996), ity to the power-law (χ /ν = 21.5/26), but the spectrum does not allow us to distinguish between thermal and non-thermal using SAS-10.0.0. It was released on 2010 July 16 and can be found at models. However, at 2 kpc, the blackbody model implies an http://xmm2.esac.esa.int/external/xmm_sw_cal/calib/index.shtml. emitting region with a radius of only 1.15 × 103 cm, which is 6 Tomsick et al. too small to be physical. The Bremsstrahlung model requires a relatively high temperature (>11 keV), which makes it very 0.08 a similar to a power-law over the XMM-Newton bandpass. Al- 0.06 though we are not able to distinguish between models, it is important to note that the power-law fit above shows that the 0.04 spectrum is hard, as expected for an IGR source and espe- 0.02

cially for an HMXB. Source Rate (c/s) 0.00 3.3.2. Long- and Short-Term X-ray Variability -0.02 0.08 b Comparing the two panels of Figure 4 shows the difference between ﬂuxes in 2009 (the XMM-Newton observation) and 0.06 2007 (the Chandra observation). We did not detect the Be 0.04 star during a 5 ks Chandra/ACIS observation (Tomsick et al. 0.02

2008), and 1 count is detected by ACIS within the XMM- 0.00 Newton error circle. Thus, using Poisson statistics, this cor- Background Rate (c/s) -0.02 responds to a 90% confidence upper limit of <3.9 counts 0 5000 10000 15000 20000 25000 (Gehrels 1986). For the power-law spectral model mea- Time (s) 22 −2 sured by XMM-Newton (Γ =1.4 and NH =1.5 × 10 cm ), FIG. 6.— Source (a) and background (b) XMM-Newton light curves of this corresponds to an absorbed 0.3–10 keV flux of <1.6 × − − − − XMMU J013549.5+661243. The rates are a combination of counts from all 10 14 ergs cm 2 s 1 and an unabsorbed flux of <2.9 × 10 14 three instruments (pn, MOS1, and MOS2), and the background rates have ergs cm−2 s−1, which is >6.6 times lower than the flux seen been normalized to the rates predicted in the source region. The energy range during the XMM-Newton observation. During the original in both panels is 0.2–12 keV, and the time resolution is 250 s. The source shows variability on this time scale. outburst detected by INTEGRAL in 2004, an 8–15 keV flux of 9 mcrab (9.1 × 10−11 ergs cm−2 s−1) was measured using the JEM-X instrument. In the 8–15 keV band, the flux measured × −14 −2 −1 with XMM-Newton in 2009 was 7.9 10 ergs cm s , spectra. For example, there are several γ Cas-like systems that which is a factor of 1,150 lower than the outburst level. Thus, may harbor white dwarfs and have power-law spectral com- the quiescent upper limit obtained with Chandra represents a ponents with photon indices of 1.4–1.7 (Motch et al. 2007). flux >7,600 times lower than the level seen during outburst. However, based on our study of the X-ray sources in the We also find evidence for shorter term variability during IGR J01363+6610 field, we argue that it is very likely that the XMM-Newton observation. Figure 6 shows the source and XMMU J013549.5+661243 is the quiescent counterpart to background rates combined for all three XMM-Newton detec- IGR J01363+6610. In this case, the outburst detected by IN- tors in the 0.2–12 keV band with 250 s time bins. The light TEGRAL indicates that the source can produce an X-ray lu- curve shows that the 0.2–12 keV rates can change from values minosity as high as ∼1035 ergs s−1 along with variations in ± ± 2 as low as –0.003 0.007to 0.061 0.018. We used a χ test X-ray luminosity by a factor of >7,600 (when the Chandra to determine the significance of the variability. The best fit observation is considered as well). In contrast, the γ Cas-like 2 that can be obtained with a constant rate has a χ = 142.5 for systems have luminosities typically near 1033 ergs s−1, and 106 degrees of freedom, which indicates that the variability is while they show variability in flux by factors of a few, they significant at the 99% confidence level. do not have large outbursts (Motch et al. 2007). Thus, the most likely interpretation of the available information is that 4. DISCUSSION IGR J01363+6610 is a Be X-ray binary with a neutron star 4.1. The Nature of IGR J01363+6610 (or perhaps a black hole as discussed below) accretor. In the Our results show that XMMU J013549.5+661243 is a following, we discuss IGR J01363+6610 in this context. hard X-ray source coincident with a Be star at a distance of ∼2 kpc. At this distance, the inferred X-ray luminos- 4.2. Accretion Regimes and the X-ray Luminosity ity measured by XMM-Newton is 9.1 × 1031 ergs s−1. While Models for X-ray emission from Be X-ray binaries this would be a relatively high X-ray luminosity for an iso- are based on the picture of a neutron star with a rel- lated B1 star, it is not too far above the values reported by atively strong magnetic field (B ∼ 1012 G) and a rela- Cohen, Cassinelli & Macfarlane (1997) for such stars. How- tively slow rotation speed (Pspin ∼ 0.1–1,000 s) periodically ever, in nearly all cases of isolated B-type stars, the X- accreting from the circumstellar disk around the Be star ray emission is thermal with temperatures that are typically (Rappaport & van den Heuvel 1982; Corbet 1986). In most not above a few ×106 K (Cohen, Cassinelli & Macfarlane systems, the orbits have non-zero eccentricity, and “Type 1” 1997; Cohen 2000), and their spectra fall very steeply above X-ray outbursts occur when the neutron star makes its closest ∼0.25 keV. The very hard X-ray spectrum that we measure approach to the star so that the outburst periodicity is equal to for XMMU J013549.5+661243 is inconsistent with a thermal the orbital period, Porb. spectrum with such a low temperature, indicating that the Be Within this physical picture, there are two accretion regimes star must have a compact binary companion that is emitting that are usually considered. One regime is at higher mass ac- most of the X-ray emission. cretion rates when the accretion pressure overcomes the cen- While the X-ray luminosity and spectrum are not consis- trifugal magnetic field barrier, and matter is accreted directly tent with an isolated star, the XMM-Newton data alone are onto the poles of the neutron star (Corbet 1996). The sec- not sufficient to determine the nature of the compact binary ond regime is during times when the accretion is inhibited by companion. There are known Be systems with both white the magnetic field. This is often called the propeller regime dwarf and neutron star companions that can produce hard (Illarionov & Sunyaev 1975), and the X-ray emission is much TheBeX-raybinaryIGRJ01363+6610inquiescence 7

TABLE 3 SPECTRAL RESULTS FOR XMMU J013549.5+661243

Parameter Valuea Power-law model +0.7 × 22 −2 NH (1.5−0.5) 10 cm Γ 1.4 ± 0.3 +0.3 × −13 −2 −1 FPL (unabsorbed, 0.2–12 keV) (1.9−0.2) 10 ergs cm s χ2/ν 23.1/26 Blackbody model +0.3 × 22 −2 NH (0.3−0.2) 10 cm +0.16 kT 1.33−0.15 keV 2 2 b +1.6 × −3 Rkm/d10 (3.3−1.1) 10 χ2/ν 21.5/26 Bremsstrahlung model +0.5 × 22 −2 NH (1.5−0.4) 10 cm kT >12 keV c +1.5 × −5 −5 Nbremss (2.9−0.4) 10 cm χ2/ν 22.6/26 aThe errors on the parameters are for ∆χ2 = 2.7, corresponding to 90% confidence for one parameter of interest. b The normalization for this model is parameterized in terms of the radius of a spherical emitting region in units of km (Rkm) and the distance in units of 10 kpc (d10). c −16 −2 −5 The normalization for this model is 2.4 × 10 d R nenidV cm where d is the source distance, and ne and ni are, respectively, the electron and ion number densities within the volume V. lower in this regime because the matter being accreted only neutron star would be thermal in origin and would be ex- reaches the neutron star magnetosphere rather than falling tremely soft. Such components are seen at temperatures near onto the neutron star surface. A gap is predicted between the 50–150 eV in LMXBs (Brown, Bildsten & Rutledge 1998; minimum luminosity in which the system can be in the di- Campana et al. 1998; Degenaar et al. 2009) and although it is rect accretion regime and the maximum luminosity in which unclear whether this regime occurs for Be X-ray binaries, it the system can be in the propeller regime (Corbet 1996), and has probably been seen for the HMXB Supergiant Fast X-ray we define ∆ to be equal to this luminosity ratio. For typical Transient (SFXT) IGR J17544–2619 (in’t Zand 2005). Dur- 12 neutron star properties (1.4 M⊙ mass, 10 km radius, 10 G ing the XMM-Newton observation of IGR J01363+6610, we 2/3 magnetic field), ∆ = 170(Pspin/1 s) (Campana et al. 2002), are not seeing this regime since we measure a hard spectrum, and values of ∆ range from ∼100 for fast rotators to ∼10,000 which is either non-thermal or has a much higher temperature for slow rotators (Corbet 1996). (kT =1.3 keV) than would be expected from a non-accreting Assuming a distance of 2 kpc and a spectral shape consis- neutron star. Since we did not detect the source during the tent with the power-law model measured using XMM-Newton, Chandra observation, the spectrum is not constrained. As- 22 suming the same NH as seen during outburst (NH =1.6 × 10 IGR J01363+6610 has been observed by INTEGRAL at Lx = − 1.04 × 1035 ergs s−1 in 2004, by Chandra at <1.4 × 1031 cm 2), we can put limits on the flux and the temperature of ergs s−1 in 2007, and by XMM-Newton at 9.1 × 1031 ergs s−1 the neutron star by requiring an ACIS count rate below the in 2009. These luminosities are unabsorbed values measured upper limit given above. We use the magnetic neutron star in or extrapolated into the 0.2–12 keV band. The 2004 out- model nsa (Pavlov et al. 1995), assuming emission from the burst luminosity is at the lower end of typical values for nor- entire neutron star, a mass of 1.4M⊙, a radius of 10 km, a mal (“type I”) Be X-ray binary outbursts, which are in the distance of 2 kpc, and a magnetic field of B = 1012 G. The 1035 to 1037 ergs s−1 range (Stella, White & Rosner 1986; upper limit on the (unredshifted) surface temperature is then Wilson et al. 2005). Still, it is likely (but perhaps not cer- 80 eV, and the unabsorbed 0.2–12 keV luminosity upper limit − tain if the neutron star spin or magnetic field strength have is <2.4×1032 ergs s 1. Changes in the assumptions will affect extreme values) that the source reached a high enough ac- this, particularly changes in the assumed NH or distance. For 22 −2 cretion rate to enter the direct accretion regime. Following instance, a choice of NH =2 × 10 cm gives kT < 85 eV 32 −1 Corbet (1996) and Campana et al. (2002), we derive an upper and Lx < 3.2 × 10 ergs s , or a distance of 3 kpc gives 32 −1 limit on the luminosity produced by magnetospheric emission kT < 88 eV and Lx < 3.7 × 10 ergs s . If the X-ray emis- 35 32 −2/3 of Lm < 1.04 × 10 /∆ =6.1 × 10 (Pspin/1 s) . Thus, un- sion in this quiescent state is dominated by hot spots at the der the assumptions that the neutron star properties (mass, ra- polar caps, their temperature may be higher, but the tempera- dius, and magnetic field strength) are typical and the system ture of the rest of the neutron star surface must be even lower. reached the direct accretion regime during the 2004 outburst, 4.3. The Quiescent Luminosity and Possible Implications the luminosity measured by XMM-Newton in 2009 implies a neutron star spin period less than 17 s for IGR J01363+6610. The above analysis suggests that the IGR J01363+6610 While the propeller regime may extend to very low lumi- luminosities can be explained within the standard picture nosities, a third regime to consider is when the mass accre- for Be-neutron star X-ray binaries, but the luminosity up- tion rate onto the neutron star drops to zero. In this case, per limit that we infer from the Chandra observation is any X-ray emission produced could have contributions from lower than has been previously reported from these sys- the Be star and from the neutron star. Emission from the tems. Campanaetal. (2002) present a study that fo- cuses on quiescent X-ray observations of 3 Be X-ray bina- 8 Tomsick et al. ries, and their luminosities are (1–3)×1035 ergs s−1 (0.1– cases, the Hα lines were significantly weaker than measured 10 keV) for A 0538–66, (0.8–2)×1033 ergs s−1 (0.5–10 keV) for IGR J01363+6610. There is no precedent that we are for 4U 0115+63, and 5 × 1032 ergs s−1 (0.5–10 keV) for aware of for dissipation and re-formation of a circumstellar V 0332+53. In addition, quiescent luminosities of (3– disk in a system with an Hα EW as large as the one mea- 9)×1033 ergs s−1 and (2–4.5)×1033 ergs s−1 have been re- sured for IGR J01363+6610. While there are possible expla- ported for GRO J2058+42 and A 0535+26, respectively nations for the luminosity upper limit measured over a short (Wilson et al. 2005; Negueruela et al. 2000). The candidate period of time by Chandra, these explanations do not apply Be X-ray binary XTE J1829–098 has been observed at a lu- to the XMM-Newton observation. The luminosity measured minosity of 3×1032 (d/10 kpc)2 ergs s−1 (Halpern & Gotthelf by XMM-Newton would also make it one of, if not the, lowest 2007), which is low for a Be X-ray binary, but still not as luminosity Be X-ray binaries. Thus, it is likely that the low low as we find for IGR J01363+6610. It should, however, luminosities measured for IGR J01363+6610 require that this be noted that, for A 0535+26, Negueruela et al. (2000) found system has some unusual physical properties. the source at this luminosity after it had previously been re- The Be X-ray binaries with known compact object type ported to have a quiescent luminosity two orders of magnitude have neutron stars, but binary evolution models predict higher (Motch et al. 1991). Thus, it is clear that even when that between zero and two of the 64 known Be X-ray bi- they are not in outburst, these sources exhibit a large amount nary systems harbor a black hole instead of a neutron star of variability, and many observationsmay be requiredfor each (Belczynski & Ziolkowski 2009). The observational signa- source to define a quiescent luminosity. tures for Be X-ray binaries with black holes are not as clear as While the Chandra upper limit of <1.4 × 1031 (d/2 kpc)2 neutron star signatures, such as pulsations. At high mass ac- ergs s−1 obtained with a hard spectrum as seen in outburst is cretion rates, one would expect the X-ray emission from black significantly lower than has been found for other Be X-ray holes to be softer and brighter due to the presence of an inner binaries, this conclusion does depend on the validity of the accretion disk that cannot form in systems with highly mag- distance determination. The 2 kpc distance from Reig et al. netized neutron stars due to the magnetosphere. At moderate (2005) depends on the measurement of the V-band magni- to low mass accretion rates, accreting black holes typically tude, which is given in Reig et al. (2005) to high precision have hard spectra (Remillard & McClintock 2006). At very (V = 13.29 ± 0.02, based on three measurements), the extinc- low mass accretion rates, black hole systems would likely tion, and the spectral type. Reig et al. (2005) give a value of be fainter than neutron star systems because the latter exhibit E(B −V)=1.6, and our determination is consistent with this magnetospheric and surface emission, which are both absent suggesting that the uncertainty in the distance due to the mea- from black holes. The expected hard spectrum and the faint surementof extinction is not large. Probably the largest uncer- quiescent luminosity are both consistent with the observed tainty is related to the spectral type. Reig et al. (2005) use a properties of IGR J01363+6610. spectral type of B1V and an absolute magnitude of MV = –3.2 There may, however, be a less exotic explanation if the for the 2 kpc estimate, but they indicate that the star could be binary orbit of IGR J01363+6610 causes a neutron star to a more luminous sub-giant (i.e., B1IV). The difference in the sample an unusually low density part of the circumstellar absolute magnitude of a main sequence B-type star compared disk. This could occur if the system has a long orbital period to a B-type sub-giant is 0.6–0.7 magnitudes (Cox 2000), so with relatively low eccentricity. As pointed out by Reig et al. it is unlikely that the star is more luminous than MV ∼ –4.0, (2005), the EW of the Hα line for IGR J01363+6610 is which corresponds to a distance upper limit near 3 kpc. This, one of the strongest of any Be X-ray binary, suggesting the in turn, indicates an upper limit on the quiescent X-ray lu- presence of a large circumstellar disk, which might favor a minosity of <3.2 × 1031 ergs s−1, which is still significantly large and nearly circular orbit over a highly eccentric or- lower than the X-ray luminosities given above for compari- bit. While the system may have a large circumstellar disk, son. Negueruela & Okazaki (2001) have shown that tidal intera- Although it is not clear why the quiescent luminosity is tions of the neutron star can lead to truncation of the circum- so much lower for IGR J01363+6610, there are several po- stellar disk in low eccentricity systems. In addition to provid- tentially interesting possibilities. Among the more trivial is ing a possible explanation for the low quiescent luminosity of the possibility that the short (5 ks) Chandra observation oc- IGR J01363+6610, the fact that the disk can be truncated to a curred during an X-ray eclipse. Although eclipses are not size smaller than the Be star’s Roche lobe means that the sys- common for Be X-ray binaries due to their typically wide or- tem will not show luminous outbursts (Okazaki & Negueruela bits, we currently cannot rule out this possibility since the 2001), which is consistent with what has been seen so far for orbit and binary inclination for IGR J01363+6610 are un- IGR J01363+6610. constrained. It is also possible that Chandra did not detect Like IGR J01363+6610, the orbital parameters for the the source either because the compact object was highly ob- other Be X-ray binaries with the strongest Hα emission scured during the Chandra observation or that Chandra hap- lines, A 1118–616 (Coe et al. 1994) and IGR J01583+6713 pened to catch the source during a low point in its normal (Kaur et al. 2008), are not known. It has been suggested that short-term variability. Finally, there is also the possibility A 1118–616 has a nearly circular orbit based on its having a that the circumstellar disk, which was strongly present in small number of outbursts with no clear periodicity(Coe et al. 2004 and 2009 based on the detection of the Hα emission 1994). It has been argued that both of these systems have rel- line with an EW of ∼–50Å, had dissipated during the 2007 atively long orbital periods (∼0.6–2.2 years) based on mea- Chandra observation. Although dissipation and re-formation sured neutron star spin periods of 405.6 s and 469.2 s and of a circumstellar disk on this time scale did occur in the the Porb–Pspin relationship for Be X-ray binaries (Corbet 1986; case of GRO J1008–57 (Coe et al. 2007), and significant disk Coe et al. 1994; Kaur et al. 2008; Doroshenko et al. 2010). loss was observed for IGR J06074+2205 on a time scale of However, for IGR J01583+6713, the detection of X-ray pul- ∼3 years (Reig, Zezas & Gkouvelis 2010), in both of these sations was noted as being marginal in Kaur et al. (2008). Al- though there is no direct indication of the orbital or spin pe- TheBeX-raybinaryIGRJ01363+6610inquiescence 9 riod for IGR J01363+6610,it should be noted that the limit of possibility is that this could be a Be-black hole system. The Pspin < 17 s mentioned above could be taken as evidence for a former suggestion could be confirmed by a measurement of smaller Porb; however, we view the evidence for a limit on the the orbital period (although this will be challenging since the spin period as being relatively weak. source does not seem to produce regular outbursts), and the latter suggestion could be refuted with the detection of pulsa- 5. SUMMARY AND CONCLUSIONS tions during another outburst from the source. In summary, the detection of the variable hard X-ray source XMMU J013549.5+661243 at the location of a Be star provides confirmation that IGR J01363+6610 is a Be X-ray bi- JAT acknowledges partial support from NASA XMM- nary. Although some of its properties suggest that it is a fairly Newton Guest Observer award number NNX09AP91G. COH typical HMXB, its 2007 non-detection by Chandra indicate a acknowledges support from an NSERC Discovery Grant. We quiescent X-ray luminosity that is significantly lower than has thank Ms. Jia Liu for assistance with the optical observations been measured previously for other Be X-ray binaries. While and spectral reductions. We acknowledge helpful comments some possible explanations for the low luminosity (an eclipse, from the referee, Ignacio Negueruela. This publication makes a large change in the column density or spectrum, an extended use of data products from the Two Micron All Sky Survey, drop in mass accretion rate, or dissipation of the Be star’s which is a joint project of the University of Massachusetts circumstellar disk) may be consistent with relatively normal and the Infrared Processing and Analysis Center/California Be X-ray binary properties, other possibilities would require Institute of Technology, funded by NASA and the National that IGR J01363+6610 has unusual properties. One possibil- Science Foundation. This research makes use of the US- ity is that the system has a large and relatively circular orbit, NOFS Image and Catalog Archive operated by the United which could explain the low duty cycle for outbursts and, per- States Naval Observatory, Flagstaff Station and the SIMBAD haps, the low quiescent luminosity. Another very interesting database, operated at CDS, Strasbourg, France.

REFERENCES Balucinska-Church, M., & McCammon, D., 1992, ApJ, 400, 699 Kalberla, P. M. W., Burton, W. B., Hartmann, D., Arnal, E. M., Bajaja, E., Belczynski, K., & Ziolkowski, J., 2009, ApJ, 707, 870 Morras, R., & Pöppel, W. G. L., 2005, A&A, 440, 775 Bird, A. J., et al., 2010, ApJS, 186, 1 Kaur, R., Paul, B., Kumar, B., & Sagar, R., 2008, MNRAS, 386, 2253 Bird, A. J., et al., 2007, ApJS, 170, 175 Krivonos, R., Revnivtsev, M., Lutovinov, A., Sazonov, S., Churazov, E., & Bodaghee, A., et al., 2007, A&A, 467, 585 Sunyaev, R., 2007, A&A, 475, 775 Brown, E. F., Bildsten, L., & Rutledge, R. E., 1998, ApJ, 504, L95 Moon, D.-S., Kaplan, D. L., Reach, W. T., Harrison, F. A., Lee, J.-E., & Campana, S., Colpi, M., Mereghetti, S., Stella, L., & Tavani, M., 1998, Martin, P. G., 2007, ApJ, 671, L53 A&A Rev., 8, 279 Motch, C., Lopes de Oliveira, R., Negueruela, I., Haberl, F., & Janot-Pacheco, Campana, S., Stella, L., Israel, G. L., Moretti, A., Parmar, A. N., & Orlandini, E., 2007, in Active OB-Stars: Laboratories for Stellare and Circumstellar M., 2002, ApJ, 580, 389 Physics, ed. A. T. Okazaki, S. P. Owocki, & S. Steﬂ, Vol. 361, 117 Chakrabarty, D., & Roche, P., 1997, ApJ, 489, 254 Motch, C., Stella, L., Janot-Pacheco, E., & Mouchet, M., 1991, ApJ, 369, 490 Charles, P. A. & Coe, M. J., 2006, Optical, ultraviolet and infrared Negueruela, I., & Okazaki, A. T., 2001, A&A, 369, 108 observations of X-ray binaries, in Compact stellar X-ray sources. Edited Negueruela, I., Reig, P., Finger, M. H., & Roche, P., 2000, A&A, 356, 1003 by Walter Lewin & Michiel van der Klis: Cambridge University Press, 215 Okazaki, A. T., & Negueruela, I., 2001, A&A, 377, 161 Chaty, S., Rahoui, F., Foellmi, C., Tomsick, J. A., Rodriguez, J., & Walter, Oke, J. B., & Gunn, J. E., 1983, ApJ, 266, 713 R., 2008, A&A, 484, 783 Orosz, J. A., et al., 2002, ApJ, 568, 845 Coe, M. J., et al., 2007, MNRAS, 378, 1427 Patel, S. K., et al., 2007, ApJ, 657, 994 Coe, M. J., et al., 1994, A&A, 289, 784 Pavlov, G. G., Shibanov, Y. A., Zavlin, V. E., & Meyer, R. D., 1995, in The Cohen, D. H., 2000, in IAU Colloq. 175: The Be Phenomenon in Early-Type Lives of the Neutron Stars, ed. M. A. Alpar, U. Kiziloglu, & J. van Paradijs, Stars, ed. M. A. Smith, H. F. Henrichs, & J. Fabregat, Vol. 214, 156 71 Cohen, D. H., Cassinelli, J. P., & Macfarlane, J. J., 1997, ApJ, 487, 867 Porter, J. M., & Rivinius, T., 2003, PASP, 115, 1153 Corbet, R. H. D., 1986, MNRAS, 220, 1047 Rappaport, S., & van den Heuvel, E. P. J., 1982, in Be Stars, ed. M. Jaschek Corbet, R. H. D., 1996, ApJ, 457, L31 & H.-G. Groth, Vol. 98, 327 Cox, A. N., 2000, Allen’s astrophysical quantities, 4th ed. Publisher: New Reig, P., Negueruela, I., Manousakis, A., Marco, A., & Papamastorakis, G., York: AIP Press; Springer, Edited by Arthur N. Cox. 2004, The Astronomer’s Telegram, 343 Degenaar, N., et al., 2009, MNRAS, 396, L26 Reig, P., Negueruela, I., Papamastorakis, G., Manousakis, A., & Kougentakis, den Hartog, P. R., Hermsen, W., Kuiper, L., Vink, J., in’t Zand, J. J. M., & T., 2005, A&A, 440, 637 Collmar, W., 2006, A&A, 451, 587 Reig, P., Zezas, A., & Gkouvelis, L., 2010, arXiv:1006.4935 Doroshenko, V., Suchy, S., Santangelo, A., Staubert, R., Kreykenbohm, I., Remillard, R. A., & McClintock, J. E., 2006, ARA&A, 44, 49 Rothschild, R., Pottschmidt, K., & Wilms, J., 2010, A&A, 515, L1 Rodriguez, J., Tomsick, J. A., & Chaty, S., 2009, A&A, 494, 417 Fender, R. P., Russell, D. M., Knigge, C., Soria, R., Hynes, R. I., & Goad, Ryter, C. E., 1996, Ap&SS, 236, 285 M., 2009, MNRAS, 393, 1608 Smith, D. M., Heindl, W. A., Markwardt, C. B., Swank, J. H., Negueruela, I., Filliatre, P., & Chaty, S., 2004, ApJ, 616, 469 Harrison, T. E., & Huss, L., 2006, ApJ, 638, 974 Garmire, G. P., Bautz, M. W., Ford, P. G., Nousek, J. A., & Ricker, Stella, L., White, N. E., & Rosner, R., 1986, ApJ, 308, 669 G. R., 2003, in X-Ray and Gamma-Ray Telescopes and Instruments Strüder, L., et al., 2001, A&A, 365, L18 for Astronomy. Edited by Joachim E. Truemper, Harvey D. Tananbaum. Tomsick, J. A., Chaty, S., Rodriguez, J., Walter, R., & Kaaret, P., 2008, ApJ, Proceedings of the SPIE, 4851, 28 685, 1143 Gehrels, N., 1986, ApJ, 303, 336 Turner, M. J. L., et al., 2001, A&A, 365, L27 González, G., & González, G., 1956, Boletin de los Observatorios Walter, R., et al., 2006, A&A, 453, 133 Tonantzintla y Tacubaya, 2, 16 Wilms, J., Allen, A., & McCray, R., 2000, ApJ, 542, 914 Grebenev, S. A., Ubertini, P., Chenevez, J., Orr, A., & Sunyaev, R. A., 2004, Wilson, C. A., Weisskopf, M. C., Finger, M. H., Coe, M. J., Greiner, J., Reig, The Astronomer’s Telegram, 275 P., & Papamastorakis, G., 2005, ApJ, 622, 1024 Halpern, J. P., & Gotthelf, E. V., 2007, ApJ, 669, 579 Winkler, C., et al., 2003, A&A, 411, L1 Herbig, G. H., 1975, ApJ, 196, 129 Illarionov, A. F., & Sunyaev, R. A., 1975, A&A, 39, 185 in’t Zand, J. J. M., 2005, A&A, 441, L1