The discovery of weak coherent pulsations in the ultraluminous X-ray source NGC 1313 X-2

TL;DR

This paper reports the discovery of weak, sinusoidal X-ray pulsations from the ultraluminous X-ray source NGC 1313 X-2, providing evidence for a neutron star accretor and insights into its magnetic field and accretion processes.

Contribution

First detection of weak pulsations in NGC 1313 X-2, supporting neutron star models for ULXs and suggesting strong winds or jets are involved in accretion.

Findings

Pulsation period of ~1.5 seconds with 5% pulsed fraction.

Detections unlikely due to false noise, confirmed by simulations.

Implications for neutron star magnetic field and accretion history.

Abstract

We report the detection of weak pulsations from the archetypal ultraluminous X-ray source (ULX) NGC 1313 X-2. Acceleration searches reveal sinusoidal pulsations in segments of two out of six new deep observations of this object, with a period of $\sim$ 1.5 s and a pulsed fraction of $\sim$ 5%. We use Monte Carlo simulations to demonstrate that the individual detections are unlikely to originate in false Poisson noise detections given their very close frequencies; their strong similarity to other pulsations detected from ULXs also argues they are real. The presence of a large bubble nebula surrounding NGC 1313 X-2 implies an age of order 1 Myr for the accreting phase of the ULX, which implies that the neutron star's magnetic field has not been suppressed over time by accreted material, nor has the neutron star collapsed into a black hole, despite an average energy output into the nebula two orders of magnitude above Eddington. This argues that most of the accreted material has been expelled over the lifetime of the ULX, favouring physical models including strong winds and/or jets for neutron star ULXs.

Tables

Table 1. Table 1: Characteristics of the detected pulsations.

Data set	Date${}^{\text{a}}$	$f^{\text{b}}$	$\dot{f^{\text{c}}}$	PRESTO significance	Pulsed fraction	$L_{\text{X}}^{\text{d}}$
Obs.ID	(MJD)	(Hz)	(Hz s-1)	($\sigma$)	(%)	(1040 erg s-1)
0803990401	57998.0	0.686041(1)	6.5(5) $\times {10}^{-9}$	3.00	5.3 $\pm$ 2	1.99
0803990601	58096.0	0.6843642(6)	1.5(2) $\times {10}^{-8}$	4.12	6.4 $\pm$ 2	1.44

Full text

††pubyear: 2019††pagerange: The discovery of weak coherent pulsations in the ultraluminous X-ray source NGC 1313 X-2–The discovery of weak coherent pulsations in the ultraluminous X-ray source NGC 1313 X-2

1 Introduction

Ultraluminous X-ray sources (ULXs) are extragalactic off-nuclear X-ray sources with apparent X-ray luminosities that match or exceed the Eddington limit for a 10 $M_{\odot}$ black hole (see Kaaret et al. 2017 for a recent review). For many years the debate over ULXs focussed on the mass of the black hole powering this phenomenon (e.g. Colbert & Mushotzky 1999; King et al. 2001; Zampieri & Roberts 2009). However the detection of coherent pulsations in four ULXs implies that each must host a neutron star (NS) primary radiating at up to several hundred times its Eddington limit (Bachetti et al. 2014; Israel et al. 2017a, b; Fürst et al. 2016; Carpano et al. 2018). Hence, the key questions for ULXs now revolve around what fraction of the population host NSs, and the physics of super-Eddington accretion in strong magnetic fields. The discovery of further NS-ULXs is crucial to answering both questions, with signatures such as cyclotron resonant scattering features and possible transitions in and out of a propeller regime revealing more NS-ULX candidates (e.g. Brightman et al. 2018; Earnshaw et al. 2018). Indeed, broad-band spectroscopic similarities in a small sample of ULXs suggest NSs may be ubiquitous in the ULX population (Walton et al. 2018). This can only be confirmed by more evidence for the presence of NSs in individual objects. This evidence is growing with, for example, the first detection of a Galactic pulsating ULX (PULX, e.g. Wilson-Hodge et al. 2018). Here we report a new detection of pulsations from a hitherto archetypal extragalactic ULX, NGC 1313 X-2. This PULX is located on the outskirts of the barred spiral galaxy NGC 1313, and is historically amongst the best studied ULXs in terms of both its X-ray properties (e.g Feng & Kaaret 2007; Pintore & Zampieri 2012) and its optical counterpart and environment (e.g. Pakull & Mirioni 2002; Grisé et al. 2008a). In particular, a lack of radial velocity variations in its optical spectrum (Roberts et al. 2011) and the patterns of change in its hard X-ray spectrum and flux variability could suggest it is viewed at a low inclination angle (Sutton et al. 2013; Middleton et al. 2015), affording a direct view of the central regions of its accretion flow.

2 Observations and data analysis

The barred spiral galaxy NGC 1313 was observed by XMM-Newton seven times between 2017 June 16 and 2017 September 24, as part of a Large Programme observing both ULXs in the galaxy: NGC 1313 X-1 and X-2 (see Pinto et al. in prep. for a summary of the data obtained). The latter is the target of this study, and it lies $\sim$ 7 arcminutes off-axis from the nominal aim-point of the observations. The observations were carried out with the EPIC-pn and MOS detectors in full-frame imaging mode. The raw observation data files (ODF) were processed with the Science Analysis Software (sas)v.16. We excluded observation 0803990701 from further analysis at this stage as it contained only $\sim 15$ ks of exposure; all other observations were exposed for $>120$ ks. Source events were extracted in a circular aperture of 25 arcsecond radius centred on the target, while the background region was chosen to be adjacent to the source on the same chip. We then used the the sas task evselect to derive the background-subtracted 0.3-10 keV light curves in 73 ms time bins, which corresponds to the maximal time resolution of the EPIC-pn camera (in full-frame mode). We did not utilise data from the MOS detectors, since they do not have sufficient time resolution to be able to detect pulsations with frequencies larger than $0.2$ Hz. We implemented the standard filtering of events as recommended in the XMM-Newton science threads. Namely, we selected only single or double patterned events (PATTERN $<=4$) and used the standard flagging criteria (i.e. #XMMEA_EP) that ensures the exclusion of hot pixels and events falling outside the field of view. In some observations the source emission was spread over more than one chip, and in such circumstances, we combined the data from both chips to construct the light curves.

2.1 Discovery of coherent pulsations

As the starting point of our timing analysis, we converted the event times from the EPIC-pn detector to the barycentre of the Solar System with the sas task barycen, using the Chandra source position (RA = 49.594 degrees, DEC = -66.601 degrees) and the DE-200 Solar System ephemeris (Standish 1990). We computed the power spectra of the six remaining observations using the heasarc (v6.22) tool powspec. Initially, these did not indicate any signals peaking above the $1\sigma$ white-noise threshold (van der Klis 1989). However, if the source has a relatively strong period derivative $\dot{f}$, this causes the signal to drift across multiple Fourier bins in the power spectrum, with the number of bins scaling as $\dot{f}T_{\text{obs}}^{2}$ for observations of length $T_{\text{obs}}$, thus reducing the amplitude of the signal power and the sensitivity of any given search. In order to correct for such an effect, we implemented the Fourier domain acceleration search method by using the presto software (Ransom, Eikenberry & Middleditch 2002). This method is optimised if the acceleration (or the period derivative) of the source is constant throughout the observation, which is best satisfied if the exposure time is smaller than a tenth of the orbital period of the system. Therefore, we split the light curve of each observation, first into (i) contiguous 30 ks segments, which accounts for the strongest period derivative measured so far in PULXs (i.e. $\sim 10^{-7}$ Hz s*-1*; Carpano et al. 2018), and later into (ii) larger 60 ks segments, which allows for a stronger detection of candidates with smaller period derivatives. We show the full light curves, and their segmentation, in Figure 1; in all cases we ran the search separately on each segment.

This resulted in the possible detection of coherent pulsations in a 30 ks and a 60 ks segment of two separate observations, out of the total of six, indicating that NGC 1313 X-2 is a new addition to the class of PULXs. Both of these detections were found to occur at similar Fourier frequencies and with a statistical significance larger than 3$\sigma$ (see Figures 2 and 3)111Here, we refer to the significance reported by the accelsearch routine in presto, which accounts for the total number of trials. We folded the light curves of these two segments on the best-fit values of the pulse frequency and its derivative (determined using presto) in order to obtain the pulse profiles (Larsson 1996; Figure 2). These appear sinusoidal, as per other PULXs, but the pulsed fractions at $\sim 5\%$ are notably smaller than previously seen (e.g. Bachetti et al. 2014; Fürst et al. 2016). We further verified our timing analysis through an epoch folding search on the two light curve segments shown in Figure 1 by using the hendrics tool HENzsearch (Bachetti 2018). Here, a test for periodicity was performed by evaluating the Z2 statistic over a grid of $f-\dot{f}$ values (c.f. Buccheri et al. 1983 for further details). The results of the test are shown in Figure 3, with the peak Z2 values in both observations being larger than the $3\sigma$ significance level, consistent with the results obtained through presto.

We calculated the unabsorbed 0.3-10 keV X-ray luminosity for each observation by fitting a simple power-law plus multi-colour disc model to the EPIC-pn spectra in xspec and correcting for both the Galactic and intrinsic absorption columns. The X-ray spectra were extracted by following the standard procedure outlined in Walton et al. (2017). We tabulate these luminosities alongside the characteristics of the detected pulses in Table 1. Finally, we note that we do not detect pulsations in any archival observations with similar quality of data and exposure as in the six observations from our programme.

2.2 Verifying the detection significance with Monte Carlo simulations

In van der Klis (1989), the significance level for the detection of a coherent signal (above a given threshold) is evaluated by considering the distribution of the noise powers and the number of Fourier bins in the power spectrum. The significance given by presto is calculated in the same manner. However, it does not consider the fact that we split each light curve into several segments, and also searched for signals dispersed over an initially unknown number of Fourier bins $z_{\text{max}}$. Searching over a range of different $z_{\text{max}}$ values for each segmented light curve vastly increases the parameter space, making it more likely to obtain a false positive. Therefore, Monte Carlo simulations are required in order to account for these factors and derive an estimate of the true significance of our detections.

We simulated 10000 fake light curves containing pure Poisson noise with exposure times set to match those of the real observations, and mean count rates sampled randomly from the range of count rates observed in archival XMM-Newton and Swift observations of NGC 1313 X-2 (with the latter converted to equivalent count rates in XMM-Newton using pimms). Each light curve was split in the same manner as the real data (including the addition of GTI timing gaps), and was searched for pulsations using presto as described in Section 2.1. We then picked 12 simulated light curves at random to match the 12 long observations in the actual data (including archival observations). We looked for two detections of pulsations within these simulated observations with statistical significances equal to or larger than those measured for the real data, and detected pulse frequencies differing by no more than the difference between our two detections (i.e. $<1.5\times 10^{-3}$ Hz). We repeated this random sampling of our simulated light curves 10000 times. We found a false positive rate of 1 in 10000, implying that our detections are significant to more than $3.5\sigma$ even after accounting for the "look elsewhere" effect.

2.3 Constraints on orbital parameters

In order to set preliminary constraints on the orbital parameters of the system, we will conservatively assume that any variation in the pulsation frequency is dominated by an orbital modulation (although see Section 3). In this case, the line-of-sight acceleration of the source $\alpha$ can be expressed in terms of the spin parameters as $\alpha=c(\dot{f}/f$). The value of $\dot{f}$ measured in the sixth observation is $1.6(2)\times 10^{-8}$ Hz s*-1*, which corresponds to $\alpha=7(1)$ m s*-2*. In order to produce such a large acceleration, the companion mass $M_{c}$ and the orbital period $P_{\text{orb}}$ must be restricted to the range of values illustrated in the top panel of Figure 4. Moreover, if the observed pulse frequency shift between the two observations is solely due to the orbital motion of the source, this imposes further limits on the orbital parameters (as illustrated in the bottom panel of Figure 4). In deriving these estimates we have assumed that the pulsations indicate the presence of a NS primary in the binary system with a mass $\sim 1.4M_{\odot}$. We also assume a circular orbit, which appears to be the case for several other PULXs likely due to the action of strong tidal torques (Bachetti et al. 2014, Israel et al. 2017a, Fürst et al. 2016). We consider the constraints on the orbital parameters for two separate inclination angles, and find that it is particularly difficult to explain the observed frequency shift if the system is viewed moderately face-on (i.e. $<$ 30 degrees; second figure in the bottom panel). However, since the source is classified to be hard ultraluminous from its X-ray spectrum, its central regions are likely to be seen close to face-on rather than edge-on (Sutton, Roberts & Middleton 2013). This apparent contradiction can be reconciled if the spin and orbital angular momenta of the system are misaligned. Further, Grisé et al. (2008b) have derived an upper limit on the mass of the companion star from studies of the young OB association in its vicinity. By fitting stellar isochrones to the HST optical magnitudes of these stars, they were able to measure the cluster age to be $\sim$ 20 Myrs. This limits the maximum mass of the companion star of the ULX to be $\lesssim 12M_{\odot}$, provided that it is associated with the star cluster. A typical B-type companion consistent with this mass limit has a radius of $\sim 4R_{\odot}$ and must overflow its Roche-Lobe during the ULX phase. In such a scenario, the period-density relation (Frank et al. 2002) independently constrains the orbital period, which we find to be consistent with the range of values inferred from X-ray timing (provided that there is some spin-orbit misalignment). Finally, we can use the constraints derived from this analysis to verify that our definition of a false positive in the previous section (i.e. any two detections with pulse frequencies differing by $\sim 1.5\times 10^{-3}$ Hz or less) is not physically unreasonable. Indeed, for the allowed set of parameters shown in Figure 4, the expected change in the pulse frequency in a single orbital cycle will not exceed $5\times 10^{-3}$ Hz. If we re-run the Monte Carlo simulations by allowing the pulse frequencies to differ by as much as the latter value, we find that this does not change the false positive rate quoted in the previous section.

3 Discussion & conclusions

We report the detection of coherent pulsations in individual segments of two XMM-Newton observations of NGC 1313 X-2 (out of a total of twelve). The detected pulsations are notably weaker than for previous PULXs, with the measured pulsed fractions ($\sim 5\%$) being roughly a factor of two smaller. As a result, each individual detection is only marginally significant at $\sim 3-4\sigma$; however the two detections have remarkably similar pulsation frequencies, and we confirm using Monte Carlo simulations that the chance of detecting two false positives with such similar frequencies is extremely low, so we are confident that we are observing a bona fide new member of the PULX class. This is supported by the detected frequencies (that can be associated with the spin of the neutron star) and their derivatives, in addition to the sinusoidal pulse profiles, all of which are very similar to those of previously discovered sources (Bachetti et al. 2014; Israel et al. 2017a, b). The pulsations are transient on intra-observational timescales, which is reminiscent of the behaviour seen in M82 X-2 (Bachetti et al. 2014), although the physical mechanism for this behaviour remains poorly understood.

We were unable to measure an orbital modulation of the rotational frequency $f$ given the sparse number of pulsation detections. However, given the constraints on orbital acceleration provided by the frequency derivatives we were able to calculate the dependency of orbital period on mass of the companion to the NS; it is notable that we rule out periods as large as the 6.1 day period reported by Liu et al. (2009) on the basis of HST data for the optical counterpart to NGC 1313 X-2 (although the veracity of this period was questioned by Zampieri et al. 2012 on the basis of a wider dataset). The lack of constraints on the orbital modulation of the pulsations means that we cannot estimate secular (i.e. accretion-driven) changes in the pulsation frequency with any great certainty. Although, we can at least state that the observed frequency reduction of $\sim 1$ mHz over 100 days between the two pulse detections suggests a spin down rate $\dot{f}\sim-1.2\times 10^{-10}$ Hz s*-1*, which is notable as other PULXs all show secular spin increases. However, the instantaneous values of $\dot{f}$ measured during both epochs have the opposite sign and are significantly larger (see Table 1). Therefore, we caution that the observed frequency variation may well be dominated by orbital motion. We note that if the source is indeed spinning-down between the two observations, then this is probably not due to a transition to the propeller regime (Illarionov & Sunyaev 1975), as this is naturally accompanied by a dramatic decline in the X-ray luminosity by several orders of magnitude. Unlike other PULXs (Tsygankov et al. 2016; see also Earnshaw et al. 2018) no such diminution of luminosity has been seen for NGC 1313 X-2 (cf. Figure 1; also Weng & Feng 2018).

NGC 1313 X-2 is unusual amongst ULXs as the presence of a large ($\sim 400$ pc diameter), shock-ionised bubble nebula around it means that we can infer an age for the active phase of the ULX. Pakull & Grisé (2008) calculate this age as $\sim 1$ Myr for NGC 1313 X-2. This appears to contradict the assertion of King, Lasota & Kluźniak (2017) that PULXs are rare because the accretion of material onto the surface of the neutron star acts to rapidly diminish its magnetic field; NGC 1313 X-2 demonstrates that this is not the case. Indeed, current binary evolution models show that there are only a small number of scenarios in which stable Roche Lobe overflow persists for more than $\sim 1$ Myr onto a NS (e.g. Wiktorowicz et al. 2017). This likely implies that NGC 1313 X-2 is close to the end of its lifetime as a ULX, and yet it is still pulsating. As the bubble nebula requires an average input energy of $\sim 10^{40}$ erg s*-1* over its entire lifetime to inflate, from the kinetic output of the ULX wind/jets, this also suggests that NGC 1313 X-2 has sustained extreme ($\sim 100$ times Eddington) accretion rates from its secondary stellar companion. It is likely that the majority of this matter was expelled in radiatively driven outflows, as if this amount of matter were retained by the NS it would have accreted $\gtrsim 1M_{\odot}$ of material and potentially collapsed to a black hole, resulting in a cessation of pulsations. This necessity for both a strong outflow and a pulsating compact object supports the recent hybrid models for PULXs, such as described by Mushtukov et al. (2019) and is consistent with the recent detection of outflowing material from a PULX (Kosec et al., 2018).

Finally, we note that NGC 1313 X-2 is an archetype for the ’hard ultraluminous’ X-ray spectral class for ULXs, in which the object is thought to be viewed directly down the wind-delimited funnel to its central engine (Sutton et al. 2013). This observational configuration is obviously well-suited to detecting pulsations as the emission of the central regions is visible and/or beamed into an observer’s line of sight; notwithstanding the transient nature of the detected pulsations, or their faintness, the advanced age of NGC 1313 X-2 as a ULX offers some encouragement that observations with sufficient depth and baseline should detect pulsations from all those ULXs with a face-on aspect, if they contain NSs.

Acknowledgements

The authors wish to thank the referee for their useful comments that helped improve this paper. RS is grateful for the studentship from the STFC (ST/N50404X/1). TPR was funded as part of the STFC consolidated grant ST/K000861/1. HPE acknowledges support under NASA contract NNG08FD60C. CP is thankful for an ESA Research Fellowship, and DJW acknowledges support from an STFC ERF. This work is based on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA.