skip to main content

Title: 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/(Myr−1) over the last 50 Myr and 150–199 HMXBs/(Myr−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.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Page Range / eLocation ID:
Article No. 114
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    We present an analysis of the two-point spatial correlation functions of high-mass X-ray binary (HMXB) and young star cluster (YSC) populations in M31 and M33. We find evidence that HMXBs are spatially correlated with YSCs to a higher degree than would be expected from random chance in both galaxies. When supplemented with similar studies in the Milky Way, Small Magellanic Cloud, and NGC 4449, we find that the peak value of the spatial correlation function correlates strongly with the specific star formation rate of the host galaxy. We additionally perform an X-ray stacking analysis of 211 non-X-ray detected YSCs in M31 and 463 YSCs in M33. We do not detect excess X-ray emission at the stacked cluster locations down to 3σ upper limits of ∼1033 erg s−1 (0.35–8 keV) in both galaxies, which strongly suggests that dynamical formation within YSCs is not a major HMXB formation channel. We interpret our results in the context of (1) the recent star formation histories of the galaxies, which may produce differences in the demographics of compact objects powering the HMXBs, and (2) the differences in natal kicks experienced by compact objects during formation, which can eject newly formed HMXBs from their birth clusters.

    more » « less
  2. Abstract

    We measure the spatially resolved recent star formation history (SFH) of M33 using optical images taken with the Hubble Space Telescope as part of the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) survey. The area covered by the observations used in this analysis covers a de-projected area of ∼38 kpc2and extends to ∼3.5 and ∼2 kpc from the center of M33 along the major and semimajor axes, respectively. We divide the PHATTER optical survey into 2005 regions that measure 24 arcsec, ∼100 pc, on a side and fit color–magnitude diagrams for each region individually to measure the spatially resolved SFH of M33 within the PHATTER footprint. There are significant fluctuations in the SFH on small spatial scales and also galaxy-wide scales that we measure back to about 630 Myr ago. We observe a more flocculent spiral structure in stellar populations younger than about 80 Myr, while the structure of the older stellar populations is dominated by two spiral arms. We also observe a bar in the center of M33, which dominates at ages older than about 80 Myr. Finally, we find that the mean star formation rate (SFR) over the last 100 Myr within the PHATTER footprint is 0.32 ± 0.02 Myr−1. We measure a current SFR (over the last 10 Myr) of 0.20 ± 0.03 Myr−1. This SFR is slightly higher than previous measurements from broadband estimates, when scaled to account for the fraction of the D25 area covered by the PHATTER survey footprint.

    more » « less
  3. Abstract

    We measure the star cluster mass function (CMF) for the Local Group galaxy M33. We use the catalog of stellar clusters selected from the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region survey. We analyze 711 clusters in M33 with7.0<log(Age/yr)<8.5, and log(M/M) > 3.0 as determined from color–magnitude diagram fits to individual stars. The M33 CMF is best described by a Schechter function with power-law slopeα= −2.060.13+0.14, and truncation mass log(Mc/M)=4.240.13+0.16. The data show strong evidence for a high-mass truncation, thus strongly favoring a Schechter function fit over a pure power law. M33's truncation mass is consistent with the previously identified linear trend betweenMc, and star formation rate surface density, ΣSFR. We also explore the effect that individual cluster mass uncertainties have on derived mass function parameters, and find evidence to suggest that large cluster mass uncertainties have the potential to bias the truncation mass of fitted mass functions at the 1σlevel.

    more » « less
  4. Abstract The component black holes (BHs) observed in gravitational-wave (GW) binary black hole (BBH) events tend to be more massive and slower spinning than those observed in black hole X-ray binaries (BH-XRBs). Without modeling their evolutionary histories, we investigate whether these apparent tensions in the BH populations can be explained by GW observational selection effects alone. We find that this is indeed the case for the discrepancy between BH masses in BBHs and the observed high-mass X-ray binaries (HMXBs), when we account for statistical uncertainty from the small sample size of just three HMXBs. On the other hand, the BHs in observed low-mass X-ray binaries (LMXBs) are significantly lighter than the astrophysical BBH population, but this may just be due to a correlation between component masses in a binary system. Given their light stellar companions, we expect light BHs in LMXBs. The observed spins in HMXBs and LMXBs, however, are in tension with the inferred BBH spin distribution at the >99.9% level. We discuss possible scenarios behind the significantly larger spins in observed BH-XRBs. One possibility is that a small subpopulation (conservatively <30%) of BBHs have rapidly spinning primary components, indicating that they may have followed a similar evolutionary pathway to the observed HMXBs. In LMXBs, it has been suggested that BHs can spin up by accretion. If LMXB natal spins follow the BBH spin distribution, we find LMXBs must gain an average dimensionless spin of 0.47 − 0.11 + 0.10 , but if their natal spins follow the observed HMXB spins, the average spin-up must be <0.03. 
    more » « less
  5. Abstract

    Many core-collapse supernovae (SNe) with hydrogen-poor and low-mass ejecta, such as ultra-stripped SNe and type Ibn SNe, are observed to interact with dense circumstellar material (CSM). These events likely arise from the core collapse of helium stars that have been heavily stripped by a binary companion and have ejected significant mass during the last weeks to years of their lives. In helium star models run to days before core collapse we identify a range of helium core masses ≈2.5–3Mwhose envelopes expand substantially due to the helium shell burning while the core undergoes neon and oxygen burning. When modeled in binary systems, the rapid expansion of these helium stars induces extremely high rates of late-stage mass transfer (Ṁ102Myr1) beginning weeks to decades before core collapse. We consider two scenarios for producing CSM in these systems: either mass transfer remains stable and mass loss is driven from the system in the vicinity of the accreting companion, or mass transfer becomes unstable and causes a common envelope event (CEE) through which the helium envelope is unbound. The ensuing CSM properties are consistent with the CSM masses (∼10−2–1M) and radii (∼1013–1016cm) inferred for ultra-stripped SNe and several type Ibn SNe. Furthermore, systems that undergo a CEE could produce short-period neutron star binaries that merge in less than 100 Myr.

    more » « less