Search

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., 21 pages, 19 figures, 6 tables. Accepted for publication in MNRAS

Modern and future surveys effectively provide a panchromatic view for large numbers of extragalactic objects. Consistently modeling these multiwavelength survey data is a critical but challenging task for extragalactic studies. The Code Investigating GALaxy Emission (cigale) is an efficient python code for spectral energy distribution (SED) fitting of galaxies and active galactic nuclei (AGNs). Recently, a major extension of cigale (named x-cigale) has been developed to account for AGN/galaxy X-ray emission and improve AGN modeling at UV-to-IR wavelengths. Here, we apply x-cigale to different samples, including Cosmological Evolution Survey (COSMOS) spectroscopic type 2 AGNs, Chandra Deep Field-South X-ray detected normal galaxies, Sloan Digital Sky Survey quasars, and COSMOS radio objects. From these tests, we identify several weaknesses of x-cigale and improve the code accordingly. These improvements are mainly related to AGN intrinsic X-ray anisotropy, X-ray binary emission, AGN accretion-disk SED shape, and AGN radio emission. These updates improve the fit quality and allow for new interpretation of the results, based on which we discuss physical implications. For example, we find that AGN intrinsic X-ray anisotropy is moderate, and can be modeled as ${L}_{X}(\theta )\propto 1+\cos \theta $, where θ is the viewing angle measured from the AGN axis. We merge the new code into the major branch of cigale, and publicly release this new version as cigale v2022.0 on https://cigale.lam.fr., G.Y. and C.P. acknowledge support from the George P. and Cynthia Woods Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, from the National Science Foundation through grants AST1614668 and AST-2009442, and from the NASA/ESA/CSA James Webb Space Telescope through the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS5-03127. M.B. acknowledges support by the ANID BASAL project FB21000 and the FONDECYT regular grants 1170618 and 1211000. W.N.B. acknowledges the support from NASA grant 80NSSC19K0961. K.M. has been supported by the Polish National Science Centre (UMO2018/30/E/ST9/00082). G.M. acknowledges support by the Agencia Estatal de Investigación, Unidad de Excelencia María de Maeztu, ref. MDM-2017-0765. M.S. acknowledges support by the Ministry of Education, Science and Technological Development of the Republic of Serbia through contract No. 451-03-9/2021-14/200002 and by the Science Fund of the Republic of Serbia, PROMIS 6060916, BOWIE.

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

High-energy emission associated with star formation has been proposed as a significant source of interstellar medium (ISM) ionization in low-metallicity starbursts and an important contributor to the heating of the intergalactic medium (IGM) in the high-redshift ($z > 8$) Universe. Using Chandra observations of a sample of 30 galaxies at $D \approx$~200--450 Mpc that have high specific star-formation rates of 3--9 Gyr$^{-1}$ and metallicities near $Z \approx 0.3 Z_\odot$, we provide new measurements of the average 0.5--8 keV spectral shape and normalization per unit star-formation rate (SFR). We model the sample-combined X-ray spectrum as a combination of hot gas and high-mass X-ray binary (HMXB) populations and constrain their relative contributions. We derive scaling relations of $\log L_{\rm 0.5-8 keV}^{\rm HMXB}$/SFR $= 40.19 \pm 0.06$ and $\log L_{\rm 0.5-2 keV}^{\rm gas}$/SFR $= 39.58^{+0.17}_{-0.28}$; significantly elevated compared to local relations. The HMXB scaling is also somewhat higher than $L_{\rm 0.5-8 keV}^{\rm HMXB}$-SFR-$Z$ relations presented in the literature, potentially due to our galaxies having relatively low HMXB obscuration and young and X-ray luminous stellar populations. The elevation of the hot gas scaling relation is at the level expected for diminished attenuation due to a reduction of metals; however, we cannot conclude that an $L_{\rm 0.5-2 keV}^{\rm gas}$-SFR-$Z$ relation is driven solely by changes in ISM metal content. Finally, we present SFR-scaled spectral models (both emergent and intrinsic) that span the X-ray--to--IR band, providing new benchmarks for studies of the impact of ISM ionization and IGM heating in the early Universe., Accepted for publication in ApJ (25 pages, 9 figures, and 4 tables)

The Chandra Deep Field-South and North surveys (CDFs) provide unique windows into the cosmic history of X-ray emission from normal (nonactive) galaxies. Scaling relations of normal-galaxy X-ray luminosity (L X) with star formation rate (SFR) and stellar mass (M ⋆) have been used to show that the formation rates of low-mass and high-mass X-ray binaries (LMXBs and HMXBs, respectively) evolve with redshift across z ≈ 0–2 following L HMXB/SFR ∝ (1 + z) and L LMXB/M ⋆ ∝ (1 + z)2−3. However, these measurements alone do not directly reveal the physical mechanisms behind the redshift evolution of X-ray binaries (XRBs). We derive star formation histories for a sample of 344 normal galaxies in the CDFs, using spectral energy distribution (SED) fitting of FUV-to-FIR photometric data, and construct a self-consistent, age-dependent model of the X-ray emission from the galaxies. Our model quantifies how X-ray emission from hot gas and XRB populations vary as functions of host stellar-population age. We find that (1) the ratio L X/M ⋆ declines by a factor of ∼1000 from 0 to 10 Gyr and (2) the X-ray SED becomes harder with increasing age, consistent with a scenario in which the hot gas contribution to the X-ray SED declines quickly for ages above 10 Myr. When dividing our sample into subsets based on metallicity, we find some indication that L X/M ⋆ is elevated for low-metallicity galaxies, consistent with recent studies of X-ray scaling relations. However, additional statistical constraints are required to quantify both the age and metallicity dependence of X-ray emission from star-forming galaxies.

We develop and implement an inclination-dependent attenuation prescription for spectral energy distribution (SED) fitting and study its impact on derived star-formation histories. We apply our prescription within the SED fitting code Lightning to a clean sample of 82, z=0.21-1.35 disk-dominated galaxies in the Great Observatories Origins Deep Survey North and South fields. To compare our inclination-dependent attenuation prescription with more traditional fitting prescriptions, we also fit the SEDs with the inclination-independent Calzetti et al. (2000) attenuation curve. From this comparison, we find that fits to a subset of 58, z < 0.7 galaxies in our sample, utilizing the Calzetti et al. (2000) prescription, recover similar trends with inclination as the inclination-dependent fits for the far-UV-band attenuation and recent star-formation rates. However, we find a difference between prescriptions in the optical attenuation (AV) that is strongly correlated with inclination (p-value < 10^-11). For more face-on galaxies, with i < 50 deg, (edge-on, i = 90 deg), the average derived AV is 0.31 +\- 0.11 magnitudes lower (0.56 +\- 0.16 magnitudes higher) for the inclination-dependent model compared to traditional methods. Further, the ratio of stellar masses between prescriptions also has a significant (p-value < 10^-2) trend with inclination. For i = 0-65 deg, stellar masses are systematically consistent between fits, with log(Mstar_inc/Mstar_Calzetti) = -0.05 +/- 0.03 dex and scatter of 0.11 dex. However, for i = 80-90 deg, derived stellar masses are lower for the Calzetti et al. (2000) fits by an average factor of 0.17 +\- 0.03 dex and scatter of 0.13 dex. Therefore, these results suggest that SED fitting assuming the Calzetti et al. (2000) attenuation law potentially underestimates stellar masses in highly inclined disk-dominated galaxies., 38 pages, 22 figures. 30 page body and 8 page appendix. Updated with published version

We examine possible environmental sources of the enhanced star formation and active galactic nucleus (AGN) activity in the $z = 3.09$ SSA22 protocluster using Hubble WFC3 F160W ($\sim1.6\ \rm \mu m$) observations of the SSA22 field, including new observations centered on eight X-ray selected protocluster AGN. To investigate the role of mergers in the observed AGN and star formation enhancement, we apply both quantitative (S\'ersic-fit and Gini-$M_{20}$) and visual morphological classifications to F160W images of protocluster Lyman break galaxies (LBGs) in the fields of the X-ray AGN and $z \sim 3$ field LBGs in SSA22 and GOODS-N. We find no statistically significant differences between the morphologies and merger fractions of protocluster and field LBGs, though we are limited by small number statistics in the protocluster. We also fit the UV-to-near-IR spectral energy distributions (SED) of F160W-detected protocluster and field LBGs to characterize their stellar masses and star formation histories (SFH). We find that the mean protocluster LBG is by a factor of $\sim2$ times more massive and more attenuated than the mean $z \sim 3$ field LBG. We take our results to suggest that ongoing mergers are not more common among protocluster LBGs than field LBGs, though protocluster LBGs appear to be more massive. We speculate that the larger mass of the protocluster LBGs contributes to the enhancement of SMBH mass and accretion rate in the protocluster, which in turn drives the observed protocluster AGN enhancement., Comment: 43 pages, 18 figures. 32 page body and 11 page appendix. Updated with published version

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 (D50% 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., 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

We examine and quantify how hybrid (e.g., UV+IR) star formation rate (SFR) estimators and the $A_{\rm FUV}$-$\beta$ relation depend on inclination for disk-dominated galaxies using spectral energy distribution modeling that utilizes the inclination-dependent attenuation curves described in Doore et al. We perform this analysis on a sample of 133 disk-dominated galaxies from the CANDELS fields and 18 disk galaxies from the Spitzer Infrared Nearby Galaxies Survey and Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel samples. We find that both the hybrid SFR estimators and the $A_{\rm FUV}$-$\beta$ relation present clear dependencies on inclination. To quantify this dependence in the hybrid SFR estimators, we derive an inclination and a far-UV and near-IR color-dependent parametric relation for converting observed UV and IR luminosities into SFRs. For the $A_{\rm FUV}$-$\beta$ relation, we introduce an inclination-dependent component that accounts for the majority of the inclination dependence with the scatter of the relation increasing with inclination. We then compare both of these inclination-dependent relations to similar inclination-independent relations found in the literature. From this comparison, we find that the UV+IR correction factor and $A_{\rm FUV}$ for our hybrid and $A_{\rm FUV}$-$\beta$ relations, respectively, result in a reduction in the residual scatter of our sample by approximately a factor of 2. Therefore, we demonstrate that inclination must be considered in hybrid SFR estimators and the $A_{\rm FUV}$-$\beta$ relation to produce more accurate SFR estimates in disk-dominated galaxies., Comment: 24 pages, 15 figures. 23 page body and 1 page appendix. Updated with published version

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., Accepted for publication in ApJ

We investigate the population of high-redshift ($3\leq z < 6$) AGN selected in the two deepest X-ray surveys, the 7 Ms \textit{Chandra} Deep Field-South and 2 Ms \textit{Chandra} Deep Field-North. Their outstanding sensitivity and spectral characterization of faint sources allow us to focus on the sub-$L_*$ regime (log$L_{\mathrm{X}}\lesssim44$), poorly sampled by previous works using shallower data, and the obscured population. Taking fully into account the individual photometric-redshift probability distribution functions, the final sample consists of $\approx102$ X-ray selected AGN at $3\leq z < 6$. The fraction of AGN obscured by column densities log$N_{\mathrm{H}}>23$ is $\sim0.6-0.8$, once incompleteness effects are taken into account, with no strong dependence on redshift or luminosity. We derived the high-redshift AGN number counts down to $F_{\mathrm{0.5-2\,keV}}=7\times10^{-18}\,\mathrm{erg\,cm^{-2}\,s^{-1}}$, extending previous results to fainter fluxes, especially at $z>4$. We put the tightest constraints to date on the low-luminosity end of AGN luminosity function at high redshift. The space-density, in particular, declines at $z>3$ at all luminosities, with only a marginally steeper slope for low-luminosity AGN. By comparing the evolution of the AGN and galaxy densities, we suggest that such a decline at high luminosities is mainly driven by the underlying galaxy population, while at low luminosities there are hints of an intrinsic evolution of the parameters driving nuclear activity. Also, the black-hole accretion rate density and star-formation rate density, which are usually found to evolve similarly at $z\lesssim3$, appear to diverge at higher redshifts., Accepted 2017 September 15 (MNRAS). 31 Pages, 25 figures

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

We present 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., Accepted for publication in ApJS (37 pages, 16 figures, and 7 tables). Data products available (https://lehmer.uark.edu/downloads/)

Neutron star-neutron star mergers are known to be associated with short gamma-ray bursts. If the neutron star equation of state is sufficiently stiff, at least some of such mergers will leave behind a supramassive or even a stable neutron star that spins rapidly with a strong magnetic field (i.e., a magnetar). Such a magnetar signature may have been observed as the X-ray plateau following a good fraction (up to 50%) of short gamma-ray bursts, and it has been expected that one may observe short gamma-ray burst-less X-ray transients powered by double neutron star mergers. A fast X-ray transient (CDF-S XT1) was recently found to be associated with a faint host galaxy whose redshift is unknown. Its X-ray and host-galaxy properties allow several possibleexplanations including a short gamma-ray burst seen off axis, a low-luminosity gamma-ray burst at high redshift, or a tidal disruption event involving an intermediate mass black hole and a white dwarf. Here we report a second X-ray transient, CDF-S XT2, that is associated with a galaxy at redshift z = 0.738. The light curve is fully consistent with being powered by a millisecond magnetar. More intriguingly, CDF-S XT2 lies in the outskirts of its star-forming host galaxy with a moderate offset from the galaxy center, as short bursts often do. The estimated event rate density of similar X-ray transients, when corrected to the local value, is consistent with the double neutron star merger rate density inferred from the detection of GW170817., Comment: 29 pages, 4 figures, 3 tables, published in Nature on 11 April 2019

Galaxy clusters in the local universe descend from high-redshift overdense regions known as protoclusters. The large gas reservoirs and high rate of galaxy interaction in protoclusters are expected to trigger star-formation activity and luminous SMBH accretion in the host galaxies. We investigated the AGN content of a gas-rich and starbursting protocluster at $z=4$, known as the Distant Red Core (DRC). We observed with Chandra (139 ks) the 13 identified members of the structure, and searched for luminous and possibly obscured AGN among them. We also tested whether a hidden AGN can power the $Ly\alpha$ blob (LAB) detected with VLT/MUSE in the DRC. We detected obscured X-ray emission from the two most gas-rich members of the DRC, named DRC-1 and DRC-2. Both of them are resolved into multiple interacting clumps in high-resolution ALMA and HST observations. In particular, DRC-2 is found to host a luminous ($L_{2-10\,\mathrm{keV}}\approx3\times10^{45}\,\mathrm{erg\,s^{-1}}$) Compton-thick ($N_H\gtrsim10^{24}\,\mathrm{cm^{-2}}$) QSO, comparable to the most luminous QSOs known at all cosmic times. The AGN fraction among DRC members is consistent with results found for lower redshift protoclusters. However, X-ray stacking analysis reveals that SMBH accretion is likely also taking place in other DRC galaxies that are not detected individually by Chandra. Our results point toward the presence of a strong link between large gas reservoirs, galaxy interactions, and luminous and obscured nuclear activity in protocluster members. The powerful and obscured QSO detected in DRC-2 is likely powering the nearby LAB detected with VLT/MUSE, possibly through photoionization; however, we propose that the diffuse $Ly\alpha$ emission may be due to gas shocked by a massive outflow launched by DRC-2 over a $\approx10$ kpc scale., Comment: 13 pages, accepted for publication in A&A

We investigate X-ray binary (XRB) luminosity function (XLF) scaling relations for Chandra detected populations of low-mass XRBs (LMXBs) within the footprints of 24 early-type galaxies. Our sample includes Chandra and HST observed galaxies at D < 25 Mpc that have estimates of the globular cluster (GC) specific frequency (SN) reported in the literature. As such, we are able to directly classify X-ray-detected sources as being either coincident with unrelated background/foreground objects, GCs, or sources that are within the fields of the galaxy targets. We model the GC and field LMXB population XLFs for all galaxies separately, and then construct global models characterizing how the LMXB XLFs vary with galaxy stellar mass and SN. We find that our field LMXB XLF models require a component that scales with SN, and has a shape consistent with that found for the GC LMXB XLF. We take this to indicate that GCs are "seeding" the galactic field LMXB population, through the ejection of GC-LMXBs and/or the diffusion of the GCs in the galactic fields themselves. However, we also find that an important LMXB XLF component is required for all galaxies that scales with stellar mass, implying that a substantial population of LMXBs are formed "in situ," which dominates the LMXB population emission for galaxies with SN < 2. For the first time, we provide a framework quantifying how directly-associated GC LMXBs, GC-seeded LMXBs, and in-situ LMXBs contribute to LMXB XLFs in the broader early-type galaxy population., Accepted for publication in ApJS. Data products and catalogs available at https://lehmer.uark.edu/downloads/

We present a broadband X-ray spectral analysis of the M51 system, including the dual active galactic nuclei (AGN) and several off-nuclear point sources. Using a deep observation by NuSTAR, new high-resolution coverage of M51b by Chandra, and the latest X-ray torus models, we measure the intrinsic X-ray luminosities of the AGN in these galaxies. The AGN of M51a is found to be Compton thick, and both AGN have very low accretion rates ($��_{\rm Edd} 10^{39}$ erg s$^{-1}$ and find that a spectral turnover is present below 10 keV in most such sources, in line with recent results on ultraluminous X-ray sources., Accepted for publication in ApJ

We present CO J=4-3 line and 3 mm dust continuum observations of a 100 kpc-scale filamentary Ly{\alpha} nebula (SSA22 LAB18) at z=3.1 using the Atacama Large Millimeter/submillimeter Array (ALMA). We detected the CO J=4-3 line at a systemic z(CO)=3.093 {\pm} 0.001 at 11 {\sigma} from one of the ALMA continuum sources associated with the Ly{\alpha} filament. We estimated the CO J=4-3 luminosity of L'CO(4-3)=(2.3\pm0.2)x10^9 K km s^{-1} pc^2 for this CO source, which is one order of magnitude smaller than those of typical z>1 dusty star-forming galaxies (DSFGs) of similar far-infrared luminosity L(IR)~10^{12} Lsun. We derived a molecular gas mass of Mgas=(4.4^{+0.9}_{-0.6})x10^9 Msun and a star-formation rate of SFR=270\pm160 Msun yr^{-1}. We also estimated a gas depletion time of {\tau}(dep)=17\pm10 Myr, being shorter than those of typical DSFGs. It is suggested that this source is in a transition phase from DSFG to a gas-poor, early-type galaxy. From ALMA to Herschel multi-band dust continuum observations, we measured a dust emissivity index {\beta}=2.3\pm0.2, which is similar to those of local gas-poor, early-type galaxies. Such a high {\beta} can be reproduced by specific chemical compositions for interstellar dust at the submillimeter wavelengths from recent laboratory experiments. ALMA CO and multi-band dust continuum observations can constrain the evolutionary stage of high-redshift galaxies through {\tau}(dep) and {\beta}, and thus we can investigate dust chemical compositions even in the early Universe.

We present an X-ray point-source catalog from the XMM-Large Scale Structure survey region (XMM-LSS), one of the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS) fields. We target the XMM-LSS region with $1.3$ Ms of new XMM-Newton AO-15 observations, transforming the archival X-ray coverage in this region into a 5.3 deg$^2$ contiguous field with uniform X-ray coverage totaling $2.7$ Ms of flare-filtered exposure, with a $46$ ks median PN exposure time. We provide an X-ray catalog of 5242 sources detected in the soft (0.5-2 keV), hard (2-10 keV), and/or full (0.5-10 keV) bands with a 1% expected spurious fraction determined from simulations. A total of 2381 new X-ray sources are detected compared to previous source catalogs in the same area. Our survey has flux limits of $1.7\times10^{-15}$, $1.3\times10^{-14}$, and $6.5\times10^{-15}$ erg cm$^{-2}$ s$^{-1}$ over 90% of its area in the soft, hard, and full bands, respectively, which is comparable to those of the XMM-COSMOS survey. We identify multiwavelength counterpart candidates for 99.9% of the X-ray sources, of which 93% are considered as reliable based on their matching likelihood ratios. The reliabilities of these high-likelihood-ratio counterparts are further confirmed to be $\approx 97$% reliable based on deep Chandra coverage over $\approx 5$% of the XMM-LSS region. Results of multiwavelength identifications are also included in the source catalog, along with basic optical-to-infrared photometry and spectroscopic redshifts from publicly available surveys. We compute photometric redshifts for X-ray sources in 4.5 deg$^2$ of our field where forced-aperture multi-band photometry is available; $>70$% of the X-ray sources in this subfield have either spectroscopic or high-quality photometric redshifts., MNRAS, accepted. 34 pages, 25 figures, and 8 tables. The data products are available at this http url: http://personal.psu.edu/wnb3/xmmservs/xmmservs.html

We present 15 high mass X-ray binary (HMXB) candidates in the disk of M31 for which we are able to infer compact object type, spectral type of the donor star, and age using multiwavelength observations from NuSTAR, Chandra, and the Hubble Space Telescope (HST). The hard X-ray colors and luminosities from NuSTAR permit the tentative classification of accreting X-ray binary systems by compact object type, distinguishing black hole from neutron star systems. We find hard state black holes, pulsars, and non-magnetized neutron stars associated with optical point source counterparts with similar frequency. We also find nine non-magnetized neutron stars coincident with globular clusters and an equal number of pulsars with and without point source optical counterparts. We perform spectral energy distribution (SED) fitting for the most likely optical counterparts to the HMXB candidates, finding 7 likely high mass stars and 1 possible red Helium burning star. The remaining 7 HMXB optical counterparts have poor SED fits, so their companion stars remain unclassified. Using published star formation histories, we find that the majority of HMXB candidates --- X-ray sources with UV-bright point source optical counterpart candidates --- are found in regions with star formation bursts less than 50 Myr ago, with 3 associated with young stellar ages (, 16 pages, 9 figures, 4 tables

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

In deep X-ray surveys, active galactic nuclei (AGNs) with a broad range of luminosities have been identified. However, cosmologically distant low-luminosity AGN (LLAGN, $L_{\mathrm{X}} \lesssim 10^{42}$ erg s$^{-1}$) identification still poses a challenge due to significant contamination from host galaxies. Based on the 7 Ms Chandra Deep Field-South (CDF-S) survey, the longest timescale ($\sim 17$ years) deep X-ray survey to date, we utilize an X-ray variability selection technique to search for LLAGNs that remain unidentified among the CDF-S X-ray sources. We find 13 variable sources from 110 unclassified CDF-S X-ray sources. Except for one source which could be an ultraluminous X-ray source, the variability of the remaining 12 sources is most likely due to accreting supermassive black holes. These 12 AGN candidates have low intrinsic X-ray luminosities, with a median value of $7 \times10^{40}$ erg s$^{-1}$. They are generally not heavily obscured, with an average effective power-law photon index of 1.8. The fraction of variable AGNs in the CDF-S is independent of X-ray luminosity and is only restricted by the total number of observed net counts, confirming previous findings that X-ray variability is a near-ubiquitous property of AGNs over a wide range of luminosities. There is an anti-correlation between X-ray luminosity and variability amplitude for high-luminosity AGNs, but as the luminosity drops to $\lesssim 10^{42}$ erg s$^{-1}$, the variability amplitude no longer appears dependent on the luminosity. The entire observed luminosity-variability trend can be roughly reproduced by an empirical AGN variability model based on a broken power-law power spectral density function., Comment: 18 pages, 11 figures, accepted for publication in ApJ

We present the X-ray point source population of NGC 7457 based on 124 ks of Chandra observations. Previous deep Chandra observations of low mass X-ray binaries (LMXBs) in early-type galaxies have typically targeted the large populations of massive galaxies. NGC 7457 is a nearby, early-type galaxy with a stellar luminosity of $1.7\times10^{10} L_{K\odot}$, allowing us to investigate the populations in a relatively low mass galaxy. We classify the detected X-ray sources into field LMXBs, globular cluster LMXBs, and background AGN based on identifying optical counterparts in new HST/ACS images. We detect 10 field LMXBs within the $r_{ext}$ ellipse of NGC 7457 (with semi-major axis $\sim$ 9.1 kpc, ellipticity = 0.55). The corresponding number of LMXBs with $L_{x}>2\times10^{37}erg/s$ per stellar luminosity is consistent with that observed in more massive galaxies, $\sim 7$ per $10^{10} L_{K\odot}$. We detect a small globular cluster population in these HST data and show that its colour distribution is likely bimodal and that its specific frequency is similar to that of other early type galaxies. However, no X-ray emission is detected from any of these clusters. Using published data for other galaxies, we show that this non-detection is consistent with the small stellar mass of these clusters. We estimate that 0.11 (and 0.03) LMXBs are expected per $10^{6}M_{\odot}$ in metal-rich (and metal-poor) globular clusters. This corresponds to 1100 (and 330) LMXBs per $10^{10} L_{K\odot}$, highlighting the enhanced formation efficiency of LMXBs in globular clusters. A nuclear X-ray source is detected with $L_{x}$ varying from $2.8-6.8\times10^{38}erg/s$. Combining this $L_{x}$ with a published dynamical mass estimate for the central SMBH in NGC 7457, we find that $L_{x}/L_{Edd}$ varies from $0.5-1.3\times10^{-6}$., 9 pages, 7 figures, 1 table, accepted for publication in MNRAS

The stellar initial mass function (IMF), which is often assumed to be universal across unresolved stellar populations, has recently been suggested to be "bottom-heavy" for massive ellipticals. In these galaxies, the prevalence of gravity-sensitive absorption lines (e.g. Na I and Ca II) in their near-IR spectra implies an excess of low-mass ($m = 8$ $M_\odot$) would lead to a corresponding deficit of neutron stars and black holes, and therefore of low-mass X-ray binaries (LMXBs), per unit near-IR luminosity in these galaxies. Peacock et al. (2014) searched for evidence of this trend and found that the observed number of LMXBs per unit $K$-band luminosity ($N/L_K$) was nearly constant. We extend this work using new and archival Chandra X-ray Observatory (Chandra) and Hubble Space Telescope (HST) observations of seven low-mass ellipticals where $N/L_K$ is expected to be the largest and compare these data with a variety of IMF models to test which are consistent with the observed $N/L_K$. We reproduce the result of Peacock et al. (2014), strengthening the constraint that the slope of the IMF at $m >= 8$ $M_\odot$ must be consistent with a Kroupa-like IMF. We construct an IMF model that is a linear combination of a Milky Way-like IMF and a broken power-law IMF, with a steep slope ($��_1=$ $3.84$) for stars < 0.5 $M_\odot$ (as suggested by near-IR indices), and that flattens out ($��_2=$ $2.14$) for stars > 0.5 $M_\odot$, and discuss its wider ramifications and limitations., Accepted for publication in ApJ; 7 pages, 2 figures, 1 table

We present X-ray source catalogs for the $\approx7$ Ms exposure of the Chandra Deep Field-South (CDF-S), which covers a total area of 484.2 arcmin$^2$. Utilizing WAVDETECT for initial source detection and ACIS Extract for photometric extraction and significance assessment, we create a main source catalog containing 1008 sources that are detected in up to three X-ray bands: 0.5-7.0 keV, 0.5-2.0 keV, and 2-7 keV. A supplementary source catalog is also provided including 47 lower-significance sources that have bright ($K_s\le23$) near-infrared counterparts. We identify multiwavelength counterparts for 992 (98.4%) of the main-catalog sources, and we collect redshifts for 986 of these sources, including 653 spectroscopic redshifts and 333 photometric redshifts. Based on the X-ray and multiwavelength properties, we identify 711 active galactic nuclei (AGNs) from the main-catalog sources. Compared to the previous $\approx4$ Ms CDF-S catalogs, 291 of the main-catalog sources are new detections. We have achieved unprecedented X-ray sensitivity with average flux limits over the central $\approx1$ arcmin$^2$ region of $\approx1.9\times10^{-17}$, $6.4\times10^{-18}$, and $2.7\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$ in the three X-ray bands, respectively. We provide cumulative number-count measurements observing, for the first time, that normal galaxies start to dominate the X-ray source population at the faintest 0.5-2.0 keV flux levels. The highest X-ray source density reaches $\approx50\,500$ deg$^{-2}$, and $47\%\pm4\%$ of these sources are AGNs ($\approx23\,900$ deg$^{-2}$)., Comment: 31 pages, 31 figures, 8 tables, ApJS in press, minor text updates; full-resolution PDF version and data products available at http://www.astro.psu.edu/users/niel/cdfs/cdfs-chandra.html

Galaxy-scale bars are expected to provide an effective means for driving material towards the central region in spiral galaxies, and possibly feeding supermassive black holes (BHs). Here we present a statistically-complete study of the effect of bars on average BH accretion. From a well-selected sample of 50,794 spiral galaxies (with M* ~ 0.2-30 x 10^10 Msun) extracted from the Sloan Digital Sky Survey Galaxy Zoo 2 project, we separate those sources considered to contain galaxy-scale bars from those that do not. Using archival data taken by the Chandra X-ray Observatory, we identify X-ray luminous (L_X >~ 10^41 erg/s) active galactic nuclei (AGN) and perform an X-ray stacking analysis on the remaining X-ray undetected sources. Through X-ray stacking, we derive a time-averaged look at accretion for galaxies at fixed stellar mass and star formation rate, finding that the average nuclear accretion rates of galaxies with bar structures are fully consistent with those lacking bars (Mdot_acc ~ 3 x 10^-5 Msun/yr). Hence, we robustly conclude that large-scale bars have little or no effect on the average growth of BHs in nearby (z < 0.15) galaxies over gigayear timescales., 10 pages, 5 figures, accepted to ApJ

We present a new technique for empirically calibrating how the X-ray luminosity function (XLF) of X-ray binary (XRB) populations evolves following a star-formation event. We first utilize detailed stellar population synthesis modeling of far-UV to far-IR photometry of the nearby face-on spiral galaxy M51 to construct maps of the star-formation histories (SFHs) on subgalactic (~400 pc) scales. Next, we use the ~850 ks cumulative Chandra exposure of M51 to identify and isolate 2-7 keV detected point sources within the galaxy, and we use our SFH maps to recover the local properties of the stellar populations in which each X-ray source is located. We then divide the galaxy into various subregions based on their SFH properties (e.g., star-formation rate [SFR] per stellar mass [M*] and mass-weighted stellar age) and group the X-ray point sources according to the characteristics of the regions in which they are found. Finally, we construct and fit a parameterized XLF model that quantifies how the XLF shape and normalization evolves as a function of the XRB population age. Our best-fit model indicates the XRB XLF per unit stellar mass declines in normalization, by ~3-3.5 dex, and steepens in slope from ~10 Myr to ~10 Gyr. We find that our technique recovers results from past studies of how XRB XLFs and XRB luminosity scaling relations vary with age and provides a self-consistent picture for how the XRB XLF evolves with age., Comment: Accepted for publication in ApJ (19 pages, 15 figures, and 2 tables)

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., Accepted for publication in ApJ (15 pages, 4 tables, 6 figures)

We exploit ALMA 870um observations to measure the star-formation rates (SFRs) of eight X-ray detected Active Galactic Nuclei (AGNs) in a z~3.1 protocluster, four of which reside in extended Ly-alpha haloes (often termed Ly-alpha blobs: LABs). Three of the AGNs are detected by ALMA and have implied SFRs of ~220-410~M_sun/yr; the non detection of the other five AGNs places SFR upper limits of 100 kpc) do not host more luminous star formation than the smaller LABs, despite being an order of magnitude brighter in Ly-alpha emission. We use these results to discuss star formation as the power source of LABs., MNRAS, in press: 10 pages, 3 figures, 2 tables

Understanding galaxy formation and evolution requires studying the interplay between the growth of galaxies and the growth of their black holes across cosmic time. Here we explore a sample of Ha-selected star-forming galaxies from the HiZELS survey and use the wealth of multi-wavelength data in the COSMOS field (X-rays, far-infrared and radio) to study the relative growth rates between typical galaxies and their central supermassive black holes, from z=2.23 to z=0. Typical star-forming galaxies at z~1-2 have black hole accretion rates (BHARs) of 0.001-0.01 Msun/yr and star formation rates (SFRs) of ~10-40 Msun/yr, and thus grow their stellar mass much quicker than their black hole mass (~3.3 orders of magnitude faster). However, ~3% of the sample (the sources detected directly in the X-rays) show a significantly quicker growth of the black hole mass (up to 1.5 orders of magnitude quicker growth than the typical sources). BHARs fall from z=2.23 to z=0, with the decline resembling that of star formation rate density or the typical SFR. We find that the average black hole to galaxy growth (BHAR/SFR) is approximately constant for star-forming galaxies in the last 11 Gyrs. The relatively constant BHAR/SFR suggests that these two quantities evolve equivalently through cosmic time and with practically no delay between the two., 10 pages, 8 figures, accepted to MNRAS

Ly-alpha blobs (LABs) offer insight into the complex interface between galaxies and their circumgalactic medium. Whilst some LABs have been found to contain luminous star-forming galaxies and active galactic nuclei that could potentially power the Ly-alpha emission, others appear not to be associated with obvious luminous galaxy counterparts. It has been speculated that LABs may be powered by cold gas streaming on to a central galaxy, providing an opportunity to directly observe the `cold accretion' mode of galaxy growth. Star-forming galaxies in LABs could be dust obscured and therefore detectable only at longer wavelengths. We stack deep SCUBA-2 observations of the SSA22 field to determine the average 850um flux density of 34 LABs. We measure S_850 = 0.6 +/- 0.2mJy for all LABs, but stacking the LABs by size indicates that only the largest third (area > 1794 kpc^2) have a mean detection, at 4.5 sigma, with S_850 = 1.4 +/- 0.3mJy. Only two LABs (1 and 18) have individual SCUBA-2 > 3.5 sigma detections at a depth of 1.1mJy/beam. We consider two possible mechanisms for powering the LABs and find that central star formation is likely to dominate the emission of Ly-alpha, with cold accretion playing a secondary role., Accepted MNRAS 16 May 2016

We present the results of a joint Swift-NuSTAR spectroscopy campaign on M31. We focus on the five brightest globular cluster X-ray sources in our fields. Two of these had previously been argued to be black hole candidates on the basis of apparent hard-state spectra at luminosities above those for which neutron stars are in hard states. We show that these two sources are likely to be Z-sources (i.e. low magnetic field neutron stars accreting near their Eddington limits), or perhaps bright atoll sources (low magnetic field neutron stars which are just a bit fainter than this level) on the basis of simultaneous Swift and NuSTAR spectra which cover a broader range of energies. These new observations reveal spectral curvature above 6-8 keV that would be hard to detect without the broader energy coverage the NuSTAR data provide relative to Chandra and XMM-Newton. We show that the other three sources are also likely to be bright neutron star X-ray binaries, rather than black hole X-ray binaries. We discuss why it should already have been realized that it was unlikely that these objects were black holes on the basis of their being persistent sources, and we re-examine past work which suggested that tidal capture products would be persistently bright X-ray emitters. We discuss how this problem is likely due to neglecting disk winds in older work that predict which systems will be persistent and which will be transient., 12 pages, 3 figures, accepted to MNRAS

We present improved point-source catalogs for the 2 Ms Chandra Deep Field-North (CDF-N) and the 250 ks Extended Chandra Deep Field-South (E-CDF-S), implementing a number of recent improvements in Chandra source-cataloging methodology. For the CDF-N/E-CDF-S, we provide a main catalog that contains 683/1003 X-ray sources detected with wavdetect at a false-positive probability threshold of $10^{-5}$ that also satisfy a binomial-probability source-selection criterion of $P, 45 pages, 57 figures, 14 tables, ApJS in press; Data products and images available at http://www2.astro.psu.edu/users/niel/hdf/hdf-chandra.html and http://www2.astro.psu.edu/users/niel/ecdfs/ecdfs-chandra.html arXiv admin note: text overlap with arXiv:1105.5643

We present results from sensitive, multi-epoch NuSTAR observations of the late-type star-forming galaxy M83 (d=4.6 Mpc), which is the first investigation to spatially resolve the hard (E>10 keV) X-ray emission of this galaxy. The nuclear region and ~ 20 off-nuclear point sources, including a previously discovered ultraluminous X-ray (ULX) source, are detected in our NuSTAR observations. The X-ray hardnesses and luminosities of the majority of the point sources are consistent with hard X-ray sources resolved in the starburst galaxy NGC 253. We infer that the hard X-ray emission is most likely dominated by intermediate accretion state black hole binaries and neutron star low-mass X-ray binaries (Z-sources). We construct the X-ray binary luminosity function (XLF) in the NuSTAR band for an extragalactic environment for the first time. The M83 XLF has a steeper XLF than the X-ray binary XLF in NGC 253, consistent with previous measurements by Chandra at softer X-ray energies. The NuSTAR integrated galaxy spectrum of M83 drops quickly above 10 keV, which is also seen in the starburst galaxies NGC253, NGC 3310 and NGC 3256. The NuSTAR observations constrain any AGN to be either highly obscured or to have an extremely low luminosity of $_{\sim}^, Comment: Accepted for publication in ApJ (25 pages, 17 figures)

We present measurements of the evolution of normal-galaxy X-ray emission from $z \approx$ 0-7 using local galaxies and galaxy samples in the 6 Ms Chandra Deep Field-South (CDF-S) survey. The majority of the CDF-S galaxies are observed at rest-frame energies above 2 keV, where the emission is expected to be dominated by X-ray binary (XRB) populations; however, hot gas is expected to provide small contributions to the observed- frame < 1 keV emission at $z < 1$. We show that a single scaling relation between X-ray luminosity ($L_{\rm X}$) and star-formation rate (SFR) is insufficient for characterizing the average X-ray emission at all redshifts. We establish that scaling relations involving not only SFR, but also stellar mass ($M_\star$) and redshift, provide significantly improved characterizations of the average X-ray emission from normal galaxy populations at $z \approx$ 0-7. We further provide the first empirical constraints on the redshift evolution of X-ray emission from both low-mass XRB (LMXB) and high-mass XRB (HMXB) populations and their scalings with $M_\star$ and SFR, respectively. We find $L_{\rm 2-10 keV}({\rm LMXB})/M_\star \propto (1 + z)^{2-3}$ and $L_{\rm 2-10 keV}({\rm HMXB})$/SFR $\propto (1 + z)$, and show that these relations are consistent with XRB population-synthesis model predictions, which attribute the increase in LMXB and HMXB scaling relations with redshift as being due to declining host galaxy stellar ages and metallicities, respectively. We discuss how emission from XRBs could provide an important source of heating to the intergalactic medium in the early Universe, exceeding that of active galactic nuclei., Accepted for publication in ApJ (23 pages, 13 figures, 5 tables)

We study the X-ray properties of 393 optically selected early-type galaxies (ETGs) over the redshift range of z≈ 0.0–1.2 in the Chandra Deep Fields (CDFs). To measure the average X-ray properties of the ETG population, we use X-ray stacking analyses with a subset of 158 passive ETGs (148 of which were individually undetected in X-ray). This ETG subset was constructed to span the redshift ranges of z= 0.1–1.2 in the ≈4 Ms CDF-South and ≈2 Ms CDF-North and z= 0.1–0.6 in the ≈250 ks Extended-CDF-South where the contribution from individually undetected active galactic nuclei (AGN) is expected to be negligible in our stacking. We find that 55 of the ETGs are detected individually in X-ray, and 12 of these galaxies have properties consistent with being passive hot-gas-dominated systems (i.e. systems not dominated by an X-ray bright AGN). On the basis of our analyses, we find little evolution in the mean 0.5–2 keV to B-band luminosity ratio (LX/LB∝ [1 +z]1.2) since z≈ 1.2, implying that some heating mechanism prevents the gas from cooling in these systems. We consider that feedback from radio-mode AGN activity could be responsible for heating the gas. We select radio AGN in the ETG population using their far-infrared/radio flux ratio. Our radio observations allow us to constrain the duty cycle history of radio AGN activity in our ETG sample. We estimate that if scaling relations between radio and mechanical power hold out to z≈ 1.2 for the ETG population being studied here, the average mechanical power from AGN activity is a factor of ≈1.4–2.6 times larger than the average radiative cooling power from hot gas over the redshift range z≈ 0–1.2. The excess of inferred AGN mechanical power from these ETGs is consistent with that found in the local Universe for similar types of galaxies.

We investigate the local supermassive black hole (SMBH) density function and relative mass accretion rates of all active galactic nuclei (AGNs) identified in a volume-limited sample of infrared (IR) bright galaxies (L_IR > 3 x 10^9 L_sun) to D

We present the results of a 400ks Chandra survey of 29 extended Ly-alpha emitting nebulae (Ly-alpha Blobs, LABs) in the z=3.09 proto-cluster in the SSA22 field. We detect luminous X-ray counterparts in five LABs, implying a large fraction of active galactic nuclei (AGN) in LABs, f_AGN = 17% down to L_2-32keV ~ 10^44 erg/s. All of the AGN appear to be heavily obscured, with spectral indices implying obscuring column densities of N_H > 10^23 cm^-2. The AGN fraction should be considered a lower limit, since several more LABs not detected with Chandra show AGN signatures in their mid-infrared emission. We show that the UV luminosities of the AGN are easily capable of powering the extended Ly-alpha emission via photo-ionization alone. When combined with the UV flux from a starburst component, and energy deposited by mechanical feedback, we demonstrate that 'heating' by a central source, rather than gravitational cooling is the most likely power source of LABs. We argue that all LABs could be powered in this manner, but that the luminous host galaxies are often just below the sensitivity limits of current instrumentation, or are heavily obscured. No individual LABs show evidence for extended X-ray emission, and a stack equivalent to a >9Ms exposure of an average LAB also yields no statistical detection of a diffuse X-ray component. The resulting diffuse X-ray/Ly-alpha luminosity limit implies there is no hot (T>10^7 K) gas component in these halos, and also rules out inverse Compton scattering of cosmic microwave background photons, or local far-infrared photons, as a viable power source for LABs., 9 pages, 4 figures. Accepted for publication in the Astrophysical Journal. This version updated to accepted version

We present point-source catalogs for the ~2 Ms exposure of the Chandra Deep Field-South (CDF-S); this is one of the two most-sensitive X-ray surveys ever performed. The survey covers an area of ~436 arcmin^2 and reaches on-axis sensitivity limits of ~1.9x10^{-17} and ~1.3x10^{-16} ergs/cm^2/s for the 0.5-2.0 and 2-8 keV bands, respectively. Four hundred and sixty-two X-ray point sources are detected in at least one of three X-ray bands that were searched; 135 of these sources are new compared to the previous ~1 Ms CDF-S detections. Source positions are determined using centroid and matched-filter techniques; the median positional uncertainty is ~0.36". The X-ray-to-optical flux ratios of the newly detected sources indicate a variety of source types; ~55% of them appear to be active galactic nuclei while ~45% appear to be starburst and normal galaxies. In addition to the main Chandra catalog, we provide a supplementary catalog of 86 X-ray sources in the ~2 Ms CDF-S footprint that was created by merging the ~250 ks Extended Chandra Deep Field-South with the CDF-S; this approach provides additional sensitivity in the outer portions of the CDF-S. A second supplementary catalog that contains 30 X-ray sources was constructed by matching lower significance X-ray sources to bright optical counterparts (R, 18 pages, 15 figures, ApJS in press. Data and images available at http://www.astro.psu.edu/users/niel/cdfs/cdfs-chandra.html

We present VLT VIMOS, Keck DEIMOS and Keck LRIS multi-object spectra of 367 sources in the field of the z ~ 3.09 protocluster SSA22. Sources are spectroscopically classified via template matching, allowing new identifications for 206 extragalactic sources, including 36 z > 2 Lyman-break galaxies (LBGs) and Lyman-�� emitters (LAEs), 8 protocluster members, and 94 X-ray sources from the ~ 400 ks Chandra deep survey of SSA22. Additionally, in the area covered by our study, we have increased by ~ 4, 13, and 6 times the number of reliable redshifts of sources at 1.0 < z < 2.0, at z > 3.4, and with X-Ray emission, respectively. We compare our results with past spectroscopic surveys of SSA22 to investigate the completeness of the LBGs and the X-Ray properties of the new spectroscopically-classified sources in the SSA22 field., 27 pages, 14 figures, 6 tables. Paper has been accepted for publication in MNRAS

The over-dense environments of protoclusters of galaxies in the early Universe (z>2) are expected to accelerate the evolution of galaxies, with an increased rate of stellar mass assembly and black hole accretion compared to co-eval galaxies in the average density `field'. These galaxies are destined to form the passive population of massive systems that dominate the cores of rich clusters today. While signatures of accelerated growth of galaxies in the SSA22 protocluster (z=3.1) have been observed, the mechanism driving this remains unclear. Here we show an enhanced rate of galaxy-galaxy mergers could be responsible. We morphologically classify Lyman-break Galaxies (LBGs) in the SSA22 protocluster and compare these to those of galaxies in the field at z=3.1 as either active mergers or non-merging using Hubble Space Telescope ACS/F814W imaging, probing the rest-frame ultraviolet stellar light. We measure a merger fraction of 48+/-10 per cent for LBGs in the protocluster compared to 30+/-6 per cent for the field. Although the excess is marginal, an enhanced rate of mergers in SSA22 hints that galaxy-galaxy mergers are one of the key channels driving accelerated star formation and AGN growth in protocluster environments., Accepted MNRAS 2015 October 20

We present nearly simultaneous Chandra and NuSTAR observations of two actively star-forming galaxies within 50 Mpc: NGC 3256 and NGC 3310. Both galaxies are detected by both Chandra and NuSTAR, which together provide the first-ever spectra of these two galaxies spanning 0.3-30 keV. The X-ray emission from both galaxies is spatially resolved by Chandra; we find that hot gas dominates the E < 1-3 keV emission while ultraluminous X-ray sources (ULXs) dominate at E > 1-3 keV. The NuSTAR galaxy-wide spectra of both galaxies follow steep power-law distributions with Gamma ~ 2.6 at E > 5-7 keV, similar to the spectra of bright individual ULXs and other galaxies that have been studied by NuSTAR. We find that both NGC 3256 and NGC 3310 have X-ray detected sources coincident with nuclear regions; however, the steep NuSTAR spectra of both galaxies restricts these sources to be either low luminosity AGN or non-AGN in nature (e.g., ULXs or crowded X-ray sources that reach L2-10 keV ~ 10^40 erg/s cannot be ruled out). Combining our constraints on the 0.3-30 keV spectra of NGC 3256 and NGC 3310 with equivalent measurements for nearby star-forming galaxies M83 and NGC 253, we analyze the SFR-normalized spectra of these starburst galaxies. The spectra of all four galaxies show sharply declining power-law slopes above 3-6 keV due to ULX populations. Our observations therefore constrain the average spectra of luminous accreting binaries (i.e., ULXs). This result is similar to the super-Eddington accreting ULXs that have been studied individually in a targeted NuSTAR ULX program. We also find that NGC 3310 exhibits a factor of ~3-10 elevation of X-ray emission over the other star-forming galaxies. We argue that the excess is most likely explained by the relatively low metallicity of the young stellar population in NGC 3310., 13 pages, 6 figures, accepted for publication in ApJ

We report the results of $1^{\prime}.5 \times3^{\prime}$ mapping at 1.1~mm with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the central region of the $z=3.09$ SSA22 protocluster. By combining our source catalog with archival spectroscopic redshifts, we find that eight submillimeter galaxies (SMGs) with flux densities, $S_{\rm 1.1~mm}=0.7-6.4$~mJy ($L_{\rm IR}\sim10^{12.1}-10^{13.1}L_\odot$) are at $z=3.08-3.10$. Not only are these SMGs members of the protocluster but they in fact reside within the node at the junction of the 50 Mpc-scale filamentary three-dimensional structure traced by Lyman-$\alpha$ emitters (LAEs) in this field. The eight SMGs account for a star formation rate density (SFRD) $\sim$10 $M_\odot$ yr$^{-1}$ Mpc$^{-3}$ in the node, which is two orders of magnitudes higher than the global SFRD at this redshift. We find that four of the eight SMGs host a X-ray luminous active galactic nuclei (AGN). Our results suggest that the vigorous star formation activity and the growth of super massive black holes (SMBHs) occurred simultaneously in the densest regions at $z\sim3$, which may correspond to the most active historical phase of the massive galaxy population found in the core of the clusters in the present universe. Two SMGs are associated with Lyman-$\alpha$ blobs (LABs), implying that the two populations coexist in high density environments for a few cases., Comment: Accepted for Publication in ApJL. 6 pages, 3 figures

We present partial-correlation analyses that examine the strengths of the relationships between L_UV, L_X, Alpha_OX, and redshift for optically-selected AGNs. We extend the work of Strateva et al. (2005), that analyzed optically-selected AGNs from the Sloan Digital Sky Survey (SDSS), by including 52 moderate-luminosity, optically-selected AGNs from the COMBO-17 survey with corresponding deep (~250 ks to 1 Ms) X-ray observations from the Extended Chandra Deep Field-South. The COMBO-17 survey extends ~3 magnitudes deeper than the SDSS and probes the moderate-luminosity AGNs that numerically dominate the AGN population in the Universe. We also include recently published observations of 19 high-redshift, optically-selected AGNs, and 46 luminous, low-redshift AGNs from the Bright Quasar Survey. The full sample used in our analysis consists of 333 AGNs, extending out to z~6, with 293 (88%) having X-ray detections. The sample spans five decades in UV luminosity and four decades in X-ray luminosity. We confirm that Alpha_OX is strongly anti-correlated with L_UV (13.6 sigma), the highest significance found for this relation to date, and find evidence suggesting that the slope of this relation may be dependent on L_UV. We find that no significant correlation exists between Alpha_OX and redshift (1.3 sigma), and constrain the maximum evolution of AGN UV-to-X-ray flux ratios to be less than 30% (1 sigma) out to z=5. Using our sample's high X-ray detection fraction, we also find a significant anti-correlation (3.0 sigma) between Alpha_OX and L_X . We make comparisons to earlier studies on this topic and discuss implications for X-ray vs. optical luminosity functions., Comment: 16 pages, 10 figures, emulateapj. Accepted for publication in AJ

We report the identification of a bright hard X-ray source dominating the M31 bulge above 25 keV from a simultaneous NuSTAR-Swift observation. We find that this source is the counterpart to Swift J0042.6+4112, which was previously detected in the Swift BAT All-sky Hard X-ray Survey. This Swift BAT source had been suggested to be the combined emission from a number of point sources; our new observations have identified a single X-ray source from 0.5 to 50 keV as the counterpart for the first time. In the 0.5-10 keV band, the source had been classified as an X-ray binary candidate in various Chandra and XMM studies; however, since it was not clearly associated with Swift J0042.6+4112, the previous E < 10 keV observations did not generate much attention. This source has a spectrum with a soft X-ray excess (kT~ 0.2 keV) plus a hard spectrum with a power law of Gamma ~ 1 and a cutoff around 15-20 keV, typical of the spectral characteristics of accreting pulsars. Unfortunately, any potential pulsation was undetected in the NuSTAR data, possibly due to insufficient photon statistics. The existing deep HST images exclude high-mass (>3 Msun) donors at the location of this source. The best interpretation for the nature of this source is an X-ray pulsar with an intermediate-mass (, Comment: 10 pages, 6 figures, accepted for publication in ApJ

We have discovered a persistent, but highly variable X-ray source in the nearby starburst galaxy NGC 253. The source varies at the level of a factor of about 5 in count rate on timescales of a few hours. Two long observations of the source with Chandra and XMM-Newton show suggestive evidence for the source having a period of about 14-15 hours, but the time sampling in existing data is insufficient to allow a firm determination that the source is periodic. Given the amplitude of variation and the location in a nuclear starburst, the source is likely to be a Wolf-Rayet X-ray binary, with the tentative period being the orbital period of the system. In light of the fact that we have demonstrated that careful examination of the variability of moderately bright X-ray sources in nearby galaxies can turn up candidate Wolf-Rayet X-ray binaries, we discuss the implications of Wolf-Rayet X-ray binaries for predictions of the gravitational wave source event rate, and, potentially, interpretations of the events., 9 pages, 3 figures, accepted to MNRAS

We present direct constraints on how the formation of low-mass X-ray binary (LMXB) populations in galactic fields depends on stellar age. In this pilot study, we utilize Chandra and Hubble Space Telescope (HST) data to detect and characterize the X-ray point source populations of three nearby early-type galaxies: NGC 3115, 3379, and 3384. The luminosity-weighted stellar ages of our sample span 3-10 Gyr. X-ray binary population synthesis models predict that the field LMXBs associated with younger stellar populations should be more numerous and luminous per unit stellar mass than older populations due to the evolution of LMXB donor star masses. Crucially, the combination of deep Chandra and HST observations allows us to test directly this prediction by identifying and removing counterparts to X-ray point sources that are unrelated to the field LMXB populations, including LMXBs that are formed dynamically in globular clusters, Galactic stars, and background AGN/galaxies. We find that the "young" early-type galaxy NGC 3384 (~2-5 Gyr) has an excess of luminous field LMXBs (L_X > (5-10) x 10^37 erg/s) per unit K-band luminosity (L_K; a proxy for stellar mass) than the "old" early-type galaxies NGC 3115 and 3379 (~8-10 Gyr), which results in a factor of ~2-3 excess of LX/LK for NGC 3384. This result is consistent with the X-ray binary population synthesis model predictions; however, our small galaxy sample size does not allow us to draw definitive conclusions on the evolution field LMXBs in general. We discuss how future surveys of larger galaxy samples that combine deep Chandra and HST data could provide a powerful new benchmark for calibrating X-ray binary population synthesis models., Accepted for publication in ApJ (16 pages, 9 figures, 6 tables)

We report on simultaneous observations of the local starburst system Arp 299 with NuSTAR and Chandra, which provides the first resolved images of this galaxy up to energies of ~ 45 keV. Fitting the 3-40 keV spectrum reveals a column density of $N_{\rm H}$ ~ 4 x10^{24} cm^{-2}, characteristic of a Compton-thick AGN, and a 10-30 keV luminosity of 1.2x 10^{43} ergs s^{-1}. The hard X-rays detected by NuSTAR above 10 keV are centered on the western nucleus, Arp 299-B, which previous X-ray observations have shown to be the primary source of neutral Fe-K emission. Other X-ray sources, including Arp 299-A, the eastern nucleus which is also thought to harbor an AGN, as well as X-ray binaries, contribute $\lesssim 10%$ to the 10-20 keV emission from the Arp 299 system. The lack of significant emission above 10 keV other than that attributed to Arp 299-B suggests that: a) any AGN in Arp 299-A must be heavily obscured ($N_{\rm H}$ > 10^{24} cm^{-2}) or have a much lower luminosity than Arp 299-B and b) the extranuclear X-ray binaries have spectra that cut-off above ~10 keV. Such soft spectra are characteristic of ultraluminous X-ray (ULX) sources observed to date by NuSTAR., 9 pages; accepted for publication in Astrophysical Journal

We present results from a Chandra study of ultraluminous X-ray sources (ULXs) in a sample of 17 nearby (D_L, Accepted for publication in MNRAS; 26 pages, 6 figures, 12 tables