An Empirical Framework Characterizing the Metallicity and Star-Formation History Dependence of X-ray Binary Population Formation and Emission in Galaxies

We present a new empirical framework modeling the metallicity and star-formation history (SFH) dependence of X-ray luminous ($L>{10}^{36}$ ergs s${}^{-1}$) point-source population luminosity functions (XLFs) in normal galaxies. We expect the X-ray point-source populations are dominated by X-ray binaries (XRBs), with contributions from supernova remnants near the low luminosity end of our observations. Our framework is calibrated using the collective statistical power of 3,731 X-ray detected point-sources within 88 Chandra-observed galaxies at $D<$ 40 Mpc that span broad ranges of metallicity ($Z\approx$ 0.03-2 $Z_{\odot }$), SFH, and morphology (dwarf irregulars, late-types, and early-types). Our best-fitting models indicate that the XLF normalization per unit stellar mass declines by $\approx$2-3 dex from 10 Myr to 10 Gyr, with a slower age decline for low-metallicity populations. The shape of the XLF for luminous X-ray sources ($L<{10}^{38}$ ergs s${}^{-1}$) significantly steepens with increasing age and metallicity, while the lower-luminosity XLF appears to flatten with increasing age. Integration of our models provide predictions for X-ray scaling relations that agree very well with past results presented in the literature, including, e.g., the $L_{\mathrm{X}}$-SFR-$Z$ relation for high-mass XRBs (HMXBs) in young stellar populations as well as the $L_{\mathrm{X}}\mathrm{/}{M}_{\star }$ ratio observed in early-type galaxies that harbor old populations of low-mass XRBs (LMXBs). The model framework and data sets presented in this paper further provide unique benchmarks that can be used for calibrating binary population synthesis models.