The Spectral and Temporal Properties of Transient Sources in Early-Type Galaxies
We report the spectral and temporal variability properties of 18 candidate transient and potential transient (TC and PTC) sources detected in deep multi-epoch Chandra observation of the nearby elliptical galaxies, NGC 3379, NGC 4278 and NGC 4697. Only one source can be identified with a background counterpart, leaving 17 TCs + PTCs in the galaxies. Of these, 14 are in the galaxy field, supporting the theoretical picture that the majority of field X-ray binaries (XRBs) will exhibit transient accretion for 75% of their lifetime. Three sources are coincident with globular clusters (GCs), including two high-luminosity candidate black hole (BH) XRBs, with =5.4, and =2.8 , respectively.
The spectra, luminosities and temporal behavior of these 17 sources suggest that the transient population is heterogeneous, including neutron star (NS) and BH XRBs in both normal and high-rate accretion modes, and super soft sources containing white dwarf binaries. Our TC and PTC detections are noticeably fewer that the number expected from the populations synthesis (PS) models of Fragos et al. (2009), tailored to our new Chandra pointings of NGC 4278. We attribute this discrepancy to the PS assumption that the transient population is composed of NS XRBs, as well as differences between the statistical analysis and error estimates used in the model and our observations.
Subject headings:galaxies: individual (NGC 3379) — galaxies: individual (NGC 4278) — galaxies: individual (NGC 4697) — X-rays: galaxies — X-rays: binaries
Chandra observations of elliptical and S0 galaxies have led to the widespread detection of low-mass X-ray binary (LMXB) populations (see review Fabbiano, 2006). The existence of these populations had been inferred from the first X-ray observations of early-type galaxies with the Einstein Observatory in the 1980s (Trinchieri & Fabbiano 1985; see Fabbiano 1989 and references therein), but their first detection required the sub-arcsecond resolution of Chandra (Sarazin et al., 2000). These populations of extra-galactic LMXBs provide a wider observational base for the understanding of LMXB formation and evolution, which has been debated since their discovery in the Milky Way (see e.g. Giacconi 1974; Grindlay 1984; review by Verbunt & van den Heuvel 1995; Piro & Bildsten 2002; King 2002; Bildsten & Deloye 2004). A key observational result that can constrain the nature of these LMXBs is their time variability, as luminous LMXBs ( 10 ) in relatively wide orbits are expected to exhibit transient behavior, according to the ‘standard’ formation scenario of LMXBs in galactic fields (Verbunt & van den Heuvel, 1995; Shahbaz et al., 1997; Piro & Bildsten, 2002; King, 2002; Wu et al., 2010). For sources originating as ultra-compact binaries formed in Globular Clusters (GC) instead, transient behavior is only expected at lower luminosities (Bildsten & Deloye, 2004). Until recently, only the most luminous extra-galactic LMXBs ( a few 10 in the 038.0 keV band) have been observed with Chandra, with limited time coverage, making the detection of transients problematic. Recent deep monitoring observations of the nearby elliptical galaxies NGC 3379 (Brassington et al. 2008; from hereafter B08), NGC 4278 (Brassington et al. 2009; hereafter B09) and NGC 4697 (Sivakoff et al. 2008a; Sivakoff et al. 2008b; hereafter S08) provide the means necessary to pursue this line of investigation.
Each of these galaxies have been observed with Chandra as part of monitoring campaigns; NGC 3379 and NGC 4278 were observed as part of a legacy study of LMXB populations, with two new observations of NGC 4278 performed in 2010 also included (PI: Fabbiano). Multi-epoch observations of NGC 4697 were obtained from the archive. All three galaxies have been determined to be old 10 Gyr111Terlevich & Forbes (2002) estimated ages for the central galactic region for all of three galaxies, based on spectral line indices calibrated to the Lick system, via interpolation with the results of the Worthey (1994) and Worthey & Ottaviani (1997) single stellar population evolutionary tracks. Further, for NGC 3379 and NGC 4278 SAURON data within the central region, again calibrated with the Lick indices and single stellar population synthesis models of Schiavon (2007), were used to determine age estimates (Kuntschner et al., 2010). For NGC 4697 Rogers et al. (2010) used spectral energy distributions of the stellar light, as well as absorption lines, which were then both compared to four different population models. This results in age ranges of: 9.314.7 Gyr, 10.714.1 Gyr and 8.111.0 Gyr for NGC 3379, NGC 4278 and NGC 4697 respectively., ensuring that the X-ray point source population is not contaminated by younger sources (e.g. high-mass X-ray binaries or supernova remnants). In this paper distances of 10.6 Mpc (NGC 3379), 16.1 Mpc (NGC 4278), and 11.8 Mpc (NGC 4697) are adopted. These are based on the surface brightness fluctuation analysis by (Tonry et al., 2001). At these distances, 1 corresponds to 3.1 kpc, 4.7 kpc, and 3.4 kpc, respectively. A summary of the properties of the galaxies and their Chandra observations is presented in Table 1, where column (1) gives the galaxy name, (2) distance, (3) the B-band luminosity, (4) the average stellar age, derived from single stellar population models, (5) globular cluster specific frequencies, (6) the Chandra observation ID, (7) the date of each pointing, (8) the cleaned exposure time and column (9) the number of point sources detected in each observation within the overlapping area covered by all pointings. In this table, in addition to the single pointings, information from the coadded (All) observations are also provided.
The catalogs of point sources detected in NGC 3379 and NGC 4278 are presented in B08 and B09 respectively. The point sources of NGC 4697 presented in this paper were detected and analyzed using the techniques detailed in B08, all of these sources are included in the catalog presented in S08. A study of the X-ray luminosity functions (XLFs) of these galaxies in both the field and in GCs is presented in Kim et al. (2009). Brassington et al. (2010) and Fabbiano et al. (2010) provide a detailed analysis of the bright persistent point source population of NGC 3379 and NGC 4278 respectively. PS modeling of the XLFs and transient population were pursued by Fragos et al. (2008) and Fragos et al. (2009). The new point sources detected in the 2010 observations of NGC 4278 have been analyzed using the techniques presented in B09, and will be presented in a forthcoming paper.
This paper is structured as follows: In Section 2 we describe the definition of transient candidate (TC) and potential transient candidate (PTC) sources, provide an overview of the properties of these sources from the B08 & B09 catalogs, and our new analysis of NGC 4697, and present our results. In Section 3 we present our discussion of these sources along with comparisons to population synthesis modeling. Our conclusions are summarized in Section 4.
2. Analysis and Results
2.1. Transient Candidate Selection
A transient candidate is usually defined to be a source that either appears or disappears, or is visible for only a limited amount of ‘contiguous’ time during observations, where the flux ratio between the peak ‘on-state’ emission and the non-detection upper limit, or ‘off-state’ is greater than a certain value (usually between 510, e.g. Williams et al. 2008). However, using only this ratio as a discriminator can lead to overestimating the number of transients by including lower luminosity sources, which have poorly constrained values. Following B08 and B09, we use instead the Bayesian model developed by Park et al. (2006)(see Section 2.4 of B08 for details), which estimates the uncertainties in the ratio, as well as lower-bounds222This lower-bound value accounts for the uncertainties in the source counts and upper limit values from the individual observations and thereby provides the lowest measure of ‘on’ to ‘off’ values for each source. We classify sources with a lower bound ratio ¿10 as TCs and sources with ratios between 5 and 10 as PTCs. We note here that a small number of Galactic sources have been observed to vary by ratios 10 but do not go into periods of quiescence, e.g. 4U 1705-440 (Homan et al., 2009) and 4U 0513-40 (Maccarone et al., 2010b), where we particularly note that this latter source has been confirmed to be an ultra-compact binary (Zurek et al., 2009).
From the multi-epoch observations of the galaxies presented, the Bayesian analysis leads to the detection of 18 TCs and PTCs, listed in Table 2. In this table column (1) indicates the source numbering used in this paper, and (2) indicates the original source number as presented in each catalog paper (B08, B09, S08 or ‘new’ for sources detected in the NGC 4278 cycle 11 observations). Columns (3) and (4) provide RA and Dec of each source, (5) indicates the observation ID(s) the source was detected in. (6) provides net counts for that pointing (0.38.0 keV), (7) the mode ratio and (8) the lower-bound ratio from the Bayesian analysis, (9) marks the source as a TC or PTC, (10) the long-term light curve behavior (see Section 2.2 for an explanation). In column (11) additional notes for each source are provided and (12) provides band magnitude (Vegamag) optical upperlimits for sources that have no detected optical counterpart. These upperlimit values are the 3 flux values centered on the RA and Dec values in columns (3) and (4). In B08 there was no optical coverage for source A8 but from recent HST observations it can now be confirmed that this source is coincident with a GC (with a separation 0.1), with and band (Vegamag) magnitudes of 23.4 mag and 21.9 mag respectively. This corresponds to a color of 0.8 (ABMAGs), which indicates that the cluster is blue333These values and the upper limit values presented in column 12 of Table 2 have been derived using the methods presented in Kundu et al. (2007).. The long-term light curves are shown in Figure 1 where the top, middle and bottom panels present the TCs and PTCs in NGC 3379, NGC 4278 and NGC 4697, respectively (In Figure 2 individual long-term light curves are presented for the four additional sources that are not presented in B08 or B09).
In summary, we observe 5 TCs and 3 PTCs in NGC 3379, 3 TCs and 5 PTCs in NGC 4278, and 1 TC and 1 PTC in NGC 4697. Using the HST WFPC2 and ACS observations of these galaxies we find that three transients reside in GCs, and 14 in the galaxy stellar field. The source A3 in the NGC 3379 field, is found to be coincident with a background object; we only list this source and report its properties for completeness.
2.2. Flux Variability
The long-term light curves (Figure 1) show three TCs with ‘on states’ of at least three months in NGC 3379. Of these sources, one has a maximum outburst time of 9 months (A1), while the other two (A6 & A8) could have maximum outburst times up to 4 years and 5 months, respectively. In NGC 4278, only source B7 was observed in two epochs (both in March 2010), with a minimum outburst time of 5 days and unconstrained maximum outburst time; the other TCs and PTCs were detected in single epochs, giving minimum outburst times of hours. Upper limits to the duration of the outbursts ranging from months to a few years can be placed on seven TC/PTCs, while the remaining seven were detected at the beginning or end of the monitoring observations and therefore the duration of the outburst is unconstrained. In Table 2, column 10 lists our estimated outburst durations.
Source A5 is denoted as only being in outburst in one observation in Table 2, although in B08 it was reported as being detected in three out of five observations. However, A5 was also flagged as a possible double source in B08, a conclusion supported by our further investigation (see Figure 3). The counts from observation 1 arise almost exclusively from within the top white circular region shown in this figure, whereas in the two subsequent observations the counts are centered in the lower green circle. The luminosity from observations 2 and 3 is more than a factor of 15 lower, suggesting that this secondary source is a variable lower flux object.
Short-term variability was investigated for all sources with net counts in a single observation, using both the Kolmogorov-Smirnov test (K-S test), and the Bayesian blocks method (BB) (Scargle 1998; see Section 2.4 of B08 for more details). Of the 18 transient sources presented in this paper, only A5 and A8 exhibited detectable short-term variability during their outburst.
The intraobservation variability of these two sources was investigated in more detail by extracting the lightcurves from the ‘on’ observations (and the results are discussed in sections 3.5.1 and 3.3 respectively). This analysis was performed with the CXC CIAO software suite (v4.2)444http://asc.harvard.edu/ciao and HEASOFT (v5.3.1), where barycenter corrected binned lightcurves were extracted with the CIAO tool dmextract. The binning for these two sources was selected to provide well constrained values, allowing variability to be identified, while remaining fine enough to ensure that any changes in flux were not masked by the selected binning scheme. Counts were extracted from circular source region files with radii as in B08. Background counts were extracted from a large elliptical source-free region, located within the overlapping area covered in all observations.
Hardness ratios (HR) and colors were derived for each source, as described in Section 2.3 of B08. They are listed in Table 3 along with log . These luminosity values were derived from the net count rate, assuming a power law shape of =1.7 and Galactic (see Section 2.1 in B08 for further details).
2.4. Spectral Analysis
Spectral extraction was performed for the ‘on’ observations for the TC/PTC sources. using the CIAO tool psextract. We used the same circular extraction regions as in B08 or B09; for the TC/PTC discovered in our new analysis, the radii are: 2.5″ (B7), 3.0″(B8), 2.2″ (C1) and 3.0″(C2); background counts were extracted from surrounding annuli, with outer radii 2-3 times larger, depending on the presence of nearby sources. In cases where the source region was found to overlap with a nearby source, the extractions region radius was reduced (to a minimum of 1.5″) and the area of the overlapping source was excluded.
The source spectra were fitted in XSPEC (v12.5.1). We used only data between 0.38.0 keV, to avoid the calibration uncertainties at the low energies, and the high cosmic ray background above 8.0 keV. All the spectra were fitted to two models that describe the properties of X-ray binary spectra well (see e.g. review by Remillard & McClintock 2006): the multicolor disc blackbody (DISKBB in XSPEC; hereon in referred to as DBB), and the power law (PO) models. For each model all parameters were allowed to vary freely, with the exception of the absorption column , described by the X-ray absorption model tbabs, which in cases where the best-fit value was below that of the Galactic absorption, was frozen to the Galactic value555Galactic was calculated with the tool COLDEN:
http://cxc.harvard.edu/toolkit/colden.jsp. This tool provides values of (2.79, 1.79 & 2.14) for NGC 3379, NGC 4278 and NGC 4697 respectively.. In instances where there were sufficient counts to bin the data to at least 20 counts per bin (allowing Gaussian error approximation to be used), the minimum method was used to fit the data. The Cash statistic (Cash, 1979) was otherwise used, however this statistical method has the disadvantage that it does not provide a goodness of fit measure like .
Table 4 summarizes the best-fit models, with the three galaxies listed separately. Column (1) contains the source number, (2) the observations used in each fit, (3) the net source counts from the spectral extraction, (4) indicates which spectral model was used, (5) the fit statistic ( and number of degrees of freedom , or the C statistics - indicated by ), (6) the null hypothesis probability (or goodness when using Cash). Columns (7), (8) & (9) present the best-fit values of the fit parameters with 1 errors for each interesting parameter (F denotes that the value was frozen): , for the power law model and , the temperature of the innermost stable orbit of an accretion disc, for the multicolor disc model. Column (10) indicates the intrinsic value of and (11) the luminosity range from the lowest (1) lower bound to the highest (1) upper bound of each observation included in the joint fit666Throughout this paper luminosities quoted are the 0.38.0 keV Chandra bandpass. corrections for the HS sources will be between 515 (Portegies Zwart et al., 2004) and corrections for SSS/QSS will be between 1.53 (e.g. McGowan et al. 2005). Bolometric corrections for sources in a TD state are small, with 90% of the emission occurring between 0.38.0 keV (Portegies Zwart et al., 2004)..
Figure 4 shows that the majority of TC/PTCs have both HR and colors consistent with those of the persistent LMXBs in the galaxies. Five of these sources have insufficient counts to allow even simple spectral modeling (A2, A7, B5, C2 & the background source A3, which will not be discussed further; see Table 4) and therefore single-component models with canonical values have been applied ( or =1.0 keV Remillard & McClintock 2006). Six more sources (A1, A4, B1, B4, B6 & B7) have spectra well described by single-component models. Four of these (A4, B1, B4 & B7), are well fitted with a PO model with 1.7 and Galactic . A1, instead requires a large value of in the PO fit. Based on the simulations of Brassington et al. (2010), this result may favor disk emission. The single-component PO model of source B6 also shows an elevated value of , however, the two-dimensional errors on this value indicate that the best-fit is consistent with Galactic absorption (Figure 5).
Sources B2 & B3, have colors and HR values softer than the persistent LMXB population (Figure 4). Their spectra are well fitted by both a single PO and single DBB model. However, the PO model in both cases requires an elevated value of , indicative of the source containing a significant disc component (Brassington et al., 2010). In the case of B2, the best-fit value of from the DBB model tends to zero when left free to vary, indicating that there is likely a non-thermal component as well as disc emission in this source (Brassington et al., 2010). When freezing the absorption column to the Galactic value, the best-fit disc temperature is 0.41 keV, with an X-ray luminosity of 1.0. The spectrum of source B3 is well described by a DBB model with an absorption value more than twice that of the Galactic absorption value (albeit poorly constrained), indicating that more than 75% of the emission is arising from the disc component, with a temperature of 0.55 keV.
Four further sources (A5, A6, B8 & C1; see Figure 4 to compare with the persistent sources), all have little emission above 2 keV, and are therefore super soft, or quasi-soft, sources (SSS and QSS respectively). In all cases, the PO fits yielded extremely steep power-law indices (¿6), reflecting the very soft nature of the emission. Their spectra were also fitted with black-body models (BB), as reported in Table 4.
3.1. The Field Transient Population
One of the clear results from this work is that TC/PTC sources are predominantly found in the field with only three sources out of the 17 TC/PTCs (excluding the background source) coincident with a GC777In comparison observations of the total X-ray binary population suggest that between 30%70% of LMXBs are coincident with a GC (e.g. Angelini et al. 2001, B08, B09). We also note that the optical coverage of all three galaxies provided a representative view of the GC population of each of the systems (see B08, B09 and Kim et al. 2009 for further details).. These numbers support the theory which suggests that the LMXB field population largely arises from relatively detached systems evolved from native field binaries and as such will exhibit transient accretion for 75% of their lifetime (Piro & Bildsten, 2002). This is as opposed to transient GC-LMXBs, which have a more difficult formation channel (as discussed in Section 3.2).
From the detailed spectral analysis presented in this paper we have been able to determine that out of the 14 confirmed field sources, 10 exhibit colors and spectra consistent with typical emission from LMXBs, with power-law ( 1.32.0) or disk ( 0.41.0 keV) emission and on-state intrinsic luminosities ranging from 4 to 2.32 (see Table 4). The spectra of these sources are consistent with those of NS or BH binaries in either a hard state (HS) or thermally dominant (TD) state (Remillard & McClintock, 2006).
Only one of these sources has a luminosity significantly above the Eddington limit for a 1.4 neutron star (B6) and is a strong BH candidate. A further three sources that have been determined to be in a TD state (A1, B2, B3) exhibit inner-disc temperatures that are softer than the spectra typically observed in Galactic NS-LMXBs ( 12 keV), indicating that they could be BH candidates. A further three field sources (A4, B4 & B7) have been determined to be in a HS, but have X-ray luminosities 10% of the Eddington limit. This high accretion is not observed in NS-LMXBs in a HS and suggests that the primary of these binaries may be BHs. However, due to the quality of our data, while a HS is the preferred interpretation, we cannot rule out any of these sources emitting in a TD state, and therefore cannot exclude any of these sources from being NS-LMXBs. The remaining three sources have low-count data and have values derived from canonical models, therefore these sources could be NS or BH-LMXBs. A further soft source (B8) has been classified as a QSS but is also consistent with a NS-LMXB (or WD-LMXB; see section 3.5). This results in a ratio of BH to NS-LMXBs in transient field sources ranging from at least 10% up to 90% of sources. Such an unconstrained ratio is consistent with a study of the transient population of M31 Williams et al. (2006) which suggest that the majority of the sources that they identify to be LMXBs contain BHs.
3.1.1 Comparison with PS models of field LMXBs
Population synthesis modeling of the LMXB populations of NGC 3379 and NGC 4278 were presented in Fragos et al. (2008) and a study comparing the theoretical and observed field transient populations of these galaxies was provided in Fragos et al. (2009). From this work it was determined that a constant duty cycle (DC) for all transient systems did not match the data and instead a variable DC for each system provided results that were consistent with the observations. Furthermore Fragos et al. (2009) suggested that the observed number of TCs and PTCs is proportional to the total number of field LMXBs (or the stellar mass of the galaxy) and not the total observed LMXB population (which will be enhanced by GC-LMXB formation; Kim et al. 2006).
These models were compared to the transient population of NGC 4278 prior to the 2010 observations and as such enabled us to estimate the number of new TC and PTC sources that would be observed with repeat pointings. Based on the models’ prediction the additional 2010 observations of NGC 4278 would reveal an additional 6 TCs and 6 PTCs. However no TCs and only two new PTCs have been determined from these recent observations. The discrepancy between the model predictions and the observed number of new sources can be attributed to three main factors. Fragos et al. (2009) assumed in their modeling calibration that all TCs and PTCs determined from B08 & B09 were NS-LMXBs. This assumption was based on the fact that NS-LMXBs were the dominant population in their population synthesis models, and the lack of detailed characterization of the observed transient population. From the spectral modeling that has been presented in this paper, it can now be determined that only 2 TCs and 1 PTC, and 1 TC and 3 PTCs in NGC 3379 and NGC 4278 respectively are classified as NS-LMXBs (compared to 5 TC and 3 PTC, and 3 TC and 3 PTC previously). Combining this fact with the low number statistics used in the model-observation comparison by Fragos et al. (2009) can possibly explain the over-prediction of their models.
Furthermore, Fragos et al. (2009) used the ‘on-to-off’ flux ratio of a source in order to characterize it as a TC or PTC, instead of the more conservative Bayesian analysis employed here. Finally, we should also consider the uncertainties of the bolometric corrections used by Fragos et al. (2008, 2009) in order to convert mass-transfer rates to X-ray luminosities, which would significantly increase the error bars in the predicted number of transients. These last two factor can explain the discrepancy between the model predictions and observations based solely on statistical arguments.
Three of the TC/PTCs presented in this work are in GCs. It has been suggested that there is a trend for GCs to have a higher fraction of ultra-compact binaries, compared to the field LMXB population (Deutsch et al., 2000; Heinke, 2010). However for sources brighter than such objects are expected to be persistent (Lasota et al., 2008; Bildsten & Deloye, 2004) and are therefore not the transient sources presented here. Instead these sources could be NS with main sequence or red giant donors, where population synthesis models of the formation of LMXBs in GCs predict that these systems, in addition to ultra-compact binaries, also contribute to the GC-LMXB population (Ivanova et al., 2008). Such sources were shown to exhibit transient emission from 60% to all of their mass-transfer lifetimes. We therefore suggest that A2, a source with low counts which has therefore been described with a canonical PO spectrum, is likely to be a NS-main sequence or NS-red giant system with (although NS-red giants are less favored as Ivanova et al. (2008) predicted that their numbers are low when considering their formation rates and short lifetimes).
The second GC-LMXB, B1, has a spectrum that is well described by both a PO and DBB model, however the best-fit from the PO model indicates that the source is emitting greatly above the 24 percent Eddington luminosity that is typically seen from a LMXBs in a HS (Maccarone, 2003) and the DBB model is therefore preferred888It is also noted that while the absorption column of the PO model indicates a value consistent with the Galactic value, and therefore indicative of a HS (c.f. Brassington et al. 2010), this value is poorly constrained and consequently does not rule out the presence of a thermal component.. The best-fit values from the DBB model determines an inner disc temperature of =1.3 keV and =5.4. Therefore, because of this high luminosity, which is above the Eddington limit for a NS binary (although could be close to the limit of a heavy (23 ; Kalogera & Baym 1996) NS, or a 1.4 NS with a He or C/O donor), it is possible that B1 is a BH binary instead of a NS-main sequence system.
The third GC-LMXB in this sample is A8 (presented in section 3.3), which could be a BH-ULX undergoing a period of super-Eddington accretion.
The detection of a BH binary in a GC is thought to be rare, since in tidally captured persistent sources BHs are likely to be expelled from the GC (Spitzer, 1969); they could be transients with very low duty cycles from an exchange interaction (Kalogera et al., 2004). More recent work has indicated that a significant number of BH can be retained in GCs (Mackey et al., 2007; Moody & Sigurdsson, 2009). In particular, Maccarone et al. (2011) point out that once a number of BHs are ejected from their host GC the ratio of the masses between the heavy (BHs) and light (non-BH) components falls below the critical value and the Spitzer instability criterion is no longer met, therefore not all BHs will be expunged.
Observationally, the existence of accreting BHs in GCs has been supported by the detection of five luminous (above the NS Eddington limit), variable sources (Maccarone et al., 2007; Brassington et al., 2010; Shih et al., 2010; Irwin et al., 2010; Maccarone et al., 2011). Using the TC/PTC criteria defined in this paper, four of these binary systems are persistent and could be ultra-compact binaries with either a BH (Gnedin et al., 2009) or NS with super-Eddington accretion and rather mild beaming (King, 2011)999Although Peacock et al. (2012) noted that the model presented in King (2011) is inconsistent with the observed emission of bright X-ray sources in GCs, where the observed X-ray emission is harder. They also suggested that many of the sources considered in King (2011) will be dominated by carbon and oxygen edges instead of Thomson scattering, and therefore, as Thomson scattering is assumed in King (2011), this model cannot be applied.. The BH-GC candidates A8 and B1 exhibit large flux variability, with , therefore ruling out genuine disc instability transients. They are also unlikely to be the mildly beamed super-Eddington NS-LMXBs suggested by King (2011).
Instead, A8 and B1 could represent BHs in a GC that have been formed through an exchange interaction, resulting in transient behavior (Kalogera et al., 2004). In Barnard et al. (2011) this formation channel has been advocated to explain a recurring transient in M31, which has been observed with an outburst luminosity and therefore could be a BH-GC candidate (although a NS binary cannot be ruled out). In their work they also suggest that the source could have formed through a complex interaction of a triple system resulting in a BH-WD binary, as first theorised by Ivanova et al. (2010). Further, additional observations of the first strong BH-GC candidate (first presented in Maccarone et al. 2007), have lead to the suggestion that this source is also a triple system, with the inner binary comprising a BH and WD (Maccarone et al., 2010a).
The properties of the host clusters for the five BH-GC candidates are presented in Maccarone et al. (2011), where they demonstrate that this small population suggests that BH binaries are favored in massive ‘red’ GCs. The properties of the GCs hosting A8 and B1 are also massive with =23.4 and =21.6 respectively. However, both of these sources have color values indicating that these GCs are blue, (although B6 has a () of 1.000.02 which is only marginally classified as a ‘blue’ low metallicity cluster: Fabbiano et al. 2010; blue clusters have ). This further indicates that the mass of the GC is important in the formation of GC-LMXBs in these environments. However, we do not comment on the influence of GC colors in the formation of BH-GCs, as the metallicity correlation presented in Maccarone et al. (2011) was only suggestive.
3.3. A8: ULX in Outflow in a GC
The GC source A8 in NGC 3379 is very luminous (Table 4) and variable (Section 2.2). The photometric parameters are consistent with those of the persistent LMXB population (Figure 4), and a single-component PO model provides a statistically acceptable description of the spectrum, with 1.1 (Table 4). However, there is a statistically significant deficiency in the spectrum at 1.0 keV as well as excess emission 0.7 keV (Figure 6). These features cannot be modeled with either a composite DBB plus PO model, or additional neutral absorption components, or a thermal plasma model. A single-component power law model, modified by an ionized intrinsic absorber (absori), was fitted to the spectrum. The results of this fit are presented in Table 5, where the best-fit values from an XSPEC model with tbabstbabsabsoripo (column 2) are shown with one of the neutral absorption components being frozen to Galactic .
A8 is one of the two TC/PTC sources exhibiting short-term variability, where the count-rate falls from 0.01 cnt s to 0.005 cnt s 42 ks into the observation (see Figure 7). The spectrum of the high flux state requires an ionized absorption component. The best-fit values of the absori power law model are presented in column 3 of Table 5, and the spectrum is presented in the left hand panel of Figure 8. These values are similar to those determined for the spectrum from the full observation, except that the best-fit is 1.89 and is larger (although both parameters are consistent with the values derived from the whole observation when considering the uncertainties). The derived X-ray luminosity for the source in this higher flux state is 2.8 , in the Ultraluminous X-ray Source (ULX) range.
The spectrum of the lower state is instead adequately described by a single-component PO model (right hand panel Figure 8), with consistent with Galactic (column 4 of Table 5). However the best-fit photon index is very flat, similar to the value derived when fitting the whole spectrum to a single-component model. This could indicate that ionized absorption may still be present, although there are too few counts to be able to investigate it. The hardness ratios, binned to 5000s over this period, are constant. From this model a luminosity of 8.9 is derived.
This source was also detected in the 3rd observation of NGC 3379 (OBS ID 7074), although it only has 9 net counts, too few to perform any meaningful spectral analysis. B08 estimated an X-ray luminosity of =2.5 . This luminosity then fell to an upper limit of 2 in the subsequent pointing, taken three months later.
The parameters determined from the model during the period of high flux emission are similar to those reported for the NGC 1365 flaring ULX by Soria et al. (2007). In A8, the properties of the spectrum (see Table 5) suggest that the source is in some sort of outflowing phase, similar to the decline period of the ULX in NGC 1365. During this phase the spectrum is fitted with an ionized absorption column of and an ionization parameter of 150 erg cm s, largely consistent with the parameters determined in Soria et al. (2007), of and 111 respectively, where they also determine 1.9 and neutral absorption column . The lower flux spectrum, following the outflow phase reported here, indicates that there has been a three-fold reduction in over the period of 12 hours (assuming a PO model), similar to the factor of 2 seen after three days in NGC 1365.
Soria et al. (2007) suggest that such properties could be analogous to flares observed in Galactic X-ray binaries in a steep-power law state, where the outflow commences after the source exceeds the Eddington limit. Alternatively, some Galactic binaries have been observed to remain in a power-law spectrum HS during short outbursts (Yu & Yan, 2009). It has been suggested that some sources with such HS events could reach flux levels which classify them as ULX during their short outburst events (Yu & Yan, 2009), as we find in the case of A8. This further strengthens the suggestion of Yu & Yan (2009), that some ULXs which stay in the HS during a flaring event harbor stellar-mass compact stars. We note that Galactic BH binaries have been observed to typically proceed through hysteresis cycles (Remillard & McClintock, 2006) meaning that the HS event should be followed by a TD state, which is not observed in the subsequent pointing three months later. However, these canonical state transitions observed for Galactic binaries are often not observed in ULXs (e.g. Feng & Soria 2011 and references therein).
3.4. Super Soft Sources
SSS are sources with little or no emission above 1 keV, with typical X-ray luminosities between 1 to 1 . Following the discovery of SSS in the Large Magellanic Cloud (Long et al., 1981), van den Heuvel et al. (1992) proposed a model of quasi-steady nuclear burning on the surface of a white dwarf (WD) accreting matter from a Roche lobe-filling companion with high accretion rates, with emission described by a simple absorbed blackbody (BB). Further observations have revealed a heterogeneous class of objects with sources detected in both early- and late-type galaxies, in the field and in GCs, and with a range of temporal behavior, and luminosities as high as 1 , in the ULX range (e.g. Carpano et al. 2007; Di Stefano & Kong 2003; Fabbiano et al. 2003). It has been suggested that these sources could be WDs in super-Eddington outbursts, stellar-mass BHs, or even intermediate-mass black holes (e.g. Di Stefano et al. 2004).
A6 was observed in the second and third observations of NGC 3379 (Obs. ID 7073 & 7074). The HR and colors do not change between pointings; consequently the spectra from these two observations were jointly fitted with all parameters apart from the normalization tied between observations. Of the three models, the PO does not provide a physically realistic description of the spectrum, determining a best-fit photon index of , and can be discounted. The best-fit parameters between the DBB and BB model are very similar, with an inner disc temperature of 91 15 eV, and an absorption column consistent with the Galactic value for the BB model, and an of more than twice the Galactic value for the DBB model, although within errors this is consistent with Galactic . These best-fit models yield =1.1 (BB) and =1.5 (DBB), consistent within the uncertainties.
The spectrum is that of a ‘classical’ SSS. Due to its transient nature it is possible that the observed X-ray emission arises from a classical nova (CNe) where episodic thermonuclear explosions on the surface of the WD lasting from months to years can appear as a transient SSS (Pietsch et al., 2007). From X-ray/optical studies of M31 and M33 Pietsch et al. (2005) suggested that these sources are a major class of SSS.
3.4.1 The extremely luminous SSS C1
The second SSS, C1, is even softer than A6, with an inner disc temperature of 6010 eV and a large absorption column (), yielding an intrinsic extremely high X-ray luminosity of 5.8 . The spectrum is shown in Figure 9, along with 1 and 2 contour parameters for the BB model. Even using the 1 lower limits of the inner disc temperature and absorption column of this fit, the luminosity would still be extreme, 5.6 . A neutron star atmosphere model applied to this spectrum also resulted in a large best-fit luminosity .
We further investigated if this source could be a foreground object (e.g. a flaring M dwarf) by fitting a thermal plasma model (APEC) to the spectrum. This resulted in a best-fit temperature 0.3 keV and a flux of (corresponding to at the distance of NGC 4697). By assuming that C1 is a star it is possible to calculate the / ratio by combining the photometric g-band upper limit for this source with the semi-empirical stellar SEDs of Kraus & Hillenbrand (2007) to provide , with values scaled from the best-fit X-ray luminosity derived from the APEC model. The resulting value, /, is much higher than values that are typically seen in stars, where in a recent large survey (Wright et al., 2011) / values were 1 with a dispersion of an order of magnitude. Therefore, C1 is unlikely to be a foreground object unless it was observed during an incredibly rare flaring event for which the X-ray luminosity increased by more than an order of magnitude. To determine the likelihood of observing such a flare we compare to the studies of Feigelson et al. (2004) and Wright et al. (2010), who both studied the properties of stellar X-ray sources in deep, high Galactic latitude fields. From a total of 71 sources and a cumulative exposure of 26.68 Ms they observed only 2 stars with flaring events with peak amplitudes greater than a factor of 10. From their detection thresholds and the relative exposure times we estimate that we should detect 5 stars in the observation of NGC 4697. This equates to a cumulative exposure of 660 ks on these stars, and therefore a probability of observing a significantly large flare in one of these stars of 0.05. This seems highly unlikely and therefore we conclude that the emission from C1 does not arise from a flare star.
C1 was previously discussed in S08 (Section 4.4), who, after determining that a background AGN is unlikely, suggest an intermediate mass black hole with a mass range of , although they note that the formation of such a large black hole outside of the nucleus or a GC is a challenge for current black hole formation theory.
X-ray luminosities in excess of 10 have only been reported for two previous SSS; ULX-1 in M101 (Kong & Di Stefano, 2005) and a luminous SSS transient in NGC 4631 (Carpano et al., 2007). Both of these galaxies are late-type star forming systems, unlike NGC 4697. These sources, when described by a single-component BB model, exhibited similar properties to C1, with soft X-ray spectra (150 eV), high intrinsic absorption, and resulting extremely large . However, in both cases, alternate models have been suggested resulting in less extreme luminosities (Mukai et al. (2005) ; (Liu, 2009); Soria & Ghosh (2009); Carpano et al. (2007). For C1, even though a single-component model is statistically acceptable, there is a possible excess in the residuals around 0.8 keV (see Figure 9). Following Carpano et al. (2007) an additional Gaussian line was included with the BB model in the fit (we also attempted to use an absorption edge to model the spectrum, as presented in Soria & Ghosh (2009) for the NGC 4631 source, but were able to constrain the model parameters). The values obtained from the BB+Gaussian model are presented in Table 6 and the spectrum and the 1 and 2 contours are shown in Figure 10. The best-fit parameters indicate a lower intrinsic absorption, consistent with the Galactic value within errors, and a temperature of 70 eV (compared to 60 eV in the single-component BB model). The peak energy of the Gaussian component, 0.7 keV, is likely to arise from a blend of unresolved emission lines (e.g. O viii; Ness et al. 2005), or there could alternatively be absorption edges, as seen in Soria & Ghosh (2009). The resulting 0.38.0 keV luminosity, 3.4 , although poorly constrained, is three orders of magnitude lower than that from the simple BB model. Although in excess of the emission from a hydrogen burning 1.35 WD system, this luminosity could arise from an extreme super Eddington event (a fireball scenario Soria & Ghosh 2009).
The 2XMM survey reports a source within 2″of C1. Given its soft hardness ratios (e.g. HR2=-0.95101010Derivation of Hardness Ratios are provided in section 3.1.2 in the XMM-Newton Serendipitous Source Catalogue Users Guide), and the lower spatial resolution of XMM-Newton, it is possible that this is emission arising from C1. The XMM-Newton observation was taken in July 2003 (compared to August 2004 for the ‘on’ Chandra observation) yielding an X-ray luminosity of 1.30 , from the 2XMM catalogue (compared to 4.72 for the model derived from the Chandra observation over the 0.212 keV energy range). If this observed emission is from C1, it suggests that the source underwent two outbursts, with a recurrence time of 1 year. From simulations Starrfield et al. (2004) estimated that the hydrogen layer involved in the surface nuclear burning would have a mass of and that steady surface hydrogen burning (for a 1.35 WD) only occurs for accretion rates below 1 . If the outbursts of C1 are a consequence of this fireball scenario this implies accretion rates and therefore the hydrogen shell could be replenished within 1 year. In fact, Soria & Ghosh (2009) estimates that the SSS in NGC 4631 has an accretion rate of around 1 , which could therefore replenish the shell within a month.
However, this interpretation should be treated with caution due to the large uncertainties in the best-fit parameters. Furthermore, the inclusion of the Gaussian component, possibly related to the photospheric expansion period, is speculative, although it serves the purpose of lowering the intrinsic luminosity from the extreme value of to a more typical SSS value of 3.4 .
3.5. Quasi-Soft Sources
QSS are systems with little or no emission above energies of 2 keV, with temperatures between 100 eV to 350 eV (Di Stefano & Kong, 2004). QSS are too hot to be WDs, unless there is significant upscattering of photons emitted by the WD, or the emission emanates from a limited portion of the surface. Alternative scenarios include QSS as NS or BH binary systems (Di Stefano et al., 2010) (and SNRs; Orio 2006). The results of the spectral analysis put B8 and A5 in this category (see Table 4).
Source B8, was observed in one of the cycle 11 Chandra observations and has a maximum outburst time of 3 years. Spectral modeling determines a temperature of between 280360 eV and X-ray luminosity (depending on if a DBB or BB model is preferred). Such values are consistent with the properties of previously observed QSS.
3.5.1 The flaring QSS A5
In the case of the flaring source A5, the best-fit DBB model yields inner disc temperature of 220 eV, with intrinsic absorption of 3, Both the temperature and absorption column values from the BB model are lower than the DBB model, although are consistent within errors, and =2.5 . The outburst data (of 3 ks; Figure 11) yield similar best-fit parameters to those determined from the full observation. However, due to the reduced statistics it was necessary to fix the absorption column; this was frozen to the Galactic value as well as the best-fit column density of 31 and 19 for the DBB and BB models respectively.
Although the average best-fit luminosity of A5 is =2.5 , during the flare (see Table 4), the luminosity reaches 2.9 for the high absorption column, and 9.5 for the Galactic absorption value for the DBB model; the corresponding values for the BB model are 1.7 and 9.4 .
This 3000s flaring event is not typical of QSS (e.g. Di Stefano et al. 2004; Liu 2011). Instead it could arise from a flare star and, to investigate this, a thermal plasma model was also fitted to the spectrum (e.g. Wargelin et al. 2008). This model resulted in a best-fit temperature of 1 keV and a flux of (corresponding to = at the distance of NGC 3379). Using the g-band upper limit derived for A5 (see section 2.1), and following the same method described in section 3.4.1, /. The most luminous stellar flares ever observed reach only /0.1 (Güdel et al., 2004) making this is an unfeasibly large value for a stellar source and making the possibility that this is a flaring star highly unlikely.
Instead this behavior is suggestive of a superburst from a NS. These events were first observed with BeppoSax (Cornelisse et al., 2000) and exhibit outbursts of s, with luminosities and energetics ergs. Because of the large total energy released from these events, these bursts are thought to be powered by carbon burning instead of hydrogen and helium (see Cumming & Bildsten 2001; Cooper & Narayan 2005 and references therein), although it has also been suggested that this emission could arise from runaway electron capture (Kuulkers et al., 2002).
The spectra of superbursts are well described by a BB model, with temperatures 2 keV. For A5, when only the outburst spectrum is fitted, a temperature of eV is determined. However, as we only model the high flux component of the spectrum, the luminosity of the source is = when intrinsic absorption (constrained from the spectrum of the whole observation) is included. When is fixed to the value of Galactic absorption =, with a lower limit of . From these parameters the fluence of this source could be as high as 9.5ergs (with a lower limit of 1.3ergs). This emission is more luminous and energetic than the reported superburst sources, although both of these parameters are consistent within errors. However, the significant difference between A5 and superburst sources arises from the best-fit temperature, where A5 is much cooler. Under the assumption that the source is emitting isotropically a lower limit of R cm is derived. This value indicates that this source is larger than a NS and therefore cannot be a superburst.
Instead, given the implied radius (which has a maximum limit of R cm), this emission could arise from a WD. If such a high fluence event is powered by a thermo-nuclear processes this source could be a helium nova explosion. Such an object was first considered by Kato et al. (1989), where the case of a degenerate WD accreting helium from its helium-rich companion was presented. In this scenario, for accretion rates yr, when sufficient helium has been accumulated along a shell on the surface of the WD a thermo-nuclear reaction, or helium shell flash, occurs resulting in a nova outburst. In 2000 Kato et al. (2000) reported the discovery of such an object; V445 Puppis. This type of outburst event has also been predicted to emit in the X-ray, although possibly not until the post-outburst stage, before which the wind phase may lead to absorption of the soft X-rays (Iben & Tutukov, 1994; Ashok & Banerjee, 2003). However, the period of X-ray emission is expected to be similar to that of classical novae, which lasts for periods months (Pietsch et al., 2007), therefore, an X-ray outburst of 3-ks indicates that this emission does not arise from a Helium Nova.
Another possibility is that this emission arises from a very short period system such as a double WD binary, where these short periods imply that mass transfer is driven by angular momentum loss via gravitational radiation. This would allow for the high X-ray luminosity observed and suggests a mass of 0.40.9 for the accreting white dwarf in A5 (=(49.4)). Assuming that the mass transfer is driven by gravitational radiation gives donor masses slightly smaller than the accretors. Taking slightly smaller WD mass pairs gives a somewhat smaller average luminosity, allowing for some intrinsic variability, although the cause of this is unclear. A possible complication is that the Eddington luminosity may be depressed below the usual value for this type of accretion, since for CO-dominated compositions the Kramers opacity can exceed Thomson at high temperatures (Peacock et al., 2012). Although ‘opacity’ here refers to the Rosseland mean, and is therefore not strictly applicable for considering radiation pressure, a fuller exploration of this type of model would be needed.
In this paper we have identified the transient and potential transient candidates within three elliptical galaxies using deep multi-epoch Chandra data. From spectral analysis of the seventeen sources that have been confirmed as a TC or PTCs (excluding the background source) a large variety of properties have been determined, revealing the heterogeneous nature of transient populations in elliptical galaxies. These sources are summarized in Table 7.
From Table 7 it can be seen that the majority of these sources have been revealed to be normally accreting LMXBs, where 12 of the objects have spectra consistent with binaries in either a hard or thermally dominant state. Of the remaining five sources, one (A8, which has also been determined to be a GC-LMXB) has been observed as a ULX, with a peak luminosity of 2.8 over a 12 hour period, after which the source drops to =8. During this bright period of emission the source has been determined to be in a hard state with evidence of enhanced ionized absorption with a column density of and an ionization parameter of erg cm s. Following this 12 hour period the lower flux emission is consistent with an LMXB in a typically HS. We suggest that this source has undergone a large flaring event prior to our observation and we have observed the subsequent period of decline, where an outflow phase has commenced after the super-Eddington accretion. This behavior is similar to the properties of the ULX observed in NGC 1365 (Soria et al., 2007) and could indicate that A8 is a stellar-mass BH binary undergoing a HS flaring event.
The two SSS transients presented in this paper have properties indicating that they are both likely to arise from WD binaries. The first, A6, has spectra consistent with ‘typical’ SSS, with a temperature of 90 eV and, given the transient nature of this source, is likely to be a classical nova. The second SSS, when fitted with a single-component DBB model has . This value is much greater than any of the previously confirmed X-ray luminosities of ultra-luminous SSS and is difficult to explain. Instead, we include a Gaussian component to the model, which we suggest could be related to an outflow or expansion, as seen in NGC 4361 (Soria & Ghosh, 2009). Our tentative interpretation of this source is a WD binary with photospheric expansion taking place, resulting in an extreme super Eddington event with an X-ray luminosity .
In addition to the SSS, two of the transients presented here have been classified as QSS. One of these sources (B8) has spectral properties that are consistent with previously observed QSS, with a temperature 300 eV and a luminosity . Source A5 on the other hand has rather unusual temporal behavior, where 70% of the emission from this source arises within 3 ks, with a peak luminosity of 9.4 (with a lower limit of 4) and a temperature of eV. This cool emission indicates that the source is a WD and, as a consequence of this short-lived highly luminous event, we tentatively suggest that this emission arises from super-Eddington accretion from a short period double WD binary.
From cross-correlations of HST data (presented in B08, B09, S08 and section 2.1 of this paper) we have been able to quantify the number of field and GC transients and determine the nature of binaries in these objects. Only three sources (A2, A8 & B1) have been determined to be GC-LMXBs. From the modeling of B1 this source could be emitting in either a hard or thermally dominant state but, due to the high luminosity of , we favor the interpretation of thermal emission, providing a best-fit inner disc temperature of 1.3 keV. Due to the high luminosity of not only this object, which is above the Eddington limit for a NS, but also A8, both sources are BH-GC candidates. Such objects are expected to be rare and these provide only the sixth and seventh unambiguous example of such a system. Further, these are only the second and third BH-GC that has been shown to exhibit transient behavior (although a BH-GC candidate emitting has recently been presented in Barnard et al. 2011). Due to the transient nature of these systems we suggest a formation channel of an exchange interaction (Kalogera et al., 2004), instead of a tidally captured source (or an ultra-compact binary; Gnedin et al. 2009; King 2011) which would lead to the formation of a persistent BH-GC.
From the optical matching we have determined that 14 sources arise from the field compared to only three from GCs (the remaining source is a background object). The large number of field sources in the transient population is expected as theoretical work predicts that the majority of the field population will arise from relatively detached systems, which will be transient 75% of their lifetime.
From population synthesis modeling (Fragos et al., 2009) predictions were made about the number of TCs and PTCs that would be detected in our new observations of NGC 4278. However, from comparing the newly determined numbers of TC/PTC sources to those predicted we have much fewer sources, in fact no new TCs were observed and only 2 PTCs, compared to the prediction of 6 TCs and 6 PTCs. We attribute this discrepancy to the original comparison of the modeling results to the observed transient population, where all transient sources included in B08 and B09 were assumed to be NS-LMXBs, and significant difference in our and their statistical analysis and error estimates.
- Angelini et al. (2001) Angelini, L., Loewenstein, M., & Mushotzky, R. F. 2001, ApJ, 557, L35
- Ashok & Banerjee (2003) Ashok, N. M. & Banerjee, D. P. K. 2003, A&A, 409, 1007
- Barnard et al. (2011) Barnard, R., Garcia, M., Li, Z., Primini, F., & Murray, S. S. 2011, ApJ, 734, 79
- Bildsten & Deloye (2004) Bildsten, L. & Deloye, C. J. 2004, ApJL, 607, L119
- Brassington et al. (2010) Brassington, N. J., Fabbiano, G., Blake, S., Zezas, A., Angelini, L., Davies, R. L., Gallagher, J., Kalogera, V., Kim, D.-W., King, A. R., Kundu, A., Trinchieri, G., & Zepf, S. 2010, ApJ, 725, 1805
- Brassington et al. (2008) Brassington, N. J., Fabbiano, G., Kim, D.-W., Zezas, A., Zepf, S., Kundu, A., Angelini, L., Davies, R. L., Gallagher, J., Kalogera, V., Fragos, T., King, A. R., Pellegrini, S., & Trinchieri, G. 2008, ApJS, 179, 142
- Brassington et al. (2009) —. 2009, ApJS, 181, 605
- Carpano et al. (2007) Carpano, S., Pollock, A. M. T., King, A. R., Wilms, J., & Ehle, M. 2007, A&A, 471, L55
- Cash (1979) Cash, W. 1979, ApJ, 228, 939
- Cooper & Narayan (2005) Cooper, R. L. & Narayan, R. 2005, ApJ, 629, 422
- Cornelisse et al. (2000) Cornelisse, R., Heise, J., Kuulkers, E., Verbunt, F., & in’t Zand, J. J. M. 2000, A&A, 357, L21
- Cumming & Bildsten (2001) Cumming, A. & Bildsten, L. 2001, ApJ, 559, L127
- Deutsch et al. (2000) Deutsch, E. W., Margon, B., & Anderson, S. F. 2000, ApJ, 530, L21
- Di Stefano et al. (2010) Di Stefano, R., Kong, A., & Primini, F. A. 2010, NewAR, 54, 72
- Di Stefano & Kong (2003) Di Stefano, R. & Kong, A. K. H. 2003, ApJ, 592, 884
- Di Stefano & Kong (2004) —. 2004, ApJ, 609, 710
- Di Stefano et al. (2004) Di Stefano, R., Kong, A. K. H., Greiner, J., Primini, F. A., Garcia, M. R., Barmby, P., Massey, P., Hodge, P. W., Williams, B. F., Murray, S. S., Curry, S., & Russo, T. A. 2004, ApJ, 610, 247
- Fabbiano (1989) Fabbiano, G. 1989, ARA&A, 27, 87
- Fabbiano (2006) Fabbiano, G. 2006, ARA&A, 44, 323
- Fabbiano et al. (2010) Fabbiano, G., Brassington, N. J., Lentati, L., Angelini, L., Davies, R. L., Gallagher, J., Kalogera, V., Kim, D.-W., King, A. R., Kundu, A., Pellegrini, S., Richings, A. J., Trinchieri, G., Zezas, A., & Zepf, S. 2010, ApJ, 725, 1824
- Fabbiano et al. (2003) Fabbiano, G., King, A. R., Zezas, A., Ponman, T. J., Rots, A., & Schweizer, F. 2003, ApJ, 591, 843
- Feigelson et al. (2004) Feigelson, E. D., Hornschemeier, A. E., Micela, G., Bauer, F. E., Alexander, D. M., Brandt, W. N., Favata, F., Sciortino, S., & Garmire, G. P. 2004, ApJ, 611, 1107
- Feng & Soria (2011) Feng, H. & Soria, R. 2011, NewAR, 55, 166
- Fragos et al. (2008) Fragos, T., Kalogera, V., Belczynski, K., Fabbiano, G., Kim, D.-W., Brassington, N. J., Angelini, L., Davies, R. L., Gallagher, J. S., King, A. R., Pellegrini, S., Trinchieri, G., Zepf, S. E., Kundu, A., & Zezas, A. 2008, ApJ, 683, 346
- Fragos et al. (2009) Fragos, T., Kalogera, V., Willems, B., Belczynski, K., Fabbiano, G., Brassington, N. J., Kim, D.-W., Angelini, L., Davies, R. L., Gallagher, J. S., King, A. R., Pellegrini, S., Trinchieri, G., Zepf, S. E., & Zezas, A. 2009, ApJ, 702, L143
- Giacconi (1974) Giacconi, R. 1974, In X-ray Astronomy, Giacconi, R. & Gursky, H. eds., p. 155, Dordrecht: Reidel
- Gnedin et al. (2009) Gnedin, O. Y., Maccarone, T. J., Psaltis, D., & Zepf, S. E. 2009, ApJ, 705, L168
- Grindlay (1984) Grindlay, J. E. 1984, In Cataclysmic Variables and Low-Mass X-ray Binaries, Bignami, G. and Sunyaev, R. A. eds, p19, Oxford: Pergamon
- Güdel et al. (2004) Güdel, M., Audard, M., Reale, F., Skinner, S. L., & Linsky, J. L. 2004, A&A, 416, 713
- Heinke (2010) Heinke, C. O. 2010, in American Institute of Physics Conference Series, Vol. 1314, American Institute of Physics Conference Series, ed. V. Kologera & M. van der Sluys, 135–142
- Homan et al. (2009) Homan, J., Kaplan, D. L., van den Berg, M., & Young, A. J. 2009, ApJ, 692, 73
- Iben & Tutukov (1994) Iben, Jr., I. & Tutukov, A. V. 1994, ApJ, 431, 264
- Irwin et al. (2010) Irwin, J. A., Brink, T. G., Bregman, J. N., & Roberts, T. P. 2010, ApJ, 712, L1
- Ivanova et al. (2010) Ivanova, N., Chaichenets, S., Fregeau, J., Heinke, C. O., Lombardi, Jr., J. C., & Woods, T. E. 2010, ApJ, 717, 948
- Ivanova et al. (2008) Ivanova, N., Heinke, C. O., Rasio, F. A., Belczynski, K., & Fregeau, J. M. 2008, MNRAS, 386, 553
- Kalogera & Baym (1996) Kalogera, V. & Baym, G. 1996, ApJ, 470, L61+
- Kalogera et al. (2004) Kalogera, V., King, A. R., & Rasio, F. A. 2004, ApJ, 601, L171
- Kato et al. (1989) Kato, M., Saio, H., & Hachisu, I. 1989, ApJ, 340, 509
- Kato et al. (2000) Kato, T., Kanatsu, K., Takamizawa, K., Takao, A., & Stubbings, R. 2000, IAU Circ, 7552, 1
- Kim et al. (2009) Kim, D.-W., Fabbiano, G., Brassington, N. J., Fragos, T., Kalogera, V., Zezas, A., Jordán, A., Sivakoff, G. R., Kundu, A., Zepf, S. E., Angelini, L., Davies, R. L., Gallagher, J. S., Juett, A. M., King, A. R., Pellegrini, S., Sarazin, C. L., & Trinchieri, G. 2009, ApJ, 703, 829
- Kim et al. (2006) Kim, E., Kim, D.-W., Fabbiano, G., Lee, M. G., Park, H. S., Geisler, D., & Dirsch, B. 2006, ApJ, 647, 276
- King (2002) King, A. R. 2002, MNRAS, 335, L13
- King (2011) —. 2011, ApJ, 732, L28+
- Kong & Di Stefano (2005) Kong, A. K. H. & Di Stefano, R. 2005, ApJ, 632, L107
- Kraus & Hillenbrand (2007) Kraus, A. L. & Hillenbrand, L. A. 2007, AJ, 134, 2340
- Kundu et al. (2007) Kundu, A., Maccarone, T. J., & Zepf, S. E. 2007, ApJ, 662, 525
- Kuntschner et al. (2010) Kuntschner, H., Emsellem, E., Bacon, R., Cappellari, M., Davies, R. L., de Zeeuw, P. T., Falcón-Barroso, J., Krajnović, D., McDermid, R. M., Peletier, R. F., Sarzi, M., Shapiro, K. L., van den Bosch, R. C. E., & van de Ven, G. 2010, MNRAS, 408, 97
- Kuulkers et al. (2002) Kuulkers, E., in’t Zand, J. J. M., van Kerkwijk, M. H., Cornelisse, R., Smith, D. A., Heise, J., Bazzano, A., Cocchi, M., Natalucci, L., & Ubertini, P. 2002, A&A, 382, 503
- Lasota et al. (2008) Lasota, J.-P., Dubus, G., & Kruk, K. 2008, A&A, 486, 523
- Liu (2009) Liu, J. 2009, ApJ, 704, 1628
- Liu (2011) —. 2011, ApJS, 192, 10
- Long et al. (1981) Long, K. S., Helfand, D. J., & Grabelsky, D. A. 1981, ApJ, 248, 925
- Maccarone (2003) Maccarone, T. J. 2003, A&A, 409, 697
- Maccarone et al. (2007) Maccarone, T. J., Kundu, A., Zepf, S. E., & Rhode, K. L. 2007, Nature, 445, 183
- Maccarone et al. (2010a) —. 2010a, MNRAS, 409, L84
- Maccarone et al. (2011) —. 2011, MNRAS, 410, 1655
- Maccarone et al. (2010b) Maccarone, T. J., Long, K. S., Knigge, C., Dieball, A., & Zurek, D. R. 2010b, MNRAS, 406, 2087
- Mackey et al. (2007) Mackey, A. D., Wilkinson, M. I., Davies, M. B., & Gilmore, G. F. 2007, MNRAS, 379, L40
- McGowan et al. (2005) McGowan, K. E., Charles, P. A., Blustin, A. J., Livio, M., O’Donoghue, D., & Heathcote, B. 2005, MNRAS, 364, 462
- Moody & Sigurdsson (2009) Moody, K. & Sigurdsson, S. 2009, ApJ, 690, 1370
- Mukai et al. (2005) Mukai, K., Still, M., Corbet, R. H. D., Kuntz, K. D., & Barnard, R. 2005, ApJ, 634, 1085
- Ness et al. (2005) Ness, J.-U., Starrfield, S., Jordan, C., Krautter, J., & Schmitt, J. H. M. M. 2005, MNRAS, 364, 1015
- Orio (2006) Orio, M. 2006, ApJ, 643, 844
- Park et al. (2006) Park, T., Kashyap, V. L., Siemiginowska, A., van Dyk, D. A., Zezas, A., Heinke, C., & Wargelin, B. J. 2006, ApJ, 652, 610
- Peacock et al. (2012) Peacock, M. B., Zepf, S. E., & Maccarone, T. J. 2012, astro-ph/1204.5800
- Pietsch et al. (2005) Pietsch, W., Fliri, J., Freyberg, M. J., Greiner, J., Haberl, F., Riffeser, A., & Sala, G. 2005, A&A, 442, 879
- Pietsch et al. (2007) Pietsch, W., Haberl, F., Sala, G., Stiele, H., Hornoch, K., Riffeser, A., Fliri, J., Bender, R., Bühler, S., Burwitz, V., Greiner, J., & Seitz, S. 2007, A&A, 465, 375
- Piro & Bildsten (2002) Piro, A. L. & Bildsten, L. 2002, ApJL, 571, L103
- Portegies Zwart et al. (2004) Portegies Zwart, S. F., Dewi, J., & Maccarone, T. 2004, MNRAS, 355, 413
- Remillard & McClintock (2006) Remillard, R. A. & McClintock, J. E. 2006, ARA&A, 44, 49
- Rogers et al. (2010) Rogers, B., Ferreras, I., Peletier, R., & Silk, J. 2010, MNRAS, 402, 447
- Sarazin et al. (2000) Sarazin, C. L., Irwin, J. A., & Bregman, J. N. 2000, ApJ, 544, L101
- Scargle (1998) Scargle, J. D. 1998, ApJ, 504, 405
- Schiavon (2007) Schiavon, R. P. 2007, ApJS, 171, 146
- Shahbaz et al. (1997) Shahbaz, T., Naylor, T., & Charles, P. A. 1997, MNRAS, 285, 607
- Shih et al. (2010) Shih, I. C., Kundu, A., Maccarone, T. J., Zepf, S. E., & Joseph, T. D. 2010, ApJ, 721, 323
- Sivakoff et al. (2008a) Sivakoff, G. R., Jordán, A., Juett, A. M., Sarazin, C. L., & Irwin, J. A. 2008a, astro-ph/0806.0626
- Sivakoff et al. (2008b) —. 2008b, astro-ph/0806.0627
- Soria et al. (2007) Soria, R., Baldi, A., Risaliti, G., Fabbiano, G., King, A., La Parola, V., & Zezas, A. 2007, MNRAS, 379, 1313
- Soria & Ghosh (2009) Soria, R. & Ghosh, K. K. 2009, ApJ, 696, 287
- Spitzer (1969) Spitzer, Jr., L. 1969, ApJ, 158, L139
- Starrfield et al. (2004) Starrfield, S., Timmes, F. X., Hix, W. R., Sion, E. M., Sparks, W. M., & Dwyer, S. J. 2004, ApJ, 612, L53
- Terlevich & Forbes (2002) Terlevich, A. I. & Forbes, D. A. 2002, MNRAS, 330, 547
- Tonry et al. (2001) Tonry, J. L., Dressler, A., Blakeslee, J. P., Ajhar, E. A., Fletcher, A. B., Luppino, G. A., Metzger, M. R., & Moore, C. B. 2001, ApJ, 546, 681
- van den Heuvel et al. (1992) van den Heuvel, E. P. J., Bhattacharya, D., Nomoto, K., & Rappaport, S. A. 1992, A&A, 262, 97
- Verbunt & van den Heuvel (1995) Verbunt, F. & van den Heuvel, E. P. J. 1995, in X-ray Binaries, Lewin, W. H. G., van Paradijs, J., van den Heuvel, E. P. J., eds., p. 457, Cambridge, UK: Cambridge Univeristy Press
- Wargelin et al. (2008) Wargelin, B. J., Kashyap, V. L., Drake, J. J., García-Alvarez, D., & Ratzlaff, P. W. 2008, ApJ, 676, 610
- Williams et al. (2008) Williams, B. F., Gaetz, T. J., Haberl, F., Pietsch, W., Shporer, A., Ghavamian, P., Plucinsky, P. P., Mazeh, T., Sasaki, M., & Pannuti, T. G. 2008, ApJ, 680, 1120
- Williams et al. (2006) Williams, B. F., Naik, S., Garcia, M. R., & Callanan, P. J. 2006, ApJ, 643, 356
- Worthey (1994) Worthey, G. 1994, ApJS, 95, 107
- Worthey & Ottaviani (1997) Worthey, G. & Ottaviani, D. L. 1997, ApJS, 111, 377
- Wright et al. (2010) Wright, N. J., Drake, J. J., & Civano, F. 2010, ApJ, 725, 480
- Wright et al. (2011) Wright, N. J., Drake, J. J., Mamajek, E. E., & Henry, G. W. 2011, ApJ, 743, 48
- Wu et al. (2010) Wu, Y. X., Yu, W., Li, T. P., Maccarone, T. J., & Li, X. D. 2010, ApJ, 718, 620
- Yu & Yan (2009) Yu, W. & Yan, Z. 2009, ApJ, 701, 1940
- Zurek et al. (2009) Zurek, D. R., Knigge, C., Maccarone, T. J., Dieball, A., & Long, K. S. 2009, ApJ, 699, 1113
|Galaxy||Age||Obs ID||Obs Date||Expo||No. Src|
Note. – No. src denotes the number of sources detected within each observation in the overlapping region of the S3 chip in all observations. Also note there that in NGC 4278 obs IDs 7081 was taken prior to obs ID 7080, we therefore retain the chronological sequence with obs ID 7081 labeled as pointing 5 in subsequent figures. Age estimates from () Terlevich & Forbes (2002), () Kuntschner et al. (2010), () Rogers et al. (2010).
|Src No.||Catalog No.||RA||Dec||OBS ID||Net Counts||Mode Ratio||Lower-bound value||TC/PTC||Light Curve||Notes||Optical Limits|
|A1||25 (B08)||10:47:54.5||+12:35:30.9||7074||174.9||161.4||43.5||TC||Mid (2)||Field||27.4|
|A2||50 (B08)||10:47:51.6||+12:35:35.9||7073||44.1||6.8||5.9||PTC||LMid (1)||GC||-|
|A3||85 (B08)||10:47:49.2||+12:34:31.1||7074||38.4||30.2||9.0||PTC||Mid (1)||Background||-|
|A4||89 (B08)||10:47:48.9||+12:34:59.1||7076||97.6||18.4||16.4||TC||Beg (1)||Field||25.0|
|A5||94 (B08)||10:47:48.0||+12:34:45.7||1587||58.3||19.2||13.2||TC||End (1)||Field/Double||26.2|
|A6||100 (B08)||10:47:47.2||+12:34:59.8||7073||54.8||20.2||14.1||TC||LMid (2)||Field||26.4|
|A7||115 (B08)||10:47:43.6||+12:34:10.8||1587||18.3||12.9||6.8||PTC||End (1)||Field||27.0|
|A8||128 (B08)||10:47:38.2||+12:33:07.6||7073||548.3||85.1||50.7||TC||LMid (2)||GC||-|
|B1||83 (B09)||12:20:08.138||+29:17:17.1||7079||249.7||57.8||35.5||TC||Mid (1)||GC||-|
|B2||91 (B09)||12:20:07.930||+29:16:57.2||7079||93.5||9.0||8.3||PTC||Mid (1)||Field||26.2|
|B3||139 (B09)||12:20:06.218||+29:16:57.5||7079||170.8||26.1||18.9||TC||Mid (1)||Field||25.8|
|B4||156 (B09)||12:20:05.669||+29:17:15.9||7077||50.5||6.2||5.5||PTC||LMid (1)||Field||26.5|
|B5||159 (B09)||12:20:05.388||+29:17:17.0||7080||23.2||13.6||5.8||PTC||LMid (1)||Field||26.6|
|B6||204 (B09)||12:20:01.582||+29:16:26.2||4741||339.4||871.7||202.9||TC||End (1)||Field||27.0|
|B7||- (new)||12:20:09.001||+29:17:23.7||11269||70.6||11.2||8.4||PTC||Mid (2)||Field||27.0|
|B8||- (new)||12:20:06.361||+29:16:54.6||11269||66.1||11.6||8.7||PTC||Beg (1)||Field||24.9|
|C1||110 (S08)||12:48:40.507||-05:50:20.7||4730||87.4||30.9||22.9||TC||Beg (1)||Field||27.0|
|C2||104 (S08)||12:48:27.484||-05:49:00.8||4730||16.5||8.0||5.4||PTC||Beg (1)||Field||27.8|
Note. – All net count values are derived from the standard processing techniques detailed in both B08 and B09. Optical limits are provided for sources without an identified optical counterpart, where 3 g-band (Vega) magnitudes are quoted. In column 10 sources that could be in outburst for over a year are labeled as ‘LMid’ and those with a shorter outburst upperlimit are ‘Mid’. Sources that were observed in the first observation are labeled as ‘End’ and sources labeled as ‘Beg’ are seen at the beginning of their ‘on state’ in the last epoch of each galaxy’s monitoring campaign. The number in brackets in this column indicates the number of observations the source was detected in.
Note. – The definition of these values are explained in Sections 2.1 & 2.3 of B08.
|Source||Obs. ID.||Net Counts||Model||/||Parameters||Range|
|Parameters||Whole Spectrum||High Flux Spectrum||Low Flux Spectrum|
|(erg cm s)||133||149||-|
|Range (10 )||1.391.93||1.973.14||0.801.03|
Note. – The difference in counts between the ‘whole’ observation and the high and low flux spectra arises from the short time gap between the two extracted spectra. This model is discussed in detail in Section 3.3.
|Range (10 )||0.003.90|
Note. – This spectral model and subsequent interpretation of the source is presented in Section 3.4.1.
|Src No.||Optical Correlation||Spectral State||Suggested Nature|
|A6||Field||SSS||1.12||Classical Nova (CNe)|
|A8||GC||Ionized HS||27.50||GC-ULX in outflow|
|B6||Field||HS (+ disc contribution)||23.22||BH-LMXB|
|C1||Field||SSS+Gauss||3.40||WD binary with|
Note. – Column (3): TD indicates a source in a thermally dominant state, HS a source in a hard state and Canonical a source with insufficient counts to allow spectral modeling to be carried out. Details of the spectral modeling and interpretation of these sources are presented in Sections 2.4 & 3 of this paper.