Identifying X-ray binary (XRB) candidates in nearby galaxies requires distinguishing them from possible contaminants including foreground stars and background active galactic nuclei. This work investigates the use of supervised machine learning algorithms to identify high-probability X-ray binary candidates. Using a catalogue of 943 Chandra X-ray sources in the Andromeda galaxy, we trained and tested several classification algorithms using the X-ray properties of 163 sources with previously known types. Amongst the algorithms tested, we find that random forest classifiers give the best performance and work better in a binary classification (XRB/non-XRB) context compared to the use of multiple classes. Evaluating our method by comparing with classifications from visible-light and hard X-ray observations as part of the Panchromatic Hubble Andromeda Treasury, we find compatibility at the 90% level, although we caution that the number of sources in common is rather small. The estimated probability that an object is an X-ray binary agrees well between the random forest binary and multiclass approaches and we find that the classifications with the highest confidence are in the X-ray binary class. The most discriminating X-ray bands for classification are the 1.7-2.8, 0.5-1.0, 2.0-4.0, and 2.0-7.0 keV photon flux ratios. Of the 780 unclassified sources in the Andromeda catalogue, we identify 16 new high-probability X-ray binary candidates and tabulate their properties for follow-up.

We present a source catalog from the first deep hard X-ray ($E>10$ keV) survey of the Small Magellanic Cloud (SMC), the NuSTAR Legacy Survey of the SMC. We observed three fields, for a total exposure time of 1 Ms, along the bar of this nearby star-forming galaxy. Fields were chosen for their young stellar and accreting binary populations. We detected 10 sources above a 3$\sigma$ significance level (4$-$25 keV) and obtained upper limits on an additional 40 sources. We reached a 3$\sigma$ limiting luminosity in the 4$-$25 keV band of $\sim$ $10^{35}$ erg s$^{-1}$, allowing us to probe fainter X-ray binary (XRB) populations than has been possible with other extragalactic NuSTAR surveys. We used hard X-ray colors and luminosities to constrain the compact-object type, exploiting the spectral differences between accreting black holes and neutron stars at $E>10$ keV. Several of our sources demonstrate variability consistent with previously observed behavior. We confirmed pulsations for seven pulsars in our 3$\sigma$ sample. We present the first detection of pulsations from a Be-XRB, SXP305 (CXO J005215.4$-$73191), with an X-ray pulse period of $305.69\pm0.16$ seconds and a likely orbital period of $\sim$1160-1180 days. Bright sources ($\gtrsim 5\times 10^{36}$ erg s$^{-1}$) in our sample have compact-object classifications consistent with their previously reported types in the literature. Lower luminosity sources ($\lesssim 5\times 10^{36}$ erg s$^{-1}$) have X-ray colors and luminosities consistent with multiple classifications. We raise questions about possible spectral differences at low luminosity between SMC pulsars and the Galactic pulsars used to create the diagnostic diagrams.

What are the most important conditions and processes governing the growth of stellar-origin compact objects? The identification of compact object type as either black hole (BH) or neutron star (NS) is fundamental to understanding their formation and evolution. To date, time-domain determination of compact object type remains a relatively untapped tool. Measurement of orbital periods, pulsations, and bursts will lead to a revolution in the study of the demographics of NS and BH populations, linking source phenomena to accretion and galaxy parameters (e.g., star formation, metallicity). To perform these measurements over sufficient parameter space, a combination of a wide-field (>5000 deg^2) transient X-ray monitor over a dynamic energy range (~1-100 keV) and an X-ray telescope for deep surveys with <5 arcsec PSF half-energy width (HEW) angular resolution are required. Synergy with multiwavelength data for characterizing the underlying stellar population will transform our understanding of the time domain properties of transient sources, helping to explain details of supernova explosions and gravitational wave event rates.

We present nearly simultaneous NuSTAR and XMM-Newton observations of the nearby (832 kpc) ultraluminous X-ray source (ULX) M33 X-8. M33 X-8 has a 0.3-10 keV luminosity of LX ~ 1.4 x 10^39 erg/s, near the boundary of the "ultraluminous" classification, making it an important source for understanding the link between typical Galactic X-ray binaries and ULXs. Past studies have shown that the 0.3-10 keV spectrum of X-8 can be characterized using an advection-dominated accretion disk model. We find that when fitting to our NuSTAR and XMM-Newton observations, an additional high-energy (>10 keV) Comptonization component is required, which allows us to rule out single advection-dominated disk and classical sub-Eddington models. With our new constraints, we analyze XMM-Newton data taken over the last 17 years to show that small (~30%) variations in the 0.3-10 keV flux of M33 X-8 result in spectral changes similar to those observed for other ULXs. The two most likely phenomenological scenarios suggested by the data are degenerate in terms of constraining the nature of the accreting compact object (i.e., black hole versus neutron star). We further present a search for pulsations using our suite of data; however, no clear pulsations are detected. Future observations designed to observe M33 X-8 at different flux levels across the full 0.3-30 keV range would significantly improve our constraints on the nature of this important source.

The X-ray source populations within galaxies are typically difficult to identify and classify from X-ray data alone. We are able to break through this barrier by combining deep new Chandra ACIS-I observations with extensive Hubble Space Telescope imaging from the PHAT survey of the M31 disk. We detect 373 X-ray sources down to 0.35-8.0 keV flux of 10$^{-15}$ erg cm$^{-2}$ s$^{-1}$ over 0.4 square degrees, 170 of which are reported for the first time. We identify optical counterpart candidates for 188 of the 373 sources, after using the HST data to correct the absolute astrometry of our Chandra imaging to 0.1$"$. While 58 of these 188 are associated with point sources potentially in M31, over half (107) of the counterpart candidates are extended background galaxies, 5 are star clusters, 12 are foreground stars, and 6 are supernova remnants. Sources with no clear counterpart candidate are most likely to be undetected background galaxies and low-mass X-ray binaries in M31. The 58 point sources that are not consistent with foreground stars are bright enough that they could be high mass stars in M31; however, all but 8 have optical colors inconsistent with single stars, suggesting that many could be background galaxies or binary counterparts. For point-like counterparts, we examine the star formation history of the surrounding stellar populations to look for a young component that could be associated with a high mass X-ray binary. For the 40 point-like counterpart candidates associated with young populations, we find that their age distribution has two peaks at 15-20 Myr and 40-50 Myr. If we only consider the 8 counterpart candidates with typical high-mass main sequence optical star colors, their age distribution peaks mimic those of the sample of 40. Finally, we find that intrinsic faintness, and not extinction, is the main limitation for finding further counterpart candidates.

Nearby galaxy surveys have long classified X-ray binaries (XRBs) by the mass category of their donor stars (high-mass and low-mass). The NuSTAR observatory, which provides imaging data at E $>10$ keV, has enabled the classification of extragalactic XRBs by their compact object type: neutron star (NS) or black hole (BH). We analyzed NuSTAR/Chandra/XMM-Newton observations from a NuSTAR-selected sample of 12 galaxies within 5 Mpc having stellar masses ($M_{\star}$) $10^{7-11}$ $M_{\odot}$ and star formation rates (SFR) $\approx0.01-15$ $M_{\odot}$ yr$^{-1}$. We detect 128 NuSTAR sources to a sensitivity of $\approx10^{38}$ erg s$^{-1}$. Using NuSTAR color-intensity and color-color diagrams we classify 43 of these sources as candidate NS and 47 as candidate BH. We further subdivide BH by accretion states (soft, intermediate, and hard) and NS by weak (Z/Atoll) and strong (accreting pulsar) magnetic field. Using 8 normal (Milky Way-type) galaxies in the sample, we confirm the relation between SFR and galaxy X-ray point source luminosity in the 4-25 and 12-25 keV energy bands. We also constrain galaxy X-ray point source luminosity using the relation $L_{\rm{X}}=\alpha M_{\star}+\beta\text{SFR}$, finding agreement with previous work. The XLF of all sources in the 4-25 and 12-25 keV energy bands matches with the $\alpha=1.6$ slope for high-mass XRBs. We find that NS XLFs suggest a decline beginning at the Eddington limit for a 1.4 $M_{\odot}$ NS, whereas the BH fraction shows an approximate monotonic increase in the 4-25 and 12-25keV energy bands. We calculate the overall ratio of BH to NS to be $\approx1$ for 4-25 keV and $\approx2$ for 12-25 keV.

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 (HST). 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 7 likely high mass stars and 1 possible red Helium burning star. The remaining 7 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, with 3 associated with young stellar ages (<10 Myr). This is consistent with similar studies of HMXB populations in the Magellanic Clouds, M33, NGC 300, and NGC 2403.

Galactic and extragalactic studies have shown that metal-rich globular clusters (GCs) are approximately three times more likely to host bright low-mass X-ray binaries (LMXBs) than metal-poor GCs. There is no satisfactory explanation for this metallicity effect. We tested the hypothesis that the number density of red giant branch (RGB) stars is larger in metal-rich GCs, and thus potentially the cause of the metallicity effect. Using Hubble Space Telescope photometry for 109 unique Milky Way GCs, we investigated whether RGB star density was correlated with GC metallicity. Isochrone fitting was used to calculate the number of RGB stars, which were normalized by the GC mass and fraction of observed GC luminosity, and determined density using the volume at the half-light radius $r_{h}$. The RGB star number density was weakly correlated with metallicity [Fe/H], giving Spearman and Kendall Rank test $p$-values of 0.00016 and 0.00021 and coefficients $r_{s} = 0.35$ and $\tau = 0.24$ respectively. This correlation may be biased by a possible dependence of $r_{h}$ on [Fe/H], although studies have shown that $r_{h}$ is correlated with Galactocentric distance and independent of [Fe/H]. The dynamical origin of the $r_{h}$-metallicity correlation (tidal stripping) suggests that metal-rich GCs may have had more active dynamical histories, which would promote LMXB formation. No correlation between the RGB star number density and metallicity was found when using only the GCs that hosted quiescent LMXBs. A complete census of quiescent LMXBs in our Galaxy is needed to further probe the metallicity effect, which will be possible with the upcoming launch of eROSITA.

This study represents the most sensitive Chandra X-ray point source catalogue of M31. Using 133 publicly available Chandra ACIS-I/S observations totalling ~1 Ms, we detected 795 X-ray sources in the bulge, northeast, and southwest fields of M31, covering an area of approximately 0.6 deg$^{2}$, to a limiting unabsorbed 0.5-8.0 keV luminosity of $10^{34}$ erg/s. In the inner bulge, where exposure is approximately constant, X-ray fluxes represent average values because they were determined from many observations over a long period of time. Similarly, our catalogue is more complete in the bulge fields since monitoring allowed more transient sources to be detected. The catalogue was cross-correlated with a previous XMM-Newton catalogue of M31's $D_{25}$ isophote consisting of 1948 X-ray sources, with only 979 within the field of view of our survey. We found 387 (49%) of our Chandra sources (352 or 44% unique sources) matched to within 5 arcsec of 352 XMM-Newton sources. Combining this result with matching done to previous Chandra X-ray sources we detected 259 new sources in our catalogue. We created X-ray luminosity functions (XLFs) in the soft (0.5-2.0 keV) and hard (2.0-8.0 keV) bands that are the most sensitive for any large galaxy based on our detection limits. Completeness-corrected XLFs show a break around $1.3\times10^{37}$ erg/s, consistent with previous work. As in past surveys, we find the bulge XLFs are flatter than the disk, indicating a lack of bright high-mass X-ray binaries in the disk and an aging population of low-mass X-ray binaries in the bulge.

X-ray binaries (XRBs) are probes of both star formation and stellar mass, but more importantly remain one of the only direct tracers of the compact object population. To investigate the XRB population in M31, we utilized all 121 publicly available observations of M31 totalling over 1 Ms from $\it{Chandra's}$ ACIS instrument. We studied 83 star clusters in the bulge using the year 1 star cluster catalogue from the Panchromatic Hubble Andromeda Treasury Survey. We found 15 unique star clusters that matched to 17 X-ray point sources within 1'' (3.8 pc). This population is composed predominantly of globular cluster low-mass XRBs, with one previously unidentified star cluster X-ray source. Star clusters that were brighter and more compact preferentially hosted an X-ray source. Specifically, logistic regression showed that the F475W magnitude was the most important predictor followed by the effective radius, while color (F475W$-$F814W) was not statistically significant. We also completed a matching analysis of 1566 HII regions and found 10 unique matches to 9 X-ray point sources within 3'' (11 pc). The HII regions hosting X-ray point sources were on average more compact than unmatched HII regions, but logistic regression concluded that neither the radius nor H$\alpha$ luminosity was a significant predictor. Four matches have no previous classification and thus are high-mass XRB candidates. A stacking analysis of both star clusters and HII regions resulted in non-detections, giving typical upper limits of $\approx10^{32}$ erg s$^{-1}$, which probes the quiescent XRB regime.

The population of low-luminosity (< 10^35 erg/s) X-Ray Binaries (XRBs) has been investigated in our Galaxy and M31 but not further. To address this problem, we have used data from the Chandra X-Ray Observatory and the Hubble Space Telescope to investigate the faint population of XRBs in the grand-design spiral galaxy M51. A matching analysis found 25 star clusters coincident with 20 X-ray point sources within 1.5" (60 pc). From X-ray and optical color-color plots we determine that this population is dominated by high-mass XRBs. A stacking analysis of the X-ray data at the positions of optically-identified star clusters was completed to probe low-luminosity X-ray sources. No cluster type had a significant detection in any X-ray energy band. An average globular cluster had the largest upper limit, 9.23 x 10^34 erg/s, in the full-band (0.3 - 8 keV) while on average the complete sample of clusters had the lowest upper limit, 6.46 x 10^33 erg/s in the hard-band (2 - 8 keV). We determined average luminosities of the young and old star cluster populations and compared the results to those from the Milky Way. We conclude that deeper X-ray data is required to identify faint sources with a stacking analysis.