BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Andreas Zezas

Andreas Zezas appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-ph.GAastro-ph.COastro-ph.IMastro-phastro-ph.SRApplicationsastro-ph.EPphysics.data-an

Catalog footprint

What is connected

46works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

46 published item(s)

preprint2026arXiv

The LIRA-Ising Model: Estimating the boundaries of irregularly shaped X-ray sources

Mapping the boundary of an extended source is a key step in the study of its morphology. The background contamination and statistical fluctuations of typical astronomical images make this a challenging statistical task, particularly for X-ray images with low surface brightness. We develop a three-step Bayesian procedure to identify the boundaries of irregularly shaped sources. We first apply a Bayesian multiscale reconstruction algorithm known as LIRA to obtain posterior pixelwise probability distributions of the source intensity that properly account for known structures, astrophysical background, and the effect of the telescope point spread function. Next, we adopt an Ising model to group pixels with similar intensities into cohesive regions corresponding to background and source. Finally, the boundary is derived on the basis of the most likely aggregation of pixels into the source region. Because the overall model combines LIRA and the Ising model, we call it LIRA-Ising. We verify the proposed method using a set of simulation studies. We then apply it to the Chandra X-ray Observatory images of two high redshift quasars, PKS J1421-0643 and 0730+257, to determine the extent and morphology of X-ray jets. Our method shows a uniform X-ray surface brightness of PKS J1421-0643 jet, and identifies knotty structure in the X-ray jet of 0730+257.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Spatial Structures in the Globular Cluster Distribution of Fornax Cluster Galaxies

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Next Generation X-ray Galaxy Survey with eROSITA

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Clarifying the population of HMXBs in the Small Magellanic Cloud

Almost all confirmed optical counterparts of HMXBs in the SMC are OB stars with equatorial decretion disks (OBe). These sources emit strongly in Balmer lines and standout when imaged through narrow-band Hα imaging. The lack of secure counterparts for a significant fraction of the HMXBs motivated us to search for more. Using the catalogs for OB/OBe stars (Maravelias et al. 2017) and for HMXBs (Haberl & Sturm 2016) we detect 70 optical counterparts (out of 104 covered by our survey). We provide the first identification of the optical counterpart to the source XTEJ0050-731. We verify that 17 previously uncertain optical counterparts are indeed the proper matches. Regarding 52 confirmed HMXBs (known optical counterparts with Hα emission), we detect 39 as OBe and another 13 as OB stars. This allows a direct estimation of the fraction of active OBe stars in HMXBs that show Hα emission at a given epoch to be at least $\sim75\%$ of their total HMXB population.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Demonstrating the likely neutron star nature of five M31 globular cluster sources with Swift-NuSTAR spectroscopy

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Exploring the Overabundance of ULXs in Metal- and Dust-poor Local Lyman Break Analogs

We have studied high mass X-ray binary (HMXB) populations within two low-metallicity, starburst galaxies, Haro 11 and VV 114. These galaxies serve as analogs to high-redshift (z>2) Lyman break galaxies, and within the larger sample of Lyman break analogs (LBAs) are sufficiently nearby (<87 Mpc) to be spatially-resolved by Chandra. Previous studies of the X-ray emission in LBAs have found that the 2-10 keV luminosity per star formation rate (SFR) in these galaxies is elevated, potentially because of their low metallicities (12+log[O/H]= 8.3-8.4). Theoretically, the progenitors of XRBs forming in lower metallicity environments lose less mass from stellar winds over their lifetimes, producing more massive compact objects (i.e., neutron stars and black holes), and thus resulting in more numerous and luminous HMXBs per SFR. We have performed an in-depth study of the only two LBAs that have spatially-resolved 2-10 keV emission with Chandra to present the bright end of the X-ray luminosity distribution of HMXBs (L$_X>10^{39}$ erg/s; ultraluminous X-ray sources, ULXs) in these low-metallicity galaxies, based on 8 detected ULXs. Comparing with the star-forming galaxy X-ray luminosity function (XLF) presented by Mineo et al. (2012), Haro 11 and VV 114 host ~4 times more L$_X>10^{40}$ erg/s sources than expected given their SFRs. We simulate the effects of source blending from crowded lower luminosity HMXBs using the star-forming galaxy XLF and then vary the XLF shapes until we reproduce the observed point source luminosity distributions. We find that these LBAs have a shallower bright end slope than the standard XLF. If we conservatively assume that the brightest X-ray source from each galaxy is powered by an AGN rather than a HMXB and eliminate these sources from consideration, the luminosity distribution becomes poorly constrained but does appear to be consistent with a standard XLF. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Star-formation history and X-ray binary populations: the case of the Large Magellanic Cloud

In the present work we investigate the link between high-mass X-ray binaries (HMXBs) and star formation in the Large Magellanic Cloud (LMC), our nearest star-forming galaxy. Using optical photometric data, we identify the most likely counterpart of 44 X-ray sources. Among the 40 HMXBs classified in this work, we find 33 Be/X-ray binaries, and 4 supergiant XRBs. Using this census and the published spatially resolved star-formation history map of the LMC, we find that the HMXBs (and as expected the X-ray pulsars) are present in regions with star-formation bursts $\sim$6-25 Myr ago, in contrast to the Small Magellanic Cloud (SMC), for which this population peaks at later ages ($\sim$25-60 Myr ago). We also estimate the HMXB production rate to be equal to 1 system per $\sim43.5\times10^{-3}$ Mo/yr, or 1 system per $\sim$143 Mo of stars formed during the associated star-formation episode. Therefore, the formation efficiency of HMXBs in the LMC is $\sim$17 times lower than that in the SMC. We attribute this difference primarily in the different ages and metallicity of the HMXB populations in the two galaxies. We also set limits on the kicks imparted on the neutron star during the supernova explosion. We find that the time elapsed since the supernova kick is $\sim$3 times shorter in the LMC than the SMC. This in combination with the average offsets of the HMXBs from their nearest star clusters results in $\sim$4 times faster transverse velocities for HMXBs in the LMC than in the SMC.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

SXP214, an X-ray Pulsar in the Small Magellanic Cloud, Crossing the Circumstellar Disk of the Companion

Located in the Small Magellanic Cloud (SMC), SXP214 is an X-ray pulsar in a high mass X-ray binary system with a Be-star companion. A recent survey of the SMC under a Chandra X-ray Visionary program found the source in a transition when the X-ray flux was on a steady rise. The Lomb-Scargle periodogram revealed a pulse period of 211.49 +/- 0.42 s, which is significantly (>5sigma) shorter than the previous measurements with XMM-Newton and RXTE. This implies that the system has gone through sudden spin-up episodes recently. The pulse profile shows a sharp eclipse-like feature with a modulation amplitude of >95%. The linear rise of the observed X-ray luminosity from <~2x to 7x10^35 erg s^-1 is correlated with steady softening of the X-ray spectrum, which can be described by the changes in the local absorption from N_H ~ 10^24 to <~10^20 cm^-2 for an absorbed power-law model. The soft X-ray emission below 2 keV was absent in the early part of the observation when only the pulsating hard X-ray component was observed, whereas at later times both soft and hard X-ray components were observed pulsating. A likely explanation is that the neutron star was initially hidden in the circumstellar disk of the companion, and later came out of the disk with the accreted material that continued fueling the observed pulsation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Merger Signatures in the Dynamics of Star-forming Gas

Spatially resolved kinematics have been used to determine the dynamical status of star-forming galaxies with ambiguous morphologies, and constrain the importance of galaxy interactions during the assembly of galaxies. However, measuring the importance of interactions or galaxy merger rates requires knowledge of the systematics in kinematic diagnostics and the visible time with merger indicators. We analyze the dynamics of star-forming gas in a set of binary merger hydrodynamic simulations with stellar mass ratios of 1:1 and 1:4. We find that the evolution of kinematic asymmetries traced by star-forming gas mirrors morphological asymmetries derived from mock optical images, in which both merger indicators show the largest deviation from isolated disks during strong interaction phases. Based on a series of simulations with various initial disk orientations, orbital parameters, gas fractions, and mass ratios, we find that the merger signatures are visible for ~0.2-0.4 Gyr with kinematic merger indicators but can be approximately twice as long for equal-mass mergers of massive gas-rich disk galaxies designed to be analogs of z~2-3 submillimeter galaxies. Merger signatures are most apparent after the second passage and before the black holes coalescence, but in some cases they persist up to several hundred Myr after coalescence. About 20-60% of the simulated galaxies are not identified as mergers during the strong interaction phase, implying that galaxies undergoing violent merging process do not necessarily exhibit highly asymmetric kinematics in their star-forming gas. The lack of identifiable merger signatures in this population can lead to an underestimation of merger abundances in star-forming galaxies, and including them in samples of star-forming disks may bias the measurements of disk properties such as intrinsic velocity dispersion.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Simulating the sensitivity to stellar point sources of Chandra X-ray observations

The Chandra Cygnus OB2 Legacy Survey is a wide and deep X-ray survey of the nearby and massive Cygnus OB2 association. The survey has detected ~8,000 X-ray sources, the majority of which are pre-main sequence X-ray emitting young stars in the association itself. To facilitate quantitative scientific studies of these sources as well as the underlying OB association it is important to understand the sensitivity of the observations and the level of completeness the observations have obtained. Here we describe the use of a hierarchical Monte Carlo simulation to achieve this goal by combining the empirical properties of the observations, analytic estimates of the source verification process, and an extensive set of source detection simulations. We find that our survey reaches a 90% completeness level for a pre-main-sequence population at the distance of Cyg OB2 at an X-ray luminosity of 4 x 10^30 ergs/s and a stellar mass of 1.3 Msun for a randomly distributed population. For a spatially clustered population such as Cyg~OB2 the 90% completeness level is reached at 1.1 Msun instead, as the sources are more concentrated in areas of our survey with a high exposure. These simulations can easily be adapted for use with other X-ray observations and surveys, and we provide X-ray detection efficiency curves for a very wide array of source and background properties to allow these simulations to be easily exploited by other users.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spectral and temporal properties of the ultra-luminous X-ray pulsar in M82 from 15 years of Chandra observations and analysis of the pulsed emission using NuSTAR

The recent discovery by Bachetti et al. (2014) of a pulsar in M82 that can reach luminosities of up to 10^40 ergs s^-1, a factor of ~100 the Eddington luminosity for a 1.4 Msol compact object, poses a challenge for accretion physics. In order to better understand the nature of this source and its duty cycle, and in the light of several physical models that have been subsequently published, we conduct a spectral and temporal analysis of the 0.5-8 keV X-ray emission from this source from 15 years of Chandra observations. We fit the Chandra spectra of the pulsar with a power-law model and a disk black body model, subjected to interstellar absorption in M82. We carefully assess for the effect of pile-up in our observations, where 4/19 observations have a pile-up fraction >10%, which we account for during spectral modeling with a convolution model. When fitted with a power-law model, the average photon index when the source is at high luminosity (L_X>10^39 ergs s^-1) is Gamma=1.33+/-0.15. For the disk black body model, the average temperature is T=3.24+/-0.65 keV, consistent with other luminous X-ray pulsars. We also investigated the inclusion of a soft excess component and spectral break, finding that the spectra are also consistent with these features common to luminous X-ray pulsars. In addition, we present spectral analysis from NuSTAR over the 3-50 keV range where we have isolated the pulsed component. We find that the pulsed emission in this band is best fit by a power-law with a high-energy cut-off, where Gamma=0.6+/-0.3 and E_C=14^{+5}_{-3} keV. While the pulsar has previously been identified as a transient, we find from our longer-baseline study that it has been remarkably active over the 15-year period, where for 9/19 (47%) observations that we analyzed, the pulsar appears to be emitting at a luminosity in excess of 10^39 ergs s^-1, greater than 10 times its Eddington limit.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Variations of the ISM Compactness Across the Main Sequence of Star-Forming Galaxies: Observations and Simulations

(abridged) The majority of star-forming galaxies follow a simple empirical correlation in the star formation rate (SFR) versus stellar mass ($M_*$) plane, usually referred to as the star formation Main Sequence (MS). Here we combine a set of hydro-dynamical simulations of interacting galactic disks with state-of-the-art radiative transfer codes to analyze how the evolution of mergers is reflected upon the properties of the MS. We present \textsc{Chiburst}, a Markov Chain Monte Carlo (MCMC) Spectral Energy Distribution (SED) code that fits the multi-wavelength, broad-band photometry of galaxies and derives stellar masses, star formation rates, and geometrical properties of the dust distribution. We apply this tool to the SEDs of simulated mergers and compare the derived results with the reference output from the simulations. Our results indicate that changes in the SEDs of mergers as they approach coalescence and depart from the MS are related to an evolution of dust geometry in scales larger than a few hundred parsecs. This is reflected in a correlation between the specific star formation rate (sSFR), and the compactness parameter $\mathcal{C}$, that parametrizes this geometry and hence the evolution of dust temperature ($T_{\rm{dust}}$) with time. As mergers approach coalescence, they depart from the MS and increase their compactness, which implies that moderate outliers of the MS are consistent with late-type mergers. By further applying our method to real observations of Luminous Infrared Galaxies (LIRGs), we show that the merger scenario is unable to explain these extreme outliers of the MS. Only by significantly increasing the gas fraction in the simulations are we able to reproduce the SEDs of LIRGs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A new candidate Wolf-Rayet X-ray binary in NGC 253

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Exploring the Small Magellanic Cloud to the Faintest X-ray Fluxes: Source Catalog, Timing and Spectral Analysis

We present the results of a pair of 100 ksec Chandra observations in the Small Magellanic Cloud to survey high mass X-ray binaries (HMXBs), stars and LMXBs/CVs down to Lx = 4.3 x 10^32 erg/s The two SMC Deep Fields are located in the most active star forming region of the bar, with Deep Field-1 positioned at the most pulsar-rich location identified from previous surveys. Two new pulsars were discovered in outburst: CXOU J004929.7-731058 (P=892s), CXOU J005252.2-721715 (P=326s), and 3 new HMXB candidates were identified. Of 15 Be-pulsars now known in the field, 13 were detected, with pulsations seen in 9 of them. Ephemerides demonstrate that 6 of the 10 pulsars known to exhibit regular outbursts were seen outside their periastron phase, and quiescent X-ray emission at Lx=10^33 - 10^34 is shown to be common. Comparison with ROSAT, ASCA, XMM-Newton catalogs resulted in positive identification of several previously ambiguous sources. Bright optical counterparts exist for 40 of the X-ray sources, of which 33 are consistent with early-type stars Mv<-2, B-V<0.2), and are the subject of a companion paper. The results point to an underlying HMXB population-density up to double that of active systems. The full catalog of 394 point-sources is presented along with detailed analyses of timing and spectral properties.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulated Galaxy Interactions as Probes of Merger Spectral Energy Distributions

We present the first systematic comparison of ultraviolet-millimeter spectral energy distributions (SEDs) of observed and simulated interacting galaxies. Our sample is drawn from the Spitzer Interacting Galaxy Survey, and probes a range of galaxy interaction parameters. We use 31 galaxies in 14 systems which have been observed with Herschel, Spitzer, GALEX, and 2MASS. We create a suite of GADGET-3 hydrodynamic simulations of isolated and interacting galaxies with stellar masses comparable to those in our sample of interacting galaxies. Photometry for the simulated systems is then calculated with the SUNRISE radiative transfer code for comparison with the observed systems. For most of the observed systems, one or more of the simulated SEDs match reasonably well. The best matches recover the infrared luminosity and the star formation rate of the observed systems, and the more massive systems preferentially match SEDs from simulations of more massive galaxies. The most morphologically distorted systems in our sample are best matched to simulated SEDs close to coalescence, while less evolved systems match well with SEDs over a wide range of interaction stages, suggesting that an SED alone is insufficient to identify interaction stage except during the most active phases in strongly interacting systems. This result is supported by our finding that the SEDs calculated for simulated systems vary little over the interaction sequence.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spatially Resolving a Starburst Galaxy at Hard X-ray Energies: NuSTAR, Chandra, AND VLBA Observations of NGC 253

Prior to the launch of NuSTAR, it was not feasible to spatially resolve the hard (E > 10 keV) emission from galaxies beyond the Local Group. The combined NuSTAR dataset, comprised of three ~165 ks observations, allows spatial characterization of the hard X-ray emission in the galaxy NGC 253 for the first time. As a follow up to our initial study of its nuclear region, we present the first results concerning the full galaxy from simultaneous NuSTAR, Chandra, and VLBA monitoring of the local starburst galaxy NGC 253. Above ~10 keV, nearly all the emission is concentrated within 100" of the galactic center, produced almost exclusively by three nuclear sources, an off-nuclear ultraluminous X-ray source (ULX), and a pulsar candidate that we identify for the first time in these observations. We detect 21 distinct sources in energy bands up to 25 keV, mostly consisting of intermediate state black hole X-ray binaries. The global X-ray emission of the galaxy - dominated by the off-nuclear ULX and nuclear sources, which are also likely ULXs - falls steeply (photon index >~ 3) above 10 keV, consistent with other NuSTAR-observed ULXs, and no significant excess above the background is detected at E > 40 keV. We report upper limits on diffuse inverse Compton emission for a range of spatial models. For the most extended morphologies considered, these hard X-ray constraints disfavor a dominant inverse Compton component to explain the γ-ray emission detected with Fermi and H.E.S.S. If NGC 253 is typical of starburst galaxies at higher redshift, their contribution to the E > 10 keV cosmic X-ray background is < 1%.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The total infrared luminosity may significantly overestimate the star formation rate of recently quenched galaxies

The total infrared (IR) luminosity is very useful for estimating the star formation rate (SFR) of galaxies, but converting the IR luminosity into an SFR relies on assumptions that do not hold for all galaxies. We test the effectiveness of the IR luminosity as an SFR indicator by applying it to synthetic spectral energy distributions generated from three-dimensional hydrodynamical simulations of isolated disc galaxies and galaxy mergers. In general, the SFR inferred from the IR luminosity agrees well with the true instantaneous SFR of the simulated galaxies. However, for the major mergers in which a strong starburst is induced, the SFR inferred from the IR luminosity can overestimate the instantaneous SFR during the post-starburst phase by greater than two orders of magnitude. Even though the instantaneous SFR decreases rapidly after the starburst, the stars that were formed in the starburst remain dust-obscured and thus produce significant IR luminosity. Consequently, use of the IR luminosity as an SFR indicator may cause one to conclude that post-starburst galaxies are still star-forming, whereas in reality, star formation was recently quenched.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Energy Feedback from X-ray Binaries in the Early Universe

X-ray photons, because of their long mean-free paths, can easily escape the galactic environments where they are produced, and interact at long distances with the inter-galactic medium, potentially having a significant contribution to the heating and reionization of the early Universe. The two most important sources of X-ray photons in the Universe are active galactic nuclei (AGN) and X-ray binaries (XRBs). In this Letter we use results from detailed, large scale population synthesis simulations to study the energy feedback of XRBs, from the first galaxies (z~ 20) until today. We estimate that X-ray emission from XRBs dominates over AGN at z>6-8. The shape of the spectral energy distribution of the emission from XRBs shows little change with redshift, in contrast to its normalization which evolves by ~4 orders of magnitude, primarily due to the evolution of the cosmic star-formation rate. However, the metallicity and the mean stellar age of a given XRB population affect significantly its X-ray output. Specifically, the X-ray luminosity from high-mass XRBs per unit of star-formation rate varies an order of magnitude going from solar metallicity to less than 10% solar, and the X-ray luminosity from low-mass XRBs per unit of stellar mass peaks at an age of ~300 Myr and then decreases gradually at later times, showing little variation for mean stellar ages > 3 Gyr. Finally, we provide analytical and tabulated prescriptions for the energy output of XRBs, that can be directly incorporated in cosmological simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fast and Furious: Shock Heated Gas as the Origin of Spatially Resolved Hard X-ray Emission in the Central 5 kpc of the Galaxy Merger NGC 6240

We have obtained a deep, sub-arcsecond resolution X-ray image of the nuclear region of the luminous galaxy merger NGC 6240 with Chandra, which resolves the X-ray emission from the pair of active nuclei and the diffuse hot gas in great detail. We detect extended hard X-ray emission from kT~6 keV (~70 million K) hot gas over a spatial scale of 5 kpc, indicating the presence of fast shocks with velocity of ~2200 km/s. For the first time we obtain the spatial distribution of this highly ionized gas emitting FeXXV, which shows a remarkable correspondence to the large scale morphology of H_2(1-0) S(1) line emission and Hαfilaments. Propagation of fast shocks originated in the starburst driven wind into the ambient dense gas can account for this morphological correspondence. With an observed L(0.5-8 keV)=5.3E+41 erg/s, the diffuse hard X-ray emission is 100 times more luminous than that observed in the classic starburst galaxy M82. Assuming a filling factor of 1% for the 70 MK temperature gas, we estimate its total mass (M_{hot}=1.8E+8 Msun) and thermal energy (E_{th}=6.5E+57 ergs). The total iron mass in the highly ionized plasma is M_{Fe}=4.6E+5 Msun. Both the energetics and the iron mass in the hot gas are consistent with the expected injection from the supernovae explosion during the starburst that is commensurate with its high star formation rate. No evidence for fluorescent Fe I emission is found in the CO filament connecting the two nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Global Star Formation Rates and Dust Emission Over the Galaxy Interaction Sequence

We measured and modeled the spectral energy distributions (SEDs) in 28 bands from the ultraviolet to the far-infrared (FIR) for 31 interacting galaxies in 14 systems. The sample is drawn from the Spitzer Interacting Galaxy Survey, which probes a range of galaxy interaction parameters at multiple wavelengths with an emphasis on the infrared bands. The subset presented in this paper consists of all galaxies for which FIR Herschel SPIRE observations are publicly available. Our SEDs combine the Herschel photometry with multi-wavelength data from Spitzer, GALEX, Swift UVOT, and 2MASS. While the shapes of the SEDs are broadly similar across our sample, strongly interacting galaxies typically have more mid-infrared emission relative to their near-infrared and FIR emission than weakly or moderately interacting galaxies. We modeled the full SEDs to derive host galaxy star formation rates (SFR), specific star formation rates (sSFR), stellar masses, dust temperatures, dust luminosities, and dust masses. We find increases in the dust luminosity and mass, SFR, and cold (15-25 K) dust temperature as the interaction progresses from moderately to strongly interacting and between non- interacting and strongly interacting galaxies. We also find increases in the SFR between weakly and strongly interacting galaxies. In contrast, the sSFR remains unchanged across all the interaction stages. The SFR derived from the SED modeling agrees well with rates estimated by proportionality relations that depend on infrared emission. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Optical spectra of 5 new Be/X-ray Binaries in the Small Magellanic Cloud and the link of the supergiant B[e] star LHA 115-S 18 with an X-ray source

The Small Magellanic Cloud (SMC) is well known to harbor a large number of High-Mass X-ray Binaries (HMXBs). The identification of their optical counterparts provides information on the nature of the donor stars and can help to constrain the parameters of these systems and their evolution. We obtained optical spectra for a number of HMXBs identified in previous Chandra and XMM-Newton surveys of the SMC using the AAOmega/2dF fiber-fed spectrograph at the Anglo-Australian Telescope. We find 5 new Be/X-ray binaries (BeXRBs; including a tentative one), by identifying the spectral type of their optical counterparts, and we confirm the spectral classification of an additional 15 known BeXRBs. We compared the spectral types, orbital periods, and eccentricities of the BeXRB populations in the SMC and the Milky Way and we find marginal evidence for difference between the spectral type distributions, but no statistically significant differences for the orbital periods and the eccentricities. Moreover, our search revealed that the well known supergiant B[e] star LHA 115-S 18 (or AzV 154) is associated with the weak X-ray source CXOU J005409.57-724143.5. We provide evidence that the supergiant star LHA 115-S 18 is the optical counterpart of the X-ray source, and we discuss different possibilities of the origin of its low X-ray luminosity (Lx ~ 4 x 10^33 erg/s).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Exceptional Soft X-ray Halo of the Galaxy Merger NGC 6240

We report on a recent ~150-ks long Chandra observation of the ultraluminous infrared galaxy merger NGC 6240, which allows a detailed investigation of the diffuse galactic halo. Extended soft X-ray emission is detected at the 3-sigma confidence level over a diamond-shaped region with projected physical size of ~110x80 kpc, and a single-component thermal model provides a reasonably good fit to the observed X-ray spectrum. The hot gas has a temperature of ~7.5 million K, an estimated density of 2.5x10^{-3} cm^{-3}, and a total mass of ~10^10 M_sun, resulting in an intrinsic 0.4-2.5 keV luminosity of 4x10^41 erg s^{-1}. The average temperature of 0.65 keV is quite high to be obviously related to either the binding energy in the dark-matter gravitational potential of the system or the energy dissipation and shocks following the galactic collision, yet the spatially-resolved spectral analysis reveals limited variations across the halo. The relative abundance of the main alpha-elements with respect to iron is several times the solar value, and nearly constant as well, implying a uniform enrichment by type II supernovae out to the largest scales. Taken as a whole, the observational evidence is not compatible with a superwind originated by a recent, nuclear starburst, but rather hints at widespread, enhanced star formation proceeding at steady rate over the entire dynamical timescale (~200 Myr). The preferred scenario is that of a starburst-processed gas component gently expanding into, and mixing with, a pre-existing halo medium of lower metallicity (Z ~ 0.1 solar) and temperature (kT ~ 0.25 keV). This picture cannot be probed more extensively with the present data, and the ultimate fate of the diffuse, hot gas remains uncertain. Under some favorable conditions, at least a fraction of it might be retained after the merger completion, and evolve into the hot halo of a young elliptical galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Hot and Energetic Universe: A White Paper presenting the science theme motivating the Athena+ mission

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Hot and Energetic Universe: The close environments of supermassive black holes

Most of the action in Active Galactic Nuclei (AGN) occurs within a few tens of gravitational radii from the supermassive black hole, where matter in the accretion disk may lose up to almost half of its energy with a copious production of X-rays, emitted via Comptonization of the disk photons by hot electrons in a corona and partly reflected by the accretion disk. Thanks to its large effective area and excellent energy resolution, Athena+ contributions in the understanding of the physics of accretion in AGN will be fundamental - and unique - in many respects. It will allow us to map the disk-corona system - which is crucial to understand the mechanism of energy extraction and the relation of the corona with winds and jets - by studying the time lags between reflected and primary photons. These lags have been recently discovered by XMM-Newton, but only Athena+ will have the sensitivity required to fully exploit this technique. Athena+ will also be able e.g. to determine robustly the spin of the black hole in nearby sources (and to extend these measurements beyond the local Universe), to establish the nature of the soft X-ray components, and to map the circumnuclear matter within the AGN inner parsec with unprecedented details.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

CfA4: Light Curves for 94 Type Ia Supernovae

We present multi-band optical photometry of 94 spectroscopically-confirmed Type Ia supernovae (SN Ia) in the redshift range 0.0055 to 0.073, obtained between 2006 and 2011. There are a total of 5522 light curve points. We show that our natural system SN photometry has a precision of roughly 0.03 mag or better in BVr'i', 0.06 mag in u', and 0.07 mag in U for points brighter than 17.5 mag and estimate that it has a systematic uncertainty of 0.014, 0.010, 0.012, 0.014, 0.046, and 0.073 mag in BVr'i'u'U, respectively. Comparisons of our standard system photometry with published SN Ia light curves and comparison stars reveal mean agreement across samples in the range of ~0.00-0.03 mag. We discuss the recent measurements of our telescope-plus-detector throughput by direct monochromatic illumination by Cramer et al (in prep.). This technique measures the whole optical path through the telescope, auxiliary optics, filters, and detector under the same conditions used to make SN measurements. Extremely well-characterized natural-system passbands (both in wavelength and over time) are crucial for the next generation of SN Ia photometry to reach the 0.01 mag accuracy level. The current sample of low-z SN Ia is now sufficiently large to remove most of the statistical sampling error from the dark energy error budget. But pursuing the dark-energy systematic errors by determining highly-accurate detector passbands, combining optical and near-infrared (NIR) photometry and spectra, using the nearby sample to illuminate the population properties of SN Ia, and measuring the local departures from the Hubble flow will benefit from larger, carefully measured nearby samples.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Deep Chandra Monitoring Observations of NGC 4649: II. Wide-Field Hubble Space Telescope Imaging of the Globular Clusters

We present g and z photometry and size estimates for globular clusters (GCs) in the massive Virgo elliptical NGC 4649 (M60) using a five-pointing Hubble Space Telescope/Advanced Camera for Surveys mosaic. The metal-poor GCs show a monotonic negative metallicity gradient of (-0.43 +/- 0.10) dex per dex in radius over the full radial range of the data, out to ~ 24 kpc. There is evidence for substantial color substructure among the metal-rich GCs. The metal-poor GCs have typical sizes ~ 0.4 pc larger than the metal-rich GCs out to large galactocentric distances (~> 20 kpc), favoring an intrinsic explanation for the size difference rather than projection effects. There is no clear relation between half-light radius and galactocentric distance beyond ~ 15 kpc, suggesting that the sizes of GCs are not generically set by tidal limitation. Finally, we identify ~ 20 candidate ultra-compact dwarfs that extend down to surprisingly faint absolute magnitudes (M_z ~ -8.5), and may bridge the gap between this class and "extended clusters" in the Local Group. Three of the brighter candidates have published radial velocities and can be confirmed as bona fide ultra-compact dwarfs; follow-up spectroscopy will determine the nature of the remainder of the candidates.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The X-ray Star Formation Story as Told by Lyman Break Galaxies in the 4 Ms CDF-S

We present results from deep X-ray stacking of >4000 high redshift galaxies from z~1 to 8 using the 4 Ms Chandra Deep Field South (CDF-S) data, the deepest X-ray survey of the extragalactic sky to date. The galaxy samples were selected using the Lyman break technique based primarily on recent HST ACS and WFC3 observations. Based on such high specific star formation rates (sSFRs): log SFR/M* > -8.7, we expect that the observed properties of these LBGs are dominated by young stellar populations. The X-ray emission in LBGs, eliminating individually detected X-ray sources (potential AGN), is expected to be powered by X-ray binaries and hot gas. We find, for the first time, evidence of evolution in the X-ray/SFR relation. Based on X-ray stacking analyses for z<4 LBGs (covering ~90% of the Universe's history), we find that the 2-10 keV X-ray luminosity evolves weakly with redshift (z) and SFR as log LX = 0.93 log (1+z) + 0.65 log SFR + 39.80. By comparing our observations with sophisticated X-ray binary population synthesis models, we interpret that the redshift evolution of LX/SFR is driven by metallicity evolution in HMXBs, likely the dominant population in these high sSFR galaxies. We also compare these models with our observations of X-ray luminosity density (total 2-10 keV luminosity per Mpc^3) and find excellent agreement. While there are no significant stacked detections at z>5, we use our upper limits from 5<z<8 LBGs to constrain the SMBH accretion history of the Universe around the epoch of reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

X-ray Binary Evolution Across Cosmic Time

High redshift galaxies permit the study of the formation and evolution of X-ray binary populations on cosmological timescales, probing a wide range of metallicities and star-formation rates. In this paper, we present results from a large scale population synthesis study that models the X-ray binary populations from the first galaxies of the universe until today. We use as input to our modeling the Millennium II Cosmological Simulation and the updated semi-analytic galaxy catalog by Guo et al. (2011) to self-consistently account for the star formation history and metallicity evolution of the universe. Our modeling, which is constrained by the observed X-ray properties of local galaxies, gives predictions about the global scaling of emission from X-ray binary populations with properties such as star-formation rate and stellar mass, and the evolution of these relations with redshift. Our simulations show that the X-ray luminosity density (X-ray luminosity per unit volume) from X-ray binaries in our Universe today is dominated by low-mass X-ray binaries, and it is only at z>2.5 that high-mass X-ray binaries become dominant. We also find that there is a delay of ~1.1 Gyr between the peak of X-ray emissivity from low-mass Xray binaries (at z~2.1) and the peak of star-formation rate density (at z~3.1). The peak of the X-ray luminosity from high-mass X-ray binaries (at z~3.9), happens ~0.8 Gyr before the peak of the star-formation rate density, which is due to the metallicity evolution of the Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Deep Chandra ACIS Study of NGC 4151. I. the X-ray Morphology of the 3 kpc-diameter Circum-nuclear Region and Relation to the Cold Interstellar Medium

We report on the imaging analysis of 200 ks sub-arcsecond resolution Chandra ACIS-S observations of the nearby Seyfert 1 galaxy NGC 4151. Bright, structured soft X-ray emission is observed to extend from 30 pc to 1.3 kpc in the south-west from the nucleus, much farther than seen in earlier X-ray studies. The terminus of the north-eastern X-ray emission is spatially coincident with a CO gas lane, where the outflow likely encounters dense gas in the host galactic disk. X-ray emission is also detected outside the boundaries of the ionization cone, which indicates that the gas there is not completely shielded from the nuclear continuum, as would be the case for a molecular torus collimating the bicone. In the central r<200 pc region, the subpixel processing of the ACIS data recovers the morphological details on scales of <30~pc (<0.5") first discovered in Chandra HRC images. The X-ray emission is more absorbed towards the boundaries of the ionization cone, as well as perpendicular to the bicone along the direction of a putative torus in NGC 4151. The innermost region where X-ray emission shows the highest hardness ratio, is spatially coincident with the near-infrared resolved H_2 emission and dusty spirals we find in an HST V-H color image. The agreement between the observed H_2 line flux and the value predicted from X-ray-irradiated molecular cloud models supports photo-excitation by X-rays from the active nucleus as the origin of the H_2 line, although contribution from UV fluorescence or collisional excitation cannot be fully ruled out with current data. The discrepancy between the mass of cold molecular gas inferred from recent CO and near-infrared H_2 observations may be explained by the anomalous CO abundance in this X-ray dominated region. The total H_2 mass derived from the X-ray observation agrees with measurement in Storchi-Bergmann et al.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Deep Chandra ACIS Study of NGC 4151. II. The Innermost Emission Line Region and Strong Evidence for Radio Jet-NLR Cloud Collision

We have studied the X-ray emission within the inner 150 pc radius of NGC 4151 by constructing high spatial resolution emission line images of OVII, OVIII, and NeIX. These maps show extended structures that are spatially correlated with the radio outflow and optical [OIII] emission. We find strong evidence for jet--gas cloud interaction, including morphological correspondences with regions of X-ray enhancement, peaks of near-infrared [FeII] emission, and optical clouds. In these regions, moreover, we find evidence of elevated NeIX/OVII ratios; the X-ray emission of these regions also exceeds that expected from nuclear photoionization. Spectral fitting reveals the presence of a collisionally ionized component. The thermal energy of the hot gas suggests that >0.1% of the estimated jet power is deposited into the host interstellar medium through interaction between the radio jet and the dense medium of the circum-nuclear region. We find possible pressure equilibrium between the collisionally ionized hot gas and the photoionized line-emitting cool clouds. We also obtain constraints on the extended iron and silicon fluorescent emission. Both lines are spatially unresolved. The upper limit on the contribution of an extended emission region to the Fe Kalpha emission is <5% of the total, in disagreement with a previous claim that 65% of the Fe Kalpha emission originates in the extended narrow line region.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Deep Chandra ACIS Study of NGC 4151. III. the Line Emission and Spectral Analysis of the Ionization Cone

This paper is the third in a series in which we present deep Chandra ACIS-S imaging spectroscopy of the Seyfert 1 galaxy NGC 4151, devoted to study its complex circum-nuclear X-ray emission. Emission features in the soft X-ray spectrum of the bright extended emission (L[0.3-2keV]~10^40 erg/s) at r>130 pc (2") are consistent with the brighter OVII, OVIII, and NeIX lines seen in the Chandra HETGS and XMM-Newton RGS spectra below 2 keV. We construct emission line images of these features and find good morphological correlations with the narrow line region clouds mapped in [OIII]5007A. Self-consistent photoionization models provide good descriptions of the spectra of the large scale emission, as well as resolved structures, supporting the dominant role of nuclear photoionization, although displacement of optical and X-ray features implies a more complex medium. Collisionally ionized emission is estimated to be <12% of the extended emission. Presence of both low and high ionization spectral components and extended emission in the X-ray image perpendicular to the bicone indicates leakage of nuclear ionization, likely filtered through warm absorbers, instead of being blocked by a continuous obscuring torus. The ratios of [OIII]/soft X-ray flux are approximately constant (~15) for the 1.5 kpc radius spanned by these measurements, indicating a relatively constant ionization parameter, consistent with the photoionized outflow of a wind-like density profile. Using spatially resolved features, we estimate that the mass outflow rate in NGC 4151 is ~2Msun/yr at 130 pc and the kinematic power of the ionized outflow is 1.7x10^41 erg/s, approximately 0.3% of the bolometric luminosity of NGC 4151.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Accounting for Calibration Uncertainties in X-ray Analysis: Effective Areas in Spectral Fitting

While considerable advance has been made to account for statistical uncertainties in astronomical analyses, systematic instrumental uncertainties have been generally ignored. This can be crucial to a proper interpretation of analysis results because instrumental calibration uncertainty is a form of systematic uncertainty. Ignoring it can underestimate error bars and introduce bias into the fitted values of model parameters. Accounting for such uncertainties currently requires extensive case-specific simulations if using existing analysis packages. Here we present general statistical methods that incorporate calibration uncertainties into spectral analysis of high-energy data. We first present a method based on multiple imputation that can be applied with any fitting method, but is necessarily approximate. We then describe a more exact Bayesian approach that works in conjunction with a Markov chain Monte Carlo based fitting. We explore methods for improving computational efficiency, and in particular detail a method of summarizing calibration uncertainties with a principal component analysis of samples of plausible calibration files. This method is implemented using recently codified Chandra effective area uncertainties for low-resolution spectral analysis and is verified using both simulated and actual Chandra data. Our procedure for incorporating effective area uncertainty is easily generalized to other types of calibration uncertainties.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hubble Space Telescope Photometry of Globular Clusters in M81

We perform aperture photometry and profile fitting on 419 globular cluster (GC) candidates with mV \leq 23 mag identified in Hubble Space Telescope Advanced Camera for Surveys BVI imaging, and estimate the effective radii of the clusters. We identify 85 previously known spectroscopically-confirmed clusters, and newly identify 136 objects as good cluster candidates within the 3σ color and size ranges defined by the spectroscopically confirmed clusters, yielding a total of 221 probable GCs. The luminosity function peak for the 221 probable GCs with estimated total dereddening applied is V ~(20.26 \pm 0.13) mag, corresponding to a distance of ~3.7\pm0.3 Mpc. The blue and red GC candidates, and the metal-rich (MR) and metal-poor (MP) spectroscopically confirmed clusters, are similar in half-light radius, respectively. Red confirmed clusters are about 6% larger in median half-light radius than blue confirmed clusters, and red and blue good GC candidates are nearly identical in half-light radius. The total population of confirmed and "good" candidates shows an increase in half-light radius as a function of galactocentric distance.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Interaction between the intergalactic medium and central radio source in the NGC 4261 group of galaxies

Using observations from the Chandra and XMM-Newton X-ray observatories, we examine the interaction between the intra-group medium and central radio source in the nearby NGC 4261 galaxy group. We confirm the presence of cavities associated with the radio lobes and estimate their enthalpy to be ~2.4x10^58 erg. The mechanical power output of the jets is >=10^43 erg/s, at least a factor of 60 greater than the cooling luminosity in the region the lobes inhabit. We identify rims of compressed gas enclosing the lobes, but find no statistically significant temperature difference between them and their surroundings, suggesting that the lobe expansion velocity is approximately sonic (Mach<=1.05). The apparent pressure of the radio lobes, based on the synchrotron minimum energy density argument, is a factor of 5 lower than that of the intra-group medium. Pressure balance could be achieved if entrainment of thermal gas provided additional non-radiating particles in the lobe plasma, but the energy required to heat these particles would be ~20 per cent. of the mechanical energy output of the radio source. NGC 4261 has a relatively compact cool core, which should probably be categorised as a galactic corona. The corona is capable of fuelling the active nucleus for considerably longer than the inferred source lifetime, but can be only inefficiently heated by the AGN or conduction. The expansion of the radio lobes has affected the structure of the gas in the galaxy, compressing and moving the material of the corona without causing significant shock heating, and expelling gas from the immediate neighbourhood of the jets. We discuss the possible implications of this environment for the duration of the AGN outburst, and consider mechanisms which might lead to the cessation of nuclear activity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The SEDs of Interacting Galaxies

The evolution of galaxies is greatly influenced by their interactions. As part of a program to study interacting galaxies, we have measured and modeled the spectral energy distri- butions (SEDs) from the ultraviolet (UV) to the far-infrared (FIR). We describe the constraints imposed on star formation histories by these SEDs, and the variations therein seen across the interaction sequence, and we compare the results of different star formation rate prescriptions applied to the data. The sample itself is based on the Spitzer Interacting Galaxy Survey (SIGS) of 111 galaxies in 50 systems, a project designed to probe a range of galaxy interaction parameters in the infrared. Our SEDs combine the Spitzer results with multiwavelength data from other missions, in particular GALEX and Herschel. The subset presented here is the sample for which FIR Herschel observations are currently publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A multi-wavelength study of Supernova Remnants in six nearby galaxies. I: Detection of new X-ray selected Supernova Remnants with Chandra

We present results from a study of the Supernova Remnant (SNR) population in a sample of six nearby galaxies (NGC 2403, NGC 3077, NGC 4214, NGC 4449, NGC 4395 and NGC 5204) based on Chandra archival data. We have detected 244 discrete X-ray sources down to a limiting flux of 10^{-15} erg/s. We identify 37 X-ray selected thermal SNRs based on their X-ray colors or spectra, 30 of which are new discoveries. In many cases the X-ray classification is confirmed based on counterparts with SNRs identified in other wavelengths. Three of the galaxies in our sample (NGC 4214, NGC 4395 and NGC 5204) are studied for the first time, resulting in the discovery of 13 thermal SNRs. We discuss the properties (luminosity, temperature, density) of the X-ray detected SNRs in the galaxies of our sample in order to address their dependence on their environment. We find that X-ray selected SNRs in irregular galaxies appear to be more luminous than those in spirals. We attribute this to the lower metalicities and therefore more massive progenitor stars of irregular galaxies or the higher local densities of the ISM. We also discuss the X-ray selected SNR populations in the context of the Star Formation Rate of their host galaxies. A comparison of the numbers of observed luminous X-ray selected SNRs with those expected based on the luminosity functions of X-ray SNRs in the MCs and M33 suggest different luminosity distributions between the SNRs in spiral and irregular galaxies with the latter tending to have flatter distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Extended X-ray Emission in the HI Cavity of NGC 4151: Galaxy-scale AGN Feedback?

We present the Chandra discovery of soft diffuse X-ray emission in NGC 4151 (L[0.5-2keV]~10^{39} erg s$^{-1}$), extending ~2 kpc from the active nucleus and filling in the cavity of the HI material. The best fit to the X-ray spectrum requires either a kT~0.25 keV thermal plasma or a photoionized component. In the thermal scenario, hot gas heated by the nuclear outflow would be confined by the thermal pressure of the HI gas and the dynamic pressure of inflowing neutral material in the galactic disk. In the case of photoionization, the nucleus must have experienced an Eddington limit outburst. For both scenarios, the AGN-host interaction in NGC 4151 must have occured relatively recently (some 10^4 yr ago). This very short timescale to the last episode of high activity phase may imply such outbursts occupy $\gtrsim$1% of AGN lifetime.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On Computing Upper Limits to Source Intensities

A common problem in astrophysics is determining how bright a source could be and still not be detected. Despite the simplicity with which the problem can be stated, the solution involves complex statistical issues that require careful analysis. In contrast to the confidence bound, this concept has never been formally analyzed, leading to a great variety of often ad hoc solutions. Here we formulate and describe the problem in a self-consistent manner. Detection significance is usually defined by the acceptable proportion of false positives (the TypeI error), and we invoke the complementary concept of false negatives (the TypeII error), based on the statistical power of a test, to compute an upper limit to the detectable source intensity. To determine the minimum intensity that a source must have for it to be detected, we first define a detection threshold, and then compute the probabilities of detecting sources of various intensities at the given threshold. The intensity that corresponds to the specified TypeII error probability defines that minimum intensity, and is identified as the upper limit. Thus, an upper limit is a characteristic of the detection procedure rather than the strength of any particular source and should not be confused with confidence intervals or other estimates of source intensity. This is particularly important given the large number of catalogs that are being generated from increasingly sensitive surveys. We discuss the differences between these upper limits and confidence bounds. Both measures are useful quantities that should be reported in order to extract the most science from catalogs, though they answer different statistical questions: an upper bound describes an inference range on the source intensity, while an upper limit calibrates the detection process. We provide a recipe for computing upper limits that applies to all detection algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Star formation history and X-ray binary populations: the case of the Small Magellanic Cloud

Using Chandra, XMM-Newton, and optical photometric catalogs we study the young X-ray binary (XRB) populations of the Small Magellanic Cloud. We find that the Be/X-ray binaries (Be-XRBs) are observed in regions with star formation rate bursts ~25-60 Myr ago. The similarity of this age with the age of maximum occurrence of the Be phenomenon (~40 Myr) indicates that the presence of a circumstellar decretion disk plays a significant role in the number of observed XRBs in the 10-100 Myr age range. We also find that regions with strong but more recent star formation (e.g., the Wing) are deficient in Be-XRBs. By correlating the number of observed Be-XRBs with the formation rate of their parent populations, we measure a Be-XRB production rate of ~1 system per 3 x 10^(-3) M$_{\odot}$/yr. Finally, we use the strong localization of the Be-XRB systems in order to set limits on the kicks imparted on the neutron star during the supernova explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Optical spectroscopy of 20 Be/X-ray Binaries in the Small Magellanic Cloud

We present a large sample (20 in total) of optical spectra of Small Magellanic Cloud (SMC) High-Mass X-ray Binaries obtained with the 2dF spectrograph at the Anglo-Australian Telescope. All of these sources are found to be Be/X-ray binaries (Be-XRBs), while for 5 sources we present original classifications. Several statistical tests on this expanded sample support previous findings for similar spectral-type distributions of Be-XRBs and Be field stars in the SMC, and of Be-XRBs in the Large Magellanic Cloud and the Milky Way, although this could be the result of small samples. On the other hand, we find that Be-XRBs follow a different distribution than Be stars in the Galaxy, also in agreement with previous studies. In addition, we find similar Be spectral type distributions between the Magellanic Clouds samples. These results reinforce the relation between the orbital period and the equivalent width of the Halpha line that holds for Be-XRBs. SMC Be stars have larger Halpha equivalent widths when compared to Be-XRBs, supporting the notion of circumstellar disk truncation by the compact object.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Chandra survey of the SMC "Bar": II. Optical counterparts of X-ray sources

We present the most likely optical counterparts of 113 X-ray sources detected in our Chandra survey of the central region of the Small Magellanic Cloud (SMC) based on the OGLE-II and MCPS catalogs. We estimate that the foreground contamination and chance coincidence probability are minimal for the bright optical counterparts (corresponding to OB type stars; 35 in total). We propose here for the first time 13 High-Mass X-ray Binaries (HMXBs), of which 4 are Be/X-ray binaries (Be-XRBs), and we confirm the previous classification of 18 Be-XRBs. We estimate that the new candidate Be-XRBs have an age of 15-85 Myr, consistent with the age of Be stars. We also examine the "overabundance" of Be-XRBs in the SMC fields covered by Chandra, in comparison with the Galaxy. In luminosities down to about 10^{34} erg/s, we find that SMC Be-XRBs are 1.5 times more common when compared to the Milky Way even after taking into account the difference in the formation rates of OB stars. This residual excess can be attributed to the lower metallicity of the SMC. Finally, we find that the mixing of Be-XRBs with other than their natal stellar population is not an issue in our comparisons of Be-XRBs and stellar populations in the SMC. Instead, we find indication for variation of the SMC XRB populations on kiloparsec scales, related to local variations of the formation rate of OB stars and slight variation of their age, which results in different relative numbers of Be stars and therefore XRBs.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Bayesian Estimation of Hardness Ratios: Modeling and Computations

A commonly used measure to summarize the nature of a photon spectrum is the so-called Hardness Ratio, which compares the number of counts observed in different passbands. The hardness ratio is especially useful to distinguish between and categorize weak sources as a proxy for detailed spectral fitting. However, in this regime classical methods of error propagation fail, and the estimates of spectral hardness become unreliable. Here we develop a rigorous statistical treatment of hardness ratios that properly deals with detected photons as independent Poisson random variables and correctly deals with the non-Gaussian nature of the error propagation. The method is Bayesian in nature, and thus can be generalized to carry out a multitude of source-population--based analyses. We verify our method with simulation studies, and compare it with the classical method. We apply this method to real world examples, such as the identification of candidate quiescent Low-mass X-ray binaries in globular clusters, and tracking the time evolution of a flare on a low-mass star.

0 citations0 reviewsNo code linkedOpen access