Showing 1–50 of 234 results for author: Zezas, A

Blueberry galaxies up to 200 Mpc and their optical and infrared properties

Authors: K. Kouroumpatzakis, J. Svoboda, A. Zezas, A. Borkar, E. Kyritsis, P. G. Boorman, C. Daoutis, B. Adamcová, R. Grossová

Abstract: 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… ▽ More 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~λ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. △ Less

Submitted 26 June, 2024; v1 submitted 6 May, 2024; originally announced May 2024.

Comments: 18 pages, 11 + 4 (appendix) figures. Accepted for publication in Astronomy & Astrophysics on 2024-05-02. Updated to match the published version

Report number: aa49766-24

The first all-sky survey of star-forming galaxies with eROSITA: Scaling relations and a population of X-ray luminous starbursts

Authors: E. Kyritsis, A. Zezas, F. Haberl, P. Weber, A. Basu-Zych, N. Vulic, C. Maitra, S. Hämmerich, J. Wilms, M. Sasaki, A. Hornschemeier, A. Ptak, A. Merloni, J. Comparat

Abstract: 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… ▽ More 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. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: Submitted to A&A

O Corona, where art thou? eROSITA's view of UV-optical-IR variability-selected massive black holes in low-mass galaxies

Authors: R. Arcodia, A. Merloni, J. Comparat, T. Dwelly, R. Seppi, Y. Zhang, J. Buchner, A. Georgakakis, F. Haberl, Z. Igo, E. Kyritsis, T. Liu, K. Nandra, Q. Ni, G. Ponti, M. Salvato, C. Ward, J. Wolf, A. Zezas

Abstract: Finding massive black holes (MBHs, $M_{BH}\approx10^4-10^7 M_{\odot}$) in the nuclei of low-mass galaxies ($M_{*}\lessapprox10^{10} M_{\odot}$) is crucial to constrain seeding and growth of black holes over cosmic time, but it is particularly challenging due to their low accretion luminosities. Variability selection via long-term photometric ultraviolet, optical, or infrared (UVOIR) light curves h… ▽ More Finding massive black holes (MBHs, $M_{BH}\approx10^4-10^7 M_{\odot}$) in the nuclei of low-mass galaxies ($M_{*}\lessapprox10^{10} M_{\odot}$) is crucial to constrain seeding and growth of black holes over cosmic time, but it is particularly challenging due to their low accretion luminosities. Variability selection via long-term photometric ultraviolet, optical, or infrared (UVOIR) light curves has proved effective and identifies lower-Eddington ratios compared to broad and narrow optical spectral lines searches. In the inefficient accretion regime, X-ray and radio searches are effective, but they have been limited to small samples. Therefore, differences between selection techniques have remained uncertain. Here, we present the first large systematic investigation of the X-ray properties of a sample of known MBH candidates in dwarf galaxies. We extracted X-ray photometry and spectra of a sample of $\sim200$ UVOIR variability-selected MBHs and significantly detected 17 of them in the deepest available \emph{SRG}/eROSITA image, of which four are newly discovered X-ray sources and two are new secure MBHs. This implies that tens to hundreds of LSST MBHs will have SRG/eROSITA counterparts, depending on the seeding model adopted. Surprisingly, the stacked X-ray images of the many non-detected MBHs are incompatible with standard disk-corona relations, typical of active galactic nuclei, inferred from both the optical and radio fluxes. They are instead compatible with the X-ray emission predicted for normal galaxies. After careful consideration of potential biases, we identified that this X-ray weakness needs a physical origin. A possibility is that a canonical X-ray corona might be lacking in the majority of this population of UVOIR-variability selected low-mass galaxies or that unusual accretion modes and spectral energy distributions are in place for MBHs in dwarf galaxies. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: Accepted for publication in A&A

The High Energy X-ray Probe: Resolved X-ray Populations in Extragalactic Environments

Authors: Bret D. Lehmer, Kristen Garofali, Breanna A. Binder, Francesca Fornasini, Neven Vulic, Andreas Zezas, Ann Hornschemeier, Margaret Lazzarini, Hannah Moon, Toni Venters, Daniel Wik, Mihoko Yukita, Matteo Bachetti, Javier A. García, Brian Grefenstette, Kristin Madsen, Kaya Mori, Daniel Stern

Abstract: 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 en… ▽ More 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. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: Accepted for publication in Frontiers in Astronomy and Space Sciences

A versatile classification tool for galactic activity using optical and infrared colors

Authors: Charalampos Daoutis, Elias Kyritsis, Konstantinos Kouroumpatzakis, Andreas Zezas

Abstract: (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,… ▽ More (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. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: 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

Journal ref: A&A 679, A76 (2023)

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

Authors: Lacey West, Kristen Garofali, Bret D. Lehmer, Andrea Prestwich, Rafael Eufrasio, Wasutep Luangtip, Timothy P. Roberts, Andreas Zezas

Abstract: 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 originat… ▽ More 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. △ Less

Submitted 30 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ (17 pages, 3 tables, 8 figures)

X-ray Binaries in External Galaxies

Authors: Marat Gilfanov, Giuseppina Fabbiano, Bret Lehmer, Andreas Zezas

Abstract: 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 h… ▽ More 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. △ Less

Submitted 27 April, 2023; originally announced April 2023.

Comments: Invited chapter for the Handbook of X-ray and Gamma-ray Astrophysics. Editors: Cosimo Bambi, Andrea Santangelo. Publisher: Springer Singapore, 2023

A versatile classification tool for galactic activity using optical and infrared colors

Authors: Elias Kyritsis, Charalampos Daoutis, Andreas Zezas, Konstantinos Kouroumpatzakis

Abstract: 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 a… ▽ More 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. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: 3 pages, 3 figures, the manuscript will be published on Memorie della SAIt

Star-formation rate and stellar mass calibrations based on infrared photometry and their dependence on stellar population age and extinction

Authors: Konstantinos Kouroumpatzakis, Andreas Zezas, Elias Kyritsis, Samir Salim, Jiri Svoboda

Abstract: 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… ▽ More 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. △ Less

Submitted 17 March, 2023; originally announced March 2023.

Comments: 18 pages, 10 figures. Accepted for publication in Astronomy & Astrophysics on 16 March 2023

Journal ref: A&A 673, A16 (2023)

X-ray luminosity function of high-mass X-ray binaries: Studying the signatures of different physical processes using detailed binary evolution calculations

Authors: Devina Misra, Konstantinos Kovlakas, Tassos Fragos, Margaret Lazzarini, Simone S. Bavera, Bret D. Lehmer, Andreas Zezas, Emmanouil Zapartas, Zepei Xing, Jeff J. Andrews, Aaron Dotter, Kyle A. Rocha, Philipp M. Srivastava, Meng Sun

Abstract: 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 variou… ▽ More 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. △ Less

Submitted 14 March, 2023; v1 submitted 12 September, 2022; originally announced September 2022.

Comments: 31 pages, 32 figures, Accepted by A&A. Fixed typos and updated references. Referee's comments were addressed

Journal ref: A&A 672, A99 (2023)

Identifying diffuse spatial structures in high-energy photon lists

Authors: Minjie Fan, Jue Wang, Vinay L. Kashyap, Thomas C. M. Lee, David A. van Dyk, Andreas Zezas

Abstract: 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 c… ▽ More 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. △ Less

Submitted 4 November, 2022; v1 submitted 15 August, 2022; originally announced August 2022.

Optical emission-line Luminosity Function models for populations of Supernova Remnants

Authors: M. Kopsacheili, A. Zezas, I. Leonidaki

Abstract: 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… ▽ More 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. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: 11 pages, Accepted for publication in MNRAS

Average bolometric corrections and optical to X-ray flux measurements as a function of accretion rate for X-ray binaries

Authors: K. Anastasopoulou, A. Zezas, J. F. Steiner, P. Reig

Abstract: 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 a… ▽ More 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. △ Less

Submitted 4 April, 2022; originally announced April 2022.

Comments: 16 pages, 6 figures, 14 tables, accepted for publication in MNRAS

Spatial Structures in the Globular Cluster Distribution of Fornax Cluster Galaxies

Authors: Raffaele D'Abrusco, David Zegeye, Giuseppina Fabbiano, Michele Cantiello, Maurizio Paolillo, Andreas Zezas

Abstract: 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… ▽ More 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. △ Less

Submitted 17 January, 2022; originally announced January 2022.

Comments: 40 pages, 26 figures, accepted for publication in ApJ

An expanded ultraluminous X-ray source catalogue

Authors: M. C. i Bernadich, A. D. Schwope, K. Kovlakas, A. Zezas, I. Traulsen

Abstract: 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-D… ▽ More 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. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: 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

Journal ref: A&A 659, A188 (2022)

A new automated tool for the spectral classification of OB stars

Authors: E. Kyritsis, G. Maravelias, A. Zezas, P. Bonfini, K. Kovlakas, P. Reig

Abstract: (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… ▽ More (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. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: 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

Journal ref: A&A 657, A62 (2022)

The Supernova Remnant Populations of the galaxies NGC 45, NGC 55, NGC 1313, NGC 7793: Luminosity and Excitation Functions

Authors: M. Kopsacheili, A. Zezas, I. Leonidaki, P. Boumis

Abstract: 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… ▽ More 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. △ Less

Submitted 17 August, 2021; originally announced August 2021.

Comments: Accepted in MNRAS. 19 pages, 21 figures

Evidence for low kick velocities among high-mass X-ray binaries in the Small Magellanic Cloud from the spatial correlation function

Authors: Arash Bodaghee, Vallia Antoniou, Andreas Zezas, John A. Tomsick, Zachary Jordan, Eric Frechette, Brenton Jackson, Ryan Agnew, Ann E. Hornschemeier, Jerome Rodriguez

Abstract: We present the two-point cross-correlation function between high-mass X-ray binaries (HMXBs) in the Small Magellanic Cloud (SMC) and their likely birthplaces (OB Associations: OBAs). This function compares the spatial correlation between the observed HMXB and OBA populations against mock catalogs in which the members are distributed randomly across the sky. A significant correlation (15 sigma) is… ▽ More 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. △ Less

Submitted 7 July, 2021; v1 submitted 6 July, 2021; originally announced July 2021.

Comments: 10 pages, 7 figures, accepted for publication in the Astrophysical Journal

The Star Formation Reference Survey-V: the effect of extinction, stellar mass, metallicity, and nuclear activity on star-formation rates based on H$α$ emission

Authors: Konstantinos Kouroumpatzakis, Andreas Zezas, Alexandros Maragkoudakis, Steven P. Willner, Paolo Bonfini, Matthew L. N. Ashby, Paul H. Sell, Thomas H. Jarrett

Abstract: We present new H$α$ 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$α$-based star-formation rates (SFRs) with and without correction for the contribution of [$\rm N_{^{II}}$] emission. We consider the effect of extinction c… ▽ More We present new H$α$ 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$α$-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 $μ$m + H$α$, 8 $μ$m + H$α$, FIR + FUV, and H$α$ 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$α$ emission, and we measure a dependence of the radio-to-H$α$ 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. △ Less

Submitted 2 July, 2021; originally announced July 2021.

Comments: 20 pages, 21 Figures, accepted for publication in MNRAS

The Star Formation Reference Survey. IV. Stellar mass distribution of local star-forming galaxies

Authors: Paolo Bonfini, Andreas Zezas, Matthew L. N. Ashby, Steven P. Willner, Alexandros Maragkoudakis, Konstantinos Kouroumpatzakis, Paul. H. Sell, Konstantinos Kovlakas

Abstract: 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 image… ▽ More 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 $ρ_{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. △ Less

Submitted 2 July, 2021; originally announced July 2021.

Comments: 33 pages, 21 figures, published in MNRAS, 504, 3831-3861 (2021)

Journal ref: MNRAS, 504, 3831-3861 (2021)

The eROSITA Final Equatorial-Depth Survey (eFEDS): Presenting The Demographics of X-ray Emission From Normal Galaxies

Authors: N. Vulic, A. E. Hornschemeier, F. Haberl, A. R. Basu-Zych, E. Kyritsis, A. Zezas, M. Salvato, A. Ptak, A. Bogdan, K. Kovlakas, J. Wilms, M. Sasaki, T. Liu, A. Merloni, T. Dwelly, H. Brunner, G. Lamer, C. Maitra, K. Nandra, A. Santangelo

Abstract: 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 provide… ▽ More 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. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Comments: 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

Journal ref: A&A 661, A16 (2022)

The Heraklion Extragalactic Catalogue (HECATE): a value-added galaxy catalogue for multi-messenger astrophysics

Authors: Konstantinos Kovlakas, Andreas Zezas, Jeff J. Andrews, Antara Basu-Zych, Tassos Fragos, Ann Hornschemeier, Konstantinos Kouroumpatzakis, Bret Lehmer, Andrew Ptak

Abstract: 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 an… ▽ More 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. △ Less

Submitted 22 June, 2021; originally announced June 2021.

Comments: 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

eBASCS: Disentangling Overlapping Astronomical Sources II, using Spatial, Spectral, and Temporal Information

Authors: Antoine D. Meyer, David A. van Dyk, Vinay L. Kashyap, Luis F. Campos, David E. Jones, Aneta Siemiginowska, Andreas Zezas

Abstract: 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-… ▽ More 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. △ Less

Submitted 18 May, 2021; originally announced May 2021.

Change point detection and image segmentation for time series of astrophysical images

Authors: Cong Xu, Hans Moritz Günther, Vinay L. Kashyap, Thomas C. M. Lee, Andreas Zezas

Abstract: 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)… ▽ More 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. △ Less

Submitted 26 January, 2021; originally announced January 2021.

Comments: 22 pages, 10 figures

Targeted modeling of GW150914's binary black hole source with dart_board

Authors: Jeff J. Andrews, Julianne Cronin, Vicky Kalogera, Christopher Berry, Andreas Zezas

Abstract: 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… ▽ More 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. △ Less

Submitted 16 May, 2021; v1 submitted 27 November, 2020; originally announced November 2020.

Comments: 14 pages, 5 figures, 1 table, accepted for publication in ApJL

The Metallicity Dependence of the High-Mass X-ray Binary Luminosity Function

Authors: Bret D. Lehmer, Rafael T. Eufrasio, Antara Basu-Zych, Keith Doore, Tassos Fragos, Kristen Garofali, Konstantinos Kovlakas, Benjamin J. Williams, Andreas Zezas, Luidhy Santana-Silva

Abstract: 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… ▽ More 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. △ Less

Submitted 18 November, 2020; originally announced November 2020.

Comments: Accepted for publication in ApJ

Metallicity and X-ray luminosity variations in NGC 922

Authors: Konstantinos Kouroumpatzakis, Andreas Zezas, Anna Wolter, Antonella Fruscione, Konstantina Anastasopoulou, Andrea Prestwich

Abstract: 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… ▽ More 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. △ Less

Submitted 24 October, 2020; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: 15 pages, 10 figures

A census of ultraluminous X-ray sources in the local Universe

Authors: Konstantinos Kovlakas, Andreas Zezas, Jeff J. Andrews, Antara Basu-Zych, Tassos Fragos, Ann Hornschemeier, Bret Lehmer, Andrew Ptak

Abstract: 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 w… ▽ More 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. △ Less

Submitted 24 August, 2020; originally announced August 2020.

Comments: 21 pages, 19 figures, 6 tables. Accepted for publication in MNRAS

The Next Generation X-ray Galaxy Survey with eROSITA

Authors: Antara R. Basu-Zych, Ann E. Hornschemeier, Frank Haberl, Neven Vulic, Jorn Wilms, Andreas Zezas, Konstantinos Kovlakas, Andrew Ptak, Thomas Dauser

Abstract: 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 scientifi… ▽ More 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 $σ$ significance) at 50-200 Mpc, which is ~100X more normal galaxies than detected in any X-ray survey to date. △ Less

Submitted 4 August, 2020; originally announced August 2020.

Comments: 18 pages, 15 figures. Accepted for publication in MNRAS

Constraining X-ray reflection in the low-luminosity AGN NGC 3718 using NuSTAR and XMM--Newton

Authors: Y. Diaz, P. Arévalo, L. Hernández-García, L. Bassani, A. Malizia, O. González-Martín, C. Ricci, G. Matt, D. Stern, D. May, A. Zezas, F. E. Bauer

Abstract: 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 st… ▽ More 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. △ Less

Submitted 16 June, 2020; originally announced June 2020.

Comments: 15 pages and 16 Figures

Sub-galactic scaling relations between X-ray luminosity, star-formation rate, and stellar mass

Authors: Konstantinos Kouroumpatzakis, Andreas Zezas, Paul Sell, Konstantinos Kovlakas, Paolo Bonfini, Steven Willner, Matthew Ashby, Alexandros Maragkoudakis, Thomas Jarrett

Abstract: 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… ▽ More 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$α$ shows the tightest correlation with X-ray luminosity because H$α$ emission probes stellar populations with ages similar to HMXB formation timescales, but the H$α$-based SFR is reliable only for $\rm sSFR{>}10^{-12}\,M_\odot \, yr^{-1}/M_\odot$. △ Less

Submitted 21 April, 2020; originally announced April 2020.

Comments: 19 pages, 13 figures, accepted for publication in MNRAS (2020-04-19)

Neutron Stars and Black Holes in the Small Magellanic Cloud: The SMC NuSTAR Legacy Survey

Authors: M. Lazzarini, B. F. Williams, A. E. Hornschemeier, V. Antoniou, G. Vasilopoulos, F. Haberl, N. Vulic, M. Yukita, A. Zezas, A. Bodaghee, B. D. Lehmer, T. J. Maccarone, A. Ptak, D. Wik, F. M. Fornasini, Jaesub Hong, J. A. Kennea, J. A. Tomsick, T. Venters, A. Udalski, A. Cassity

Abstract: 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$σ$ significance level (4… ▽ More 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$σ$ significance level (4$-$25 keV) and obtained upper limits on an additional 40 sources. We reached a 3$σ$ 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$σ$ 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. △ Less

Submitted 24 September, 2019; originally announced September 2019.

Comments: 8 tables, 18 figures

Discrete star formation events in the central Bar of the Small Magellanic Cloud

Authors: A. Strantzalis, D. Hatzidimitriou, A. Zezas, V. Antoniou, S. Lianou, S Tsilia

Abstract: 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.… ▽ More 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. △ Less

Submitted 19 September, 2019; originally announced September 2019.

Comments: 12 pages, 11 figures; accepted for publication in MNRAS

A diagnostic tool for the identification of Supernova Remnants

Authors: M. Kopsacheili, A. Zezas, I. Leonidaki

Abstract: We present new diagnostic tools for distinguishing supernova remnants (SNRs) from HII regions. Up to now, sources with flux ratio [S II]/H$\rmα$ 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 p… ▽ More We present new diagnostic tools for distinguishing supernova remnants (SNRs) from HII regions. Up to now, sources with flux ratio [S II]/H$\rmα$ 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α$ - [O II]/H$\rmβ$ - [O III]/H$\rmβ$. This method gives $98.95\% $ completeness in the SNR selection and $1.20\%$ contamination. We also define the [O I]/H$\rmα$ SNR selection criterion and we measure its efficiency in comparison to other selection criteria. △ Less

Submitted 12 September, 2019; originally announced September 2019.

Comments: Accepted for publication in the Monthly Notices of the Royal Astronomical Society, 18 pages, 30 figures

The Complete Evolution of a Neutron-Star Binary through a Common Envelope Phase Using 1D Hydrodynamic Simulations

Authors: Tassos Fragos, Jeff J. Andrews, Enrico Ramirez-Ruiz, Georges Meynet, Vicky Kalogera, Ronald E. Taam, Andreas Zezas

Abstract: 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 mode… ▽ More 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 $α_{\rm CE}$-equivalent efficiency of $\approx$5. △ Less

Submitted 29 July, 2019; originally announced July 2019.

Comments: 10 pages, 5 figures, submitted to The Astrophysical Journal Letters. Comments welcome

X-ray Binary Luminosity Function Scaling Relations for Local Galaxies Based on Subgalactic Modeling

Authors: Bret D. Lehmer, Rafael T. Eufrasio, Panayiotis Tzanavaris, Antara Basu-Zych, Tassos Fragos, Andrea Prestwich, Mihoko Yukita, Andreas Zezas, Ann E. Hornschemeier, Andrew Ptak

Abstract: 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)… ▽ More 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. △ Less

Submitted 13 May, 2019; originally announced May 2019.

Comments: Accepted for publication in ApJS (37 pages, 16 figures, and 7 tables). Data products available (https://lehmer.uark.edu/downloads/)

Double Neutron Star Formation: Merger Times, Systemic Velocities, and Travel Distances

Authors: Jeff J. Andrews, Andreas Zezas

Abstract: 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… ▽ More 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. △ Less

Submitted 12 April, 2019; originally announced April 2019.

Comments: 16 pages, 10 figures, accepted for publication in MNRAS

Astro 2020 Science White Paper: Time Domain Studies of Neutron Star and Black Hole Populations: X-ray Identification of Compact Object Types

Authors: N. Vulic, A. E. Hornschemeier, V. Antoniou, A. R. Basu-Zych, B. Binder, F. M. Fornasini, F. Furst, F. Haberl, M. Heida, B. D. Lehmer, T. J. Maccarone, A. F. Ptak, G. R. Sivakoff, P. Tzanavaris, D. R. Wik, B. F. Williams, J. Wilms, M. Yukita, A. Zezas

Abstract: 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, a… ▽ More 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. △ Less

Submitted 23 March, 2019; originally announced March 2019.

Comments: 9 pages, 2 figures. Submitted to the Astro2020 Decadal Survey

STROBE-X: X-ray Timing and Spectroscopy on Dynamical Timescales from Microseconds to Years

Authors: Paul S. Ray, Zaven Arzoumanian, David Ballantyne, Enrico Bozzo, Soren Brandt, Laura Brenneman, Deepto Chakrabarty, Marc Christophersen, Alessandra DeRosa, Marco Feroci, Keith Gendreau, Adam Goldstein, Dieter Hartmann, Margarita Hernanz, Peter Jenke, Erin Kara, Tom Maccarone, Michael McDonald, Michael Nowak, Bernard Phlips, Ron Remillard, Abigail Stevens, John Tomsick, Anna Watts, Colleen Wilson-Hodge , et al. (134 additional authors not shown)

Abstract: 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 techniqu… ▽ More 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. △ Less

Submitted 8 March, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: 50 pages, Probe class mission concept study report submitted to NASA for Astro2020 Decadal Survey

A ~60-day super-orbital period originating from the ultraluminous X-ray pulsar in M82

Authors: Murray Brightman, Fiona A. Harrison, Matteo Bachetti, Yanjun Xu, Felix Fürst, Dominic J. Walton, Andrew Ptak, Mihoko Yukita, Andreas Zezas

Abstract: 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$… ▽ More 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. △ Less

Submitted 29 January, 2019; originally announced January 2019.

Comments: Accepted for publication in the Astrophysical Journal

Deep Chandra survey of the Small Magellanic Cloud. III. Formation efficiency of High-Mass X-ray binaries

Authors: Vallia Antoniou, Andreas Zezas, Jeremy J. Drake, Carles Badenes, Frank Haberl, Nicholas J. Wright, Jaesub Hong, Rosanne Di Stefano, Terrance J. Gaetz, Knox S. Long, Paul P. Plucinsky, Manami Sasaki, Benjamin F. Williams, P. Frank Winkler

Abstract: 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 \cit… ▽ More 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. △ Less

Submitted 4 January, 2019; originally announced January 2019.

Comments: 9 pages, 4 figures, submitted to ApJ, comments welcome

Different generations of HMXBs: clues about their formation efficiency from Magellanic Clouds studies

Authors: Vallia Antoniou, Andreas Zezas, Jeremy J. Drake, Carles Badenes, Frank Haberl, Jaesub Hong, Paul P. Plucinsky, the SMC XVP Collaboration Team

Abstract: 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 C… ▽ More 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. △ Less

Submitted 28 December, 2018; originally announced December 2018.

Comments: 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

Clarifying the population of HMXBs in the Small Magellanic Cloud

Authors: Grigoris Maravelias, Andreas Zezas, Vallia Antoniou, Despina Hatzidimitriou, Frank Haberl

Abstract: 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α 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 (Habe… ▽ More 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α 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α 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α emission at a given epoch to be at least $\sim75\%$ of their total HMXB population. △ Less

Submitted 21 December, 2018; v1 submitted 27 November, 2018; originally announced November 2018.

Comments: 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

Journal ref: Proc. IAU 14 (2018) 350-352

Do sub-galactic regions follow the galaxy-wide X-ray scaling relations? The example of NGC 3310 and NGC 2276

Authors: K. Anastasopoulou, A. Zezas, V. Gkiokas, K. Kovlakas

Abstract: We present results from Chandra observations of the X-ray starburst galaxies NGC 3310 and NGC 2276. We detect 27 discrete sources in NGC 3310, and 19 discrete sources in NGC 2276 with luminosities above $\mathrm{1.0\times 10^{38}\ erg\ s^{-1}}$. The majority of the sources have photon indices of 1.7-2.0, typical for X-ray binaries. Both galaxies have large numbers of ultra-luminous X-ray sources (… ▽ More We present results from Chandra observations of the X-ray starburst galaxies NGC 3310 and NGC 2276. We detect 27 discrete sources in NGC 3310, and 19 discrete sources in NGC 2276 with luminosities above $\mathrm{1.0\times 10^{38}\ erg\ s^{-1}}$. The majority of the sources have photon indices of 1.7-2.0, typical for X-ray binaries. Both galaxies have large numbers of ultra-luminous X-ray sources (ULXs; sources with $\mathrm{L(0.3-10.0\ keV)>10^{39}\ erg\ s^{-1}}$), 14 for NGC 3310 concentrated on the circumnuclear star-forming ring and north spiral arm and 11 for NGC 2276 with the brighter ones on the west side of the galaxy which is compressed due to harassment by the intra-group medium it is moving into. We find for both galaxies that the ULX-hosting areas are located above the general Lx-SFR scaling relations while other areas either follow or fall below the scaling relations. This indicates that sub-galactic regions follow the galaxy-wide scaling relations but with much larger scatter resulting from the age (and possibly metallicity) of their local stellar populations in agreement with recent theoretical and observational results. Such differences in age could be the origin of the scatter we observe in the low SFR regime in the Lx-SFR scaling relations. △ Less

Submitted 15 November, 2018; originally announced November 2018.

Comments: 25 pages, 15 figures, 12 tables. Accepted for publication in MNRAS

On the Nature of the X-ray Emission from the Ultraluminous X-ray Source, M33 X-8: New Constraints from NuSTAR and XMM-Newton

Authors: Lacey A. West, Bret D. Lehmer, Daniel Wik, Jun Yang, Dominic J. Walton, Vallia Antoniou, Frank Haberl, Ann Hornschemeier, Thomas J. Maccarone, Paul P. Plucinsky, Andrew Ptak, Benjamin F. Williams, Neven Vulic, Mihoko Yukita, Andreas Zezas

Abstract: 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 k… ▽ More 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. △ Less

Submitted 26 October, 2018; originally announced October 2018.

Comments: Accepted for publication in ApJ (15 pages, 4 tables, 6 figures)

Search for AGN counterparts of unidentified Fermi-LAT sources with optical polarimetry: Demonstration of the technique

Authors: N. Mandarakas, D. Blinov, I. Liodakis, K. Kouroumpatzakis, A. Zezas, G. V. Panopoulou, I. Myserlis, E. Angelakis, T. Hovatta, S. Kiehlmann, K. Kokolakis, E. Paleologou, A. Pouliasi, R. Skalidis, V. Pavlidou

Abstract: The third Fermi-LAT catalog (3FGL) presented the data of the first four years of observations from the Fermi Gamma-ray Space Telescope mission. There are 3034 sources, 1010 of which still remain unidentified. Identifying and classifying gamma-ray emitters is of high significance with regard to studying high-energy astrophysics. We demonstrate that optical polarimetry can be an advantageous and pra… ▽ More The third Fermi-LAT catalog (3FGL) presented the data of the first four years of observations from the Fermi Gamma-ray Space Telescope mission. There are 3034 sources, 1010 of which still remain unidentified. Identifying and classifying gamma-ray emitters is of high significance with regard to studying high-energy astrophysics. We demonstrate that optical polarimetry can be an advantageous and practical tool in the hunt for counterparts of the unidentified gamma-ray sources (UGSs). Using data from the RoboPol project, we validated that a significant fraction of active galactic nuclei (AGN) associated with 3FGL sources can be identified due to their high optical polarization exceeding that of the field stars. We performed an optical polarimetric survey within $3σ$ uncertainties of four unidentified 3FGL sources. We discovered a previously unknown extragalactic object within the positional uncertainty of 3FGL J0221.2+2518. We obtained its spectrum and measured a redshift of $z=0.0609\pm0.0004$. Using these measurements and archival data we demonstrate that this source is a candidate counterpart for 3FGL J0221.2+2518 and most probably is a composite object: a star-forming galaxy accompanied by AGN. We conclude that polarimetry can be a powerful asset in the search for AGN candidate counterparts for unidentified Fermi sources. Future extensive polarimetric surveys at high galactic latitudes (e.g., PASIPHAE) will allow the association of a significant fraction of currently unidentified gamma-ray sources. △ Less

Submitted 12 February, 2019; v1 submitted 15 October, 2018; originally announced October 2018.

Comments: accepted to A&A

Journal ref: A&A 623, A61 (2019)

The integrated properties of the CALIFA galaxies: Model-derived galaxy parameters and quenching of star formation

Authors: T. Bitsakis, S. F. Sanchez, L. Ciesla, P. Bonfini, V. Charmandaris, B. Cervantes Sodi, A. Maragkoudakis, T. Diaz-Santos, A. Zezas

Abstract: We present a study of the integrated properties of the 835 galaxies in the CALIFA survey. To derive the main physical parameters of the galaxies we have fitted their UV-to-IR spectral energy distributions (SED) with sets of theoretical models using CIGALE. We perform a comparison of the integrated galaxy parameters derived from multi-band SED fitting with those obtained from modelling the Integral… ▽ More We present a study of the integrated properties of the 835 galaxies in the CALIFA survey. To derive the main physical parameters of the galaxies we have fitted their UV-to-IR spectral energy distributions (SED) with sets of theoretical models using CIGALE. We perform a comparison of the integrated galaxy parameters derived from multi-band SED fitting with those obtained from modelling the Integral Field Unit (IFU) spectra and show the clear advantage of using the SED-derived star formation rates (SFR). A detailed analysis of galaxies in the SFR/Mstar plane as a function of their properties reveals that quenching of star formation is caused by a combination of gas deficiency and the inefficiency of the existing gas to form new stars. Exploring the plausible mechanisms that could produce this effect, we find a strong correlation with galaxy morphology and the build-up of central bulge. On the other hand, the presence of AGN and/or a stellar bar, as well as the local environment have only temporal effects on the current star formation, a result also consistent with their model-derived star formation histories. △ Less

Submitted 12 October, 2018; originally announced October 2018.

Comments: Accepted for publication in MNRAS

The Star Formation Reference Survey III: A Multi-wavelength View of Star Formation in Nearby Galaxies

Authors: Smriti Mahajan, M. L. N. Ashby, S. P. Willner, P. Barmby, G. G. Fazio, A. Maragkoudakis, S. Raychaudhury, A. Zezas

Abstract: We present multi-wavelength global star formation rate (SFR) estimates for 326 galaxies from the Star Formation Reference Survey (SFRS) in order to determine the mutual scatter and range of validity of different indicators. The widely used empirical SFR recipes based on 1.4 GHz continuum, 8.0 $μ$m polycyclic aromatic hydrocarbons (PAH), and a combination of far-infrared (FIR) plus ultraviolet (UV)… ▽ More We present multi-wavelength global star formation rate (SFR) estimates for 326 galaxies from the Star Formation Reference Survey (SFRS) in order to determine the mutual scatter and range of validity of different indicators. The widely used empirical SFR recipes based on 1.4 GHz continuum, 8.0 $μ$m polycyclic aromatic hydrocarbons (PAH), and a combination of far-infrared (FIR) plus ultraviolet (UV) emission are mutually consistent with scatter of $\raise{-0.8ex}\stackrel{\textstyle <}{\sim }$0.3 dex. The scatter is even smaller, $\raise{-0.8ex}\stackrel{\textstyle <}{\sim }$0.24 dex, in the intermediate luminosity range 9.3<log(L(60 $μ$m/L$_\odot$)<10.7. The data prefer a non-linear relation between 1.4 GHz luminosity and other SFR measures. PAH luminosity underestimates SFR for galaxies with strong UV emission. A bolometric extinction correction to far-ultraviolet luminosity yields SFR within 0.2 dex of the total SFR estimate, but extinction corrections based on UV spectral slope or nuclear Balmer decrement give SFRs that may differ from the total SFR by up to 2 dex. However, for the minority of galaxies with UV luminosity ${>}5\times10^9$ L$_{\odot}$ or with implied far-UV extinction <1 mag, the UV spectral slope gives extinction corrections with 0.22~dex uncertainty. △ Less

Submitted 19 October, 2018; v1 submitted 2 October, 2018; originally announced October 2018.

Comments: To be published in Monthly Notices of the Royal Astronomical Society. Data tables available at https://www.cfa.harvard.edu/irac/SFRS/ . V2 has minor clarifications from V1; no changes in data or results

Black Holes and Neutron Stars in Nearby Galaxies: Insights from NuSTAR

Authors: Neven Vulic, Ann E. Hornschemeier, Daniel R. Wik, Mihoko Yukita, Andreas Zezas, Andrew F. Ptak, Bret D. Lehmer, Valsamo Antoniou, Thomas J. Maccarone, Benjamin F. Williams, Francesca M. Fornasini

Abstract: 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 samp… ▽ More 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}}=αM_{\star}+β\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 $α=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. △ Less

Submitted 16 August, 2018; originally announced August 2018.

Comments: 38 pages, 12 figures, 8 tables. ApJ, in press

Low-luminosity AGN and X-ray binary populations in COSMOS star-forming galaxies

Authors: Francesca M. Fornasini, Francesca Civano, Giuseppina Fabbiano, Martin Elvis, Stefano Marchesi, Takamitsu Miyaji, Andreas Zezas

Abstract: We present an X-ray stacking analysis of $\sim$75,000 star-forming galaxies between $0.1<z<5.0$ using the Chandra COSMOS Legacy survey to study the X-ray emission of low-luminosity active galactic nuclei (AGN) and its connection to host galaxy properties. The stacks at $z<0.9$ have luminosity limits as low as $10^{40}-10^{41}$ erg s$^{-1}$, a regime in which X-ray binaries (XRBs) can dominate the… ▽ More We present an X-ray stacking analysis of $\sim$75,000 star-forming galaxies between $0.1<z<5.0$ using the Chandra COSMOS Legacy survey to study the X-ray emission of low-luminosity active galactic nuclei (AGN) and its connection to host galaxy properties. The stacks at $z<0.9$ have luminosity limits as low as $10^{40}-10^{41}$ erg s$^{-1}$, a regime in which X-ray binaries (XRBs) can dominate the X-ray emission. Comparing the measured luminosities to established XRB scaling relations, we find that the redshift evolution of the luminosity per star formation rate (SFR) of XRBs depends sensitively on the assumed obscuration and may be weaker than previously found. The XRB scaling relation based on stacks from the Chandra Deep Field South overestimates the XRB contribution to the COSMOS high specific SFR (sSFR) stacks, possibly due to a bias affecting the CDF-S stacks because of their small galaxy samples. After subtracting the estimated XRB contribution from the stacks, we find that most stacks at $z>1.3$ exhibit a significant X-ray excess indicating nuclear emission. The AGN emission is strongly correlated with stellar mass but does not exhibit an additional correlation with SFR. The hardness ratios of the high-redshift stacks indicate that the AGN are substantially obscured ($N_{\mathrm{H}}\sim10^{23}$ cm$^{-2}$). These obscured AGN are not identified by IRAC color selection and have $L_X\sim10^{41}-10^{43}$ erg s$^{-1}$, consistent with accretion at an Eddington rate of $\sim10^{-3}$ onto $10^7-10^8 M_{\odot}$ black holes. Combining our results with other X-ray studies suggests that AGN obscuration depends on stellar mass and an additional variable, possibly the Eddington rate. △ Less

Submitted 2 August, 2018; originally announced August 2018.

Comments: 24 pages, 12 figures, 3 tables (2 tables available in full online or by request)