Multi-wavelength properties of IGR J05007-7047 (LXP 38.55\xspace) and identification as a Be X-ray binary pulsar in the LMC
We report on the results of a 40 d multi-wavelength monitoring of the Be X-ray binary system IGR J05007-7047 (LXP 38.55\xspace). During that period the system was monitored in the X-rays using the Swift\xspacetelescope and in the optical with multiple instruments. When the X-ray luminosity exceeded erg s we triggered an XMM-Newton\xspaceToO observation. Timing analysis of the photon events collected during the XMM-Newton\xspaceobservation reveals coherent X-ray pulsations with a period of 38.551(3) s (1), making it the 17 known high-mass X-ray binary pulsar in the LMC. During the outburst, the X-ray spectrum is fitted best with a model composed of an absorbed power law () plus a high-temperature black-body (kT keV) component. By analysing yr of available OGLE optical data we derived a 30.776(5) d optical period, confirming the previously reported X-ray period of the system as its orbital period. During our X-ray monitoring the system showed limited optical variability while its IR flux varied in phase with the X-ray luminosity, which implies the presence of a disk-like component adding cooler light to the spectral energy distribution of the system.
keywords:X-ray: binaries – individual: – (galaxies:) Magellanic Clouds
Be/X-ray binaries (BeXRBs) are a major subclass of high-mass X-ray binaries (HMXBs) that accord the majority of accreting X-ray pulsar systems. These systems are composed of a highly magnetised neutron star (NS), which accretes matter from a non super-giant Be type donor star. Since a young stellar population is needed in order to form such systems, their number within a galaxy is correlated with the recent star formation. BeXRBs are highly variable in the X-rays exhibiting moderate ( erg s) or major outbursts ( erg s). The X-ray variability is attributed to their irregular mass transfer mechanism, as accretion is fuelled by material ejected from the Be star towards its equatorial plane, forming a decretion disk around it. The spectral energy distribution (SED) of OBe stars is composed of two distinct components, a hot component originating from the stellar atmosphere and an infra-red (IR) excess produced by the equatorial decretion disk (Dougherty et al. 1994). For a comprehensive review, of the major characteristics of BeXRBs, we refer the reader to Reig (2011).
The Magellanic Clouds (MCs) offer a unique possibility to study the population of high-energy sources in a whole galaxy. Their moderate and well measured distances of 50 kpc for the Large Magellanic Cloud (LMC) (Pietrzyński et al. 2013) and 62 kpc for the Small Magellanic Cloud (SMC) (Graczyk et al. 2014) as well as their low Galactic foreground absorption (cm) makes them ideal targets for studying X-ray binary systems.
Within the last decade, numerous discoveries of new X-ray binaries have been made by monitoring the MCs (and the whole sky) in hard X-rays. The IBIS/ISGRI telescope (Ubertini et al. 2003) on board the INTEGRAL satellite (Winkler et al. 2003) and the Burst Alert Telescope (BAT, Barthelmy et al. 2005) on board the Swift\xspaceobservatory (Gehrels et al. 2004) have had a significant contribution. On the other hand, X-ray surveys of the Magellanic Clouds in the 0.3-10 keV energy band, like the XMM-Newton\xspacesurvey of the LMC (PI: F. Haberl) and the SMC (PI: F. Haberl, Haberl et al. 2012; Sturm et al. 2013) or Swift\xspace/XRT observations during the UV survey of both MCs (PI: S. Immler, Hagen et al. 2014), benefit from their better angular resolution for identifying BeXRB candidate systems. In particular the high sensitivity and sufficient energy resolution of XMM-Newton\xspaceallows the spectral characterization of the source and the identification of the optical counterpart.
The discovery of IGR J05007-7047\xspaceas an X-ray transient was made with INTEGRAL\xspacein the 17-60 keV band while its identification as a HMXB was based on a Chandra\xspaceobservation (Sazonov et al. 2005). D’Aì et al. (2011) reported a 30.77 d X-ray periodic modulation based on Swift\xspace/BAT survey data. The optical counterpart of the X-ray source was identified as USNO-B1.0 0192-0057570 having a spectral type of B2 III (Masetti et al. 2006).
In the current study we present the results of the X-ray spectral and temporal analysis of an XMM-Newton\xspaceToO observation that was triggered during the period of maximum luminosity of a moderate X-ray outburst in November 2014. The XMM-Newton\xspaceobservation revealed coherent 38.55 s pulsations. Following the nomenclature introduced by Coe et al. (2005) for the BeXRB pulsars in the SMC we suggest the alternative name LXP 38.55. Moreover, we report the results of our optical and X-ray monitoring during an orbital period of the system at the time of the X-ray outburst. Finally, we present a complete analysis of the 12 yr optical light-curve (LC) of the system as observed by the Optical Gravitational Lensing Experiment (OGLE, Udalski et al. 2015).
2 Analysis and results of X-ray observations
|#||OBS-ID||Instrument||T_start||orbital phase||exposure||count rate||UVOT filter||AB magnitudes|
|MJD [d]||[s]||cts s|
a Observation setup: For the XMM-Newton\xspaceToO, full-frame mode was used for all EPIC instruments, while medium filter was used for the EPIC/MOS cameras and thin filter for the EPIC/PN. For Swift\xspace/XRT the photon-counting mode (PC) was used.
b The zero orbital phase is selected to match Fig. 1 as d.
c Exposure time for XMM-Newton\xspaceEPIC cameras, Chandra\xspace/ACIS, or total Swift\xspace/XRT exposure time of all snapshots taken within the duration of the obs-id.
d Vignetting corrected count rate in the 0.3-10.0 keV energy band.
e UVOT filter used during the Swift\xspace/UVOT integration.
f Swift\xspace/UVOT AB magnitudes, not corrected for extinction.
g count rate of the slew survey is given tin the 0.2-12 keV band.
2.1 X-ray observations
We used Swift\xspace/XRT to monitor the most promising BeXRB candidate systems of the MCs (for most recent catalogues see: Haberl & Sturm 2016; Antoniou & Zezas 2016) to search for an X-ray outburst. We detected IGR J05007-7047\xspaceon 2014 October 26, at an X-ray Luminosity of 2 erg s in the 0.3-10.0 keV band. We continued monitoring the system with Swift\xspace/XRT for a period of 40 days to achieve an X-ray coverage of about one orbital period (30.771 d, see §3.3). During the maximum X-ray luminosity we triggered one of our granted XMM-Newton\xspaceanticipated ToO (PI: R. Sturm) that was performed on 2014 November 20. Looking at archival data, the system was detected during an XMM-Newton\xspaceslew (Saxton et al. 2008) on 2006 November 29 (XMMSL1 J050045.3-704441) at an angular distance of 4.7″(error 7.7″) from the XMM-Newton\xspaceToO position (see §2.2). Moreover, the system was also observed with Chandra\xspaceon 2005 June 16. The results of the Chandra\xspaceobservation were published in Sazonov et al. (2005), but the derived luminosities were based on a power-law model with photon index characteristic of an AGN (). Thus, we re-analysed the Chandra\xspaceobservation in order to calculate the luminosity of the system at the given time. The complete log of the Swift\xspace, XMM-Newton\xspaceand Chandra\xspaceobservations obtained and analysed for the current work is summarised in Table 1. The X-ray, optical and near IR light curves of the system during our monitoring are presented in Fig. 1.
2.2 X-ray position
The position of LXP 38.55\xspacewas determined from the XMM-Newton\xspaceToO observation that provided the best statistics. X-ray images were created from all the EPIC cameras using the XMM-Newton\xspacestandard energy sub-bands (Watson et al. 2009). Source detection was performed simultaneously on all the images using the SAS task edetect_chain. The source position was determined to R.A. = 05004606 and Dec. = –70°44′374 (J2000), with a statistical uncertainty of 0.04″. The total 1 positional error, however, is dominated by the remaining systematic uncertainty assumed to be 0.5″ (see section 4.3 of Sturm et al. 2013).
2.3 Timing analysis
We used the SAS task barycen to correct the XMM-Newton\xspaceEPIC event arrival times to the solar-system barycentre.
To search for a periodic signal we used an epoch folding technique (Davies 1990; Larsson 1996).
We initially searched for a periodic signal in the EPIC-pn event data. The detected signal was longer than the frame time of the EPIC-MOS instrument (2.6 s), thus
in order to increase the signal-to-noise ratio we used the merged event list of EPIC-pn and EPIC-MOS with the common good-time intervals for our timing analysis
By using the five standard energy bands (0.2-0.5, 0.5-1.0, 1.0-2.0, 2.0-4.5, 4.5-10 keV), we can define four hardness ratios as , with R denoting the background-subtracted count rate in energy band i. The period-folded pulse profiles in the EPIC standard energy bands, together with the hardness ratios derived from the pulse profiles in two adjacent energy bands, are plotted in Fig. 3.
2.4 Spectral analysis
We used xspec (version 12.8.2, Arnaud 1996) to perform the X-ray spectral analysis. The XMM-Newton\xspace/EPIC cameras provided a total of 35k counts enabling us to use statistics in the fitting procedure. The three spectra were fitted simultaneously with the same model with the addition of a scaling factor to account for instrumental differences.
For the EPIC-pn we fixed the scaling factor at 1, while for both EPIC-MOS, we obtained values of and , which is consistent with the expected values (see Stuhlinger
et al. 2006, or the latest version of the XMM-Newton\xspacecalibration manual
The spectra are well fitted by an absorbed power law (photon index 0.58, ), however the addition of an additional component can improve the fit at low energies and smooth out the residuals at higher energies (Fig. 4). It is not uncommon for HMXBs to show an additional component in their X-ray spectrum. Hickox et al. (2004) provided several interpretations for the origin of such emission that they refer to as soft excess. In the case of intermediate X-ray luminosity systems like LXP 38.55\xspace, this excess can be a result of reprocessing of hard X-rays from the NS by optically thick accreting material located at the inner edge of the accretion disk, by photo-ionized or collisionally heated diffuse gas, or thermal emission from the NS surface. Moreover, La Palombara et al. (2013a) showed that many X-ray pulsars show a spectral feature that can be interpreted as black body emission from a small region, most likely the polar cap of the NS. Multiple components may contribute to the spectrum of LXP 38.55\xspace. Although the simple power-law fit is already acceptable and the statistical quality of the EPIC spectra does not require additional model components, we test if emission components which are expected to exist, improve the spectral fit. We restrict the model to a maximum of two continuum components, i.e. a combination of the power law with a thermal component. We found two solutions which improve the fit at a) low energies (mainly below 1.5 keV) and b) higher energies (above 2 keV).
For case a) both a single-temperature black-body (e.g. BBodyRad in xspec) or a multi-temperature accretion disk (e.g. EzdiskBB in xspec) with temperature ( keV) provide the same improvement in fit quality.
It is however important to note that each of the models implies a different physical origin for the radiation process.
As the normalization parameter of both components translates to an emission region with radius
The best fit is provided by the combination of an absorbed power-low plus a high-temperature black-body component with temperature 2.04 keV (case b). This model improves significantly the residuals of the fit at higher energies and provides a statistical improvement above a 3.5 level compared the low temperature black-body plus power law model. There is clear evidence of an emission line present at 6.4 keV that originates from neutral Fe (K line). The width of the line ( keV) was found to be bounded by the spectral resolution of the instrument, thus we fixed the width of the line to zero since an exact value could not be resolved. The parameters for the models described above are given in Table 2.
We note that both the low and high temperature black-body components improve the residuals on different parts of the X-ray spectrum and account for different physical processes, thus in principle they could co-exist in the system. However, modelling the spectrum with a combination of three continuum components is not statistically justified. Moreover, the use of three continuum components together with the column density as free parameter adds degeneracy to the solution. While we cannot exclude the existence of more components in the spectrum, the available statistics does not justify the use of a more complicated spectral model.
The Swift\xspace/XRT observations could not provide sufficient statistics for modelling the X-ray spectrum, we were however able to derive values for the hardness of the spectrum in two adjacent energy bands (0.3-2.0 keV vs. 2.0-10 keV) and compare this quantity along the evolution of the burst. There was no significant evidence for a monotonous change in the hardness of the spectrum with luminosity along the 40 day monitoring cycle (see Fig.5), but a rather weak trend of the system becoming harder with increasing luminosity.
|Power-law plus High Energy Black body (T1 keV)|
|tbnew*tbnew*(gaussian + bbodyrad + powerlaw)|
|Flux/Norm||(7.7)||erg cm s|
|Observed Flux||1.03||erg cm s|
|Unabsorbed Luminosity||erg s|
|Power-law plus Low Energy Black body (T1 keV)|
|tbnew*tbnew*(gaussian + bbodyrad + powerlaw)|
|Flux/Norm||(10.920.15)||erg cm s|
|Observed Flux||1.05||erg cm s|
|Unabsorbed Luminosity||erg s|
|Power-law plus low Energy multi-temperature disk Black body (T1 keV)|
|tbnew*tbnew*(gaussian + EzdiskBB + powerlaw)|
|Flux/Norm||(10.920.16)||erg cm s|
|Observed Flux||1.05||erg cm s|
|Unabsorbed Luminosity||erg s|
a Unabsorbed flux of the power-law component in the 0.3-10. keV energy band.
b Absorbed flux of the fitted model in the 0.3-10. keV energy band.
c Unabsorbed Luminosity of the fitted model in the 0.3-10. keV energy band, for a distance of 50 kpc (Pietrzyński et al. 2013).
3 Analysis and results of optical data
3.1 Swift Uvot
Our Swift\xspacemonitoring enabled us to measure the UV magnitudes (using the Swift\xspaceUV/Optical Telescope, UVOT). We used the default filter of the day for the Swift\xspace/UVOT instrument setup. This resulted in nine observations performed with the U filter (central wavelength 3465 Å), seven with the UVM2 (2246 Å) and one with the UVW2 (1928 Å). To derive the systems UV magnitudes we used the uvotsource tool. A 5″ radius was used for performing aperture photometry for all filters. No significant variation was measured in any of the filters; the mean values for the derived magnitudes were: 14.230.02 for the U, 14.120.05 for the UVM2 and 14.1660.06 for the UVW2 filter.
LXP 38.55\xspacewas observed thirteen times between 2014 November 4 and December 2 with the Gamma-Ray Burst Optical/Near-Infrared Detector (GROND; Greiner et al. 2008) at the MPI/ESO 2.2 m telescope (La Silla, Chile). Seven bands were used, providing a coverage in both optical (g’, r’, i’ and z’) and near-IR (J, H and K) wavelengths. Each individual observation consists of 24 dithered exposures of 10 s in the near-IR and 4 dithered exposures in the optical with various exposure times (35 s in the eight first epochs and 66 s in the five last epochs), taken at a mean airmass of 1.3 and a mean seeing of 1.4”.
Single dithered exposures were
reduced (bias subtraction, flat-fielding, distortion correction) and stacked using standard
Having simultaneously obtained the photometry of the optical companion from NIR to optical, we can construct its SED and look for an IR excess. As BeXRBs are intrinsically variable in the optical, having both orbital and super-orbital variability, GROND data offer us the unique opportunity to study the SED of such systems. Moreover, since our Swift\xspace/UVOT fluxes showed no significant variability (compared to the NIR fluctuations) we can supplement the GROND magnitudes with the UV data from Swift\xspace, constructing an SED with (quasi-)simultaneous observations. To de-redden the optical observations, we used the extinction curve given by Misselt et al. (1999), which is most appropriate for the LMC. The foreground Galactic reddening was set to mag (Schlafly & Finkbeiner 2011). We used the Kurucz models of stellar atmospheres (Kurucz 1979; Castelli et al. 1997) to fit the observed photometry, while constraining the metallicity to that of the LMC ([Fe/H]=-0.5).
a The ratio was computed based on the best fit model and the observed flux during the first night. The uncertainties were computed accounting for both the observational uncertainties and the systems variability.
The spectral types of the optical counterparts of BeXRBs in the SMC are reported to be in the range between O5 and down to B9 (Coe et al. 2005). However, in the new compilation of spectral types by Haberl & Sturm (2016) only stars with uncertain classification may have types as late as B5. There are no B6, B7 or B8 counterparts and only one B9 which may cast doubts on the spectral classification or the correct identification. This makes the spectral type distribution more similar to that of the Galaxy which seems to end at B3. A B2 III star, like our system in study, is expected to have a temperature of 21500 K (Cox 2000), thus we limit the temperature parameter to the range of 20000-24000 K. We note that the decretion disk of a Be star can account up to 50% of its total luminosity (de Wit et al. 2006), thus stellar atmosphere models are expected to generally produce poor fits to the variable SED of those systems. A brute force fit of a stellar atmosphere model to our data without accounting for any prior knowledge of the systems properties (spectral type or distance) produces questionable results. This method delivers a best fit model of a much cooler than expected star (T10000 K) and at a much closer distance than that of the LMC. By limiting the temperature parameter search grid to the range of 20000-24000 K, and by excluding the optical and NIR GROND data from the fit we can achieve a good fit, where the normalization and the radius of the star describe our prior assumptions. The best derived model was that of a star with temperature of 22000 K and surface gravity . The fitted SED of the system is presented in Fig. 7. From the fitted model it is clear that the contribution of the disk is not restricted to the NIR magnitudes but affects the optical emission of the system as well.
To stress the importance of UV measurements to the SED fitting we repeated the above exercise by using only the GROND magnitudes. We then get a best fit model of a star with a temperature range of 8500-9500 K, with a normalization pointing to a source located much closer than the LMC distance. This suggests the presence of a cool disk that contributes to the SED of the B star. Further investigation of the contribution of the disk in the SED of the star lies beyond the scope of the present work as it requires precise modelling of the decretion disk.
The Optical Gravitational Lensing Experiment (OGLE) started its initial observations in 1992 (Udalski et al. 1992) and continues observing till today (OGLE IV, Udalski et al. 2015). Observations were made with the 1.3 m Warsaw telescope at Las Campanas Observatory, Chile. Images are taken in the V and I filter pass-bands, while the data reduction is described in Udalski et al. (2015). OGLE photometric magnitudes are calibrated to the standard VI system.
In the present work, we used the OGLE I and V band magnitudes for the counterpart of LXP 38.55\xspacethat were taken between Modified Julian dates (MJD) 52166 and 57275. The latest I and V band light curves are shown in Fig. 8. The photometric data can be downloaded from the OGLE-IV real time monitoring of X-ray sources web-page
By using the 14 year of available OGLE observations in the I band we were able to derive a refined orbital period of 30.7760.005 d (see Fig. 9).
Additionally, a super-orbital quasi-periodic modulation can be seen in the data.
We note that the super-orbital modulation should not be treated as a periodic signal like the orbital modulation, its nature is probably related to the precession of the Be disk.
Period uncertainties were computed using the bootstrap method and repeating the calculations for 10000 light curves.
Both orbital and super-orbital modulations are confirmed well above a 99% confidence level by performing 10000 white and and red noise (using REDFIT, Schulz &
By monitoring the BeXRB system IGR J05007-7047 we were able to promptly trigger an XMM-Newton\xspaceToO. Timing analysis of the photon events collected during the XMM-Newton\xspaceToO reveals coherent X-ray pulsations (see Fig. 2) with a period of 38.551(3) s (1), making LXP 38.55\xspacethe 17 known HMXB pulsar in the LMC (see Antoniou & Zezas 2016, and references therein). During the outburst, the X-ray spectrum is well represented by an absorbed power law (), however we argue for the presence of additional components that can both improve the residuals of the fit and account for known physical processes resulting in X-ray emission from BeXRB systems. Our analysis indicates that a model composed of an absorbed power-law plus a high-temperature black-body (kT 2 keV) model (Fig. 4) improves the fit. A model composed by an absorbed power law () and a multi-temperature low-temperature disk black body (kT keV) provides an acceptable fit but only improves the residuals at lower energies. The tested models can help us to study different emission mechanisms and understand the physical properties of different regions in the binary system, like the accretion column towards the magnetic poles of the neutron star, or the inner edge of the accretion disk at the outer parts of the magnetosphere.
The normalization of the best fit model with a high-energy black-body component yields an accretion region with a radius of 0.7 km, and can be interpreted as emission from a polar cap. The normalization of the multi-temperature black-body component translates to an emitting region with a size of , which could be the inner region of an accretion disk around the NS. Another possible explanation is that this low temperature multi-temperature black-body component originates from reprocessing of X-rays in the inner part of the accretion disk. Hard X-rays originating from the NS can illuminate an area at distance under a solid angle . This radius is different than the one computed from the normalization of the black-body and can be calculated from (Hickox et al. 2004), where is the total X-ray luminosity and is the temperature of the black-body. From fitted parameters we calculate a radius of .
Expressing the X-ray luminosity in terms of mass accretion rate we find M using a typical accretion efficiency value of 0.2 (Sibgatullin & Sunyaev 2000). Assuming the magnetospheric radius of the NS has the same size as the accretion disk inner radius we can make an order of magnitude estimation of the magnetic field strength of the NS. Following Frank et al. (2002), for an X-ray luminosity of 4 erg s, a magnetospheric radius of 200 km and typical values for the mass and radius (=1.4 M and =12.5 km) we derive a magnetic field strength of G. Moreover we note that the magnetospheric radius is much smaller than the corotation radius of the NS, which is , where is the spin period of the NS. If we assume that the accretion originates from a Keplerian disk with inner radius where the velocity of the disk exceeds the spin angular velocity of the NS, we conclude that the interaction between the Keplerian disk and the NS results in a spin-up of the system (accretor regime).
The 30.77 d optical modulation can be attributed to the orbital period of the NS-Be star system. This modulation is associated to the truncation of the Be disk due to gravitational interaction with the NS during an orbital revolution. The colour modulation –redder when brighter– is also a result of changes in the decretion disk around the Be star. A face-on geometry of the disk is expected to cause the above behaviour (Rajoelimanana et al. 2011). Similar behaviour has been encountered in numerous BeXRB systems, like XMMU J010743.1-715953 (Coe et al. 2012), LXP 168.8 (Maggi et al. 2013), LXP 8.04 (Vasilopoulos et al. 2014) or isolated Be Stars in the SMC (de Wit et al. 2006). The disk adds cooler light to the stellar emission of the B Starwithout obscuring the star itself (face-on geometry), while variations in the disk size or radius within the binary orbit result in a variable SED. This is also noticeable in the GROND observations where the optical part of the SED (dominated by the B star) is much less variable than the IR (larger contribution from the Be disk).
The long term optical variability (super-orbital) is believed to be related to the formation and depletion of the circumstellar disk around the Be star (Rajoelimanana et al. 2011). The orbital modulation was still present in all the available data, and the same colour modulation –redder when brighter– was detected (see Fig. 10). During the OGLE-III phase the system exhibited high long-term variability evident of a variable Be disk. There is no evidence for significant changes in the size of the disk during the OGLE-IV epoch, where the amplitude of orbital modulation in the optical band is comparable to the long term variability of the system. Multi-epoch observations in the X-ray and optical would allow for a comparison between the long term X-ray and optical variability, unfortunately all the available X-ray observations where performed when the OGLE I band magnitude of the system was between 14.55 and 14.7.
By comparing the available X-ray observations that were performed on times separated by more than one orbital period, but correspond to the same orbital phase we can advocate about the long term X-ray variability of the system. The XMM-Newton\xspaceslew detection was performed at a phase when only Swift\xspaceobservations were available (phase 0.87 of Fig. 1 & 9). By comparing their fluxes we derive that at the time of the XMM-Newton\xspaceslew detection the luminosity of the system was higher by a factor of compared to the one measured with Swift\xspace. The Chandra\xspaceobservation of the system was performed on MJD53537.1 d corresponding to phase 0.61. The luminosity of the system was found to be 2.6 erg s in the 0.3-10.0 keV band. By comparing the reported luminosity to the Swift\xspaceluminosity corresponding to the same orbital phase we found that the system was brighter by a factor of during the Swift\xspace/XRT monitoring performed on November 2014. Our data confirm the findings of D’Aì et al. (2011) that the systems long term X-ray emission is persistent within a variability factor of 3-4, which is lower than the orbital variability (variability factor ).
Our Swift\xspace/XRT monitoring was performed when the average OGLE I band magnitude of the system was 14.65. Based on the SED modelling we found that in the I band the system has an excess of 24% compared to the template flux of a B star with same temperature (see Tab. 3). Interestingly this translates to a difference of 0.23 in magnitude, which is exactly the difference between the average I magnitude during our monitoring and the lowest flux level observed in the OGLE I band; I mag of 14.88 at MJD 52930 d. This might be evident of an epoch where the Be disk was almost completely depleted.
D’Aì et al. (2011) reported the system as a possible wind fed X-ray system, but our results reveal the presence of a cooler component in the SED of the companion star suggesting that this is a NS pulsar with a Be companion. Moreover, the hot thermal excess in the X-ray spectrum adds to the population of BeXRB systems sharing this feature. Recently La Palombara et al. (2013b) reported a population of long spin (200 s), low luminosity (L erg s), X-ray persistent Galactic BeXRBs systems showing this characteristic feature. A similar hot thermal excess was also found in the X-ray spectra of other BeXRBs in the Magellanic Clouds during moderate outbursts with L erg s (Vasilopoulos et al. 2013; Bartlett et al. 2013). LXP 38.55\xspaceis likely another example of such a system, showing evidence of a hot black-body component at higher X-ray luminosity, suggesting the black-body component is more common and not limited to the population originally proposed by La Palombara et al. (2013b). It is not clear if this feature is model dependent, but the physical parameters derived from such a phenomenological spectral treatment where the degeneracy of the fitted parameters is significant should be treated as order of magnitude approximations. Yet the study of such systems in the Magellanic Clouds has the major advantage of the low column density in the line of sight, thus through the study of these systems we are able to improve the modelling of their X-ray spectrum and derive stringent constraints on our findings.
The position of LXP 38.55\xspacein the Corbet diagram (Knigge et al. 2011) lies close to the best fit line derived from the population of confirmed MC pulsars (see Fig. 11). Looking at the properties of other BeXRB systems located near LXP 38.55\xspacein the vs. diagram, we found no notable similarities between them and the newly reported system. Most of these systems show much larger long term X-ray variability, with variability factors 1000, as expected for such short-spin systems (Haberl & Sturm 2016). However since LXP 38.55\xspaceis a newly discovered system, future observations might reveal a higher variability factor. Moreover the lack of deep X-ray observations during their peak luminosity prohibits us from comparing their X-ray spectral properties. SXP 25.55 (Haberl et al. 2008a) with an orbital period of 22.5 d (Rajoelimanana et al. 2011) is an SMC BeXRB pulsar located very close to the position of LXP 38.55\xspacein the Corbet diagram. An XMM-Newton\xspacespectrum is available for this system, but its peak luminosity is about 30 times smaller than that of LXP 38.55\xspace, while its spectrum is affected by stronger absorption. LXP 28.28 (Şahiner et al. 2016) is the closest LMC system having an orbital period of d. (Coe et al. 2015).
A triggered XMM-Newton\xspaceToO observation of the candidate BeXRB system IGR J05007-7047\xspacelocated in the LMC revealed coherent X-ray pulsations with a period of 38.551(3) s. The X-ray spectrum of the system can be described by an absorbed power-law () plus a hot black body component ( keV) that can be interpreted as black body emission from a small region, most likely the polar cap of the NS. Additionally there is a statistically significant evidence for a soft excess that can originate from a region close to the magnetospheric radius and is a result of reprocessing of hard X-rays emitted by the accretion column of the NS. Analysis of optical data confirms the 30.776(5) d orbital period that was previously reported. By using simultaneously obtained optical and NIR data we constructed the SED of the binary system that reveals the presence of a decretion disk around the massive optical companion of the NS. Our analysis confirms LXP 38.55\xspaceas the 17 HMXB pulsar in this star forming galaxy.
We would like to thank the anonymous referee for his/her critical review and suggestions that helped improve the manuscript. The XMM-Newton\xspace project is supported by the Bundesministerium für Wirtschaft und Technologie / Deutsches Zentrum für Luft- und Raumfahrt (BMWi/DLR, FKZ 50 OX 0001) and the Max-Planck Society. The OGLE project has received funding from the National Science Centre, Poland, grant MAESTRO 2014/14/A/ST9/00121 to AU. We thank the Swift\xspaceteam for accepting and carefully scheduling the target of opportunity observations, and we acknowledge the use of public data from the Swift\xspacedata archive. G. V. acknowledge support from the BMWi/DLR grant FKZ 50 OR 1208.
- For the analysis the photon arrival times are randomised over the frame time of each instrument.
- , where is the inner radius of the disk in km, is the distance to the source in units of 10 kpc, is the inclination, and is the color to effective temperature ratio
- IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.
- XROM: http://ogle.astrouw.edu.pl/ogle4/xrom/xrom.html
- REDFIT: http://www.ncdc.noaa.gov/paleo/softlib/redfit/redfit.html
- Antoniou V., Zezas A., 2016, preprint, (arXiv:1603.08011)
- Arnaud K. A., 1996, in Jacoby G. H., Barnes J., eds, Astronomical Society of the Pacific Conference Series Vol. 101, Astronomical Data Analysis Software and Systems V. p. 17
- Barthelmy S. D., et al., 2005, \ssr, 120, 143
- Bartlett E. S., Coe M. J., Ho W. C. G., 2013, \mnras, 436, 2054
- Castelli F., Gratton R. G., Kurucz R. L., 1997, \aap, 318, 841
- Coe M. J., Edge W. R. T., Galache J. L., McBride V. A., 2005, \mnras, 356, 502
- Coe M. J., et al., 2012, \mnras, 424, 282
- Coe M. J., Finger M., Bartlett E. S., Udalski A., 2015, \mnras, 447, 1630
- Cox A. N., 2000, Allen’s astrophysical quantities
- D’Aì A., La Parola V., Cusumano G., Segreto A., Romano P., Vercellone S., Robba N. R., 2011, \aap, 529, A30
- Davies S. R., 1990, \mnras, 244, 93
- Dickey J. M., Lockman F. J., 1990, \araa, 28, 215
- Dougherty S. M., Waters L. B. F. M., Burki G., Cote J., Cramer N., van Kerkwijk M. H., Taylor A. R., 1994, \aap, 290
- Frank J., King A., Raine D. J., 2002, Accretion Power in Astrophysics: Third Edition
- Gehrels N., et al., 2004, \apj, 611, 1005
- Graczyk D., et al., 2014, \apj, 780, 59
- Gregory P. C., Loredo T. J., 1996, \apj, 473, 1059
- Greiner J., et al., 2008, \pasp, 120, 405
- Haberl F., Sturm R., 2016, \aap, 586, A81
- Haberl F., Eger P., Pietsch W., Corbet R. H. D., Sasaki M., 2008a, \aap, 485, 177
- Haberl F., Eger P., Pietsch W., 2008b, \aap, 489, 327
- Haberl F., et al., 2012, \aap, 545, A128
- Hagen L., Siegel M., Gronwall C., Hoversten E. A., Immler S., 2014, in American Astronomical Society Meeting Abstracts #223. p. #442.30
- Hickox R. C., Narayan R., Kallman T. R., 2004, \apj, 614, 881
- Horne J. H., Baliunas S. L., 1986, \apj, 302, 757
- Kalberla P. M. W., Burton W. B., Hartmann D., Arnal E. M., Bajaja E., Morras R., Pöppel W. G. L., 2005, \aap, 440, 775
- Knigge C., Coe M. J., Podsiadlowski P., 2011, \nat, 479, 372
- Krühler T., et al., 2008, \apj, 685, 376
- Kurucz R. L., 1979, \apjs, 40, 1
- La Palombara N., Mereghetti S., Sidoli L., Tiengo A., Esposito P., 2013a, preprint, (arXiv:1301.5120)
- La Palombara N., Mereghetti S., Sidoli L., Tiengo A., Esposito P., 2013b, \memsai, 84, 626
- Larsson S., 1996, \aaps, 117, 197
- Leahy D. A., Darbro W., Elsner R. F., Weisskopf M. C., Kahn S., Sutherland P. G., Grindlay J. E., 1983, \apj, 266, 160
- Maggi P., Haberl F., Sturm R., Pietsch W., Rau A., Greiner J., Udalski A., Sasaki M., 2013, \aap, 554, A1
- Masetti N., et al., 2006, \aap, 459, 21
- Misselt K. A., Clayton G. C., Gordon K. D., 1999, \apj, 515, 128
- Monet D. G., et al., 2003, \aj, 125, 984
- Pietrzyński G., et al., 2013, \nat, 495, 76
- Rajoelimanana A. F., Charles P. A., Udalski A., 2011, \mnras, 413, 1600
- Reig P., 2011, \apss, 332, 1
- Rolleston W. R. J., Trundle C., Dufton P. L., 2002, \aap, 396, 53
- Saxton R. D., Read A. M., Esquej P., Freyberg M. J., Altieri B., Bermejo D., 2008, \aap, 480, 611
- Sazonov S., Churazov E., Revnivtsev M., Vikhlinin A., Sunyaev R., 2005, \aap, 444, L37
- Scargle J. D., 1982, \apj, 263, 835
- Schlafly E. F., Finkbeiner D. P., 2011, \apj, 737, 103
- Schulz M., Mudelsee M., 2002, Computers and Geosciences, 28, 421
- Sibgatullin N. R., Sunyaev R. A., 2000, Astronomy Letters, 26, 699
- Skrutskie M. F., et al., 2006, \aj, 131, 1163
- Stuhlinger M., et al., 2006, in Wilson A., ed., ESA Special Publication Vol. 604, The X-ray Universe 2005. p. 937 (arXiv:astro-ph/0511395)
- Sturm R., et al., 2013, \aap, 558, A3
- Ubertini P., et al., 2003, \aap, 411, L131
- Udalski A., 2008, \actaa, 58, 187
- Udalski A., Szymanski M., Kaluzny J., Kubiak M., Mateo M., 1992, \actaa, 42, 253
- Udalski A., Szymański M. K., Szymański G., 2015, \actaa, 65, 1
- Vasilopoulos G., Maggi P., Haberl F., Sturm R., Pietsch W., Bartlett E. S., Coe M. J., 2013, \aap, 558, A74
- Vasilopoulos G., Haberl F., Sturm R., Maggi P., Udalski A., 2014, \aap, 567, A129
- Verner D. A., Ferland G. J., Korista K. T., Yakovlev D. G., 1996, \apj, 465, 487
- Watson M. G., et al., 2009, \aap, 493, 339
- Wilms J., Allen A., McCray R., 2000, \apj, 542, 914
- Winkler C., et al., 2003, \aap, 411, L1
- Yoldaş A. K., Krühler T., Greiner J., Yoldaş A., Clemens C., Szokoly G., Primak N., Klose S., 2008, in Galassi M., Palmer D., Fenimore E., eds, American Institute of Physics Conference Series Vol. 1000, American Institute of Physics Conference Series. pp 227–231, doi:10.1063/1.2943450
- York D. G., et al., 2000, \aj, 120, 1579
- Şahiner Ş., Serim M. M., Baykal A., İnam S. Ç., 2016, \mnras, 456, 845
- de Wit W. J., Lamers H. J. G. L. M., Marquette J. B., Beaulieu J. P., 2006, \aap, 456, 1027