Search

We present Chandra and XMM-Newton X-ray observations of the compact group HCG 57, and optical integral field spectroscopy of the interacting galaxy pair HCG 57A/D. These two spiral galaxies recently suffered an off-axis collision with HCG 57D passing through the disk of A. We find evidence of a gas bridge linking the galaxies, containing ~10^8 Msol of hot, ~1 keV thermal plasma and warm ionized gas radiating in H$\alpha$, H$\beta$, [OIII] and [NII] lines. The optical emission lines in the central regions of HCG 57D show excitation properties consistent with HII-regions, while the outer rim of HCG 57D, parts of the bridge and the outer regions of HCG 57A show evidence of shocked gas consistent with shock velocities of 200-300 km/s. In contrast, the X-ray emitting gas requires a collision velocity of ~650-750 km/s to explain the observed temperatures. These different shock velocities can be reconciled by considering the contributions of rotation to collision velocity in different parts of the disks, and the clumpy nature of the pre-shock medium in the galaxies, which likely lead to different shock velocities in different components of the turbulent post-shocked gas. We examine the diffuse X-ray emission in the group members and their associated point sources, identifying X-ray AGN in HCG 57A, B, and D. We also confirm the previously reported ~1 keV intra-group medium and find it to be relaxed with a low central entropy (18.0+-1.7 kev cm^2 within 20 kpc) but a long cooling time (5.9+-0.8 Gyr)., Comment: 19 pages, 11 figures, accepted for publication in ApJ

Context. A major challenge in astrophysics is classifying galaxies by their activity. Current methods often require multiple diagnostics to capture the full range of galactic activity. Furthermore, overlapping excitation sources with similar observational signatures complicate the analysis of a galaxy's activity. Aims. This study aims to create an activity diagnostic tool that overcomes the limitations of current emission line diagnostics by identifying the underlying excitation mechanisms in mixed-activity galaxies (e.g., star formation, active nucleus, or old stellar populations) and determining the dominant ones. Methods. We use the random forest machine-learning algorithm, trained on three main activity classes -- star-forming, AGN, and passive -- that represent key gas excitation mechanisms. This diagnostic employs four distinguishing features: the equivalent widths of [O iii] ${\lambda}$5007, [N ii] ${\lambda}$6584, H${\alpha}$, and the D4000 continuum break index. Results. The classifier achieves near-perfect performance, with an overall accuracy of ~ 99% and recall scores of ~ 100% for star-forming, ~ 98% for AGN, and ~ 99% for passive galaxies. These exceptional scores allow for confident decomposition of mixed activity classes into the primary gas excitation mechanisms, overcoming the limitations of current classification methods. Additionally, the classifier can be simplified to a two-dimensional diagnostic using the D4000 index and log$_{10}$(EW([O iii])$^{2}$) without significant loss of diagnostic power. Conclusions. We present a diagnostic for classifying galaxies by their primary gas excitation mechanisms and deconstructing the activity of mixed-activity galaxies into these components. This method covers the full range of galaxy activity. Aditionally, D4000 index serves as an indicator for resolving the degeneracy among various activity components., Comment: 19 pages, 16 figures, Accepted for publication in A&A

Compact star-forming galaxies were dominant galaxy types in the early Universe. Blueberry galaxies (BBs) represent their local analogues being very compact and having intensive star formation. Motivated by high X-ray emission recently found in other analogical dwarf galaxies, called Green Peas, we probe into the X-ray properties of BBs to determine if their X-ray emission is consistent with the empirical laws for star-forming galaxies. We performed the first X-ray observations of a small sample of BBs with the XMM-Newton satellite. Spectral analysis for detected sources and upper limits measured via Bayesian-based analysis for very low-count measurements were used to determine the X-ray properties of our galaxy sample. Clear detection was obtained only for 2 sources, with one source exhibiting an enhanced X-ray luminosity to the scaling relations. For the remaining 5 sources, only an upper limit was constrained, suggesting BBs to be rather underluminous as a whole. Our analysis shows that the large scatter cannot be easily explained by the stochasticity effects. While the bright source is above and inconsistent at almost the 99% confidence level, the upper limits of the two sources are below the expected distribution. These results indicate that the empirical relations between the star formation rate, metallicity, and X-ray luminosity might not hold for BBs with uniquely high specific star formation rates. One possible explanation could be that the BBs may not be old enough to have a significant X-ray binary population. The high luminosity of the only bright source can be then caused by an additional X-ray source, such as a hidden active galactic nucleus or more extreme ultraluminous X-ray sources., Comment: 10 pages, 4 + 1 (appendix) figures, accepted for publication in Astronomy & Astrophysics on 9.8.2024

Using hard (E>10 keV) X-ray observations with NuSTAR, we are able to differentiate between accretion states, and thus compact object types, of neutron stars and black holes in X-ray binaries (XRBs) in M31, our nearest Milky Way-type neighbor. Using ten moderate-depth (20-50 ks) observations of the disk of M31 covering a total of ~0.45 deg$^{2}$, we detect 20 sources at 2$\sigma$ in the 4-25 keV band pass, 14 of which we consider to be XRB candidates. This complements an existing deeper (100-400 ks) survey covering ~0.2 deg$^{2}$ of the bulge and the northeastern disk. We make tentative classifications of 9 of these sources with the use of diagnostic color-intensity and color-color diagrams, which separate sources into various neutron star and black hole regimes, identifying 3 black holes and 6 neutron stars. In addition, we create X-ray luminosity functions for both the full (4-25 keV) and hard (12-25 keV) band, as well as sub-populations of the full band based on compact object type and association with globular clusters. Our best fit globular cluster XLF is shallower than the field XLF, and preliminary BH and NS XLFs suggest a difference in shape based on compact object type. We find that the cumulative disk XLFs in the full and hard band are best fit by power laws with indices of 1.32 and 1.28 respectively. This is consistent with models of the Milky Way XLF from Grimm et al. (2002), Voss & Ajello (2010), and Doroshenko et al. (2014)., Comment: 19 pages, 6 Figures, Published in ApJ

Dwarf highly star-forming galaxies (SFGs) dominated the early Universe and are considered the main driver of its reionization. However, direct observations of these distant galaxies are mainly confined to rest-frame ultraviolet and visible light, limiting our understanding of their complete properties. Therefore, it is still paramount to study their local analogs, the green pea (GP) and blueberry (BB) galaxies. This work aims to expand our knowledge of BBs by identifying a new sample that is closer and in the southern sky. In addition to the already known BBs, this new sample will allow for a statistically significant study of their properties probed by visible and infrared (IR) light. By utilizing the HECATE catalog, which provides photometry and characterization of galaxies, along with data from Pan-STARSS and SDSS, this study selects and analyzes a new sample of BBs. We employed spectral energy distribution fitting to derive homogeneous measurements of star-formation rates and stellar masses. Additionally, we measured emission-line fluxes, including $\rm HeII~\lambda 4686$, through spectral fitting. Through this work, we identified 48 BBs, of which 40 were first recognized as such, with the nearest at 19~Mpc. 14 of the BBs are in the south sky. The BBs tend to be extremely IR red in both WISE $W1-W2$ and $W2-W3$ colors, distinguishing them from typical SFGs. Dwarf SFGs with higher specific star-formation rates tend to have redder IR colors. Blueberry galaxies stand out as the most intensely star-forming sources in the local Universe among dwarf galaxies. They exhibit unique characteristics, such as being intrinsically bluer in visible light, redder in the infrared, and less massive. They also have higher specific star-formation rates, equivalent widths, lower metallicities, and the most strongly ionized interstellar medium compared to typical SFGs and GPs., Comment: 18 pages, 11 + 4 (appendix) figures. Accepted for publication in Astronomy & Astrophysics on 2024-05-02. Updated to match the published version

We present a study of X-ray normal galaxies using data from the first all-sky scan of the eROSITA X-ray survey. eRASS1 provides the first unbiased X-ray census of normal galaxies allowing us to study the X-ray emission from XRBs and the hot ISM in the full range of stellar population parameters present in the local Universe. By combining the HECATE value-added galaxy catalogue with the eRASS1, we study the X-ray emission from normal galaxies as a function of their SFR, M$_{*}$, Metallicity, and stellar population age. After applying optical and mid-IR activity classification criteria, we constructed a sample of 18790 star-forming galaxies with measurements of their L$_{X}$. By stacking the X-ray data in SFR-M$_{*}$-distance bins we study the correlation between the average L$_{X}$ and stellar population parameters. We also present updated L$_{\rm{X}}$-SFR and L$_{\rm{X}}$/SFR-Metallicity scaling relations accounting for the scatter dependence on the SFR. We find that the integrated L$_{X}$ of the HEC-eR1 star-forming galaxies is significantly elevated with respect to that expected from the current scaling relations. The observed scatter is also significantly larger. This excess persists even when we measure the average L$_{X}$ of galaxies in SFR-M$_{*}$-distance and metallicity bins and it is stronger in lower SFRs. The excess is not the result of hot gas, LMXBs, background AGN, LLAGN (including TDEs), or stochastic sampling of the XRB XLF. We find that while the excess correlates with lower metallicity, its primary driver is the age of the stellar populations. Our analysis reveals a sub-population of X-ray luminous starbursts with high sSFRs, low metallicities, and young stellar populations. This population drives upwards the X-ray scaling relations for star-forming galaxies, and has important implications for understanding the population of XRBs in the local and high-z Universe., Comment: Submitted to A&A

We construct simulated galaxy data sets based on the High Energy X-ray Probe (HEX-P) mission concept to demonstrate the significant advances in galaxy science that will be yielded by the HEX-P observatory. The combination of high spatial resolution imaging ($<$20 arcsec FWHM), broad spectral coverage (0.2-80 keV), and sensitivity superior to current facilities (e.g., XMM-Newton and NuSTAR) will enable HEX-P to detect hard (4-25 keV) X-ray emission from resolved point-source populations within $\sim$800 galaxies and integrated emission from $\sim$6000 galaxies out to 100 Mpc. These galaxies cover wide ranges of galaxy types (e.g., normal, starburst, and passive galaxies) and properties (e.g., metallicities and star-formation histories). In such galaxies, HEX-P will: (1) provide unique information about X-ray binary populations, including accretor demographics (black hole and neutron stars), distributions of accretion states and state transition cadences; (2) place order-of-magnitude more stringent constraints on inverse Compton emission associated with particle acceleration in starburst environments; and (3) put into clear context the contributions from X-ray emitting populations to both ionizing the surrounding interstellar medium in low-metallicity galaxies and heating the intergalactic medium in the $z > 8$ Universe., Comment: Accepted for publication in Frontiers in Astronomy and Space Sciences

(abridged) The overwhelming majority of diagnostic tools for galactic activity are focused on active galaxies. Passive or dormant galaxies are often excluded from these diagnostics which usually employ emission line features. In this work, we use infrared and optical colors in order to build an all-inclusive galactic activity diagnostic tool that can discriminate between star-forming, AGN, LINER, composite, and passive galaxies, and which can be used in local and low-redshift galaxies. We explore classification criteria based on infrared colors from the 3 WISE bands supplemented with optical colors from the u, g, and r SDSS bands. From these we aim to find the minimal combination of colors for optimal results. Furthermore, to mitigate biases related to aperture effects, we introduce a new WISE photometric scheme combing different sized apertures. We develop a diagnostic tool using machine learning methods that includes both active and passive galaxies under one unified scheme using 3 colors. We find that the combination of W1-W2, W2-W3, and g-r colors offers good performance while the broad availability of these colors for a large number of galaxies ensures wide applicability on large galaxy samples. The overall accuracy is $\sim$81% while the achieved completeness for each class is $\sim$81% for star-forming, $\sim$56% for AGN, $\sim$68% for LINER, $\sim$65% for composite, and $\sim$85% for passive galaxies. Our diagnostic provides a significant improvement over existing IR diagnostics by including all types of active, as well as passive galaxies, and extending them to the local Universe. The inclusion of the optical colors improves their performance in identifying low-luminosity AGN which are generally confused with star-forming galaxies, and helps to identify cases of starbursts with extreme mid-IR colors which mimic obscured AGN galaxies, a well-known problem for most IR diagnostics., Comment: Accepted for publication in the A&A journal. The code for the application of our model can be accessed through the GitHub repository in https://github.com/BabisDaoutis/GalActivityClassifier

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 $\sim$10) compared with predictions from scaling relations of galaxy-wide high mass X-ray binary (HMXB) luminosity with star-formation rate. In the nearby ($\approx$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 ($\lesssim$ 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_{\rm{IR}}$ within these regions on timescales $\lesssim$ 5 Myr, shorter than the timescale required for the formation of HMXBs., Comment: Accepted for publication in ApJ (17 pages, 3 tables, 8 figures)

X-ray appearance of normal galaxies is mainly determined by X-ray binaries powered by accretion onto a neutron star or a stellar mass black hole. Their populations scale with the star-formation rate and stellar mass of the host galaxy and their X-ray luminosity distributions show a significant split between star-forming and passive galaxies, both facts being consequences of the dichotomy between high- and low-mass X-ray binaries. Metallicity, IMF and stellar age dependencies, and dynamical formation channels add complexity to this picture. The numbers of high-mass X-ray binaries observed in star-forming galaxies indicate quite high probability for a massive star to become an accretion powered X-ray source once upon its lifetime. This explains the unexpectedly high contribution of X-ray binaries to the Cosmic X-ray Background, of the order of $\sim 10\%$, mostly via X-ray emission of faint star-forming galaxies located at moderate redshifts which may account for the unresolved part of the CXB. Cosmological evolution of the $L_X-{\rm SFR}$ relation can make high-mass X-ray binaries a potentially significant factor in (pre)heating of intergalactic medium in the early Universe., Comment: Invited chapter for the Handbook of X-ray and Gamma-ray Astrophysics. Editors: Cosimo Bambi, Andrea Santangelo. Publisher: Springer Singapore, 2023

We use the Random Forest (RF) algorithm to develop a tool for automated activity classification of galaxies into 5 different classes: Star-forming (SF), AGN, LINER, Composite, and Passive. We train the algorithm on a combination of mid-IR (WISE) and optical photometric data while the true labels (activity classes) are based on emission line ratios. Our classifier is built to be redshift-agnostic and it is applicable to objects up to z $\sim$0.1. It reaches a completeness $>$80 % for SF and Passive galaxies, and $\sim$60 % for AGN. Applying it to an all-sky galaxy catalog (HECATE) reveals a large population of low-luminosity AGNs outside the AGN locus in the standard mid-IR diagnostics., Comment: 3 pages, 3 figures, the manuscript will be published on Memorie della SAIt

The stellar mass ($M_\star$) and the star-formation rate (SFR) are among the most important features that characterize galaxies. Measuring these fundamental properties accurately is critical for understanding the present state of galaxies, and their history. This work explores the dependence of the IR emission of galaxies on their extinction, and the age of their stellar populations (SPs). It aims at providing accurate IR SFR and $M_\star$ calibrations that account for SP age and extinction while quantifying their scatter. We use the CIGALE spectral energy distribution (SED) fitting code to create models of galaxies with a wide range of star-formation histories, dust content, and interstellar medium properties. We fit the relations between $M_\star$ and SFR with IR and optical photometry of the model-galaxy SEDs with the MCMC method, and perform a machine-learning random forest analysis on the same data set in order to validate the latter. This work provides calibrations for the SFR using a combination of the WISE bands 1 and 3, or the JWST F200W and F2100W bands. It also provides mass-to-light ratio calibrations based on the WISE band-1, or the JWST band F200W, along with the optical $u-r$ or $g-r$ colors. These calibrations account for the biases attributed to the SP age, while they are given in the form of extinction-dependent and extinction-independent relations. They show robust estimations while minimizing the scatter and biases throughout a wide range of SFRs and stellar masses. The SFR calibration offers better results, especially in dust-free or passive galaxies where the contributions of old SPs or biases from the lack of dust are significant. Similarly, the $M_\star$ calibration yields significantly better results for dusty/high-SFR galaxies where dust emission can otherwise bias the estimations., Comment: 18 pages, 10 figures. Accepted for publication in Astronomy & Astrophysics on 16 March 2023

The ever-expanding observational sample of X-ray binaries (XRBs) makes them excellent laboratories for constraining binary evolution theory. Such constraints can be obtained by studying the effects of various physical assumptions on synthetic X-ray luminosity functions (XLFs) and comparing to observed XLFs. In this work, we focus on high-mass XRBs (HMXBs) and study the effects on the XLF of various, poorly-constrained assumptions regarding physical processes such as the common-envelope phase, the core-collapse, and wind-fed accretion. We use the new binary population synthesis code POSYDON, which employs extensive pre-computed grids of detailed stellar structure and binary evolution models, to simulate the evolution of binaries. We generate 96 synthetic XRB populations corresponding to different combinations of model assumptions. The generated HMXB XLFs are feature-rich, deviating from the commonly assumed single-power law. We find a break in our synthetic XLF at luminosity $\sim 10^{38}$ erg s$^{-1}$, similar to observed XLFs. However, we find also a general overabundance of XRBs (up to a factor of $\sim$10 for certain model parameter combinations) driven primarily by XRBs with black hole accretors. Assumptions about the transient behavior of Be-XRBs, asymmetric supernova kicks, and common-envelope physics can significantly affect the shape and normalization of our synthetic XLFs. We find that less well-studied assumptions regarding the circularization of the orbit at the onset of Roche-lobe overflow and criteria for the formation of an X-ray emitting accretion disk around wind-accreting black holes can also impact our synthetic XLFs. Our study reveals the importance of large-scale parameter studies, highlighting the power of XRBs in constraining binary evolution theory., Comment: 31 pages, 32 figures, Accepted by A&A. Fixed typos and updated references. Referee's comments were addressed

Data from high-energy observations are usually obtained as lists of photon events. A common analysis task for such data is to identify whether diffuse emission exists, and to estimate its surface brightness, even in the presence of point sources that may be superposed. We have developed a novel non-parametric event list segmentation algorithm to divide up the field of view into distinct emission components. We use photon location data directly, without binning them into an image. We first construct a graph from the Voronoi tessellation of the observed photon locations and then grow segments using a new adaptation of seeded region growing, that we call Seeded Region Growing on Graph, after which the overall method is named SRGonG. Starting with a set of seed locations, this results in an over-segmented dataset, which SRGonG then coalesces using a greedy algorithm where adjacent segments are merged to minimize a model comparison statistic; we use the Bayesian Information Criterion. Using SRGonG we are able to identify point-like and diffuse extended sources in the data with equal facility. We validate SRGonG using simulations, demonstrating that it is capable of discerning irregularly shaped low surface-brightness emission structures as well as point-like sources with strengths comparable to that seen in typical X-ray data. We demonstrate SRGonG's use on the Chandra data of the Antennae galaxies, and show that it segments the complex structures appropriately.

We present a basic model for the calculation of the luminosity distribution of supernova remnant populations. We construct theoretical Ha and joint [S II] - Ha luminosity functions for supernova remnants by combining prescriptions from a basic evolution model that provides the shock velocity and radius for SNRs of different age and pre-shock density, with shock excitation models that give the gas emissivity for shocks of different physical parameters. We assume a flat age distribution, and we explore the effect of different pre-shock density distributions or different magnetic parameters. We find very good agreement between the shape of the model Ha and the joint [S II] - Ha luminosity functions and those measured from SNR surveys in nearby galaxies., Comment: 11 pages, Accepted for publication in MNRAS

In this paper we use an RXTE library of spectral models from 10 black-hole and 9 pulsar X-ray binaries, as well as model spectra available in the literature from 13 extra-galactic Ultra-luminous X-ray sources (ULXs). We compute average bolometric corrections (BC=$\mathrm{L_{band}/L_{bol}}$) for our sample as a function of different accretion rates. We notice the same behaviour between black-hole and pulsars BCs only when ULX pulsars are included. These measurements provide a picture of the energetics of the accretion flow for an X-ray binary based solely on its observed luminosity in a given band. Moreover it can be a powerful tool in X-ray binary population synthesis models. Furthermore we calculate the X-ray (2-10 keV) to optical (V-band) flux ratios at different Eddington ratios for the black-hole X-ray binaries in our sample. This provides a metric of the maximum contribution of the disk to the optical emission of a binary system and better constraints on its nature (donor type etc). We find that the optical to X-ray flux ratio shows very little variation as a function of accretion rate, but testing for different disk geometries scenarios we find that the optical contribution of the disk increases as the $p$ value decreases ($T(r)\sim r^{-p}$). Moreover observational data are in agreement with a thicker disk scenario ($p<0.65$), which could also possibly explain the lack of observed high-inclination systems., Comment: 16 pages, 6 figures, 14 tables, accepted for publication in MNRAS

We report the discovery of statistically significant spatial structures in the projected two-dimensional distributions of the Globular Cluster (GC) systems of 10 among the brightest galaxies in the Fornax Cluster. We use a catalog of GCs extracted from the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) Fornax Cluster Survey (ACSFCS) imaging data. We characterize the size, shape and location relative to the host galaxies of the GC structures and suggest a classification based on their morphology and location that is suggestive of different formation mechanisms. We also investigate the GC structures in the context of the positions of their host galaxies relative to the general spatial distributions of galaxies and intra-cluster globular clusters in the Fornax Cluster. We finally estimate the dynamical masses of the progenitors of some GC structures, under the assumption that they are the relics of past accretion events of satellite galaxies by their hosts., Comment: 40 pages, 26 figures, accepted for publication in ApJ

Context. ULXs ($L_\textrm{X}>10^{39}$ erg s$^{-1}$) are excellent probes of accretion physics, star formation and IMBHs searches. As the sample size of X-ray data from modern observatories increases, producing extensive catalogues of ULXs and studying their collective properties is a possibility and a priority. Aims. We build a ULX catalogue based on one of the latest XMM-Newton releases, 4XMM-DR9, and the galaxy catalogue HECATE and study whether the properties of the expanded XMM-Newton ULX population are consistent with previous findings. Methods. We perform cross-matching between XMM-Newton sources and HECATE objects to identify host galaxies, and flag known interlopers in external catalogues and databases. Manual inspection of image data from PanSTARRS1 and the NASA/IPAC database is occasionally performed. We use distance and luminosity arguments to identify ULX candidates. Spectral, abundance and variability properties of the candidates are studied from 4XMM-DR9, HECATE and 4XMM-DR9s parameters. Results. We identify 779 ULX candidates. In spiral galaxies the number of ULX candidates per star forming rate is consistent with previous studies, and a significant ULX population in elliptical and lenticular is also found. Candidates hosted by late-type galaxies are more abundant, present harder spectra and undergo more and more extreme inter-observation variability than the early-type hosted ones. Around 30 candidates with $L_\textrm{X}>10^{41}$ erg s$^{-1}$ are also identified. Conclusions. Our results on the spectral and abundance properties of ULXs confirm the findings made by previous studies based on XMM-Newton and Chandra data, while our population-scale study on variability properties is unprecedented. Our study provides limited insight on the properties of the brightest ULX candidates due to the small sample size., Comment: Paper accepted for publication in A&A on 2021-10-13. 18 pages, 22 figures, 8 tables. Catalogue soon to be uploaded to the CDS - VizieR database

(abridged) We develop a tool for the automated spectral classification of OB stars according to their sub-types. We use the regular Random Forest (RF) algorithm, the Probabilistic RF (PRF), and we introduce the KDE-RF method which is a combination of the Kernel-Density Estimation and the RF algorithm. We train the algorithms on the Equivalent Width (EW) of characteristic absorption lines (features) measured in high-quality spectra from large Galactic (LAMOST,GOSSS) and extragalactic surveys (2dF,VFTS) with available spectral-types and luminosity classes. We find that the overall accuracy score is $\sim$70% with similar results across all approaches. We show that the full set of 17 spectral lines is needed to reach the maximum performance per spectral class. We apply our model in other observational data sets providing examples of potential application of our classifier on real science cases. We find that it performs well for both single massive stars and for the companion massive stars in Be X-ray binaries. In addition, we propose a reduced 10-features scheme that can be applied to large data sets with lower S/N. The similarity in the performances of our models indicates the robustness and the reliability of the RF algorithm when it is used for the spectral classification of early-type stars. The score of $\sim$70% is high if we consider (a) the complexity of such multi-class classification problems, (b) the intrinsic scatter of the EW distributions within the examined spectral classes, and (c) the diversity of the training set since we use data obtained from different surveys with different observing strategies. In addition, the approach presented in this work, is applicable to data of different quality and of different format (e.g.,absolute or normalized flux) while our classifier is agnostic to the Luminosity Class of a star and, as much as possible, metallicity independent., Comment: Accepted for publication in the A&A journal. The code for the application of our models is available through the GitHub repository https://github.com/EliasKyritsis/SpectralClassifier

We present a systematic study of the Supernova Remnant (SNR) populations in the nearby galaxies NGC 45, NGC 55, NGC 1313, and NGC 7793 based on deep Ha and [S II] imaging. We find 42 candidate and 54 possible candidate SNRs based on the [S II] / Ha > 0.4 criterion, 84 of which are new identifications. We derive the Ha and the joint [S II]-Ha luminosity functions after accounting for incompleteness effects. We find that the Ha luminosity function of the overall sample is described with a skewed Gaussian with a mean equal to log(LHa / 10^(36) erg s^(-1)) = 0.07 and m(log(LHa / 10^(36) erg s^(-1)))) = 0.58. The joint [S II]-Ha function is parameterized by a skewed Gaussian along the log([S II] / 10^(36) erg s^(-1)) = 0.88 x log(LHa / 10^(36) erg s^(-1)) - 0.06 line and a truncated Gaussian with m(log(L[S II] / 10^(36))) = 0.024 and s(log(L[S II] / 10^(36))) = 0.14, on its vertical direction. We also define the excitation function as the number density of SNRs as a function of their [S II]/Ha ratios. This function is represented by a truncated Gaussian with a mean at -0.014. We find a sub-linear [S II]-Ha relation indicating lower excitation for the more luminous objects., Comment: Accepted in MNRAS. 19 pages, 21 figures

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 sigma) 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 sigma) 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. 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., Comment: 10 pages, 7 figures, accepted for publication in the Astrophysical Journal

We constrain the mass distribution in nearby, star-forming galaxies with the Star Formation Reference Survey (SFRS), a galaxy sample constructed to be representative of all known combinations of star formation rate (SFR), dust temperature, and specific star formation rate (sSFR) that exist in the Local Universe. An innovative two-dimensional bulge/disk decomposition of the 2MASS/$K_{s}$-band images of the SFRS galaxies yields global luminosity and stellar mass functions, along with separate mass functions for their bulges and disks. These accurate mass functions cover the full range from dwarf galaxies to large spirals, and are representative of star-forming galaxies selected based on their infra-red luminosity, unbiased by AGN content and environment. We measure an integrated luminosity density $j$ = 1.72 $\pm$ 0.93 $\times$ 10$^{9}$ L$_{\odot}$ $h^{-1}$ Mpc$^{-3}$ and a total stellar mass density $\rho_{M}$ = 4.61 $\pm$ 2.40 $\times$ 10$^{8}$ M$_{\odot}$ $h^{-1}$ Mpc$^{-3}$. While the stellar mass of the \emph{average} star-forming galaxy is equally distributed between its sub-components, disks globally dominate the mass density budget by a ratio 4:1 with respect to bulges. In particular, our functions suggest that recent star formation happened primarily in massive systems, where they have yielded a disk stellar mass density larger than that of bulges by more than 1 dex. Our results constitute a reference benchmark for models addressing the assembly of stellar mass on the bulges and disks of local ($z = 0$) star-forming galaxies., Comment: 33 pages, 21 figures, published in MNRAS, 504, 3831-3861 (2021)

We present new H$\alpha$ photometry for the Star-Formation Reference Survey (SFRS), a representative sample of star-forming galaxies in the local Universe. Combining these data with the panchromatic coverage of the SFRS, we provide calibrations of H$\alpha$-based star-formation rates (SFRs) with and without correction for the contribution of [$\rm N_{^{II}}$] emission. We consider the effect of extinction corrections based on the Balmer decrement, infrared excess (IRX), and spectral energy distribution (SED) fits. We compare the SFR estimates derived from SED fits, polycyclic aromatic hydrocarbons, hybrid indicators such as 24 $\mu$m + H$\alpha$, 8 $\mu$m + H$\alpha$, FIR + FUV, and H$\alpha$ emission for a sample of purely star-forming galaxies. We provide a new calibration for 1.4 GHz-based SFRs by comparing to the H$\alpha$ emission, and we measure a dependence of the radio-to-H$\alpha$ emission ratio based on galaxy stellar mass. Active galactic nuclei introduce biases in the calibrations of different SFR indicators but have only a minimal effect on the inferred SFR densities from galaxy surveys. Finally, we quantify the correlation between galaxy metallicity and extinction., Comment: 20 pages, 21 Figures, accepted for publication in MNRAS

The $\it{eROSITA}$ Final Equatorial Depth Survey (eFEDS), completed during the calibration and performance verification phase of the $\it{eROSITA}$ instrument on $\it{Spectrum\, Roentgen\, Gamma}$, delivers data at and beyond the final depth of the four-year $\it{eROSITA}$ all-sky survey (eRASS:8), $f_{0.5-2\,\text{ keV}}$ = $1.1\times10^{-14}$ erg s$^{-1}$ cm$^{2}$, over 140 deg$^{2}$. It provides the first view of normal galaxy X-ray emission from X-ray binaries (XRBs) and the hot interstellar medium at the full depth of eRASS:8. We use the Heraklion Extragalactic Catalogue (HECATE) of galaxies to correlate with eFEDS X-ray sources and identify 94 X-ray detected normal galaxies. We classify galaxies as star-forming, early-type, composite, and AGN using SDSS and 6dF optical spectroscopy. The eFEDS field harbours 37 normal galaxies: 36 late-type (star-forming) galaxies and 1 early-type galaxy. There are 1.9 times as many normal galaxies as predicted by scaling relations via SIXTE simulations, with an overabundance of late-type galaxies and a dearth of early-type galaxies. Dwarf galaxies with high specific star formation rate (SFR) have elevated L$_{\text{X}}$/SFR when compared with specific SFR and metallicity, indicating an increase in XRB emission due to low-metallicity. We expect that eRASS:8 will detect 12,500 normal galaxies, the majority of which will be star-forming, with the caveat that there are unclassified sources in eFEDS and galaxy catalogue incompleteness issues that could increase the actual number of detected galaxies over these current estimates. eFEDS observations detected a rare population of galaxies -- the metal-poor dwarf starbursts -- that do not follow known scaling relations. eRASS is expected to discover significant numbers of these high-redshift analogues, which are important for studying the heating of the intergalactic medium at high-redshift., Comment: 17 pages, 10 figures, 3 tables. Submitted to A&A for the Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission

We present the Heraklion Extragalactic Catalogue, or HECATE, an all-sky value-added galaxy catalogue, aiming to facilitate present and future multi-wavelength and multi-messenger studies in the local Universe. It contains 204,733 galaxies up to a redshift of 0.047 (D<200 Mpc), and it is >50% complete in terms of the B-band luminosity density at distances in the 0-170 Mpc range. By incorporating and homogenising data from astronomical databases and multi-wavelength surveys, the catalogue offers positions, sizes, distances, morphological classifications, star-formation rates, stellar masses, metallicities, and nuclear activity classifications. This wealth of information can enable a wide-range of applications, such as: (i) demographic studies of extragalactic sources, (ii) initial characterisation of transient events, and (iii) searches for electromagnetic counterparts of gravitational-wave events. The catalogue is publicly available to the community at a dedicated portal, which will also host future extensions in terms of the covered volume, and data products., Comment: 21 pages. 11 figures, 5 tables. This article has been accepted for publication in the Monthly Notices of the Royal Astronomical Society, published by Oxford University Press on behalf of the Royal Astronomical Society

The analysis of individual X-ray sources that appear in a crowded field can easily be compromised by the misallocation of recorded events to their originating sources. Even with a small number of sources, that nonetheless have overlapping point spread functions, the allocation of events to sources is a complex task that is subject to uncertainty. We develop a Bayesian method designed to sift high-energy photon events from multiple sources with overlapping point spread functions, leveraging the differences in their spatial, spectral, and temporal signatures. The method probabilistically assigns each event to a given source. Such a disentanglement allows more detailed spectral or temporal analysis to focus on the individual component in isolation, free of contamination from other sources or the background. We are also able to compute source parameters of interest like their locations, relative brightness, and background contamination, while accounting for the uncertainty in event assignments. Simulation studies that include event arrival time information demonstrate that the temporal component improves event disambiguation beyond using only spatial and spectral information. The proposed methods correctly allocate up to 65% more events than the corresponding algorithms that ignore event arrival time information. We apply our methods to two stellar X-ray binaries, UV Cet and HBC515 A, observed with Chandra. We demonstrate that our methods are capable of removing the contamination due to a strong flare on UV Cet B in its companion approximately 40 times weaker during that event, and that evidence for spectral variability at timescales of a few ks can be determined in HBC515 Aa and HBC515 Ab.

Many astrophysical phenomena are time-varying, in the sense that their intensity, energy spectrum, and/or the spatial distribution of the emission suddenly change. This paper develops a method for modeling a time series of images. Under the assumption that the arrival times of the photons follow a Poisson process, the data are binned into 4D grids of voxels (time, energy band, and x-y coordinates), and viewed as a time series of non-homogeneous Poisson images. The method assumes that at each time point, the corresponding multi-band image stack is an unknown 3D piecewise constant function including Poisson noise. It also assumes that all image stacks between any two adjacent change points (in time domain) share the same unknown piecewise constant function. The proposed method is designed to estimate the number and the locations of all the change points (in time domain), as well as all the unknown piecewise constant functions between any pairs of the change points. The method applies the minimum description length (MDL) principle to perform this task. A practical algorithm is also developed to solve the corresponding complicated optimization problem. Simulation experiments and applications to real datasets show that the proposed method enjoys very promising empirical properties. Applications to two real datasets, the XMM observation of a flaring star and an emerging solar coronal loop, illustrate the usage of the proposed method and the scientific insight gained from it., Comment: 22 pages, 10 figures

We present a new method to extract statistical constraints on the progenitor properties and formation channels of individual gravitational-wave sources. Although many different models have been proposed to explain the binary black holes detected by the LIGO Scientific and Virgo Collaboration (LVC), formation through isolated binary evolution remains the best explored channel. Under the assumption of formation through binary evolution, we use the statistical wrapper dart_board coupled with the rapid binary evolution code COSMIC to model the progenitor of GW150914, the first gravitational-wave signal detected by the LVC. Our Bayesian method combines the likelihood generated from the gravitational-wave signal with a prior describing the population of stellar binaries, and the Universe's star-formation and metallicity evolution. We find that the dominant evolutionary channel for GW150914 did not involve a common-envelope phase, but instead the system most probably (70%-90%) formed through stable mass transfer. This result is robust against variations of various model parameters, and it is reversed only when dynamical instability in binaries becomes more likely when a strict condition favoring common envelopes is adopted. Our analysis additionally provides a quantitative description of the progenitors relevant to each channel., Comment: 14 pages, 5 figures, 1 table, accepted for publication in ApJL

We present detailed constraints on the metallicity dependence of the high mass X-ray binary (HMXB) X-ray luminosity function (XLF). We analyze ~5 Ms of Chandra data for 55 actively star-forming galaxies at D < 30 Mpc with gas-phase metallicities spanning 12 + log(O/H) = 7-9.2. Within the galactic footprints, our sample contains a total of 1311 X-ray point sources, of which ~49% are expected to be HMXBs, with the remaining sources likely to be low-mass X-ray binaries (LMXBs; ~22%) and unrelated background sources (~29%). We construct a model that successfully characterizes the average HMXB XLF over the full metallicity range. We demonstrate that the SFR-normalized HMXB XLF shows clear trends with metallicity, with steadily increasing numbers of luminous and ultraluminous X-ray sources (logL(erg/s) = 38-40.5) with declining metallicity. However, we find that the low-luminosity (logL(erg/s) = 36-38) HMXB XLF appears to show a nearly constant SFR scaling and slope with metallicity. Our model provides a revised scaling relation of integrated LX/SFR versus 12 + log(O/H) and a new characterization of its the SFR-dependent stochastic scatter. The general trend of this relation is broadly consistent with past studies based on integrated galaxy emission; however, our model suggests that this relation is driven primarily by the high-luminosity end of the HMXB XLF. Our results have implications for binary population synthesis models, the nature of super-Eddington accreting objects (e.g., ultraluminous X-ray sources), recent efforts to identify active galactic nucleus candidates in dwarf galaxies, and the X-ray radiation fields in the early Universe during the epoch of cosmic heating at z > 10., Comment: Accepted for publication in ApJ

We present a systematic study of the metallicity variations within the collisional ring galaxy NGC 922 based on long-slit optical spectroscopic observations. We find a metallicity difference between star-forming regions in the bulge and the ring, with metallicities ranging from almost solar to significantly sub-solar ($\rm{[12+\log(O/H)]\sim 8.2}$). We detect $\rm He\,_I$ emission in the bulge and the ring star-forming regions indicating ionization from massive stars associated with recent ($<10$ Myr) star-formation, in agreement with the presence of very young star-clusters in all studied regions. We find an anti-correlation between the X-ray luminosity and metallicity of the sub-galactic regions of NGC 922. The different regions have similar stellar population ages leaving metallicity as the main driver of the anti-correlation. The dependence of the X-ray emission of the different regions in NGC 922 on metallicity is in agreement with similar studies of the integrated X-ray output of galaxies and predictions from X-ray binary population models., Comment: 15 pages, 10 figures

Using the Chandra Source Catalog 2.0 and a newly compiled catalogue of galaxies in the local Universe, we deliver a census of ultraluminous X-ray source (ULX) populations in nearby galaxies. We find 629 ULX candidates in 309 galaxies with distance smaller than 40\,Mpc. The foreground/background contamination is ${\sim}20\%$. The ULX populations in bona-fide star-forming galaxies scale on average with star-formation rate (SFR) and stellar mass ($M_\star$) such that the number of ULXs per galaxy is $0.45^{+0.06}_{-0.09}\times\frac{\rm SFR}{\rm M_\odot\,yr^{-1}}{+}3.3^{+3.8}_{-3.2}\times\frac{M_\star}{\rm M_\odot}$. The scaling depends strongly on the morphological type. This analysis shows that early spiral galaxies contain an additional population of ULXs that scales with $M_\star$. We also confirm the strong anti-correlation of the ULX rate with the host galaxy's metallicity. In the case of early-type galaxies we find that there is a non-linear dependence of the number of ULXs with $M_\star$, which is interpreted as the result of star-formation history differences. Taking into account age and metallicity effects, we find that the predictions from X-ray binary population synthesis models are consistent with the observed ULX rates in early-type galaxies, as well as, spiral/irregular galaxies., Comment: 21 pages, 19 figures, 6 tables. Accepted for publication in MNRAS

eROSITA, launched on 13 July 2019, will be completing the first all-sky survey in the soft and medium X-ray band in nearly three decades. This 4-year survey, finishing in late 2023, will present a rich legacy for the entire astrophysics community and complement upcoming multi-wavelength surveys (with, e.g. the Large Synoptic Survey Telescope and the Dark Energy Survey). Besides the major scientific aim to study active galactic nuclei (AGN) and galaxy clusters, eROSITA will contribute significantly to X-ray studies of normal (i.e., not AGN) galaxies. Starting from multi-wavelength catalogues, we measure star formation rates and stellar masses for 60 212 galaxies constrained to distances of 50-200 Mpc. We chose this distance range to focus on the relatively unexplored volume outside the local Universe, where galaxies will be largely spatially unresolved and probe a range of X-ray luminosities that overlap with the low luminosity and/or highly obscured AGN population. We use the most recent X-ray scaling relations as well as the on-orbit eROSITA instrument performance to predict the X-ray emission from XRBs and diffuse hot gas and to perform both an analytic prediction and an end-to-end simulation using the mission simulation software, SIXTE. We consider potential contributions from hidden AGN and comment on the impact of normal galaxies on the measurement of the faint end of the AGN luminosity function. We predict that the eROSITA 4-year survey, will detect $\gtrsim$ 15 000 galaxies (3 $\sigma$ significance) at 50-200 Mpc, which is ~100X more normal galaxies than detected in any X-ray survey to date., Comment: 18 pages, 15 figures. Accepted for publication in MNRAS

One distinctive feature of low-luminosity active galactic nuclei (LLAGN) is the relatively weak reflection features they may display in the X-ray spectrum, which can result from the disappearance of the torus with decreasing accretion rates. Some material, however, must surround the active nucleus, i.e., the accretion flow itself and, possibly, a flattened-out or thinned torus. In this work, we study whether reflection is indeed absent or undetectable due to its intrinsically weak features together with the low statistics inherent to LLAGN. Here we focus on NGC 3718 ($L/L_{\rm Edd}\sim10^{-5}$) combining observations from XMM--Newton and the deepest to date NuSTAR (0.5--79 keV) spectrum of a LLAGN, to constrain potential reflectors, and analyze how the fitted coronal parameters depend on the reflection model. We test models representing both an accretion disc (Relxill) and a torus-like (MYTorus and Borus) neutral reflector. From a statistical point of view, reflection is not required, but its inclusion allows to place strong constraints on the geometry and physical features of the surroundings: both neutral reflectors (torus) tested should be Compton thin ($N_H<10^{23.2}$cm$^{-2}$) and preferentially cover a large fraction of the sky. If the reflected light instead arises from an ionized reflector, a highly ionized case is preferred. These models produce an intrinsic power-law spectral index in the range [1.81--1.87], where the torus models result in steeper slopes. The cut-off energy of the power-law emission also changes with the inclusion of reflection models, resulting in constrained values for the disc reflectors and unconstrained values for torus reflectors., Comment: 15 pages and 16 Figures

X-ray luminosity ($L_X$) originating from high-mass X-ray binaries (HMXBs) is tightly correlated with the host galaxy's star-formation rate (SFR). We explore this connection at sub-galactic scales spanning ${\sim}$7 dex in SFR and ${\sim}$8 dex in specific SFR (sSFR). There is good agreement with established relations down to ${\rm SFR {\simeq} 10^{-3}\,M_\odot \, yr^{-1}}$, below which an excess of X-ray luminosity emerges. This excess likely arises from low mass X-ray binaries. The intrinsic scatter of the $L_X$-SFR relation is constant, not correlated with SFR. Different star formation indicators scale with $L_X$ in different ways, and we attribute the differences to the effect of star formation history. The SFR derived from H$\alpha$ shows the tightest correlation with X-ray luminosity because H$\alpha$ emission probes stellar populations with ages similar to HMXB formation timescales, but the H$\alpha$-based SFR is reliable only for $\rm sSFR{>}10^{-12}\,M_\odot \, yr^{-1}/M_\odot$., Comment: 19 pages, 13 figures, accepted for publication in MNRAS (2020-04-19)

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., Comment: 8 tables, 18 figures

We present the results of the photometric analysis of a large part of the main body of the Small Magellanic Cloud. Using the 6.5m Magellan Telescope at the Las Campanas Observatory in Chile, we have acquired deep B and I images in four fields (0.44 degree each in diameter), yielding accurate photometry for 1,068,893 stars down to 24th magnitude, with a spatial resolution of 0.20 arcsec per pixel. Colour-magnitude diagrams and (completeness corrected) luminosity functions have been constructed, yielding significant new results that indicate at least two discrete star formation events over a period from 2.7 to 4 Gyr ago. Also, we have derived star formation rates as a function of look back time and have found enhancements of SF between 4-6 Gyr and at younger ages., Comment: 12 pages, 11 figures; accepted for publication in MNRAS

We present new diagnostic tools for distinguishing supernova remnants (SNRs) from HII regions. Up to now, sources with flux ratio [S II]/H$\rm{\alpha}$ higher than 0.4 have been considered as SNRs. Here, we present the combinations of three or two line ratios as more effective tools for the separation of these two kinds of nebulae, depicting them as 3D surfaces or 2D lines. The diagnostics are based on photoionization and shock excitation models (MAPPINGS III) analysed with Support Vector Machine (SVM) models for classification. The line-ratio combination that gives the most efficient diagnostic is: [O I]/H$\rm{\alpha}$ - [O II]/H$\rm{\beta}$ - [O III]/H$\rm{\beta}$. This method gives $98.95\% $ completeness in the SNR selection and $1.20\%$ contamination. We also define the [O I]/H$\rm{\alpha}$ SNR selection criterion and we measure its efficiency in comparison to other selection criteria., Comment: Accepted for publication in the Monthly Notices of the Royal Astronomical Society, 18 pages, 30 figures

Over forty years of research suggests that the common envelope phase, in which an evolved star engulfs its companion upon expansion, is the critical evolutionary stage forming short-period, compact-object binary systems, such as coalescing double compact objects, X-ray binaries, and cataclysmic variables. In this work, we adapt the one-dimensional hydrodynamic stellar evolution code, MESA, to model the inspiral of a 1.4M$_{\odot}$ neutron star (NS) inside the envelope of a 12$M_{\odot}$ red supergiant star. We self-consistently calculate the drag force experienced by the NS as well as the back-reaction onto the expanding envelope as the NS spirals in. Nearly all of the hydrogen envelope escapes, expanding to large radii ($\sim$10$^2$ AU) where it forms an optically thick envelope with temperatures low enough that dust formation occurs. We simulate the NS orbit until only 0.8M$_{\odot}$ of the hydrogen envelope remains around the giant star's core. Our results suggest that the inspiral will continue until another $\approx$0.3M$_{\odot}$ are removed, at which point the remaining envelope will retract. Upon separation, a phase of dynamically stable mass transfer onto the NS accretor is likely to ensue, which may be observable as an ultraluminous X-ray source. The resulting binary, comprised of a detached 2.6M$_{\odot}$ helium-star and a NS with a separation of 3.3-5.7R$_{\odot}$, is expected to evolve into a merging double neutron-star, analogous to those recently detected by LIGO/Virgo. For our chosen combination of binary parameters, our estimated final separation (including the phase of stable mass transfer) suggests a very high $\alpha_{\rm CE}$-equivalent efficiency of $\approx$5., Comment: 10 pages, 5 figures, submitted to The Astrophysical Journal Letters. Comments welcome

We present new Chandra constraints on the X-ray luminosity functions (XLFs) of X-ray binary (XRB) populations, and their scaling relations, for a sample of 38 nearby galaxies (D = 3.4-29 Mpc). Our galaxy sample is drawn primarily from the Spitzer infrared nearby galaxy survey (SINGS), and contains a wealth of Chandra (5.8 Ms total) and multiwavelength data, allowing for star-formation rates (SFRs) and stellar masses (M*) to be measured on subgalactic scales. We divided the 2478 X-ray detected sources into 21 subsamples in bins of specific-SFR (sSFR = SFR/M*) and constructed XLFs. To model the XLF dependence on sSFR, we fit a global XLF model, containing contributions from high-mass XRBs (HMXBs), low-mass XRBs (LMXBs), and background sources from the cosmic X-ray background (CXB) that respectively scale with SFR, M*, and sky area. We find an HMXB XLF that is more complex in shape than previously reported and an LMXB XLF that likely varies with sSFR, potentially due to an age dependence. When applying our global model to XLF data for each individual galaxy, we discover a few galaxy XLFs that significantly deviate from our model beyond statistical scatter. Most notably, relatively low-metallicity galaxies have an excess of HMXBs above ~10^38 erg/s and elliptical galaxies that have relatively rich populations of globular clusters (GCs) show excesses of LMXBs compared to the global model. Additional modeling of how the XRB XLF depends on stellar age, metallicity, and GC specific frequency is required to sufficiently characterize the XLFs of galaxies., Comment: Accepted for publication in ApJS (37 pages, 16 figures, and 7 tables). Data products available (https://lehmer.uark.edu/downloads/)

The formation and evolution of double neutron stars (DNS) have traditionally been studied using binary population synthesis. In this work, we take an alternative approach by focusing only on the second supernova (SN) forming the DNS and the subsequent orbital decay and merger due to gravitational wave radiation. Using analytic and numerical methods, we explore how different NS natal kick velocity distributions, pre-SN orbital separations, and progenitor He-star masses affect the post-SN orbital periods, eccentricities, merger times, systemic velocities, and distances traveled by the system before merging. Comparison with the set of 17 known DNSs in the Milky Way shows that DNSs have pre-SN orbital separations ranging between 1 and 44 $R_{\odot}$. Those DNSs with pre-SN separations $\sim$1 $R_{\odot}$ have merger time distributions that peak $\sim$10-100 Myr after formation, regardless of the kick velocity received by the NS. These DNSs are typically formed with systemic velocities $\sim$10$^2$ km s$^{-1}$ and may travel $\sim$1-10 kpc before merging. Depending on progenitor mass of the second-born NS, the short merger time can account for the $r$-process enrichment observed in compact stellar systems such as ultra-faint dwarf galaxies. For Milky Way-mass galaxies only DNSs with the tightest pre-SN orbits and large kick velocities ($\gtrsim$10$^2$ km s$^{-1}$) can escape. However, those DNSs that do escape may travel as far as $\sim$Mpc before merging, which as previous studies have pointed out has implications for identifying the host galaxies to short gamma ray bursts and gravitational wave events., Comment: 16 pages, 10 figures, accepted for publication in MNRAS

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., Comment: 9 pages, 2 figures. Submitted to the Astro2020 Decadal Survey

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science., Comment: 50 pages, Probe class mission concept study report submitted to NASA for Astro2020 Decadal Survey

Ultraluminous X-ray (ULX) pulsars are a new class of object powered by apparent super-critical accretion onto magnetized neutron stars. Three sources in this class identified so far; M82 X-2, NGC 5907 ULX-1 and NGC 7793 P13, have been found to have two properties in common; $\sim1$-s spin periods, and for NGC 5907 ULX-1 and NGC 7793 P13 periodic X-ray flux modulations on timescales of $\sim60-80$ days. M82 X-2 resides in a crowded field that includes the ULX M82 X-1 separated from X-2 by 5", as well as other bright point sources. A 60-day modulation has been observed from the region but the origin has been difficult to identify; both M82 X-1 and X-2 have been suggested as the source. In this paper we present the analysis of a systematic monitoring campaign by Chandra, the only X-ray telescope capable of resolving the crowded field. From a simple Lomb-Scargle periodogram analysis and a more sophisticated Gaussian Process analysis we find that only X-2 exhibits a periodic signal around 60 days supporting previous claims that it is the origin. We also construct a phase-averaged flux profile of the modulations from higher cadence Swift/XRT data and find that the flux variations in the Chandra data are fully consistent with the flux profile. Since the orbit of the neutron star and its companion is known to be 2.5 days, the $\sim60$-day period must be super-orbital in origin. The flux of the modulations varies by a factor of $\sim$100 from minimum to maximum, with no evidence for spectral variations, making the origin difficult to explain., Comment: Accepted for publication in the Astrophysical Journal

We have compiled the most complete census of High-Mass X-ray Binaries (HMXBs) in the Small Magellanic Cloud with the aim to investigate the formation efficiency of young accreting binaries in its low metallicity environment. In total, we use 127 X-ray sources with detections in our \chandra X-ray Visionary Program (XVP), supplemented by 14 additional (likely and confirmed) HMXBs identified by \cite{2016A&A...586A..81H} that fall within the XVP area, but are not either detected in our survey (9 sources) or matched with any XVP source that has at least one OB counterpart in the OGLE-III catalog (5 sources). Specifically, we examine the number ratio of the HMXBs [N(HMXBs)] to {\it (a)} the number of OB stars, {\it (b)} the local star-formation rate (SFR), and {\it (c)} the stellar mass produced during the specific star-formation burst, all as a function of the age of their parent stellar populations. Each of these indicators serves a different role, but in all cases we find that the HMXB formation efficiency increases as a function of time (following a burst of star formation) up to $\sim$40--60\,Myr, and then gradually decreases. The peak formation efficiency N(HMXB)/SFR is (49 $\pm$ 14) $[10^{-5}~{\rm M_{\odot}/yr}]^{-1}$, in good agreement with previous estimates of the average formation efficiency in the broad $\sim$20--60\,Myr age range. The frequency of HMXBs is a factor of 8$\times$ higher than at $\sim$10\,Myr, and 4$\times$ higher than at $\sim$260\,Myr, i.e. at earlier and later epochs, respectively., Comment: 9 pages, 4 figures, submitted to ApJ, comments welcome

Nearby star-forming galaxies offer a unique environment to study the populations of young ($<$100 Myr) accreting binaries. These systems are tracers of past populations of massive stars that heavily affect their immediate environment and parent galaxies. Using a Chandra X-ray Visionary program, we investigate the young neutron-star binary population in the low metallicity of the Small Magellanic Cloud (SMC) by reaching quiescent X-ray luminosity levels ($\sim$few times $10^{32}$ erg/s). We present the first measurement of the formation efficiency of high-mass X-ray binaries (HMXBs) as a function of the age of their parent stellar populations by using 3 indicators: the number ratio of HMXBs to OB stars, to the SFR, and to the stellar mass produced during the specific star-formation burst they are associated with. In all cases, we find that the HMXB formation efficiency increases as a function of time up to $\sim$40-60 Myr, and then gradually decreases., Comment: 5 pages, 4 figures. To appear in Proc. IAUS 346: High-mass X-ray binaries: illuminating the passage from massive binaries to merging compact objects

Almost all confirmed optical counterparts of HMXBs in the SMC are OB stars with equatorial decretion disks (OBe). These sources emit strongly in Balmer lines and standout when imaged through narrow-band H{\alpha} imaging. The lack of secure counterparts for a significant fraction of the HMXBs motivated us to search for more. Using the catalogs for OB/OBe stars (Maravelias et al. 2017) and for HMXBs (Haberl & Sturm 2016) we detect 70 optical counterparts (out of 104 covered by our survey). We provide the first identification of the optical counterpart to the source XTEJ0050-731. We verify that 17 previously uncertain optical counterparts are indeed the proper matches. Regarding 52 confirmed HMXBs (known optical counterparts with H{\alpha} emission), we detect 39 as OBe and another 13 as OB stars. This allows a direct estimation of the fraction of active OBe stars in HMXBs that show H{\alpha} emission at a given epoch to be at least $\sim75\%$ of their total HMXB population., Comment: 3 pages, 1 figure; submitted to the proceedings of IAU Symposium 346 "High-mass X-ray binaries: illuminating the passage from massive binaries to merging compact objects" (Vienna, Aug 27 - 31, 2018) - this version: a couple of typos corrected, style of references updated