Population Of High-mass X-ray Binaries Research Articles

The High-mass X-Ray Binary Luminosity Functions of Dwarf Galaxies

Drawing from the Chandra archive and using a carefully selected set of nearby dwarf galaxies, we present a calibrated high-mass X-ray binary (HMXB) luminosity function in the low-mass galaxy regime and search for an already hinted at dependence on metallicity. Our study introduces a new sample of local dwarf galaxies (D < 12.5 Mpc and M * < 5 × 109 M ⊙), expanding the specific star formation rates (sSFR) and gas-phase metallicities probed in previous investigations. Our analysis of the observed X-ray luminosity function indicates a shallower power-law slope for the dwarf galaxy HMXB population. In our study, we focus on dwarf galaxies that are more representative in terms of sSFR compared to prior work. In this regime, the HMXB luminosity function exhibits significant stochastic sampling at high luminosities. This likely accounts for the pronounced scatter observed in the galaxy-integrated HMXB population’s L X/SFR versus metallicity for our galaxy sample. Our calibration is necessary to understand the active galactic nuclei content of low-mass galaxies identified in current and future X-ray survey fields and has implications for binary population synthesis models, as well as X-ray-driven cosmic heating in the early Universe.

Monitoring the X-Ray Variability of Bright X-Ray Sources in M33

We present a new five-epoch Chandra X-ray Observatory monitoring survey of the nearby spiral galaxy M33 which probes X-ray variability with time sampling between two weeks and four months. We characterize the X-ray variability of 55 bright point sources outside of the nucleus, many of which are expected to be high-mass X-ray binaries (HMXBs). We detect eight new candidate transients not detected in previous X-ray catalogs of M33 and discuss their possible nature. The final catalog includes 26 known HMXB candidates identified in the literature. We extend the baseline of the X-ray light curves up to 21 yr by including archival X-ray observations of these sources. We compare the detection and nondetection epochs of the sources to suites of simulated source duty cycles and infer that most of our detected sources have duty cycles >30%. We find only four sources whose detection patterns are consistent with having duty cycles below 30%. This large fraction of sources with high duty cycles is unexpected for a population of HMXBs; thus more frequent X-ray monitoring will likely reveal many more low duty cycle HMXBs in M33.

Multiwavelength Characterization of the High-mass X-Ray Binary Population of M33

We present multiwavelength characterization of 65 high-mass X-ray binary (HMXB) candidates in M33. We use the Chandra ACIS survey of M33 (ChASeM33) catalog to select hard X-ray point sources that are spatially coincident with UV-bright point-source optical counterparts in the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region catalog, which covers the inner disk of M33 at near-IR, optical, and near-UV wavelengths. We perform spectral energy distribution fitting on multiband photometry for each point-source optical counterpart to measure its physical properties including mass, temperature, luminosity, and radius. We find that the majority of the HMXB companion star candidates are likely B-type main-sequence stars, suggesting that the HMXB population of M33 is dominated by Be X-ray binaries (Be-XRBs), as is seen in other Local Group galaxies. We use spatially resolved recent star formation history maps of M33 to measure the age distribution of the HMXB candidate sample and the HMXB production rate for M33. We find a bimodal distribution for the HMXB production rate over the last 80 Myr, with a peak at ∼10 and ∼40 Myr, which match theoretical formation timescales for the most massive HMXBs and Be-XRBs, respectively. We measure an HMXB production rate of 107–136 HMXBs/(M ⊙ yr−1) over the last 50 Myr and 150–199 HMXBs/(M ⊙ yr−1) over the last 80 Myr. For sources with compact object classifications from overlapping NuSTAR observations, we find a preference for giant/supergiant companion stars in black hole HMXBs and main-sequence companion stars in neutron star HMXBs.

The Large Deficit of HMXB Emission from Luminous Infrared Galaxies: The Case of the Circumnuclear Starburst Ring in NGC 7552

Luminous infrared galaxies (LIRGs), the most extreme star-forming galaxies in the nearby (D < 30 Mpc) universe, show a notable X-ray emission deficiency (up to a factor of ∼10) compared with predictions from scaling relations of galaxy-wide high-mass X-ray binary (HMXB) luminosity with star formation rate. In the nearby (≈20 Mpc) LIRG NGC 7552, the majority of the IR emission originates in a circumnuclear starburst ring, which has been resolved into several discrete knots of star formation. We present results from recent Chandra observations of NGC 7552, which reveal significant deficits in the 2–7 keV X-ray luminosities from two of the most powerful star-forming knots. We hypothesize that the expected luminous HMXB populations in these knots are either (1) obscured by very large column densities or (2) suppressed due to the knots having relatively high metallicity and/or very young ages (≲5 Myr). We distinguish between these possibilities using data from recent NuSTAR observations, whose sensitivity above 10 keV is capable of uncovering heavily obscured HMXB populations, since emission at these energies is more immune to absorption effects. We find no evidence of a heavily obscured HMXB population in the central region of NGC 7552, suggesting suppressed HMXB formation. We further show that metallicity-dependent scaling relations cannot fully account for the observed deficit from the most powerful star-forming knots or the central region as a whole. Thus, we suggest that recent bursts in local star formation activity likely drive the high L IR within these regions on timescales ≲5 Myr, shorter than the timescale required for the formation of HMXBs.

The Missing Link between Black Holes in High-mass X-Ray Binaries and Gravitational-wave Sources: Observational Selection Effects

There are few observed high-mass X-ray binaries (HMXBs) that harbor massive black holes (BHs), and none are likely to result in a binary black hole (BBH) that merges within a Hubble time; however, we know that massive merging BBHs exist from gravitational-wave (GW) observations. We investigate the role that X-ray and GW observational selection effects play in determining the properties of their respective detected binary populations. We find that, as a result of selection effects, detectable HMXBs and detectable BBHs form at different redshifts and metallicities, with detectable HMXBs forming at much lower redshifts and higher metallicities than detectable BBHs. We also find disparities in the mass distributions of these populations, with detectable merging BBH progenitors pulling to higher component masses relative to the full detectable HMXB population. Fewer than 3% of detectable HMXBs host BHs >35M ⊙ in our simulated populations. Furthermore, we find the probability that a detectable HMXB will merge as a BBH system within a Hubble time is ≃0.6%. Thus, it is unsurprising that no currently observed HMXB is predicted to form a merging BBH with high probability.

Broadband study and the discovery of pulsations from the Be/X-ray binary eRASSU J052914.9−662446 in the Large Magellanic Cloud

Context. The Magellanic Clouds are our nearest star-forming galaxies. While the population of high-mass X-ray binaries (HMXBs) in the Small Magellanic Cloud is relatively well studied, our knowledge about the Large Magellanic Cloud (LMC) is far from complete given its large angular extent and the insufficient coverage with X-ray observations. Aims. We conducted a search for new HMXBs in the LMC using data from eROSITA, the soft X-ray instrument on board the Spektrum-Roentgen-Gamma satellite. Methods. After confirming the nature of eRASSU J052914.9−662446 as a hard X-ray source that is positionally coincident with an early-type star, we followed it up with optical spectroscopic observations from the South African Large Telescope (SALT) and a dedicated NuSTAR observation. Results. We study the broadband timing and spectral behaviour of the newly discovered HMXB eRASSU J052914.9−662446 through eROSITA, Swift, and NuSTAR data in X-rays and the Optical Gravitational Lensing Experiment (OGLE) and SALT RSS data at the optical wavelength. We report the detection of a spin period at 1412 s and suggest that the orbital period of the system is ∼151 days. We thereby establish that eRASSU J052914.9−662446 is an accreting pulsar. Furthermore, through optical spectroscopic observations and the detection of Hα emission, the source is identified as a Be X-ray binary pulsar in the LMC. We also investigated the variability of the source in the optical and X-ray regime over the past decades and provide estimates of the possible magnetic field strength of the neutron star.

A Chandra Virgo cluster survey of spiral galaxies – I. Introduction to the survey and a new ULX sample

ABSTRACT We present an analysis of the ultraluminous X-ray source (ULX) population in 75 Virgo cluster late-type galaxies, including all those with a star formation rate ≳1 M⊙ yr−1 and a representative sample of the less star forming ones. This study is based on 110 observations obtained over 20 yr with the Chandra X-ray Observatory Advanced Camera for Imaging Spectroscopy. As part of a Large Chandra Program, new observations were obtained for 52 of these 75 galaxies. The data are complete to a sensitivity of ≈1039 erg s−1, with a typical detection limit of ≈3 × 1038 erg s−1 for the majority of the sources. The catalogue contains about 80 ULXs (0.3–10 keV luminosity &amp;gt;1039 erg s−1), and provides their location, observed flux, de-absorbed luminosity, and (for the 25 most luminous ones) simple X-ray spectral properties. We discuss the ULX luminosity function in relation to the mass and star formation rate of the sample galaxies. We show that recent models of low-mass plus high-mass X-ray binary populations (scaling with stellar mass and star formation rate, respectively) are mostly consistent with our observational results. We tentatively identify the most luminous X-ray source in the sample (a source in IC 3322A with LX ≈ 6 × 1040 erg s−1) as a recent supernova or its young remnant. The properties of the sample galaxies (morphologies, stellar masses, star formation rates, total X-ray luminosities from their point-source population) are also summarized.

Case Study on the γ-algorithm: Possible Tests from the Luminosity versus Displacement Correlation of High Mass X-Ray Binaries

Using the apparent correlation of luminosity (L X) versus displacement (R) of high mass X-ray binaries (HMXBs), we aim to constrain the common envelope (CE) mechanism, which is vital in the formation and evolution of compact binaries. We find that under the assumption of the γ-algorithm, the apparent correlation can also be reconstructed generally within a reasonable range of key parameters adopted, though the population of HMXBs is distinct with that in the canonical α CE-formalism. We compare the spatial distribution of HMXBs under the two CE mechanisms, and suggest the difference in L X versus R distribution may provide an additional clue for the study of the CE phase and to discriminate between CE models.

Evidence for Low Kick Velocities among High-mass X-Ray Binaries in the Small Magellanic Cloud from the Spatial Correlation Function

Abstract We present the two-point cross-correlation function between high-mass X-ray binaries (HMXBs) in the Small Magellanic Cloud (SMC) and their likely birthplaces (OB associations: OBAs). This function compares the spatial correlation between the observed HMXB and OBA populations against mock catalogs in which the members are distributed randomly across the sky. A significant correlation (∼15σ) is found for the HMXB and OBA populations when compared with a randomized catalog in which the OBAs are distributed uniformly over the SMC. A less significant correlation (4σ) is found for a randomized catalog of OBAs built with a bootstrap method. However, no significant correlation is detected when the randomized catalogs assume the form of a Gaussian ellipsoid or a distribution that reflects the star formation history from 40 Myr ago. Based on their observed distributions and assuming a range of migration timescales, we infer that the average value of the kick velocity inherited by an HMXB during the formation of its compact object is 2–34 km s−1. This is considerably less than the value obtained for their counterparts in the Milky Way hinting that the galactic environment affecting stellar evolution plays a role in setting the average kick velocity of HMXBs.

Multiwavelength Characterization of the High-mass X-Ray Binary Population of M31

Abstract We present our analysis of high-quality high-mass X-ray binary (HMXB) candidates in M31 selected from point-source optical counterpart candidates from the Chandra-PHAT survey catalog. We fit the spectral energy distributions of optical counterpart candidates using the Bayesian Extinction and Stellar Tool. We used the best-fit luminosity, effective temperature, radius, and dust reddening for the companion stars in combination with the local star formation history, dust maps of M31, published X-ray spectral fits from XMM-Newton observations, IR colors, and Chandra X-ray hardness ratios to determine our best sample of HMXB candidates. The age distribution of the HMXB sample appears peaked between 10 and 50 Myr, consistent with findings in other nearby galaxies. Using the age distribution and mean star formation rate (SFR), we find that 80–136 HMXBs were produced per unit of SFR over the last 50 Myr and 89–163 HMXBs were produced per unit of SFR over the last 80 Myr, if we expand the assumed age limit beyond the lifetimes of single massive stars. We also calculate the HMXB production rate (HMXBs/M ⊙) over time, which ranges from 7 × 10−7 to 4 × 10−6 HMXBs/M ⊙ over the last 80 Myr, in agreement with both theoretical predictions and measured production rates in other galaxies.

Metal-THINGS: On the Metallicity and Ionization of ULX Sources in NGC 925

Abstract We present an analysis of the optical properties of three Ultra Luminous X-ray (ULX) sources identified in NGC 925. We use Integral field unit data from the George Mitchel spectrograph in the context of the Metal-THINGS survey. The optical properties for ULX-1 and ULX-3 are presented, while the spaxel associated with ULX-2 had a low S/N, which prevented its analysis. We also report the kinematics and dimensions of the optical nebula associated with each ULX using ancillary data from the PUMA Fabry–Perot spectrograph. A BPT analysis demonstrates that most spaxels in NGC 925 are dominated by star-forming regions, including those associated with ULX-1 and ULX-3. Using the resolved gas-phase metallicities, a negative metallicity gradient is found, consistent with previous results for spiral galaxies, while the ionization parameter tends to increase radially throughout the galaxy. Interestingly, ULX-1 shows a very low gas metallicity for its galactocentric distance, identified by two independent methods, while exhibiting a typical ionization. We find that such low gas metallicity is best explained in the context of the high-mass X-ray binary population, where the low-metallicity environment favors active Roche lobe overflows that can drive much higher accretion rates. An alternative scenario invoking accretion of a low-mass galaxy is not supported by the data in this region. Finally, ULX-3 shows both a high metallicity and ionization parameter, which is consistent with the progenitor being a highly accreting neutron star within an evolved stellar population region.

The M Supergiant High-mass X-Ray Binary 4U 1954+31

Abstract The X-ray binary 4U 1954+31 has been classified as a low-mass X-ray binary containing an M giant and a neutron star (NS). It has also been included in the rare class of X-ray symbiotic binaries. The Gaia parallax, infrared colors, spectral type, abundances, and orbital properties of the M star demonstrate that the cool star in this system is not a low-mass giant but a high-mass M supergiant. Thus, 4U 1954+31 is a high-mass X-ray binary (HMXB) containing a late-type supergiant. It is the only known binary system of this type. The mass of the M I is M ⊙ giving an age of this system in the range 12–50 Myr with the NS no more than 43 Myr old. The spin period of the NS is one of the longest known, 5 hr. The existence of M I plus NS binary systems is in accord with stellar evolution theory, with this system a more evolved member of the HMXB population.

Evidence from high-mass X-ray binaries that Galactic WR components of WR+O binaries end their life with a supernova explosion

Context. Theoretical population number studies of binaries with at least one black hole (BH) component obviously depend on whether or not BHs receive a (natal) kick during their formation. Aims. Several observational facts seem to indicate that BHs do indeed receive a kick during their formation. In the present paper, we discuss additional evidence of this. Methods. The progenitors of wind-fed high-mass X-ray binaries (HMXB) with a BH component (BH HMXB) are WR+OB binaries where the Wolf–Rayet (WR) star will finally collapse and form the BH. Starting from the observed population of WR+OB binaries in the solar neighborhood, we predict the population of wind-fed BH HMXBs as a function of the BH-natal kick. Results. The simulations reveal that when WR stars collapse into a BH with a zero or low kick, we should expect 100 or more wind-fed BH HMXBs in the solar neighborhood, whereas only one is observed (Cyg X-1). We consider this as evidence that either WR components in binaries end their life as a neutron star or that they collapse into BHs, both accompanied by a supernova explosion imparting significant (natal) kicks.

High-mass X-ray binaries in nearby metal-poor galaxies: on the contribution to nebular He ii emission

ABSTRACT Despite significant progress both observationally and theoretically, the origin of high-ionization nebular He ii emission in galaxies dominated by stellar photoionization remains unclear. Accretion-powered radiation from high-mass X-ray binaries (HMXBs) is still one of the leading proposed explanations for the missing He+-ionizing photons, but this scenario has yet to be conclusively tested. In this paper, we present nebular line predictions from a grid of photoionization models with input spectral energy distributions containing the joint contribution of both stellar atmospheres and a multicolour disc model for HMXBs. This grid demonstrates that HMXBs are inefficient producers of the photons necessary to power He ii, and can only boost this line substantially in galaxies with HMXB populations large enough to power X-ray luminosities of 1042 erg s−1 per unit star formation rate (SFR). To test this, we assemble a sample of 11 low-redshift star-forming galaxies with high-quality constraints on both X-ray emission from Chandra and He ii emission from deep optical spectra, including new observations with the MMT. These data reveal that the HMXB populations of these nearby systems are insufficient to account for the observed He ii strengths, with typical X-ray luminosities or upper limits thereon of only 1040–1041 erg s−1 per SFR. This indicates that HMXBs are not the dominant source of He+ ionization in these metal-poor star-forming galaxies. We suggest that the solution may instead reside in revisions to stellar wind predictions, softer X-ray sources, or very hot products of binary evolution at low metallicity.

NuSTAR Observations of the Unidentified INTEGRAL Sources: Constraints on the Galactic Population of HMXBs

Abstract The NuSTAR Legacy program titled Unidentified INTEGRAL Sources targeted faint hard X-ray sources revealed by INTEGRAL in the Galactic plane in order to provide conclusive identification of their nature and insights on the population of faint hard X-ray sources. The NuSTAR and Swift X-Ray Telescope observations obtained in 2015–2017 contributed to the successful identification of five persistent sources. Here, we report on the spectral and variability analyses that helped to consolidate the classifications of IGR J10447–6027, IGR J16181–5407, and IGR J20569+4940 as active galactic nuclei and IGR J17402–3656 as an intermediate polar. An optical spectrum of the blazar IGR J20569+4940 is also presented. Combining these results with successful identifications of other such faint and persistent INTEGRAL sources reported in the literature, we investigate possible implications for the population of persistent high-mass X-ray binaries (HMXBs) below the identification completion limit of the INTEGRAL survey. The current trend hints at a deficit of persistent HMXBs below F 17–60 keV = 10−11 erg cm−2 s−1, but additional efforts dedicated to classifying faint hard X-ray sources are needed before we can draw solid conclusions.

Young Accreting Compact Objects in M31: The Combined Power of NuSTAR, Chandra, and Hubble.

We present 15 high-mass X-ray binary (HMXB) candidates in the disk of M31 for which we are able to infer compact object type, spectral type of the donor star, and age using multiwavelength observations from NuSTAR, Chandra, and the Hubble Space Telescope. The hard X-ray colors and luminosities from NuSTAR permit the tentative classification of accreting X-ray binary systems by compact object type, distinguishing black hole from neutron star systems. We find hard-state black holes, pulsars, and non-magnetized neutron stars associated with optical point-source counterparts with similar frequency. We also find nine non-magnetized neutron stars coincident with globular clusters and an equal number of pulsars with and without point-source optical counterparts. We perform spectral energy distribution (SED) fitting for the most likely optical counterparts to the HMXB candidates, finding seven likely high-mass stars and one possible red helium-burning star. The remaining seven HMXB optical counterparts have poor SED fits, so their companion stars remain unclassified. Using published star formation histories, we find that the majority of HMXB candidates-X-ray sources with UV-bright point-source optical counterpart candidates-are found in regions with star formation bursts less than 50 Myr ago, and three are associated with young stellar ages (<10Myr). This is consistent with similar studies of HMXB populations in the Magellanic Clouds, M33, NGC 300, and NGC 2403.

Dart_board: Binary Population Synthesis with Markov Chain Monte Carlo

Abstract By employing Monte Carlo random sampling, traditional binary population synthesis (BPS) offers a substantial improvement in efficiency over brute force, grid-based studies. Even so, BPS models typically require a large number of simulation realizations, a computationally expensive endeavor, to generate statistically robust results. Recent advances in statistical methods have led us to revisit the traditional approach to BPS. In this work we describe our publicly available code dart_board, which combines rapid binary evolution codes, typically used in traditional BPS, with modern Markov chain Monte Carlo methods. dart_board takes a novel approach that treats the initial binary parameters and the supernova kick vector as model parameters. This formulation has several advantages, including the ability to model either populations of systems or individual binaries, the natural inclusion of observational uncertainties, and the flexible addition of new constraints, which are problematic to include using traditional BPS. After testing our code with mock systems, we demonstrate the flexibility of dart_board by applying it to three examples: (i) a generic population of high-mass X-ray binaries (HMXBs); (ii) the population of HMXBs in the Large Magellanic Cloud (LMC), in which the spatially resolved star formation history is used as a prior; and (iii) one particular HMXB in the LMC, Swift J0513.4−6547, in which we include observations of the system’s component masses and orbital period. Although this work focuses on HMXBs, dart_board can be applied to a variety of stellar binaries, including the recent detections by gravitational wave observatories of merging compact object binaries.

Formation time-scales for high-mass X-ray binaries in M33

We have identified 55 candidate high-mass X-ray binaries (HMXBs) in M33 using available archival {\it HST} and {\it Chandra} imaging to find blue stars associated with X-ray positions. We use the {\it HST} photometric data to model the color-magnitude diagrams in the vicinity of each candidate HMXB to measure a resolved recent star formation history (SFH), and thus a formation timescale, or age for the source. Taken together, the SFHs for all candidate HMXBs in M33 yield an age distribution that suggests preferred formation timescales for HMXBs in M33 of $<$ 5 Myr and $\sim$ 40 Myr after the initial star formation episode. The population at 40 Myr is seen in other Local Group galaxies, and can be attributed to a peak in formation efficiency of HMXBs with neutron stars as compact objects and B star secondary companions. This timescale is preferred as neutron stars should form in abundance from $\sim$ 8 M$_{\odot}$ core-collapse progenitors on these timescales, and B stars are shown observationally to be most actively losing mass around this time. The young population at $<$ 5 Myr has not be observed in other Local Group HMXB population studies, but may be attributed to a population of very massive progenitors forming black holes very early on. We discuss these results in the context of massive binary evolution, and the implications for compact object binaries and gravitational wave sources.

Impact of ultraluminous X-ray sources on photoabsorption in the first galaxies

In the local Universe, integrated X-ray emission from high-mass X-ray binaries (HMXBs) is dominated by the brightest ultraluminous X-ray sources (ULXs) with luminosity >10^40 erg/s. Such rare objects probably also dominated the production of X-rays in the early Universe. We demonstrate that a ULX with Lx~10^40-10^41 erg/s (isotropic-equivalent luminosity in the 0.1-10 keV energy band) shining for ~10^5 years (the expected duration of a supercritically accreting phase in HMXBs) can significantly ionise the ISM in its host dwarf galaxy of total mass M~10^7-10^8 Msun and thereby reduce its opacity to soft X-rays. As a result, the fraction of the soft X-ray (below 1 keV) radiation from the ULX escaping into the intergalactic medium (IGM) can increase from ~20-50% to ~30-80% over its lifetime. This implies that HMXBs can induce a stronger heating of the IGM at z>10 compared to estimates neglecting the ULX feedback on the ISM. However, larger galaxies with M>3 x 10^8 Msun could not be significantly ionised even by the brightest ULXs in the early Universe. Since such galaxies probably started to dominate the global star-formation rate at z<10, the overall escape fraction of soft X-rays from the HMXB population probably remained low, <30%, at these epochs.

The Orbit of the Gamma-Ray Binary 1FGL J1018.6−5856

Abstract Gamma-ray binaries are a small subclass of the high mass X-ray binary population that exhibit emission across the whole electromagnetic spectrum. We present the radial velocities of 1FGL J1018.6−5856 based on the observations obtained with the Southern African Large Telescope. We combine our measurements with those published in the literature to get a broad phase coverage. The mass function obtained supports a neutron star compact object, although a black hole mass is possible for the very low inclination angles. The improved phase coverage allows constraints to be placed on the orbital eccentricity (e = 0.31 ± 0.16), which agrees with the estimates from the high-energy data.