Arxiv:2001.00642V2 [Astro-Ph.HE] 19 Feb 2020
Total Page:16
File Type:pdf, Size:1020Kb
Draft version February 21, 2020 Typeset using LATEX twocolumn style in AASTeX62 The (re)appearance of NGC 925 ULX-3, a new transient ULX Hannah P. Earnshaw,1 Marianne Heida,2, 1 Murray Brightman,1 Felix Furst,¨ 3 Fiona A. Harrison,1 Amruta Jaodand,1 Matthew J. Middleton,4 Timothy P. Roberts,5 Rajath Sathyaprakash,5 Daniel Stern,6 and Dominic J. Walton7 1Cahill Center for Astronomy and Astrophysics, California Institute of Technology, Pasadena, CA 91125, USA 2European Southern Observatory, Karl-Schwarzschild-Strasse 2, 85748, Garching bei M¨unchen,Germany 3Quasar Ltd for ESA, European Space Astronomy Centre (ESAC), Science Operations Department, E-28692, Villanueva de la Caada, Madrid, Spain 4Department of Physics and Astronomy, University of Southampton, Highfield, Southampton SO17 1BJ, UK 5Department of Physics, Centre for Extragalactic Astronomy, Durham University, South Road, Durham DH1 3LE, UK 6Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, CA 91109, USA 7Institute of Astronomy, Madingley Road, Cambridge CB3 0HA, UK (Received; Revised; Accepted) Submitted to ApJ ABSTRACT We report the discovery of a third ULX in NGC 925 (ULX-3), detected in November 2017 by 39 −1 Chandra at a luminosity of LX = (7:8 0:8) 10 erg s . Examination of archival data for NGC 925 reveals that ULX-3 was detected by±Swift×at a similarly high luminosity in 2011, as well as by 38 −1 XMM-Newton in January 2017 at a much lower luminosity of LX = (3:8 0:5) 10 erg s . With an additional Chandra non-detection in 2005, this object demonstrates a± high dynamic× range of flux of factor & 26. In its high-luminosity detections, ULX-3 exhibits a hard power-law spectrum with +0:2 Γ = 1:6 0:1, whereas the XMM-Newton detection is slightly softer, with Γ = 1:8−0:1 and also well-fitted± with a broadened disc model. The long-term light curve is sparsely covered and could be consistent either with the propeller effect or with a large-amplitude superorbital period, both of which are seen in ULXs, in particular those with neutron star accretors. Further systematic monitoring of ULX-3 will allow us to determine the mechanism by which ULX-3 undergoes its extreme variability and to better understand the accretion processes of ULXs. 1. INTRODUCTION 10 keV established the presence of a turnover and steep Ultraluminous X-ray sources (ULXs) are defined as power-law drop-off in the spectra of most ULXs above 5 keV (e.g. Bachetti et al. 2013; Mukherjee et al. 2015; extragalactic non-nuclear X-ray point sources with lumi- ∼ nosities in excess of 1039 erg s−1, the majority of which Walton et al. 2014, 2015a). are widely thought to be stellar-mass compact objects Explicit confirmation of the stellar-mass nature of at accreting at super-Eddington rates (for a recent review least a proportion of ULXs came from the detection of see Kaaret et al. 2017). One key piece of evidence that coherent pulsations in M82 X-2 (Bachetti et al. 2014), arXiv:2001.00642v2 [astro-ph.HE] 19 Feb 2020 these sources are accreting in a regime different to the allowing the compact object to be unambiguously iden- canonical sub-Eddington states was the detection of a tified as a neutron star. Since then, the population of high-energy turnover in their 5{10 keV spectrum, first known neutron star ULXs has risen to at least 7, both seen in high-quality XMM-Newton data (e.g. Stobbart through the detection of pulsations (F¨urst et al. 2016; et al. 2006; Gladstone et al. 2009) and later confirmed Israel et al. 2017a,b; Carpano et al. 2018; Rodr´ıguez by observations with NuSTAR, whose coverage above Castillo et al. 2019; Sathyaprakash et al. 2019) and through the detection of probable cyclotron resonance scattering spectral features (Brightman et al. 2018). Corresponding author: Hannah Earnshaw Their extreme luminosities imply accretion rates of hun- [email protected] dreds of times the Eddington rate for neutron stars, 2 Earnshaw et al. and they tend to exhibit pulsations with period 1 s, The transient nature of ULXs and the sparse temporal with sinusoidal pulse profiles suggesting that they∼ are X-ray coverage of most galaxies means that more ULXs not highly beamed. The low count rates due to being are being discovered all the time, without necessarily be- extragalactic sources, the presence of spin-up and or- ing `new' sources. Some discoveries, having never been bital modulation, and possible transience in pulsations detected before, have previous upper limits establish- make them very challenging to detect, however, and the ing that large increases in flux must have taken place proportion of ULXs that contain neutron star accretors (e.g. Pintore et al. 2018a; Earnshaw et al. 2019; Wang as opposed to black holes is still an open question (e.g. et al. 2019), although limited coverage generally means King & Lasota 2016; Middleton & King 2017), espe- that recent periods of high luminosity cannot be ruled cially since their broadband spectra are otherwise indis- out. Others have been detected at luminosities below tinguishable from the rest of the ULX population (e.g. 1039 erg s−1 before undergoing an increase in brightness Koliopanos et al. 2017; Pintore et al. 2017; Walton et al. that takes them into the ULX regime (e.g. Hu et al. 2018). 2018). Neutron star ULXs do share a property of extreme This paper investigates the appearance of a previ- long-term variability of various types, often over orders ously unknown ULX during the latest Chandra obser- of magnitude in flux. One example of such variability vation in 2017 of NGC 925, a spiral galaxy at 9.56 Mpc is sudden and large drops in brightness, leading to an (from Cepheids distances recorded in the NED Distances approximately bimodal distribution in flux. This may be Database1; Steer et al. 2017). The galaxy is known due to the onset of the propeller regime (e.g. Tsygankov to contain two previously studied ULXs (Pintore et al. et al. 2016b), in which the magnetospheric radius of the 2018b); we report on the discovery of a third ULX at magnetic field is larger than the corotation radius of position 02h 27m 20s.18, +33◦ 340 12.8400 (J2000) which the accretion disc, creating a centrifugal barrier to mass we designate NGC 925 ULX-3 (henceforth ULX-3 in accretion and causing the flux to drop correspondingly this paper). Upon investigation of archival data for (Illarionov & Sunyaev 1975; Stella et al. 1986). This may this source, we found that it was previously detected be a way of identifying candidate neutron star ULXs in by XMM-Newton at a sub-ULX luminosity, and was the absence of detected pulsations (Earnshaw et al. 2018; also detected at ULX luminosities by Swift in 2011. In Song et al. subm.). this paper we present the results of our analysis of the Several neutron star ULXs also exhibit superorbital spectral and long-term timing properties of ULX-3, and periods on the order of tens of days, varying by up discuss how it fits into the broader picture of extreme to factors of tens in flux, discovered through long-term variability in ULXs. monitoring using Swift (e.g. Walton et al. 2016; F¨urst et al. 2018; Brightman et al. 2019). Several mechanisms 2. DATA REDUCTION & ANALYSIS have been proposed for their cause, including preces- 2.1. Chandra sion of a warped disk obscuring the central source (e.g. Motch et al. 2014), or Lense-Thirring precession of the We observed NGC 925 using ACIS-S for 10 ks on 2017 supercritical accretion flow itself (e.g. Middleton et al. December 1, in order to obtain a precise localization 2018, 2019). Monitoring at a regular cadence for long of one of the other ULXs in the galaxy for compari- durations is important, as a high-amplitude superor- son to near-infrared data (Heida et al. in prep). We bital period with limited sampling may result in the reduced this observation (Obs. ID 20356; see Fig.1, source appearing to have a bimodal flux distribution, right) as well as an archival 2.2 ks observation taken which could then be interpreted as a propeller transi- 2005 November 23 (Obs. ID 7104) using ciao v4.10 tion (e.g. M82 X-2; Tsygankov et al. 2016b; Brightman with caldb v4.7.9. We used the wavdetect routine et al. 2019). It is also possible for ULXs to exhibit both to determine the location of the new source, which we h m s ◦ 0 00 of these effects { for example, NGC 5907 ULX-1, which found to be 02 27 20 .18, +33 34 12 .84 (J2000) with 00 has a 78 day superorbital period and also underwent statistical 1σ error of 0.03 . There were insufficient a greater∼ drop in flux potentially due to the propeller background sources for further astrometric correction effect in 2013 (Walton et al. 2015b, 2016; F¨urstet al. to be performed. We subsequently used this position 2017a; Israel et al. 2017a). It is worth noting that the to extract data products from all other observations. presence of a superorbital period in itself is not evidence Products were extracted from a circular region with ra- 00 of a neutron star accretor, as they can also be seen in dius 3 centered on the source, with an annular back- black hole systems, such as the well-studied Galactic source Cyg X-1 (e.g. Rico 2008). 1 https://ned.ipac.caltech.edu/Library/Distances/ The (re)appearance of NGC 925 ULX-3 3 Figure 1.