Search

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

We analyze the physical properties of 8 X-ray selected active galactic nuclei (AGN) and one candidate protoquasar system (ADF22A1) in the $z = 3.09$ SSA22 protocluster by fitting their X-ray-to-IR spectral energy distributions (SEDs) using our SED fitting code, Lightning. We recover star formation histories (SFH) for 7 of these systems which are well-fit by composite stellar population plus AGN models. We find indications that 4/9 of the SSA22 AGN systems we study have host galaxies below the main sequence, with $\rm SFR/SFR_{MS} \leq -0.4$. The remaining SSA22 systems, including ADF22A1, are consistent with obscured supermassive black hole (SMBH) growth in star forming galaxies. We estimate the SMBH accretion rates and masses, and compare the properties and SFH of the 9 protocluster AGN systems with X-ray detected AGN candidates in the Chandra Deep Fields (CDF), finding that the distributions of SMBH growth rates, star formation rates, SMBH masses, and stellar masses for the protocluster AGN are consistent with field AGN. We constrain the ratio between the sample-averaged SSA22 SMBH mass and CDF SMBH mass to $<1.41$. While the AGN are located near the density peaks of the protocluster, we find no statistically significant trends between the AGN or host galaxy properties and their location in the protocluster. We interpret the similarity of the protocluster and field AGN populations together with existing results as suggesting that the protocluster and field AGN co-evolve with their hosts in the same ways, while AGN-triggering events are more likely in the protocluster., Comment: 27 pages, 14 figures, 8 tables. Updated with published version

We present an updated version of Lightning, a galaxy spectral energy distribution (SED) fitting code that can model X-ray to submillimeter observations. The models in Lightning include the options to contain contributions from stellar populations, dust attenuation and emission, and active galactic nuclei (AGN). X-ray emission, when utilized, can be modeled as originating from stellar compact binary populations with the option to include emission from AGN. We have also included a variety of algorithms to fit the models to observations and sample parameter posteriors; these include an adaptive Markov-Chain Monte-Carlo (MCMC), affine-invariant MCMC, and Levenberg-Marquardt gradient decent (MPFIT) algorithms. To demonstrate some of the capabilities of Lightning, we present several examples using a variety of observational data. These examples include (1) deriving the spatially resolved stellar properties of the nearby galaxy M81, (2) demonstrating how X-ray emission can provide constrains on the properties of the supermassive black hole of a distant AGN, (3) exploring how to rectify the attenuation effects of inclination on the derived the star formation rate of the edge-on galaxy NGC 4631, (4) comparing the performance of Lightning to similar Bayesian SED fitting codes when deriving physical properties of the star-forming galaxy NGC 628, and (5) comparing the derived X-ray and UV-to-IR AGN properties from Lightning and CIGALE for a distant AGN. Lightning is an open-source application developed in the Interactive Data Language (IDL) and is available at https://github.com/rafaeleufrasio/lightning., Comment: 34 pages, 17 figures. Updated with published version

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

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., Comment: Accepted for publication in ApJ (25 pages, 9 figures, and 4 tables)

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., Comment: 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

Star-forming dwarf galaxies have properties similar to those expected in high-redshift galaxies. Hence, these local galaxies may provide insights into the evolution of the first galaxies, and the physical processes at work. We present a sample of eleven potential local analogs to high-$z$ (LAHz) galaxies. The sample consists of blue compact dwarf galaxies, selected to have spectral energy distributions that fit galaxies at $1.5

During solar flares, a large flux of energetic electrons propagate from the tops of reconnecting magnetic flux tubes toward the lower atmosphere. Over the course of the electrons' transport, a co-spatial counter-streaming return current is induced, thereby balancing the current density. In response to the return current electric field, a fraction of the ambient electrons will be accelerated into the runaway regime. However, models describing the accelerated electron beam/return-current system have generally failed to take these suprathermal runaway electrons into account self-consistently. We develop a model in which an accelerated electron beam drives a steady-state, sub-Dreicer co-spatial return-current electric field, which locally balances the direct beam current and freely accelerates a fraction of background (return-current) electrons. The model is self-consistent, i.e., the electric field induced by the co-evolution of the direct beam and the runaway current is considered. We find that (1) the return current electric field can return a significant number of suprathermal electrons to the acceleration region, where they can be further accelerated to higher energies, runaway electrons can be a few tens of percent of the return current flux returning to the nonthermal beam's acceleration region, (2) the energy gain of the suprathermal electrons can be up to $10-35$\,keV, (3) the heating rate in the corona can be reduced by an order of magnitude in comparison to models which neglect the runaway component. The results depend on the injected beam flux density, the temperature and density of the background plasma., Comment: 19 pages, 10 figures, submitted to ApJ

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

We 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., Comment: Accepted for publication in ApJS. Data products and catalogs available at https://lehmer.uark.edu/downloads/

The density-wave theory of spiral structure proposes that star formation occurs in or near a spiral-shaped region of higher density that rotates rigidly within the galactic disk at a fixed pattern speed. In most interpretations of this theory, newborn stars move downstream of this position as they come into view, forming a downstream spiral which is tighter, with a smaller pitch angle than that of the density wave itself. Rival theories, including theories which see spiral arms as essentially transient structures, may demand that pitch angle should not depend on wavelength. We measure the pitch angle of a large sample of galaxies at several wavelengths associated with star formation or very young stars (8.0 {\mu}m, H-{\alpha} line and 151 nm in the far-UV) and show that they all have the same pitch angle, which is larger than the pitch angle measured for the same galaxies at optical and near-infrared wavelengths. Our measurements in the B band and at 3.6 {\mu}m have unambiguously tighter spirals than the starforming wavelengths. In addition, we have measured in the u-band, which seems to fall midway between these two extremes. Thus, our results are consistent with a region of enhanced stellar light situated downstream of a starforming region., Comment: 12 Pages, 8 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., Comment: Accepted for publication in ApJS (37 pages, 16 figures, and 7 tables). Data products available (https://lehmer.uark.edu/downloads/)

We investigate geodesic orbits and manifolds for metrics associated with Schwarzschild geometry, considering space and time curvatures separately. For `a-temporal' space, we solve a central geodesic orbit equation in terms of elliptic integrals. The intrinsic geometry of a two-sided equatorial plane corresponds to that of a full Flamm's paraboloid. Two kinds of geodesics emerge. Both kinds may or may not encircle the hole region any number of times, crossing themselves correspondingly. Regular geodesics reach a periastron greater than the $r_S$ Schwarzschild radius, thus remaining confined to a half of Flamm's paraboloid. Singular or $s$-geodesics tangentially reach the $r_S$ circle. These $s$-geodesics must then be regarded as funneling through the `belt' of the full Flamm's paraboloid. Infinitely many geodesics can possibly be drawn between any two points, but they must be of specific regular or singular types. A precise classification can be made in terms of impact parameters. Geodesic structure and completeness is conveyed by computer-generated figures depicting either Schwarzschild equatorial plane or Flamm's paraboloid. For the `curved-time' metric, devoid of any spatial curvature, geodesic orbits have the same apsides as in Schwarzschild space-time. We focus on null geodesics in particular. For the limit of light grazing the sun, asymptotic `spatial bending' and `time bending' become essentially equal, adding up to the total light deflection of 1.75 arc-seconds predicted by general relativity. However, for a much closer approach of the periastron to $r_S$, `time bending' largely exceeds `spatial bending' of light, while their sum remains substantially below that of Schwarzschild space-time., Comment: 14 pages, 19 figures

We present constraints on variations in the initial mass function (IMF) of nine local early-type galaxies based on their low mass X-ray binary (LMXB) populations. Comprised of accreting black holes and neutron stars, these LMXBs can be used to constrain the important high mass end of the IMF. We consider the LMXB populations beyond the cores of the galaxies ($>0.2R_{e}$; covering $75-90\%$ of their stellar light) and find no evidence for systematic variations of the IMF with velocity dispersion ($\sigma$). We reject IMFs which become increasingly bottom heavy with $\sigma$, up to steep power-laws (exponent, $\alpha>2.8$) in massive galaxies ($\sigma>300$km/s), for galactocentric radii $>1/4\ R_{e}$. Previously proposed IMFs that become increasingly bottom heavy with $\sigma$ are consistent with these data if only the number of low mass stars $(<0.5M_{\odot}$) varies. We note that our results are consistent with some recent work which proposes that extreme IMFs are only present in the central regions of these galaxies. We also consider IMFs that become increasingly top-heavy with $\sigma$, resulting in significantly more LMXBs. Such a model is consistent with these observations, but additional data are required to significantly distinguish between this and an invariant IMF. For six of these galaxies, we directly compare with published IMF mismatch parameters from the Atlas3D survey, $\alpha_{dyn}$. We find good agreement with the LMXB population if galaxies with higher $\alpha_{dyn}$ have more top-heavy IMFs -- although we caution that our sample is quite small. Future LMXB observations can provide further insights into the origin of $\alpha_{dyn}$ variations., Comment: 13 pages, 7 figures, accepted for publication in ApJ

Multiwavelength data are essential in order to provide a complete picture of galaxy evolution and to inform studies of galaxies' morphological properties across cosmic time. Here we present results of a multiwavelength investigation of the morphologies of "tadpole" galaxies at intermediate redshift (0.314

NGC 6872, hereafter the Condor, is a large spiral galaxy that is interacting with its closest companion, the S0 galaxy IC 4970. The extent of the Condor provides an opportunity for detailed investigation of the impact of the interaction on the current star formation rate and its history across the galaxy, on the age and spatial distribution of its stellar population, and on the mechanism that drive the star formation activity. To address these issues we analyzed the far-ultraviolet (FUV) to near-infrared (near-IR) spectral energy distribution (SED) of 17, 10 kpc diameter, regions across the galaxy, and derived their star formation history, current star formation rate, and stellar population and mass. We find that most of the star formation takes place in the extended arms, with very little star formation in the central 5 kpc of the galaxy, in contrast to what was predicted from previous numerical simulations. There is a trend of increasing star formation activity with distance from the nucleus of the galaxy, and no evidence for a recent increase in the current star formation rate (SFR) due to the interaction. The nucleus itself shows no significant current star formation activity. The extent of the Condor also provides an opportunity to test the applicability of a single standard prescription for conversion of the FUV + IR (22 micron) intensities to a star formation rate for all regions. We find that the conversion factor differs from region to region, arising from regional differences in the stellar populations., Comment: 14 pages. 10 figures. Accepted for publication in ApJ