Next Article in Journal
Lensing Magnification Seen by Gravitational Wave Detectors
Previous Article in Journal
Position in Minimal Length Quantum Mechanics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Long-Term Study of Ultraluminous X-ray Sources in NGC 891

by
Nicholas M. Earley 
1,
Vikram V. Dwarkadas 
1,* and
Victoria Cirillo 
2
1
Department of Astronomy and Astrophysics, University of Chicago, 5640 S. Ellis Ave, Chicago, IL 60637, USA
2
Department of Physics and Engineering Physics, Fordham University, 441 E. Fordham Rd., Bronx, NY 10458, USA
*
Author to whom correspondence should be addressed.
Universe 2022, 8(1), 18; https://doi.org/10.3390/universe8010018
Submission received: 25 November 2021 / Revised: 17 December 2021 / Accepted: 21 December 2021 / Published: 29 December 2021
(This article belongs to the Section Compact Objects)

Abstract

:
We perform empirical fits to the Chandra and XMM-Newton spectra of three ultraluminous X-ray sources (ULXs) in the edge-on spiral galaxy NGC 891, monitoring the region over a 17-year time window. One of these sources was visible since the early 1990s with ROSAT and was observed multiple times with Chandra and XMM-Newton. Another was visible since 2011. We build upon prior analyses of these sources by analyzing all available data at all epochs. Where possible Chandra data is used, since its superior spatial resolution allows for more effective isolation of the emission from each individual source, thus providing a better determination of their spectral properties. We also identify a new transient ULX, CXOU J022230.1+421937, which faded from view over the course of a two month period from Nov 2016 to Jan 2017. Modeling of each source at every epoch was conducted using six different models ranging from thermal bremsstrahlung to accretion disk models. Unfortunately, but as is common with many ULXs, no single model yielded a much better fit than the others. The two known sources had unabsorbed luminosities that remained fairly consistent over five or more years. Various possibilities for the new transient ULX are explored.

1. Introduction

Ultraluminous X-ray sources (ULXs) are point-like sources found far from the nucleus of external galaxies, with X-ray luminosities that exceed 10 39 erg s 1 [1]. ULXs were first discovered in the 1980s by the Einstein Observatory, although the spectral resolution of ROSAT was needed to study their spectral properties [2,3]. Over the past 20 years, Chandra and XMM-Newton provided even greater insight into the spectral properties and potential identity of these sources [4]. ULXs were initially thought to be intermediate mass black holes (IMBHs) with masses between 10 2 10 6   M , under the assumption that their emission and accretion rates were isotropic and Eddington-limited [5,6,7]. If emission is instead arising from an accretion disk surrounding the compact object, which is presumably in a binary system, then the luminosty can exceed the Eddington limit [8,9,10]. The latter model, where a supercritical accretion disk forms, and accretion occurs in super-Eddington mode, seems to be more favored nowadays [10].
Theoretical models, and population synthesis studies, also seem to suggest that the majority of ULXs could be accreting highly magnetized NSs [11,12]. Mushtukov et al. [13] showed that pulsating ULXs were not strongly beamed, and therefore their observed luminosity must be close to their actual luminosity. King et al. [14] have argued that while only a small number of systems have shown pulsations till now, that is so because magnetic neutron stars will only emit observable pulses for a short fraction of their lifetime. They conclude that a much larger fraction of ULXs harbor neutron star accretors than the number of pulsating ULXs. It is increasingly likely that the population of ULXs as a whole comprises mainly of super-Eddington accreting compact objects, with several being confirmed as neutron stars (NSs) since 2014 [15,16,17,18,19]. Others could be super-Eddington accreting stellar-mass BHs, while the brightest could still potentially be IMBHs [10,20].
The X-ray spectra of ULXs show a curvature, with a sharp turnover above 5 keV and residuals below 1 keV [16,21]. A variety of models were used to fit the spectra and explore the physical nature of the systems [22]. The spectral features appear to be consistent with accretion onto highly magnetized NSs. Mushtukov et al. [23], for example, found that when the accretion rates were high, the resulting multicolor blackbody spectrum originating from the magnetospheric surface was suggestive of observed ULX spectra. Unfortunately, it is quite common to find degeneracy between spectral models suggestive of different physical processes [24]. The latter suggest that much deeper broadband spectra, and time resolved spectroscopy are needed to resolve discrepancies between models. Accurate spectral modeling is thus a critical first step in assessing emission mechanisms of ULXs.
Many ULXs show persistent large luminosities over several years, or even decades as we show herein. Others exhibit short timescale variability and can disappear over the course of months [25,26,27,28]. Strong variability in ULX luminosity could be attributable in some cases to NSs entering the “propeller” regime of accretion, in which accretion is cut off at the pulsar’s magnetospheric radius [29,30,31,32]. Observation of the propeller effect could allow a determination of the source magnetic field. Alternatively, transient sources have been proposed to be X-ray binaries such as low mass X-ray binaries (LMXBs) seen during outbursts [27,33].
Long-term monitoring of ULXs using a wide range of spectral models is required to fully determine and understand the nature of these sources. In this paper, we model the spectra of three ULXs in the edge-on, barred spiral galaxy NGC 891 using archival XMM-Newton and Chandra observations dating since 2000. NGC 891 is a nearby analog to the Milky Way in terms of luminosity, color, and morphology [34,35,36]. This particular region of interest was extensively observed because it hosts the supernova SN 1986J. Various ULXs appeared around the SN over time (Figure 1 and Figure 2). ULX-1 (source coordinates 02:22:31.36+42:20:24.4 [37]) was visible the longest, dating since the early 1990s when it was first observed with ROSAT [38]. ULX-2 (3XMM J022233.4+422026) was first discovered in a 2011 XMM-Newton observation by Hodges-Kluck et al. [39] (hereafter HK12). The spectral modeling presented by HK12 was based solely on XMM-Newton observations; however, contamination of the ULX-2 XMM-Newton spectrum from the nearby ULX-1 (within about 25"), as well as SN 1986J, is possible. HK12 used 6 different spectral models, characteristic of ULXs, to fit the data; we adopt this same set of models in our comprehensive study of all datasets. In this work, we also report the discovery of a new transient source detected in our 2016 Chandra observation, CXOU J022230.1+421937, whose name follows the naming convention for a newly discovered Chandra object. It is hereafter referred to as ULX-3. Using subsequent XMM-Newton observations, we show that this bright source decreased considerably in luminosity over the course of two months.
Of the 3 ULXs investigated in this work, ULX-2 appears to have been relatively well studied in the past [39,40,41]. However, spectral fits to the 2016 Chandra data were not previously reported. ULX-1, although visible for over three decades, has not been investigated as thoroughly. While flux measurements from the Chandra 2000 and 2011 data were reported [37,41,42], spectral fits using other Chandra and XMM datasets are not found in the literature. Except HK12, other authors generally used only 1 or 2 spectral models to fit the data. ULX-3 is a new transient which, so far as we can tell, was not previously reported. Altogether, about half the available data appears to have not been investigated earlier.
Irrespective of whether they were studied earlier or not, in this work we analyze all the available data on all the ULXs in this region using the same spectral models as HK12. We present a comprehensive and systematic analysis of all available XMM-Newton and Chandra data at all epochs on each of the 3 ULXs. The long-term monitoring of ULXs presented here, spanning a time baseline of almost two decades, at this level of spectral analysis, is quite rare. It represents a repository of data, all fitted using the same spectral models, that could prove extremely useful for further investigations or theoretical modeling. The Chandra data is essential in isolating the emission from each source, given the crowding of sources. With this large, rich, but difficult to analyze dataset, due to the crowding of the sources, our goal in this paper is to carry out the comprehensive spectral modeling of all three objects, at all available epochs, using both Chandra and XMM-Newton data, with spectral fitting using the 6 models as used by HK12. In a followup paper, we will present an in-depth temporal analysis of the light curve of each source at every epoch, exploring the time variability and searching for pulsations and quasiperiodic oscillations. We will also explore multiwavelength observations of the newly discovered ULX-3 in a bid to determine the nature of this transient.
In Section 2, we describe the XMM-Newton and Chandra observations of the sources; in Section 3, we perform X-ray spectral fits using thermal, power-law, and accretion disk models; and in Section 4, we discuss the results of the modeling and estimates of the fluxes which confirm these sources are ultraluminous. We also discuss the nature of the potential transient. Conclusions are given in Section 5. Throughout the paper, we use z = 0.001763 for the redshift of the galaxy [43], adopting the luminosity distance as d 9 Mpc [44], with H 0 = 71 km s 1 Mpc 1 , Ω Λ = 0.73 , Ω m = 0.27 [45].

2. Data Reduction and Analysis

The region containing the three ULXs and nearby supernova SN 1986J has been frequently observed since the advent of Chandra and XMM-Newton. It was observed four times between 2000 and 2016 with the ACIS-S instrument on Chandra and seven times between 2002 and 2017 with EPIC on XMM-Newton. The details of the observations are shown in Table 1. While the earliest Chandra observation of the region was made in 2000, ULX-1 was visible since at least 1991 when it was referred to as “XNorth” in ROSAT observations [38]. ULX-2 was visible in observations taken from 2011-2017. We identified a new source, CXOU J022230.1+421937, ULX-3, in a Chandra Nov 2016 observation.

2.1. XMM-Newton Observations

The region of interest was observed by XMM-Newton in 2002, 2011, and 2017. The data are reduced using version 18.0.0 of the XMM-SAS software following standard data reduction procedures for EPIC-pn, MOS1, and MOS2 spectra [46]. PN data are processed using epproc and MOS data using emproc. Light curves for each observation for each camera are filtered by count rate thresholds to eliminate instances of flaring. A circular 15 radius region is used to extract source spectra, while a circular 41.5 radius region, away from any source, is used for the background. Observations with the pn camera are shown for each epoch in Figure 1, with only the first of five 2017 observations being presented.
In the 2011 XMM-Newton observation, ULX-2 lies close to the edge of the chip. Only a small portion of the source extraction region encompasses the edge of the chip. Due to the high central brightness of the source, the contribution to the count rate from this edge is found to be essentially negligible.
The spatial resolution of XMM-Newton is not high enough to uniquely isolate the emission from any of the sources. ULX-1 and ULX-2 are only 24 away from each other, and ULX-3 is 24 away from SN 1986J. A 25 region in XMM-Newton observations on the other hand encircles 80% of the radiated energy. To minimize contamination from nearby sources, the centers of the 15 radius apertures are offset from the source coordinates. However, contamination of the emission of any source from a nearby source (ULX or SN) may always be present. For this reason, we rely more on the spectral fits with Chandra, with its 1 spatial resolution and a 90% encircled energy radius of 4 .

2.2. Chandra Observations

The same region was observed with Chandra ACIS-S in 2000, 2003, 2011, and 2016. The data are reduced and analyzed using the standard analysis pipeline in the CIAO software version 4.12 and CALDB 4.9.2.1 [47]. As seen in Figure 2, emission from the various sources is mostly distinct and non-overlapping in contrast to the XMM-Newton observations. However, the Chandra observations also reveal that the source termed ULX-1 actually appears to be two sources in proximity to each other. A 1.5 radius aperture is used to isolate the ultraluminous source at each epoch, which is highlighted in the 2000 observation in Figure 2. We note that it is not possible to isolate the ULX in the XMM-Newton observations in this manner. For the other ULXs, a 4 radius aperture centered on the source is used as the source region, with the background region having a 30 radius offset from each individual source. Spectra are extracted using the specextract command. While the Chandra observations in 2000, 2003, and 2016 had a large number of counts to provide high-spectral resolution for all of the sources, the 2011 epoch had <2 ks of data, not enough to provide a sufficient number of counts for all sources. For ULX-1, with three other epochs of high spectral resolution data available, we have chosen not to use the 2011 data since they are the least accurate. For ULX-2, however, with only one other epoch of Chandra data available, results from 2011 are also presented. Spectra from ULX-1 are obtained in 2000, 2003, and 2016; from ULX-2 in 2011 and 2016; and ULX-3 in 2016 only.

3. Results

The Chandra and XMM-Newton spectra are fit using Sherpa in CIAO version 4.12 [47]. Due to the high count numbers in nearly all epochs, spectra are background-subtracted and binned, and the chi2gehrels statistic is employed in each fit. The flux and errors for each model are calculated using the Sherpa command sample_flux with 1000 iterations. 1-sigma errors on best-estimate model parameter values are calculated when possible via the covar command, and with the conf command when covar fails. For XMM-Newton, to ensure the tightest constraints on the model parameters, we simultaneously fit the EPIC-pn, MOS1, and MOS2 spectra. Given the stability of the spectra of ULX-1 and ULX-2 over time, the pn, MOS1, and MOS2 spectra taken over five observations in 2017 are merged for each camera, using the SAS task epicspeccombine. The resulting merged pn, MOS1, and MOS2 spectra are then fit simultaneously as for other epochs.
ULX-2 was initially analyzed by HK12. Their empirical fits employed six different models appropriate to ULXs. In this paper, we fit the spectra of all the ULXs with the same models that they used to allow for direct comparison. HK12 solely considered the ULX-2 XMM-Newton data from 2011, which could suffer from contamination from nearby sources. Hodges-Kluck et al. [40] also fit the 2017 XMM-Newton spectra. This paper expands upon their work by investigating all ULXs in this region over all available epochs, including Chandra data which has the spatial resolution to isolate the emission from each source and reduce source contamination. All source models include a tbabs absorption component with abundances set to Wilms et al. [48] prior to fitting. Thermal bremsstrahlung (bremss) and broken power-law (bknpower) models are some of the simpler models used, describing thermal continuum emission or non-thermal power-law emission that could denote synchrotron or inverse Compton processes for example. Given that the nature of ULXs is unknown, but that they may be compact objects and their emission may arise from a disk around a stellar mass black hole or IMBH, the following accretion disk models are also used:
  • The diskpbb model. This is referred to as a “p-free” multicolor disk blackbody (MCD) model in which the accretion disk temperature T ( r ) r p . In the standard MCD model, each point of the disk is assumed to radiate like a blackbody with a temperature that scales as T ( r ) r 0.75 [49]. In this model p is allowed to vary, allowing for the possibility of advection-dominated disks ( p 0.5 ) in which photons cannot radiate away from the disk but are swept along with the accreting material [50]. In the case of super-Eddington accretion, this disk may become so optically thick that this advective process dominates.
  • The two-component diskpbb+comptt model. This is an MCD disk plus Comptonization model in which the seed photons originating from the inner disk with a Wien spectrum are Comptonized in a surrounding corona [51,52].
  • The diskir model. This is a model in which the inner regions of the disk are irradiated by Comptonized photons originating from the disk itself. Part of the resulting spectrum is due to the reflection of the Comptonized photons and the rest from the absorption and subsequent re-emission of such photons. This re-emission ultimately modifies the original disk emission had there not been any irradiation [53].
  • The kdblur2*reflionX model. This model describes emission arising from the reflection of radiated photons by an optically thick, constant density atmosphere [54,55]. The reflected spectrum is smoothed by the relativistic effects from the accretion disk surrounding a rotating black hole. The blurred reflection model is a table model that can be read into Sherpa using user-contributed scripts.
The results of the spectral fitting to the Chandra data, along with the absorbed and unabsorbed model flux in the 0.3–10.0 keV band are reported at each epoch in Table 2. XMM-Newton fits for ULX-1 and ULX-2 are in Table 3. Chandra and XMM-Newton model fits for ULX-3 are included and discussed in Section 4.3.

4. Discussion

4.1. ULX-1

ULX-1 was visible since at least the 1990s, when it was referred to as “XNorth” in observations with ROSAT [38]. However, the source still remains poorly characterized. It is difficult to distinguish between the various models, as shown comprehensively in Table 2 and demonstratively for the 2003 epoch in Figure 3. The fits appear to be statistically similar, such that there is no single best-fit model that stands out. This makes it difficult to delineate the radiative mechanism giving rise to the observed emission.
The spectral fitting with Chandra does indicate a generally consistent high column density over all models and epochs, with N H 10 22   cm 2 , in excess of the Galactic column of 0.07 × 10 22 cm 2 towards SN 1986J [56,57,58]. However, NGC 891 is viewed almost edge-on, and therefore one also needs to take into account the column density through the disk of that galaxy. This depends on the line of sight but could have a value of log N H = 20.7 [59], or could be much higher and almost comparable in value to the observed column density [60]. This also depends on whether the ULX is physically located near the SN or just happens to be viewed along a similar line-of-sight. As shown in Table 2, the ULX fits reveal a different column density for each ULX, which could mean that they are located at different depths in the galaxy, or that they have differing amounts of material surrounding them, or a combination of the two.
The column densities found for ULX-1 via fitting Chandra and XMM-Newton spectra (Table 2 and Table 3) are generally consistent in 2002, with N H 10 22 cm 2 . In 2011, the XMM N H decreases to almost half its 2002 value. In 2017, XMM gives slightly larger values for N H in comparison to 2011, but still less than the 2016 value found from the Chandra spectrum. We do not fit the Chandra spectrum of ULX-1 in 2011 due to the short exposure time of the observation. We find no decrease in N H in the Chandra results between 2000 and 2016. This discrepancy in N H between the Chandra and XMM data could possibly be due to the proximity of ULX-2 to ULX-1, given the different PSFs of the two telescopes. ULX-2 appeared in 2011 and persisted through 2017 at least. During those years, the large XMM PSF suggests that the flux from ULX-1 may have been contaminated by that from ULX-2, although we have tried to isolate the emission as best as possible. Fitting the data results in a column density which lies somewhere between that of ULX-1 and ULX-2.
The hot disk model (diskpbb) returns very high best-fit input temperatures. The value of p is often, but not always, close to 0.5, suggesting the possibility for advective flow in a disk. The temperatures are lowered by adding a Comptonization component, but are still > 1 keV and sometimes > 2 keV. Thus, the disk remains ‘hot’. The addition of the Comptonization component does not do much to improve the fit. The intrinsic absorption in a model where the inner disk is irradiated by Comptonized photons can be lower than the previous models, but both the disk and corona temperatures are high. None of the fits produce cool disks; namely, disk temperatures generally hover around 1 keV for both the Chandra and XMM fits. This model results in cooler disk temperatures, but a high corona temperature. The blurred reflection model returns a higher intrinsic absorption with a generally high emissivity index q i n , implying that disks are dominated by the emission from the inner disk, which falls off steeply with radius. The disk is generally found to have a high inclination angle, although the error bars can be large. The inner edge of the disk extends closer to the last stable orbit in the last Chandra observation, but not in the first two. The ionization parameter in all cases is extremely large, and the Fe abundance sub-solar. For these models with multiple parameters, the sample_flux command at times fails to converge. We have calculated the flux using calc_energy_flux, which unfortunately does not allow calculation of the error bars.
The Chandra models show some spectral evolution of the source from 2003 to 2016. The column density in all models is high. The unabsorbed flux (0.3–10.0 keV band) over all epochs has a value generally around 5 10 × 10 13 erg s 1 cm 2 . Using the ROSAT PSPC and HRI count rates in PIMMS, and assuming a power-law model with the same photon index found from the Chandra spectral fit, we calculate the unabsorbed fluxes from the August 1991, July 1993, and January 1995 data sets to be 5.6 × 10 13 erg s 1 cm 2 , 5.3 × 10 13 erg s 1 cm 2 , and 6.8 × 10 13 erg s 1 cm 2 respectively. These are consistent with the Chandra and XMM-Newton values.
The light curve shows a possible slight decrease in flux over time, particularly from 2000 to 2003, although the error bars are large enough that a steady flux over the entire time period is also valid. There does not seem to have been significant flux evolution over almost 30 years. At a distance of 9 Mpc, this flux corresponds to an isotropic luminosity of L X 5 10 × 10 39 erg s 1 , consistent with the luminosities of other ULXs [61]. Swartz et al. [42] analyzed the Chandra 2000 observation of ULX-1 assuming a power-law model. At the distance of 9 Mpc used in our work, their result corresponds to a luminosity of 8.4 × 10 39 erg s 1 , consistent with our derived values. Our flux values are also in agreement with those quoted by Dage et al. [41] for the 2000 observation. The long-term stability of the light curve suggests that the source is not a highly variable object over these timescales (Figure 4).

4.2. ULX-2

ULX-2 is not visible in either the 2000 or 2003 Chandra spectra. We calculate 90% confidence upper limits for non-detections, using a 4” radius source region, in the 2000 and 2003 data. Using the srcflux command in CIAO and assuming a simple redshifted power-law with N H = 0.2 × 10 22 cm 2 , we find that in the Chandra 2000 dataset, the absorbed and unabsorbed upper limits are 4.0 and 4.9 × 10 15 erg s 1 cm 2 , respectively. This is almost a factor of 200 less than the detected source flux in 2011. In 2003, the upper limits are similarly constraining, with absorbed and unabsorbed flux values near 1.3 and 1.4 × 10 15 erg s 1 cm 2 , respectively. If there was a source present at that position in these datasets, its luminosity was at least two orders of magnitude lower than when it was detected. HK12 first detected ULX-2 with XMM-Newton in 2011 and fit the spectrum. Due to XMM-Newton’s limited spatial resolution, it is possible that emission from ULX-2 could be contaminated by emission from the nearby SN 1986J and ULX-1. The much better spatial resolution of the Chandra satellite reduces or eliminates the contamination. The spectrum and fits using the various models for the 2016 epoch are shown in Figure 5. As before with ULX-1, all the fits are found to be similar, both to the eye and statistically. Excess emission near 1.8 keV is difficult to fit with any of the spectral models.
The absorption deduced for the two simpler models is comparable to that seen by HK12, but overall the similarities between the values are limited. The temperature returned by the thermal bremsstrahlung model in the Chandra fit is much higher than that found by HK12, although their value is comparable to our XMM-Newton fit. The broken power-law model is also more consistent with the XMM-Newton fit than with Chandra. The story is the same with the hot disk model, where the disk temperature is found to be higher with Chandra. Adding a Comptonized component does not decrease the temperature, although the XMM-Newton model is closer to the result of HK12. HK12 find a lower level of absorption in the diskir model, although our derived temperatures for both the disk and the corona are higher. The blurred reflection models in both the Chandra and XMM cases approach the results of HK12 in many of the parameters. Overall, even given the higher flux and the larger number of counts compared to ULX-1, it is difficult to select any model as being superior to the others. Similar to ULX-1, the temperatures returned by Chandra generally seem higher than that obtained with XMM. Our 2011 XMM-Newton values are in agreement with those derived by HK12, and the 2017 values are consistent with those obtained by Hodges-Kluck et al. [40] using a p-free blackbody disk model.
ULX-2 was observed in 2011 and 2016 with Chandra. The unabsorbed flux in the various models is about 6 12 × 10 13 erg s 1 cm 2 in the 0.3–10.0 keV band. Dage et al. [41] quote a flux of 8.8 ± 0.2 erg s 1 cm 2 using an absorbed power-law model, and 6.7 ± 2.0 erg s 1 cm 2 using an absorbed multiblackbody disk model for the 2011 observation, in agreement with our results. The flux is remarkably constant over most models in 2016 except for the last reflection model, which is higher than all the others but still within a factor of 2. The variations are larger in 2011, primarily in the flux calculated with the kdblur2*reflionX models, which is not surprising given the low spectral counts. The flux appears to be 20–50% higher than that of ULX-1, which is expected from the count rates. Again, while the best-fit values may reflect a decrease over 5 years, especially as seen by XMM, within the error bars the ULX may have had a consistent luminosity over the same time period (Figure 4). The derived column density is around 0.20 × 10 22 cm 2 with some variations. This value is lower than that found for ULX-1 by a factor of a few. This could either mean that the two are at different depths in the galaxy and not spatially close to each other, or it could mean that there is a differing amount of material around the two ULXs. The column density for ULX-2 is lower than any of the column densities calculated for the other ULXs.

4.3. ULX-3

Similar to ULX-2, ULX-3 is also not detected in 2000 or 2003 Chandra spectra. We calculate 90% confidence upper limits for non-detections in 2000 and 2003, using a redshifted power-law with N H = 2.0 × 10 22 cm 2 as found from spectral fitting in 2016. In 2000, the absorbed and unabsorbed upper limits are 2.9 and 3.5 × 10 15 erg s 1 cm 2 , respectively. In 2003, the absorbed and unabsorbed flux values are 1.1 and 1.2 × 10 15 erg s 1 cm 2 , respectively. Thus, the flux at the ULX-3 position at previous epochs is at least two orders of magnitude lower than the detection value. ULX-3 is first detected in the Chandra Nov 2016 observation. Fitting the Chandra spectrum of ULX-3 with the same set of models as used for the other sources results in an unabsorbed (0.3–10.0 keV) flux around 2 × 10 13 erg s 1 cm 2 . ULX-3 is the faintest of the three ULXs. At a distance of 9 Mpc, L X = 2 × 10 39 erg s 1 , which places the source in the lower luminosity range of detected ULXs. The absorption column on the other hand exceeds that found for ULX-1 or ULX-2 by more than a factor of two. The source temperature in the thermal model is extremely high ( k T > 20 keV). The reduced χ 2 is roughly equivalent between the thermal bremsstrahlung, power-law, and p-free multicolor disk blackbody models, but improves when more complicated disk models are considered. The fits for the Nov 2016 observation are shown in the last column of Table 2, to enable comparison with the Chandra spectra of the other sources, and the first column of Table 4, to demonstrate the source’s spectral evolution over time. Results for a redshifted power-law are shown in Table 4.
We simultaneously fit the EPIC-pn, MOS1, and MOS2 spectra, collected a couple months later. These results are shown in Table 4. We primarily consider the bremsstrahlung and power-law models, since the parameters can be well constrained, and error bars on all the parameters can be obtained. The flux and column density around the source are both found to decrease by a factor of seven from Nov 2016 to Jan 2017. At this time the source no longer qualifies as ‘ultraluminous’. This lower luminosity value is more consistent with other high-energy sources, such as X-ray binaries ( L X > 10 37 erg s 1 ) [62].
When considering the thermal model, the temperature may be declining from January to February, with k T decreasing from about 40 to 10 to 5 keV. However, large errors on the best-estimate temperature at each epoch could also suggest a constant value over the time period. In the power-law model, the power-law index remains constant within error bars, with an average value around Γ = 1.5 . For completeness, the parameters for the various disk models are also given in Table 4; however, error bars are unable to be calculated on many of the parameters since they are not well constrained. The column densities across models are consistently found to have decreased from 2016 to 2017, with the exception of the blurred reflection model which generally predicts higher N H ’s than the other models. Best-estimate values for unabsorbed model fluxes in the disk models also suggest that the source has dimmed over the duration of its observability.
The apparent dimming of ULX-3 over a period of two months raises interesting questions about its nature. ULXs are known to be highly variable and can have flux changes over a matter of hours to days [63,64,65]. The poor χ 2 in some of these fits may indicate that these models are not the best ones to use, possible if the source is not a ULX. Additional spectral modeling may help to decipher the nature of this emission.
Following Earnshaw et al. [63], we explore various possibilities that could explain the transient nature of ULX-3.

4.3.1. Supernova

ULX-3 could possibly be a new supernova (SN), in which case it exploded somewhere between 2011 and 2016. X-ray SNe can have luminosities from 10 36 to 10 42 erg s 1 [66], consistent with ULX-3’s luminosity. The Chandra spectrum of this source resembles a thermal spectrum, as has been found for the highest luminosity Type IIn SNe. Fits with both thermal and non-thermal models are however equally good. The observed decline in the flux of ULX-3 by nearly a factor of 7 over the course of two months is faster than that seen in any known SNe [66]. If the source is indeed a SN, this would suggest an extremely rapid decrease uncharacteristic of SNe. Type IIn SNe tend to decrease at a faster rate than the average SN, and SN 1986J decreases more rapidly than most SNe, at a rate L x t 1.69 [66]. On the other hand, the Type IIns tend to have luminosities an order of magnitude higher than ULX-3, at least in the first 2–3 years. Even assuming that ULX-3 had a much higher luminosity earlier, a 3- or 4-year-old source decreasing by nearly a factor of 7 in 2 months would make it the fastest decreasing SN ever. Thus, the rapid decrease makes it doubtful that it is a SN. Observations at other wavelengths could help considerably in this regard.

4.3.2. Super-Eddington Accreting Source

The majority of ULXs may be stellar mass objects accreting at super-Eddington rates. In particular, several were confirmed to be NSs, which can exhibit flux variability of an order of magnitude over short periods of time. A preliminary search for pulsations in the light curve of ULX-3, suggestive of NS super-Eddington accretion, proved negative. Given the various difficulties and issues in searching for pulsating sources (see for instance [67]), we have deferred a more in-depth study to a subsequent paper. A deeper Chandra observation could assist in the search for pulsations; with the current set of observations it is difficult to validate whether ULX-3 is a super-Eddington accretor in the ULX regime. ULX-3 lives in the lower luminosity ULX regime (<3 × 10 39 erg s 1 ). Such low luminosity sources are often “broadened-discs” undergoing accretion close to, and not much in excess of, the Eddington limit [68].

4.3.3. Transient Outbursts

Another possibility for ULX-3 is that it could be an X-ray binary. A simple redshifted power-law fit with the Chandra observation yields Γ = 1.4 ± 0.4 , typical of the sub-Eddington hard state [69,70]. Assuming a hard-only outburst in which observed luminosities generally only reach Eddington fractions of ∼10% [71], the source could be an IMBH with a mass upwards of ∼150 M . However, this would imply a much cooler disk temperature than is found in the fits with the various disk models, where T > 1 keV .
Transient outbursts such as in X-ray binaries could explain ULX-3. The suggested time scale for the dimming of the source is consistent with the time scale of transient stellar mass black hole LMXBs [27,72,73]. The initial decrease in column density could be indicative of the formation and dissipation of a disk around a compact object or the onset and aftermath of Roche-lobe overflow onto the compact object in a binary system. If an optical counterpart to the X-ray outburst is observed, it is possible that ULX-3 could be interpreted as an LMXB.

4.3.4. Microtidal Disruption Events

A microtidal disruption ( μ TDE) event occurs when a star is disrupted and accreted onto a stellar mass black hole or IMBH. This would be a one-time transit that would subsequently decay. The flux is found to decay with time as t 5 / 3 . Typical time scales of μ TDE can be of the order of a few months, consistent with the observations of ULX-3. If we assume that we happened to catch the outburst just when it happened, on day 1 (which would make the observation incredibly well timed), then a t 5 / 3 decline would mean that the source would have declined by a factor of a 1000 over 2 months, i.e., it would have completely disappeared. If the maximum flux was higher but reached earlier, say even 10 days earlier, even then the decline in 2 months would be substantial. Thus, given our current data, a μ TDE could be consistent.
Unfortunately, the paucity of data does not allow us to reach any conclusions regarding ULX-3. A second observation, with Chandra’s spatial resolution, is necessary to learn more about the nature of this intriguing source.

5. Conclusions

Three ULXs appeared over the course of nearly 30 years in the edge-on spiral galaxy NGC 891. We thus have a rare window into ULX spectral and flux evolution over an exceedingly long timescale. Previous groups identified two of these sources (ULX-1 from Houck et al. [38] in 1991 and ULX-2 from HK12 in 2011). Several groups carried out spectral modeling of some of the individual datasets; however, none of them modeled all the available data with a uniform set of models, and many datasets do not appear to have been studied previously. A comprehensive spectral modeling study of this region was not carried out prior to this work. We used the complete archival observational dataset for this region, from both XMM and Chandra. In the process, we have identified a third ULX that appeared in 2016 and whose flux decreased by a factor of seven over the course of two months. Fitting the source spectra with thermal, power-law, and accretion disk models, we are unable to distinguish between the spectral fits to designate any model as being superior to the others. This is similar to HK12’s findings for ULX-2. Additional multiwavelength observations and possible timing analyses are needed to provide more insight into their potential nature.
Our spectral fitting indicates that regardless of the origin of the emission, the two previously known sources are consistently ultraluminous and do not show large variability over the observational period. Long-term monitoring of ULX fluxes is rare, and this study provides evidence for consistent ultraluminous behavior over a time span of up to twenty years for these sources. Using Chandra’s high spatial resolution, we find the unabsorbed luminosity at 9 Mpc for ULX-1 to be approximately ( 5.8 ± 0.5 ) × 10 39 erg s 1 and for ULX-2 to be ( 9 ± 1 ) × 10 39 erg s 1 . Assuming isotropic emission and Eddington-limited accretion, these luminosities would correspond to 50 80   M black holes, placing these sources in the upper stellar black hole mass gap [74]. This would reinforce the interpretation of ULXs as simply an extension in the population of stellar mass black holes [75]. However, the high disk temperatures for these sources suggest that their spectra are not attributable to much heavier mass black holes. Multicolor disk model fitting in particular indicates that for T r p , p < 0.75 , suggesting that disk advection is important. Model fits for these two sources indicate exceedingly hot disk temperatures ( T in > 1 keV ) which could further suggest that these two sources could be super-Eddington accretors. This would be consistent with the increasingly paradigmatic view that ULXs as a population are compact objects undergoing super-Eddington accretion [9].
We also report the discovery of an additional high-energy source CXOU J022230.1 +421937 detected by Chandra in Nov 2016. The unabsorbed luminosity of this source, which we refer to as ULX-3, is near 2 × 10 39 erg s 1 , placing it on the lower boundary of the ULX luminosity range. Such a luminosity could correspond to a stellar mass black hole ∼20 M . Subsequent XMM observations show that over the course of a couple of months, ULX-3 decreased by a factor of about seven in luminosity, consistent with the behavior of a variety of astrophysical transients. Follow-up spectral modeling as well as longer exposure X-ray and multiwavelength observations could allow us to determine the nature of ULX-3 and its transitory behavior.

Author Contributions

Conceptualization, N.M.E., V.V.D. and V.C.; methodology, N.M.E. and V.V.D.; formal analysis, N.M.E., V.V.D. and V.C.; investigation, N.M.E., V.V.D. and V.C.; resources, V.V.D.; data curation, N.M.E. and V.C; writing—original draft preparation, N.M.E.; writing—review and editing, V.V.D. and N.M.E.; visualization, N.M.E.; supervision, V.V.D.; project administration, V.V.D.; funding acquisition, V.V.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work is funded by NSF grant 1911061 awarded to PI Vikram V. Dwarkadas at the University of Chicago and by the National Aeronautics and Space Administration through Chandra Award Number GO7-18066X issued by the Chandra X-ray Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf of the National Aeronautics Space Administration under contract NAS8-03060.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

We thank the anonymous referees for their careful reading of the paper, and informative comments and suggestions which have helped to substantially improve the paper. We thank Franz Bauer for helpful suggestions regarding ULXs. The scientific results reported in this article are based in part on observations made by the Chandra X-ray Observatory, and in part on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA. The paper used data obtained from the Chandra and XMM-Newton Data Archives. This research has made use of software provided by the Chandra X-ray Center (CXC) in the application packages CIAO, ChIPS, and Sherpa, and the XMM-SAS software.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Kaaret, P.; Feng, H.; Roberts, T.P. Ultraluminous X-ray Sources. Annu. Rev. Astron. Astrophys. 2017, 55, 303–341. [Google Scholar] [CrossRef] [Green Version]
  2. Long, K.S.; Dodorico, S.; Charles, P.A.; Dopita, M.A. Observations of the X-ray sources in the nearby SC galaxy M 33. Astrophys. J. Lett. 1981, 246, L61–L64. [Google Scholar] [CrossRef]
  3. Fabbiano, G. X-rays from normal galaxies. Ann. Rev. Astron. Astrophys 1989, 27, 87–138. [Google Scholar] [CrossRef]
  4. Feng, H.; Soria, R. Ultraluminous X-ray sources in the Chandra and XMM-Newton era. New Astron. Rev. 2011, 55, 166–183. [Google Scholar] [CrossRef] [Green Version]
  5. Colbert, E.J.M.; Mushotzky, R.F. The Nature of Accreting Black Holes in Nearby Galaxy Nuclei. Astrophys. J. 1999, 519, 89–107. [Google Scholar] [CrossRef] [Green Version]
  6. Kaaret, P.; Prestwich, A.H.; Zezas, A.; Murray, S.S.; Kim, D.W.; Kilgard, R.E.; Schlegel, E.M.; Ward, M.J. Chandra High-Resolution Camera observations of the luminous X-ray source in the starburst galaxy M82. Mon. Not. R. Astron. Soc. 2001, 321, L29–L32. [Google Scholar] [CrossRef] [Green Version]
  7. Gao, Y.; Wang, Q.D.; Appleton, P.N.; Lucas, R.A. Nonnuclear Hyper/Ultraluminous X-ray Sources in the Starbursting Cartwheel Ring Galaxy. Astrophys. J. Lett. 2003, 596, L171–L174. [Google Scholar] [CrossRef] [Green Version]
  8. Shakura, N.I.; Sunyaev, R.A. Reprint of 1973A&A....24..337S. Black holes in binary systems. Observational appearance. Astron. Astrophys. 1973, 500, 33–51. [Google Scholar]
  9. Poutanen, J.; Lipunova, G.; Fabrika, S.; Butkevich, A.G.; Abolmasov, P. Supercritically accreting stellar mass black holes as ultraluminous X-ray sources. Mon. Not. R. Astron. Soc. 2007, 377, 1187–1194. [Google Scholar] [CrossRef]
  10. Fabrika, S.N.; Atapin, K.E.; Vinokurov, A.S.; Sholukhova, O.N. Ultraluminous X-ray Sources. Astrophys. Bull. 2021, 76, 6–38. [Google Scholar] [CrossRef]
  11. Mushtukov, A.A.; Suleimanov, V.F.; Tsygankov, S.S.; Poutanen, J. On the maximum accretion luminosity of magnetized neutron stars: Connecting X-ray pulsars and ultraluminous X-ray sources. Mon. Not. R. Astron. Soc. 2015, 454, 2539–2548. [Google Scholar] [CrossRef] [Green Version]
  12. Kuranov, A.G.; Postnov, K.A.; Yungelson, L.R. Population Synthesis of Ultraluminous X-ray Sources with Magnetized Neutron Stars. Astron. Lett. 2020, 46, 658–676. [Google Scholar] [CrossRef]
  13. Mushtukov, A.A.; Portegies Zwart, S.; Tsygankov, S.S.; Nagirner, D.I.; Poutanen, J. Pulsating ULXs: Large pulsed fraction excludes strong beaming. Mon. Not. R. Astron. Soc. 2021, 501, 2424–2429. [Google Scholar] [CrossRef]
  14. King, A.; Lasota, J.P.; Kluźniak, W. Pulsing ULXs: Tip of the iceberg? Mon. Not. R. Astron. Soc. 2017, 468, L59–L62. [Google Scholar] [CrossRef] [Green Version]
  15. Bachetti, M.; Harrison, F.A.; Walton, D.J.; Grefenstette, B.W.; Chakrabarty, D.; Fürst, F.; Barret, D.; Beloborodov, A.; Boggs, S.E.; Christensen, F.E.; et al. An ultraluminous X-ray source powered by an accreting neutron star. Nature 2014, 514, 202–204. [Google Scholar] [CrossRef] [PubMed]
  16. Walton, D.J.; Harrison, F.A.; Grefenstette, B.W.; Miller, J.M.; Bachetti, M.; Barret, D.; Boggs, S.E.; Christensen, F.E.; Craig, W.W.; Fabian, A.C.; et al. Broadband X-ray Spectra of the Ultraluminous X-ray Source Holmberg IX X-1 Observed with NuSTAR, XMM-Newton, and Suzaku. Astrophys. J. 2014, 793, 21. [Google Scholar] [CrossRef] [Green Version]
  17. Middleton, M.J.; Roberts, T.P.; Done, C.; Jackson, F.E. Challenging times: A re-analysis of NGC 5408 X-1. Mon. Not. R. Astron. Soc. 2011, 411, 644–652. [Google Scholar] [CrossRef] [Green Version]
  18. Mukherjee, E.S.; Walton, D.J.; Bachetti, M.; Harrison, F.A.; Barret, D.; Bellm, E.; Boggs, S.E.; Christensen, F.E.; Craig, W.W.; Fabian, A.C.; et al. A Hard X-ray Study of the Ultraluminous X-ray Source NGC 5204 X-1 with NuSTAR and XMM-Newton. Astrophys. J. 2015, 808, 64. [Google Scholar] [CrossRef] [Green Version]
  19. Pinto, C.; Soria, R.; Walton, D.; D’Ai, A.; Pintore, F.; Kosec, P.; Alston, W.N.; Fuerst, F.; Middleton, M.J.; Roberts, T.P.; et al. XMM-Newton campaign on the ultraluminous X-ray source NGC 247 ULX-1: Outflows. arXiv 2021, arXiv:2104.11164. [Google Scholar] [CrossRef]
  20. Gladstone, J.C. The sub-classes of ultraluminous X-ray sources. Mem. Della Soc. Astron. Ital. 2013, 84, 629. [Google Scholar]
  21. Stobbart, A.M.; Roberts, T.P.; Wilms, J. XMM-Newton observations of the brightest ultraluminous X-ray sources. Mon. Not. R. Astron. Soc. 2006, 368, 397–413. [Google Scholar] [CrossRef] [Green Version]
  22. Gladstone, J.C.; Roberts, T.P.; Done, C. Unlocking the nature of ultraluminous X-ray sources using their X-ray spectra. Astron. Nachrichten 2011, 332, 345. [Google Scholar] [CrossRef]
  23. Mushtukov, A.A.; Suleimanov, V.F.; Tsygankov, S.S.; Ingram, A. Optically thick envelopes around ULXs powered by accreating neutron stars. Mon. Not. R. Astron. Soc. 2017, 467, 1202–1208. [Google Scholar] [CrossRef]
  24. Koliopanos, F.; Vasilopoulos, G.; Godet, O.; Bachetti, M.; Webb, N.A.; Barret, D. ULX spectra revisited: Accreting, highly magnetized neutron stars as the engines of ultraluminous X-ray sources. Astron. Astrophys. 2017, 608, A47. [Google Scholar] [CrossRef]
  25. Kaaret, P.; Feng, H. X-ray Monitoring of Ultraluminous X-ray Sources. Astrophys. J. 2009, 702, 1679–1682. [Google Scholar] [CrossRef] [Green Version]
  26. Grisé, F.; Kaaret, P.; Corbel, S.; Cseh, D.; Feng, H. A long-term X-ray monitoring of the ultraluminous X-ray source NGC 5408 X-1 with Swift reveals the presence of dips but no orbital period. Mon. Not. R. Astron. Soc. 2013, 433, 1023–1038. [Google Scholar] [CrossRef] [Green Version]
  27. Middleton, M.J.; Sutton, A.D.; Roberts, T.P.; Jackson, F.E.; Done, C. The missing link: A low-mass X-ray binary in M31 seen as an ultraluminous X-ray source. Mon. Not. R. Astron. Soc. 2012, 420, 2969–2977. [Google Scholar] [CrossRef] [Green Version]
  28. Burke, M.J.; Kraft, R.P.; Soria, R.; Maccarone, T.J.; Raychaudhury, S.; Sivakoff, G.R.; Birkinshaw, M.; Brassington, N.J.; Forman, W.R.; Hardcastle, M.J.; et al. The Fading of Two Transient Ultraluminous X-ray Sources to below the Stellar Mass Eddington Limit. Astrophys. J. 2013, 775, 21. [Google Scholar] [CrossRef]
  29. Illarionov, A.F.; Sunyaev, R.A. Why the Number of Galactic X-ray Stars Is so Small? Astron. Astrophys. 1975, 39, 185. [Google Scholar]
  30. Stella, L.; White, N.E.; Rosner, R. Intermittent Stellar Wind Acceleration and the Long-Term Activity of Population I Binary Systems Containing an X-ray Pulsar. Astrophys. J. 1986, 308, 669. [Google Scholar] [CrossRef]
  31. Tsygankov, S.S.; Mushtukov, A.A.; Suleimanov, V.F.; Poutanen, J. Propeller effect in action in the ultraluminous accreting magnetar M82 X-2. Mon. Not. R. Astron. Soc. 2016, 457, 1101–1106. [Google Scholar] [CrossRef]
  32. Erkut, M.H.; Ekşi, K.Y.; Alpar, M.A. Ultra-luminous X-ray Sources as Super-critical Propellers. Astrophys. J. 2019, 873, 105. [Google Scholar] [CrossRef] [Green Version]
  33. Earnshaw, H.P.; Roberts, T.P.; Sathyaprakash, R. Searching for propeller-phase ULXs in the XMM-Newton Serendipitous Source Catalogue. Mon. Not. R. Astron. Soc. 2018, 476, 4272–4277. [Google Scholar] [CrossRef] [Green Version]
  34. de Vaucouleurs, G.; de Vaucouleurs, A.; Corwin, Herold G.J.; Buta, R.J.; Paturel, G.; Fouque, P. Third Reference Catalogue of Bright Galaxies; Springer: New York, NY, USA, 1991. [Google Scholar]
  35. van der Kruit, P.C.; Searle, L. Surface photometry of edge-on spiral galaxies. II—The distribution of light and colour in the disk and spheroid of NGC 891. Astron. Astrophys. 1981, 95, 116–126. [Google Scholar]
  36. Garcia-Burillo, S.; Guelin, M. The distorted kinematics of molecular gas in the center of NGC 891. Astron. Astrophys. 1995, 299, 657. [Google Scholar]
  37. Berghea, C.T.; Weaver, K.A.; Colbert, E.J.M.; Roberts, T.P. Testing the Paradigm that Ultraluminous X-ray Sources as a Class Represent Accreting Intermediate-Mass Black Holes. Astrophys. J. 2008, 687, 471–487. [Google Scholar] [CrossRef]
  38. Houck, J.C.; Bregman, J.N.; Chevalier, R.A.; Tomisaka, K. Recent X-ray Observations of SN 1986J with ASCA and ROSAT. Astrophys. J. 1998, 493, 431–439. [Google Scholar] [CrossRef] [Green Version]
  39. Hodges-Kluck, E.J.; Bregman, J.N.; Miller, J.M.; Pellegrini, E. A New Ultraluminous X-ray Source in the nearby Edge-on Spiral NGC 891. Astrophys. J. Lett. 2012, 747, L39. [Google Scholar] [CrossRef] [Green Version]
  40. Hodges-Kluck, E.J.; Bregman, J.N.; Li, J. The Hot, Accreted Halo of NGC 891. Astrophys. J. 2018, 866, 126. [Google Scholar] [CrossRef] [Green Version]
  41. Dage, K.C.; Vowell, N.; Thygesen, E.; Bahramian, A.; Haggard, D.; Kovlakas, K.; Kundu, A.; Maccarone, T.J.; Strader, J.; Urquhart, R.; et al. Ultraluminous X-ray sources in seven edge-on spiral galaxies. Mon. Not. R. Astron. Soc. 2021, 508, 4008–4016. [Google Scholar] [CrossRef]
  42. Swartz, D.A.; Soria, R.; Tennant, A.F.; Yukita, M. A Complete Sample of Ultraluminous X-ray Source Host Galaxies. Astrophys. J. 2011, 741, 49. [Google Scholar] [CrossRef]
  43. Trentham, N.; Tully, R.B. Dwarf galaxies in the NGC 1023 Group. Mon. Not. R. Astron. Soc. 2009, 398, 722–734. [Google Scholar] [CrossRef] [Green Version]
  44. Temple, R.F.; Raychaudhury, S.; Stevens, I.R. X-ray observations of the edge-on star-forming galaxy NGC 891 and its supernova SN1986J. Mon. Not. R. Astron. Soc. 2005, 362, 581–591. [Google Scholar] [CrossRef] [Green Version]
  45. Spergel, D.N.; Bean, R.; Doré, O.; Nolta, M.R.; Bennett, C.L.; Dunkley, J.; Hinshaw, G.; Jarosik, N.; Komatsu, E.; Page, L.; et al. Three-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Implications for Cosmology. Astrophys. J. Suppl. Ser. 2007, 170, 377–408. [Google Scholar] [CrossRef] [Green Version]
  46. Gabriel, C. XMM-Newton Science Analysis System (SAS): Medium and long-term strategy. In Proceedings of the X-ray Universe 2017 Symposium, Rome, Italy, 6–9 June 2017; Ness, J.U., Migliari, S., Eds.; p. 84. [Google Scholar]
  47. Fruscione, A.; McDowell, J.C.; Allen, G.E.; Brickhouse, N.S.; Burke, D.J.; Davis, J.E.; Durham, N.; Elvis, M.; Galle, E.C.; Harris, D.E.; et al. CIAO: Chandra’s data analysis system. In Proceedings of the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Orlando, Florida, USA, 24–31 May 2006; Volume 6270, p. 62701V. [Google Scholar] [CrossRef]
  48. Wilms, J.; Allen, A.; McCray, R. On the Absorption of X-rays in the Interstellar Medium. Astrophys. J. 2000, 542, 914–924. [Google Scholar] [CrossRef] [Green Version]
  49. Mitsuda, K.; Inoue, H.; Koyama, K.; Makishima, K.; Matsuoka, M.; Ogawara, Y.; Shibazaki, N.; Suzuki, K.; Tanaka, Y.; Hirano, T. Energy spectra of low-mass binary X-ray sources observed from Tenma. Publ. Astron. Soc. Jpn. 1984, 36, 741–759. [Google Scholar]
  50. Abramowicz, M.A.; Czerny, B.; Lasota, J.P.; Szuszkiewicz, E. Slim Accretion Disks. Astrophys. J. 1988, 332, 646. [Google Scholar] [CrossRef]
  51. Gladstone, J.C.; Roberts, T.P.; Done, C. The ultraluminous state. Mon. Not. R. Astron. Soc. 2009, 397, 1836–1851. [Google Scholar] [CrossRef] [Green Version]
  52. Titarchuk, L. Generalized Comptonization Models and Application to the Recent High-Energy Observations. Astrophys. J. 1994, 434, 570. [Google Scholar] [CrossRef]
  53. Gierliński, M.; Done, C.; Page, K. Reprocessing of X-rays in the outer accretion disc of the black hole binary XTE J1817-330. Mon. Not. R. Astron. Soc. 2009, 392, 1106–1114. [Google Scholar] [CrossRef] [Green Version]
  54. Ross, R.R.; Fabian, A.C.; Young, A.J. X-ray reflection spectra from ionized slabs. Mon. Not. R. Astron. Soc. 1999, 306, 461–466. [Google Scholar] [CrossRef] [Green Version]
  55. Ross, R.R.; Fabian, A.C. A comprehensive range of X-ray ionized-reflection models. Mon. Not. R. Astron. Soc. 2005, 358, 211–216. [Google Scholar] [CrossRef]
  56. Dickey, J.M.; Lockman, F.J. H I in the galaxy. Ann. Rev. Astron. Astrophys 1990, 28, 215–261. [Google Scholar] [CrossRef]
  57. Kalberla, P.M.W.; Burton, W.B.; Hartmann, D.; Arnal, E.M.; Bajaja, E.; Morras, R.; Pöppel, W.G.L. The Leiden/Argentine/Bonn (LAB) Survey of Galactic HI. Final data release of the combined LDS and IAR surveys with improved stray-radiation corrections. Astron. Astrophys. 2005, 440, 775–782. [Google Scholar] [CrossRef]
  58. HI4PI Collaboration; Ben Bekhti, N.; Flöer, L.; Keller, R.; Kerp, J.; Lenz, D.; Winkel, B.; Bailin, J.; Calabretta, M.R.; Dedes, L.; et al. HI4PI: A full-sky H I survey based on EBHIS and GASS. Astron. Astrophys. 2016, 594, A116. [Google Scholar] [CrossRef] [Green Version]
  59. Das, S.; Sardone, A.; Leroy, A.K.; Mathur, S.; Gallagher, M.; Pingel, N.M.; Pisano, D.J.; Heald, G. Detection of the Diffuse H I Emission in the Circumgalactic Medium of NGC 891 and NGC 4565. Astrophys. J. 2020, 898, 15. [Google Scholar] [CrossRef]
  60. Aoki, T.E.; Hiromoto, N.; Takami, H.; Okamura, S. JHK Imaging of the Edge-on Spiral Galaxy NGC 891. Publ. Astron. Soc. Jpn. 1991, 43, 755–780. [Google Scholar]
  61. Kovlakas, K.; Zezas, A.; Andrews, J.J.; Basu-Zych, A.; Fragos, T.; Hornschemeier, A.; Lehmer, B.; Ptak, A. A census of ultraluminous X-ray sources in the local Universe. Mon. Not. R. Astron. Soc. 2020, 498, 4790–4810. [Google Scholar] [CrossRef]
  62. Fabbiano, G. Populations of X-ray Sources in Galaxies. Ann. Rev. Astron. Astrophys 2006, 44, 323–366. [Google Scholar] [CrossRef]
  63. Earnshaw, H.P.; Grefenstette, B.W.; Brightman, M.; Walton, D.J.; Barret, D.; Fürst, F.; Harrison, F.A.; Heida, M.; Pike, S.N.; Stern, D.; et al. A Broadband Look at the Old and New ULXs of NGC 6946. Astrophys. J. 2019, 881, 38. [Google Scholar] [CrossRef] [Green Version]
  64. Atapin, K.; Vinokurov, A.; Sarkisyan, A.; Solovyeva, Y.; Medvedev, A.; Fabrika, S.; Kostenkov, A. Transient Ultraluminous X-ray Source NGC5474 X-1. arXiv 2020, arXiv:2011.00795. [Google Scholar]
  65. Pintore, F.; Motta, S.; Pinto, C.; Bernardini, M.G.; Rodriguez-Castillo, G.; Salvaterra, R.; Israel, G.L.; Esposito, P.; Ambrosi, E.; Salvaggio, C.; et al. The rare X-ray flaring activity of the Ultraluminous X-ray source NGC 4559 X7. arXiv 2021, arXiv:2103.14541. [Google Scholar] [CrossRef]
  66. Dwarkadas, V.V.; Gruszko, J. What are published X-ray light curves telling us about young supernova expansion? Mon. Not. R. Astron. Soc. 2012, 419, 1515–1524. [Google Scholar] [CrossRef] [Green Version]
  67. Song, X.; Walton, D.J.; Lansbury, G.B.; Evans, P.A.; Fabian, A.C.; Earnshaw, H.; Roberts, T.P. The hunt for pulsating ultraluminous X-ray sources. Mon. Not. R. Astron. Soc. 2020, 491, 1260–1277. [Google Scholar] [CrossRef] [Green Version]
  68. Sutton, A.D.; Roberts, T.P.; Middleton, M.J. The ultraluminous state revisited: Fractional variability and spectral shape as diagnostics of super-Eddington accretion. Mon. Not. R. Astron. Soc. 2013, 435, 1758–1775. [Google Scholar] [CrossRef] [Green Version]
  69. Belloni, T.; Colombo, A.P.; Homan, J.; Campana, S.; van der Klis, M. A low/hard state outburst of XTE J1550-564. Astron. Astrophys. 2002, 390, 199–204. [Google Scholar] [CrossRef]
  70. Sidoli, L.; Paizis, A.; Mereghetti, S.; Götz, D.; Del Santo, M. XMM-Newton and INTEGRAL observations of the very faint X-ray transient IGR J17285-2922/XTE J1728-295 during the 2010 outburst. Mon. Not. R. Astron. Soc. 2011, 415, 2373–2378. [Google Scholar] [CrossRef] [Green Version]
  71. Tetarenko, B.E.; Sivakoff, G.R.; Heinke, C.O.; Gladstone, J.C. WATCHDOG: A Comprehensive All-sky Database of Galactic Black Hole X-ray Binaries. Astrophys. J. Suppl. Ser. 2016, 222, 15. [Google Scholar] [CrossRef] [Green Version]
  72. Mineshige, S.; Yamasaki, T.; Ishizaka, C. On the Exponential X-ray Decay in X-ray Novae. Publ. Astron. Soc. Jpn. 1993, 45, 707–713. [Google Scholar]
  73. Chen, W.; Shrader, C.R.; Livio, M. The Properties of X-ray and Optical Light Curves of X-ray Novae. Astrophys. J. 1997, 491, 312–338. [Google Scholar] [CrossRef] [Green Version]
  74. Woosley, S.E.; Heger, A. The Pair-Instability Mass Gap for Black Holes. arXiv 2021, arXiv:2103.07933. [Google Scholar] [CrossRef]
  75. Marchant, P.; Langer, N.; Podsiadlowski, P.; Tauris, T.M.; de Mink, S.; Mandel, I.; Moriya, T.J. Ultra-luminous X-ray sources and neutron-star-black-hole mergers from very massive close binaries at low metallicity. Astron. Astrophys. 2017, 604, A55. [Google Scholar] [CrossRef] [Green Version]
Figure 1. 2002, 2011, and 27 January 2017 EPIC-pn observations in 0.3–10.0 keV band. While only pn images are shown, pn and MOS spectra were simultaneously fit in the analysis. ULX-1 was visible since 2000, and ULX-2 appeared in 2011. SN 1986J is in the center of each image in between all sources. In the 27 Jan 2017 EPIC-pn observation, the position denoting ULX-3 does not appear to harbor a source as bright or brighter than SN 1986J.
Figure 1. 2002, 2011, and 27 January 2017 EPIC-pn observations in 0.3–10.0 keV band. While only pn images are shown, pn and MOS spectra were simultaneously fit in the analysis. ULX-1 was visible since 2000, and ULX-2 appeared in 2011. SN 1986J is in the center of each image in between all sources. In the 27 Jan 2017 EPIC-pn observation, the position denoting ULX-3 does not appear to harbor a source as bright or brighter than SN 1986J.
Universe 08 00018 g001
Figure 2. 2000, 2003, 2011, and Nov 2016 Chandra observations in the 0.3–10.0 keV band. A 1.5 radius aperture denotes ULX-1 in the 2000 observation. The short exposure time in 2011 obscures the presence of SN 1986J and ULX-1.
Figure 2. 2000, 2003, 2011, and Nov 2016 Chandra observations in the 0.3–10.0 keV band. A 1.5 radius aperture denotes ULX-1 in the 2000 observation. The short exposure time in 2011 obscures the presence of SN 1986J and ULX-1.
Universe 08 00018 g002
Figure 3. Fits to the ULX-1 Chandra spectrum in 2003, with thermal, power-law, and accretion disk models. Each model includes a thermal absorption component that was left free to vary. Spectrum is grouped by 30 counts. Best-fit parameters are given explicitly in Table 2. Residuals corresponding to each model indicate highly similar fits.
Figure 3. Fits to the ULX-1 Chandra spectrum in 2003, with thermal, power-law, and accretion disk models. Each model includes a thermal absorption component that was left free to vary. Spectrum is grouped by 30 counts. Best-fit parameters are given explicitly in Table 2. Residuals corresponding to each model indicate highly similar fits.
Universe 08 00018 g003
Figure 4. The unabsorbed model light curves for ULX-1 and ULX-2 using Chandra and XMM-Newton spectra. The derived flux for any particular model for ULX-1 is consistent in time within error bars, suggesting long-term stability. Due to the wide spread in unabsorbed flux in 2011, because of low statistics, it is not possible to determine whether ULX-2 reduced in luminosity from 2011 to 2017.
Figure 4. The unabsorbed model light curves for ULX-1 and ULX-2 using Chandra and XMM-Newton spectra. The derived flux for any particular model for ULX-1 is consistent in time within error bars, suggesting long-term stability. Due to the wide spread in unabsorbed flux in 2011, because of low statistics, it is not possible to determine whether ULX-2 reduced in luminosity from 2011 to 2017.
Universe 08 00018 g004
Figure 5. Similar to Figure 3, but for Chandra 2016 ULX-2 spectrum.
Figure 5. Similar to Figure 3, but for Chandra 2016 ULX-2 spectrum.
Universe 08 00018 g005
Table 1. 2000–2017 archival observations of NGC 891.
Table 1. 2000–2017 archival observations of NGC 891.
Obs. IDMissionDateExposure  a  (ks)Data Count Rate  a ( 10 2 Counts s 1 )
ULX-1ULX-2ULX-3
794Chandra1 November 200050.93.9
0112280101XMM-Newton22 August 20027.8/12.9/12.66.8/2.7/2.7
4613Chandra10 December 2003118.92.6
14376Chandra20 December 20111.80.78.4
0670950101XMM-Newton25 August 201194.8/112.6/113.87.9/2.7/2.839.9/13.2/11.9
19297Chandra14 November 201639.52.04.70.7
0780760101XMM-Newton27 January 201728.3/40.5/41.26.8/2.3/1.713.9/4.5/4.51.2/0.4/0.4
0780760201XMM-Newton29 January 201732.1/49.2/49.47.0/2.4/1.813.2/4.2/4.11.0/0.3/0.4
0780760401XMM-Newton19 February 201732.6/46.6/46.76.9/2.4/2.110.3/3.2/3.21.0/0.3/0.3
0780760301XMM-Newton23 February 201731.3/41.1/41.46.9/2.4/2.112.5/3.9/4.11.2/0.3/0.3
0780760501XMM-Newton25 February 201742.2/63.6/64.06.6/2.3/2.012.2/3.8/3.70.9/0.2/0.2
a For XMM-Newton, exposure times and count rates are given as those for EPIC pn/MOS1/MOS2, after filtering out background flaring.
Table 2. Chandra ULX spectral models.
Table 2. Chandra ULX spectral models.
ComponentParameterUnitsULX-1ULX-2ULX-3
2000200320162011 *20162016
TB ABS · BREMSS
TB ABS N H 10 22 cm 2 0.87 ± 0.06 0.84 ± 0.05 0.87 ± 0.16 0 . 17 0.11 + 0.15 0.16 ± 0.04 2.20 ± 0.16
BREMSS k T keV 6.62 ± 1.21 6.31 ± 0.88 10.06 ± 3.86 10.12 ± 5.66 4.56 ± 0.65 24 . 12 16.48 + a
χ 2 0.570.710.660.750.800.67
Absorbed Flux 10 13 erg s 1 cm 2 4.6 ± 0.5 2 . 8 0.3 + 0.2 3 . 5 0.6 + 0.5 7 . 4 2.9 + 2.6 5 . 0 0.6 + 0.5 1.4 ± 0.1
Unabsorbed Flux 10 13 erg s 1 cm 2 6 . 8 0.6 + 0.5 4.1 ± 0.3 4 . 9 0.6 + 0.5 8 . 5 3.3 + 3.0 6 . 0 0.5 + 0.6 2.1 ± 0.1
TB ABS · BKNPOWER
TB ABS N H 10 22 cm 2 0.90 ± 0.15 1.02 ± 0.02 0.90 ± 0.07 0 . 24 0.15 + 0.20 0.29 ± 0.08 1 . 65 0.69 + 0.92
BKNPOWER Γ 1 1.50 ± 0.37 1.89 ± 0.03 1.54 ± 0.03 1 . 62 0.33 + 0.39 1.98 ± 0.16 0 . 96 0.46 + 0.50
Γ 2 1.88 ± 0.12 1.77 ± 0.41 9 . 0 7.5 + a 6.45   a 2.14 ± 0.46 4 . 70 3.16 + a
BreakEkeV 2.04 ± 0.36 4.68   a 6.45 ± 0.40 12 . 86 12.83 + a 3.82 ± 1.51 5 . 21 2.08 + 0.88
χ 2 0.550.700.690.900.780.65
Absorbed Flux 10 13 erg s 1 cm 2 4 . 8 1.6 + 2.1 3 . 1 0.2 + 0.3 3 . 2 0.2 + 0.3 8 . 1 3.2 + 5.1 5 . 2 1.2 + 1.5 1 . 3 0.9 + 3.5
Unabsorbed Flux 10 13 erg s 1 cm 2 7 . 3 2.1 + 2.3 5 . 2 0.2 + 0.3 4 . 7 0.2 + 0.3 10 . 1 3.5 + 4.8 7 . 6 1.3 + 1.6 1 . 8 1.1 + 3.7
TB ABS · DISKPBB
TB ABS N H 10 22 cm 2 0.98 ± 0.06 0.97 ± 0.06 0 . 60 0.28 + 0.32 0.20 ± 0.18 0.28 ± 0.05 1.45 ± 0.81
DISKPBB T in keV 3 . 99 1.29 + a 4 . 32 1.26 + a 1 . 94 0.44 + 1.11 3.16 ± 1.63 3 . 87 1.22 + a 2.04 ± 0.84
p 0.533 ± 0.008 0.526 ± 0.007 0 . 73 0.14 + 0.55 0.58 ± 0.08 0.506 ± 0.008 1.0 ± 0.8
χ 2 0.570.690.630.820.760.65
Absorbed Flux 10 13 erg s 1 cm 2 5 . 8 3.5 + 4.4 3 . 6 2.3 + 2.6 3 . 2 2.5 + 8.4 5 . 8 5.2 + 31.3 5 . 4 3.4 + 4.8 2 . 0 1.8 + 8.8
Unabsorbed Flux 10 13 erg s 1 cm 2 9 . 4 5.8 + 7.1 6 . 0 3.8 + 4.4 4 . 0 3.1 + 10.9 6 . 7 6.0 + 37.1 7 . 7 4.9 + 7.2 2 . 7 2.5 + 11.4
2000200320162011 *20162016
TB ABS · ( DISKPBB + COMPTT )
TB ABS N H 10 22 cm 2 0.98 ± 0.06 1 . 45 0.64 + 0.11 1 . 10 0.33 + 0.17 0 . 51 0.38 + 0.36 0 . 27 0.02 + 0.12 1 . 82 1.79 + 1.56
DISKPBB T in keV 1.31 ± 0.53 6.6 ± 0.2 1.30   a 3.37 ± 2.20 9 . 84 9.82 + a 0.47 ± 9.4 × 10 5
p 0.53 (f) 0.53 (f) 0.53 (f) 0.58 (f) 0.51 (f)1 (f)
COMPTT T 0 keV 0.12 ± 0.07 0.10 ± 0.01 0 . 11 a + 0.10 0.10   a 0.33   a 0 . 35 a + 4.86
k T e keV 15.8   a 47.8   a 41 . 4 a + 625 33 . 9 a + 47.4 50.4   a 50.7   a
τ p 8 . 1 a + 60.2 0 . 01 a + 2.02 1 . 1 0.7 + a 0.01   a 1.25   a 2.57 ± 2.26
χ 2 0.570.740.891.090.820.96
Absorbed Flux 10 13 erg s 1 cm 2 5.2   b 3.0   b 3.8   b 7.7   b 5 . 9 3.6 + 4.6 0 . 6 0.5 + 1.6
Unabsorbed Flux 10 13 erg s 1 cm 2 5.2   b 3.0   b 3.8   b 7.7   b 8 . 3 5.0 + 6.4 0 . 7 0.5 + 1.8
TB ABS · DISKIR
TB ABS N H 10 22 cm 2 0.82 ± 0.06 0.96 ± 0.02 0.79 ± 0.02 0 . 12 0.10 + 0.32 0.22 ± 0.02 2.82 ± 0.22
DISKIR k T disk keV 0 . 82 0.31 + 0.46 0.55 ± 0.03 1.04 ± 0.23 1 . 13 0.19 + a 0.93 ± 0.22 0 . 17 0.10 + 0.12
Γ 1.73 ± 0.61 1.96 ± 0.07 2.0 (f) 1.01   a 2.02 ± 0.72 1.53 ± 0.20
L c / L d 1.54 ± 0.59 3 . 1 2.1 + a 1.42 ± 0.62 4 . 81 3.89 + a 8 . 48 7.35 + a 10.0   a
k T e keV 5.04 ± 1.17 82.8 ± 4.1 5.17 ± 0.09 5.01   a 17.85 ± 0.55 142.2 ± 33.9
f in 0.1 (f) 0.1 (f) 0.1 (f) 0.1 (f) 0.1 (f) 0.15 (f)
f out 2.0 × 10 5 a + 0.115 0.02 ± 0.01 3.0 × 10 5 a + 0.04 0.002   a 0.03 ± 0.01 4.2 × 10 6   a
r irr r in 1 . 01 0.01 + 1.00 1 . 03 a + 0.60 1.01 (f) 1.0001 (f) 1.0 ± 1.1 × 10 5 3.21   a
log r out log r in 3.84 ± 2.35 3.31 ± 0.95 3.5 (f)5 (f)5 (f)5 (f)
χ 2 0.600.750.771.180.830.995
Absorbed Flux 10 13 erg s 1 cm 2 5.0   b 3 . 2 0.8 + 1.3 4 . 1 3.2 + 6.6 8.9   b 5 . 3 3.8 + 10.0 1.4   b
Unabsorbed Flux 10 13 erg s 1 cm 2 5.0   b 5 . 1 1.3 + 2.0 5 . 6 4.4 + 8.3 8.9   b 6 . 6 4.7 + 13.3 1.4   b
2000200320162011 *20162016
TB ABS · KDBLUR 2 · REFLIONX
TB ABS N H 10 22 cm 2 1 . 30 0.12 + 0.17 1.35 ± 0.10 1 . 18 0.28 + 1.31 1.35 ± 0.32 0.63 ± 0.11 2 . 48 0.47 + 1.92
KDBLUR2 R in R G 9 . 5 a + 5.1 8.2 ± 0.2 1 . 98 a + 5.13 2.37 ± 0.47 1 . 24 a + 6.93 1 . 24 a + 8.17
R out R G 400 (f)400 (f)400 (f)400 (f)400 (f)400 (f)
ideg 85 . 3 41.2 + a 55.8 ± 8.6 84 . 5 9.2 + a 31.9 ± 15.8 52.9 ± 7.0 90 . 0 38.7 + a
q in 3 . 70 4.46 + a 6 . 31 5.87 + a 5 . 21 3.39 + a 8 . 94 6.52 + a 5.32 ± 0.59 9 . 46 16.3 + a
q out 3.0 (f) 3.0 (f) 3.0 (f) 3.0 (f) 3.0 (f) 3.0 (f)
R break R G 20.0 (f) 20.0 (f) 20.0 (f) 20.0 (f) 20.0 (f) 20.0 (f)
REFLIONX Γ 1 . 70 a + 0.54 1.85 ± 0.15 1 . 80 a + 0.48 1 . 70 a + 0.81 1 . 74 a + 0.34 1 . 70 a + 0.51
ξ erg cm s 1 1195 636 + 217 1214 ± 409 1100 98 + 202 1094 ± 547 1088 ± 558 1573  a
A Fe 0 . 41 a + 1.41 0.66 ± 0.47 0 . 1 a + 0.5 3.09   a 0 . 1 a + 2.0 1 . 01 a + 14.46
χ 2 0.600.760.751.310.870.99
Absorbed Flux 10 13 erg s 1 cm 2 4 . 5 3.8 + 8.7 2 . 7 1.6 + 2.6 3 . 2 2.1 + 2.4 4 . 1 3.8 + 9.4 5 . 7 4.3 + 7.4 1.4   b
Unabsorbed Flux 10 13 erg s 1 cm 2 12 . 0 9.5 + 18.9 6 . 4 3.8 + 5.7 8 . 5 5.1 + 5.7 14 . 7 10.7 + 22.9 12 . 1 8.4 + 11.9 1.4   b
Absorbed and Unabsorbed flux is calculated over the 0.3–10 keV band. a These parameters are not well constrained and errors are not reported. b These values were found with calc_energy_flux. * Due to the lower total number of counts in 2011 relative to the other observations, the spectrum of ULX-2 at this epoch was grouped by 10 counts.
Table 3. XMM-Newton ULX-1 and ULX-2 spectral models.
Table 3. XMM-Newton ULX-1 and ULX-2 spectral models.
ComponentParameterUnitsULX-1ULX-2
20022011201720112017
TB ABS · BREMSS
TB ABS N H 10 22 cm 2 0.71 ± 0.07 0.47 ± 0.02 0.67 ± 0.02 0.194 ± 0.004 0.178 ± 0.004
BREMSS k T keV 4.0 ± 0.7 5.9 ± 0.3 6.6 ± 0.3 3.3 ± 0.06 4.6 ± 0.1
χ 2 0.610.890.941.371.08
Absorbed Flux 10 13 erg s 1 cm 2 3.0 ± 0.5 4.3 ± 0.1 3.8 ± 0.1 16.7 ± 0.3 4.9 ± 0.1
Unabsorbed Flux 10 13 erg s 1 cm 2 4 . 8 0.6 + 0.7 5.8 ± 0.2 5.4 ± 0.1 21.5 ± 0.3 6.04 ± 0.09
TB ABS · BKNPOWER
TB ABS N H 10 22 cm 2 1.32 ± 0.31 0.36 ± 0.05 0.69 ± 0.03 0.248 ± 0.006 0.197 ± 0.008
BKNPOWER Γ 1 3.3 ± 0.9 1.1 ± 0.2 1.61 ± 0.04 1.90 ± 0.02 1.64 ± 0.03
Γ 2 2.2 ± 0.2 2.04 ± 0.07 2.6 ± 0.1 3.44 ± 0.09 2.9 ± 0.1
BreakEkeV 1.9 ± 0.3 2.0 ± 0.2 4.3 ± 0.2 3.8 ± 0.1 3.5 ± 0.2
χ 2 0.550.880.931.260.96
Absorbed Flux 10 13 erg s 1 cm 2 3 . 3 1.8 + 3.0 4 . 3 0.7 + 0.9 3 . 8 0.2 + 0.3 16.6 ± 0.4 4.9 ± 0.2
Unabsorbed Flux 10 13 erg s 1 cm 2 21 . 4 10.4 + 19.5 5.3 ± 0.1 5.6 ± 0.3 24.1 ± 0.5 6.2 ± 0.2
20022011201720112017
TB ABS · DISKPBB
TB ABS N H 10 22 cm 2 0.86   a 0.45 ± 0.02 0.61 ± 0.03 0.209 ± 0.005 0.164 ± 0.007
DISKPBB T in keV 2.4   a 2 . 1 0.1 + 0.2 2.2 ± 0.1 1.56 ± 0.02 1.72 ± 0.02
p 0.5   a 0.587 ± 0.005 0.61 ± 0.01 0.552 ± 0.003 0.592 ± 0.005
χ 2 0.590.910.881.210.95
Absorbed Flux 10 13 erg s 1 cm 2 3.1/3.7/ 3.6  * 4 . 2 1.6 + 1.8 3 . 7 0.9 + 1.1 16 . 4 1.2 + 1.3 4.8 ± 0.4
Unabsorbed Flux 10 13 erg s 1 cm 2 3.1/3.7/ 3.6  * 5 . 7 2.1 + 2.4 5 . 0 1.3 + 1.5 22 . 0 1.6 + 1.8 5.8 ± 0.5
TB ABS · ( DISKPBB + COMPTT )
TB ABS N H 10 22 cm 2 0 . 91 0.12 + 0.11 0.48   a 0.63 ± 0.01 0.22   a 0.26 ± 0.03
DISKPBB T in keV 1.59   a 2.13   a 2.1   a 1.33   a 1.70 ± 0.03
p 0.5 (f) 0.59 (f) 0.61 (f) 0.55 (f) 0.59 (f)
COMPTT T 0 keV 0.07   a 0.06   a 0.1   a 0.01   a 0.07 ± 0.01
k T e keV 23.4   a 44.3   a 51.1   a 48.5   a 42.9   a
τ p 1.04   a 0.65   a 1.00   a 0.81   a 0.018 ± 0.017
χ 2 0.610.910.901.410.92
Absorbed Flux 10 13 erg s 1 cm 2 2.3/3.4/ 4.0  * 4.5/4.5/ 4.5  * 3.8/4.1/ 3.6  * 16.5/16.1/ 16.0  * 4.8/5.0/ 5.0  *
Unabsorbed Flux 10 13 erg s 1 cm 2 2.3/3.4/ 4.0  * 4.5/4.5/ 4.5  * 3.8/4.1/ 3.6  * 16.5/16.1/ 16.0  * 4.8/5.0/ 5.0  *
TB ABS · DISKIR
TB ABS N H 10 22 cm 2 0.92   a 0.37   a 0.53 ± 0.10 0.175   a 0.17   a
DISKIR k T disk keV 0.42   a 1.0   a 1.0   a 0.89   a 1.16   a
Γ 5.0   a 2.1   a 2.15 ± 0.31 5.0   a 4.09   a
L c / L d 0.90   a 0.78   a 1.24   a 0.39   a 0.36   a
k T e keV 37.8   a 5.1   a 6.5 ± 4.5 5.0   a 407  a
f in 0.1 (f) 0.1 (f) 0.1 (f) 0.1 (f) 0.1 (f)
f out 0.1   a 0.02   a 0 . 02 a + 0.05 0.10   a 0.1   a
r irr r in 1.002   a 1.02 1.2 ± 0.1 1.01   a 1.02   a
log r out log r in 5 (f)5 (f)5 (f)5 (f)5 (f)
χ 2 0.570.860.921.500.97
Absorbed Flux 10 13 erg s 1 cm 2 3.1/3.7/ 3.7  * 4.3/4.4/ 4.5  * 3.8/4.0/ 3.7  * 16.6/16.0/ 16.1  * 4.9/5.1/ 5.0  *
Unabsorbed Flux 10 13 erg s 1 cm 2 3.1/3.7/ 3.7  * 4.3/4.4/ 4.5  * 3.8/4.0/ 3.7  * 16.6/16.0/ 16.1  * 4.9/5.1/ 5.0  *
20022011201720112017
TB ABS · KDBLUR 2 · REFLIONX
TB ABS N H 10 22 cm 2 1 . 52 0.31 + 0.26 2.27 ± 0.32 1 . 09 0.23 + 7.07 2.12 ± 0.07 1.57 ± 0.05
KDBLUR2 R in R G 8 . 2 a + 2.9 1 . 24 a + 0.13 4 . 43 a + 0.75 1.38 ± 4.2 × 10 6 1 . 38 a + 0.23
R out R G 400 (f)400 (f)400 (f)400 (f)400 (f)
ideg 38 . 4 a + 12.3 61.2 ± 5.0 61 . 3 a + 26.6 35.3 ± 2.4 30.8 ± 2.3
q in 10 . 0 18.4 + a 7.6 ± 1.2 8 . 9 2.6 + a 7.9 ± 0.4 6.6 ± 0.2
q out 3.0 (f) 3.0 (f) 3.0 (f) 3.0 (f) 3.0 (f)
R break R G 20.0 (f) 20.0 (f) 20.0 (f) 20.0 (f) 20.0 (f)
REFLIONX Γ 2 . 3 0.15 + 0.36 2.03 ± 0.07 1.90 ± 0.06 2.19 ± 0.02 1.8   a
ξ erg cm s 1 1229 1047 + a 1365 ± 257 1141 12 + 81 1492 ± 123 1143 ± 37
A Fe 0 . 90 0.47 + 0.88 0 . 11 a + 0.20 0 . 29 0.10 + 0.08 1.08 ± 0.21 1 . 06 0.08 + 2.20
χ 2 0.660.801.271.621.21
Absorbed Flux 10 13 erg s 1 cm 2 3 . 7 2.5 + 3.0 4 . 4 1.3 + 1.6 4 . 2 2.3 + 2.6 16 . 6 2.5 + 2.6 4.9 ± 0.5
Unabsorbed Flux 10 13 erg s 1 cm 2 9 . 0 6.2 + 6.8 10 . 8 3.1 + 3.9 9 . 9 5.5 + 6.0 37.8 ± 5.5 8.8 ± 1.1
Absorbed and Unabsorbed flux is calculated over the 0.3–10 keV band. a These parameters are not well constrained and errors are not reported. * Sample flux does not converge, values found via calc_energy_flux for pn/MOS1/MOS2.
Table 4. 2016 Chandra and 2017 XMM-Newton ULX-3 spectral models.
Table 4. 2016 Chandra and 2017 XMM-Newton ULX-3 spectral models.
ComponentParameterUnits14 Nov 201627 Jan 201729 Jan 201719 Feb 201723 Feb 201725 Feb 2017
TB ABS · BREMSS
TB ABS N H 10 22 cm 2 2.2 ± 0.2 0.30 ± 0.15 0 . 24 0.07 + 0.09 0.20 ± 0.07 0.29 ± 0.12 0.22 ± 0.08
BREMSS k T keV 24 . 1 16.5 + a 38 . 5 a + 82.4 12 . 7 5.1 + 13.5 9.49 ± 4.44 4.68 ± 2.05 7.08 ± 3.01
χ 2 0.670.890.540.670.490.92
Absorbed Flux 10 13 erg s 1 cm 2 1.4 ± 0.1 0 . 31 * 0 . 36 0.05 + 0.04 0 . 27 0.06 + 0.05 0 . 23 0.09 + 0.07 0.21 ± 0.05
Unabsorbed Flux 10 13 erg s 1 cm 2 2.1 ± 0.1 0 . 31 * 0 . 42 0.06 + 0.05 0.32 ± 0.06 0 . 30 0.10 + 0.08 0 . 26 0.06 + 0.05
TB ABS · ZPOWERLW
TB ABS N H 10 22 cm 2 2 . 2 0.7 + 1.0 0.29 ± 0.20 0.28 ± 0.10 0.26 ± 0.11 0.40 ± 0.19 0.28 ± 0.11
ZPOWERLW Γ 1.4 ± 0.4 1.28 ± 0.25 1.52 ± 0.18 1.63 ± 0.23 1.93 ± 0.32 1.71 ± 0.22
χ 2 0.680.900.590.680.530.97
Absorbed Flux 10 13 erg s 1 cm 2 1 . 4 0.9 + 1.7 0 . 32 0.13 + 0.18 0 . 39 0.11 + 0.14 0 . 29 0.09 + 0.15 0 . 26 0.11 + 0.15 0 . 24 0.08 + 0.10
Unabsorbed Flux 10 13 erg s 1 cm 2 2 . 2 1.2 + 2.0 0 . 36 0.14 + 0.19 0 . 47 0.12 + 0.14 0 . 36 0.11 + 0.14 0 . 40 0.14 + 0.17 0 . 31 0.08 + 0.09
TB ABS · DISKPBB
TB ABS N H 10 22 cm 2 1.45 ± 0.81 0.09   a 0 . 08 0.08 + 0.13 0.18 ± 0.10 0.07   a 0 . 07 a + 0.16
DISKPBB T in keV 2.0 ± 0.8 1.7   a 1.8 ± 0.4 2 . 5 1.0 + a 1.1   a 1.4 ± 0.4
p 1.0 ± 0.8 1.0   a 0.82 ± 0.21 0.61 ± 0.02 1.0 a 0.83 ± 0.26
χ 2 0.650.830.450.680.460.89
Absorbed Flux 10 13 erg s 1 cm 2 2 . 0 1.8 + 8.8 0.26 * 0 . 36 0.28 + 0.89 0.28 * 0.22 * 0 . 21 0.17 + 0.62
Unabsorbed Flux 10 13 erg s 1 cm 2 2 . 7 2.5 + 11.4 0.26 * 0 . 37 0.28 + 0.94 0.28 * 0.22 * 0 . 23 0.19 + 0.69
TB ABS · ( DISKPBB + COMPTT )
TB ABS N H 10 22 cm 2 1 . 82 1.79 + 1.56 0.34   a 0 . 13 a + 0.07 0 . 24 a + 0.13 0.11   a 0 . 18 a + 0.46
DISKPBB T in keV 0.47 ± 9.4 × 10 5 1.1   a 1.5   a 0.73   a 1.04   a 1.35   a
p 1 (f)1 (f) 0.82 (f) 0.61 (f)1 (f) 0.83 (f)
COMPTT T 0 keV 0 . 35 a + 4.86 0.1   a 0.05   a 0.1   a 0.08   a 0.06   a
k T e keV 50.7   a 59.0   a 41.5   a 31.5   a 26.6   a 38.8   a
τ p 2.6 ± 2.3 4.2   a 4.3   a 1.9   a 3.1   a 1.3   a
χ 2 0.960.930.510.820.541.04
Absorbed Flux 10 13 erg s 1 cm 2 0 . 6 0.5 + 1.6 0.31 * 0.34 0.21 + 0.36 0.38 * 0.25 * 0.27 *
Unabsorbed Flux 10 13 erg s 1 cm 2 0 . 7 0.5 + 1.8 0.31 * 0 . 50 0.28 + 0.51 0.38 * 0.25 * 0.27 *
TB ABS · DISKIR
TB ABS N H 10 22 cm 2 2.82 ± 0.22 0.99   a 0.09   a 0.17   a 0.09   a 0.07   a
DISKIR k T disk keV 0 . 17 0.10 + 0.12 1.0   a 0.92   a 0.90   a 0.68   a 0.94   a
Γ 1.5 ± 0.2 3.1   a 5.0 a 5.0   a 4.1   a 5.0   a
L c / L d 10.0   a 0.2   a 2.0 a 1.3   a 2.2 a 1.2   a
k T e keV 142.2 ± 33.9 435.1   a 7.3   a 6.7   a 12.8   a 5.1   a
f in 0.15 (f) 0.1 (f) 0.1 (f) 0.1 (f) 0.1 (f) 0.1 (f)
f out 4.2 × 10 6   a 7.7 × 10 4   a 0.004   a 0.03   a 1.3 × 10 5 a 0.003   a
r irr r in 3.21   a 1.03   a 1.17   a 1.02   a 1.59   a 1.68   a
log r out log r in 5 (f)5 (f)5 (f)5 (f)5 (f)5 (f)
χ 2 0.9951.320.500.800.541.02
Absorbed Flux 10 13 erg s 1 cm 2 1.4 * 0.19 * 0.35 * 0.28 * 0.22 * 0.20 *
Unabsorbed Flux 10 13 erg s 1 cm 2 1.4 * 0.19 * 0.35 * 0.28 * 0.22 * 0.20 *
ComponentParameterUnitsNov 14 2016Jan 27 2017Jan 29 2017Feb 19 2017Feb 25 2017
TB ABS · KDBLUR 2 · REFLIONX
TB ABS N H 10 22 cm 2 2 . 48 0.47 + 1.92 1.02 ± 0.13 2 . 24 0.33 + 0.93 1.37 ± 0.22 2 . 41 0.51 + 0.90
KDBLUR2 R in R G 1 . 24 a + 8.17 2.60 ± 0.68 1 . 78 a + 0.72 2.09 ± 0.12 1 . 78 0.45 + 0.41
R out R G 400 (f)400 (f)400 (f)400 (f)400 (f)
ideg 90 . 0 38.7 + a 85.4 ± 1.1 29 . 4 a + 9.2 45.8 ± 5.5 0 a + 28
q in 9 . 5 16.3 + a 9.4 ± 0.3 6 . 6 2.1 + 3.2 10 . 0 7.2 + a 6 . 6 1.0 + a
q out 3.0 (f) 3.0 (f) 3.0 (f) 3.0 (f) 3.0 (f)
R break R G 20.0 (f) 20.0 (f) 20.0 (f) 20.0 (f) 20.0 (f)
REFLIONX Γ 1 . 70 a + 0.51 1 . 70 0.19 + 0.40 1 . 78 a + 0.59 1.69 ± 0.13 1 . 60 a + 0.67
ξ erg cm s 1 1573  a 1004 473 + 325 1018 657 + 1317 1060 ± 124 1152 577 208
A Fe 1 . 01 a + 14.46 0 . 10 a + 0.45 19 . 3 14.6 + a 5.8 ± 2.5 7 . 8 5.5 + a
χ 2 0.991.050.400.780.93
Absorbed Flux 10 13 erg s 1 cm 2 1.4 * 0 . 29 0.16 + 0.26 0 . 20 0.17 + 0.71 0 . 29 0.10 + 0.12 0 . 16 0.14 + 0.25
Unabsorbed Flux 10 13 erg s 1 cm 2 1.4 * 0 . 72 0.41 + 0.51 0 . 72 0.56 + 1.38 0 . 76 0.27 + 0.34 0 . 62 0.46 + 0.96
Absorbed and Unabsorbed flux is calculated over the 0.3–10 keV band. a These parameters are not well constrained and errors are not reported. * Sample flux does not converge, values found via calc_energy_flux. XMM-Newton fluxes are reported for pn spectra. Spectrum on 23 Feb 2017 was not able to be fitted with the blurred reflection model.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Earley , N.M.; Dwarkadas , V.V.; Cirillo , V. A Long-Term Study of Ultraluminous X-ray Sources in NGC 891. Universe 2022, 8, 18. https://doi.org/10.3390/universe8010018

AMA Style

Earley  NM, Dwarkadas  VV, Cirillo  V. A Long-Term Study of Ultraluminous X-ray Sources in NGC 891. Universe. 2022; 8(1):18. https://doi.org/10.3390/universe8010018

Chicago/Turabian Style

Earley , Nicholas M., Vikram V. Dwarkadas , and Victoria Cirillo . 2022. "A Long-Term Study of Ultraluminous X-ray Sources in NGC 891" Universe 8, no. 1: 18. https://doi.org/10.3390/universe8010018

APA Style

Earley , N. M., Dwarkadas , V. V., & Cirillo , V. (2022). A Long-Term Study of Ultraluminous X-ray Sources in NGC 891. Universe, 8(1), 18. https://doi.org/10.3390/universe8010018

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop