Source author record

Joel N. Bregman

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

53works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

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

53 published item(s)

preprint2022arXiv

Absorption Line Search Through Three Local Group Dwarf Galaxy Halos

Dwarf galaxies are missing nearly all of their baryons and metals from the stellar disk, presumed to be in a bound halo or expelled beyond the virial radius. The virial temperature for galaxies with $M_{\rm h} \sim 10^9 - 10^{10}$ $M_{\odot}$ is similar to the collisional ionization equilibrium temperature for the C IV ion. We searched for UV absorption from C IV in six sightlines toward three dwarf galaxies in the anti-M31 direction and at the periphery of the Local Group ($D \approx$ 1.3 Mpc; Sextans A, Sextans B, and NGC 3109). The C IV doublet is detected in only one of six sightlines, toward Sextans A, with $\log N({\rm C IV})$ = $13.06 \pm 0.08$. This is consistent with our gaseous halo models, where the halo gas mass is determined by the cooling rate, feedback, and the star formation rate; the inclusion of photoionization is an essential ingredient. This model can also reproduce the higher detection rate of O VI absorption in other dwarf samples (beyond the Local Group), and with C IV only detectable within $\sim 0.5R_{\rm vir}$.

preprint2022arXiv

CHANG-ES XXIX: The Sub-kpc Nuclear Bubble of NGC 4438

AGN bubbles could play an important role in accelerating high-energy CRs and galactic feedback. Only in nearby galaxies could we have high enough angular resolution in multi-wavelengths to study the sub-kpc environment of the AGN, where the bubbles are produced and strongly interact with the surrounding ISM. In this paper, we present the latest Chandra observations of the Virgo cluster galaxy NGC 4438, which hosts multi-scale bubbles detected in various bands. The galaxy also has low current star formation activity, so these bubbles are evidently produced by the AGN rather than a starburst. We present spatially resolved spectral analysis of the Chandra data of the $\sim3^{\prime\prime}\times5^{\prime\prime}$ ($\sim200{\rm~pc}\times350\rm~pc$) nuclear bubble of NGC 4438. The power law tail in the X-ray spectra can be most naturally explained as synchrotron emission from high-energy CR leptons. The hot gas temperature increases, while the overall contribution of the non-thermal X-ray emission decreases with the vertical distance from the galactic plane. We calculate the synchrotron cooling timescale of the CR leptons responsible for the non-thermal hard X-ray emission to be only a few tens to a few hundreds of years. The thermal pressure of the hot gas is about three times the magnetic pressure, but the current data cannot rule out the possibility that they are still in pressure balance. The spatially resolved spectroscopy presented in this paper may have important constraints on how the AGN accelerates CRs and drives outflows. We also discover a transient X-ray source only $\sim5^{\prime\prime}$ from the nucleus of NGC 4438. The source was not detected in 2002 and 2008, but became quite X-ray bright in March 2020, with an average 0.5-7 keV luminosity of $\sim10^{39}\rm~ergs~s^{-1}$.

preprint2022arXiv

Characterizing the Variable X-ray and UV-Optical Flux Behavior of Blazars

The variability of blazars in the X-ray and optical regions informs both the physics of their emitting region and places demands on the observer if a program requires that the object be bright or faint. The extensive simultaneous X-ray and optical observation by Swift provides the best insight into the variable nature of these objects. This program uses \textit{Swift} data for 19 X-ray-bright blazars, generally at $z > 0.1$, to determine their variability properties. The analysis is based on structure functions and provides insight into the nature of the variability and how it depends on time, luminosity, and redshift. We also consider strategies for observing blazars at or above average brightness, given a time delay between planning an observation and obtaining the data. This is critical to observations with orbiting X-ray telescopes, current or future. The variability in the soft X-ray band is typically three to eight times larger than at UV-optical wavelengths, at fixed time differences (i.e., 30 or 100 days). There is almost no difference in the amplitude of variation (X-ray or UV-optical) as a function of redshift (time delay of 30 days) and a modest positive correlation with luminosity. In the X-ray band, blazars that become brighter than normal typically remain bright for at least 2-3 months, although with significant flickering. One can avoid observing objects that are significantly below the average X-ray flux by scheduling the observation when the $F_X > 0.9 F_{X,avg}$, which requires monitoring observations near the time of the scheduling activity.

preprint2022arXiv

eDIG-CHANGES I: Extended Hα Emission from the Extraplanar Diffuse Ionized Gas (eDIG) around CHANG-ES Galaxies

The extraplanar diffuse ionized gas (eDIG) represents the cool/warm ionized gas reservoir around galaxies. We present a spatial analysis of H$α$ images of 22 nearby edge-on spiral galaxies from the CHANG-ES sample (the eDIG-CHANGES project), taken with the APO 3.5m telescope, in order to study their eDIG. We conduct an exponential fit to the vertical intensity profiles of the sample galaxies, of which 16 can be decomposed into a thin disk plus an extended thick disk component. The median value of the scale height (h) of the extended component is $1.13\pm 0.14$ kpc. We find a tight sublinear correlation between h and the SFR. Moreover, the offset of individual galaxies from the best-fit SFR-h relation shows significant anti-correlation with SFR_SD. This indicates that galaxies with more intense star formation tend to have disproportionately extended eDIG. Combined with data from the literature, we find that the correlations between the eDIG properties and the galaxies' properties extend to broader ranges. We further compare the vertical extension of the eDIG to multi-wavelength measurements of other CGM phases. We find the eDIG to be slightly more extended than the neutral gas (HI 21-cm line), indicating the existence of some extended ionizing sources. Most galaxies have an X-ray scale height smaller than the h, suggesting that the majority of the X-ray emission detected in shallow observations are actually from the thick disk. The h is comparable to the L-band radio continuum scale height, both slightly larger than that at higher frequencies (C-band), where the cooling is stronger and the thermal contribution may be larger. The comparable H$α$ and L-band scale height indicates that the thermal and non-thermal electrons have similar spatial distributions. This further indicates that the thermal gas, the cosmics rays, and the magnetic field may be close to energy equipartition.

preprint2022arXiv

Hot Extended Galaxy Halos Around Local L* Galaxies From Sunyaev-Zeldovich Measurements

Most of the baryons in L* galaxies are unaccounted for and are predicted to lie in hot gaseous halos (T ~ 3E6 K) that may extend beyond R200. A hot gaseous halo will produce a thermal Sunyaev-Zeldovich signal that is proportional to the product of the gas mass and the mass-weighted temperature. To best detect this signal, we used a Needlet Independent Linear Combination all-sky Planck map that we produced from the most recent Planck data release, also incorporating WMAP data. The sample is 12 L* spiral galaxies with distances of 3-10 Mpc, which are spatially resolved so that contamination from the optical galaxy can be excluded. One galaxy, NGC 891, has a particularly strong SZ signal, and when excluding it, the stack of 11 galaxies is detected at about 4sigma (declining with radius) and is extended to at least 250 kpc (~R_{200}) at > 99% confidence. The gas mass within a spherical volume to a radius of 250 kpc is 9.8 +/- 2.8 E10 Msun, for Tavg = 3E6 K. This is about 30% of the cosmic baryon content of the average galaxy (3.1E11 Msun), and about equal to the mass of stars, disk gas, and warm halo gas. The remaining missing baryons (~ 1.4E11 Msun, 40-50% of the total baryon content) are likely to be hot and extend to the 400-500 kpc volume, if not beyond. The result is higher than predictions, but within the uncertainties.

preprint2022arXiv

The Solar-Cycle Temporal Variation of the Solar Wind Charge Exchange X-ray Lines

Solar wind charge exchange (SWCX) is the primary contamination to soft X-ray emission lines from the Milky Way (MW) hot gas. We report a solar-cycle ($\approx 10$ yr) temporal variation of observed \ion{O}{7} and \ion{O}{8} emission line measurements in the {\it XMM-Newton} archive, which is tightly correlated with the solar cycle traced by the sunspot number (SSN). This temporal variation is expected to be associated with the heliospheric SWCX. Another observed correlation is that higher solar wind (SW) fluxes lead to higher O VII or O VIII fluxes, which is due to the magnetospheric SWCX. We construct an empirical model to reproduce the observed correlation between the line measurements and the solar activity (i.e., the SW flux and the SSN). With this model we discovered a lag of $0.91_{-0.22}^{+0.20}$ yr between the O VII flux and the SSN. This time lag is a combination of the SW transit time within the heliosphere, the lag of the neutral gas distribution responding to solar activity, and the intrinsic lag between the SSN and the launch of a high-energy SW (i.e., $\rm O^{7+}$ and $\rm O^{8+}$). MW O VII and O VIII fluxes have mean values of 5.4 L.U. and 1.7 L.U., which are reduced by $50\%$ and $30\%$, compared to studies where the SWCX contamination is not removed. This correction also changes the determination of the density distribution and the temperature profile of the MW hot gas.

preprint2021arXiv

Optical Confirmation of X-ray selected Galaxy clusters from the Swift AGN and Cluster survey with MDM and Pan-STARRS Data (Paper III)

To understand structure formation in the universe and impose stronger constraints on the cluster mass function and cosmological models, it is important to have large galaxy cluster catalogs. The Swift AGN and Cluster Survey is a serendipitous X-ray survey aimed at building a large statistically selected X-ray cluster catalog with 442 cluster candidates in its first release. Our initial SDSS follow-up study confirmed $50\%$ of clusters in the SDSS footprint as z $<$ 0.5 clusters. Here, we present further optical follow-up analysis of 248 (out of 442) cluster candidates from the Swift cluster catalog using multi-band imaging from the MDM $2.4m$ telescope and the Pan-STARRS survey. We report the optical confirmation of 55 clusters with $> 3σ$ galaxy overdensities and detectable red sequences in the color-magnitude space. The majority of these confirmed clusters have redshifts z $<$ 0.6. The remaining candidates are potentially higher redshift clusters that are excellent targets for infrared observations. We report the X-ray luminosity and the optical richness for these confirmed clusters. We also discuss the distinction between X-ray and optical observables for the detected and non-detected cluster candidates.

preprint2021arXiv

The Hosts of X-ray Absorption Lines Toward AGNs

Most baryonic matter in the universe exists in gaseous form and can be found in structures such as galactic halos and the low-density intergalactic medium. proposed-ray spectroscopy missions such as Athena, Arcus, and Lynx will have the capability to identify absorption lines in spectra toward bright active galactic nuclei (AGNs), which can be used as a tool to probe this missing matter. In this study, we examine the optical fields surrounding 15 primary observing targets and identify the foreground galaxies and galaxy groups that are potential hosts of absorption. We record the basic properties of the potential host and their angular and physical separation from the AGN line of sight. This process is done by marking the location of various galaxies and groups in optical images of the field surrounding the target and plotting their angular separation vs. redshift to gauge physical proximity to the background source. We identify the surrounding objects according to those which have measured redshifts and those that require them.

preprint2021arXiv

The Warm Gas in the MW: the Kinematical Model of C IV and its Connection with Si IV

We compose a 265-sight line MW C IV line shape sample using the HST/COS archive, which is complementary to the existing Si IV samples. C IV has a higher ionization potential ($47 - 64$ eV) than Si IV ($33 - 45$ eV), so it also traces warm gas, which is roughly cospatial with Si IV. The spatial density distribution and kinematics of C IV is identical with Si IV within $\approx 2 σ$. C IV is more sensitive to the warm gas density distribution at large radii with a higher element abundance. Applying the kinematical model to the C IV sample, we find two possible solutions of the density distribution, which are distinguished by the relative extension along the disk mid-plane and the normal-line direction. Both two solutions can reproduce the existing sample, and suggest a warm gas disk mass of $\log M(M_\odot) \approx 8$ and an upper limit of $\log M(M_\odot) < 9.3$ within 250 kpc, which is consistent with Si IV. There is a decrease of the C IV/Si IV column density ratio from the Galactic center to outskirts by $0.2-0.3$ dex, which may suggest a phase transition or different ionization mechanisms for C IV and Si IV. Also, we find that the difference between C IV and Si IV is an excellent tracer of small-scale features, and we find a typical size of $5^\circ-10^\circ$ for possible turbulence within individual clouds ($\approx 1\rm~kpc$).

preprint2020arXiv

An X-ray and SZ bright diffuse source toward M31: a Local Hot Bridge

We report a large-scale ($r\approx 20^\circ$) X-ray and Sunyaev-Zeldovich (SZ)-bright diffuse enhancement toward M31, which might be a Local Hot Bridge connecting the Milky Way (MW) with M31. We subtract the Galactic emission from the all-sky O VII and O VIII emission line measurement survey, and find that the emission of these two ions is enhanced within $r\approx20^\circ$ around M31. The mean emission enhancements are $5.6\pm 1.3$ L.U., and $2.8\pm0.6$ L.U. for O VII and O VIII, respectively ($>4σ$ for both ions). We also extract the SZ signal around M31, which suggests a surface brightness $y$ of $2-4\times10^{-7}$, an enhancement $>2.5σ$ (and a best fit of $5.9σ$). These three measurements trace the hot gas with a temperature $\log~T({\rm K})> 6$, showing similar plateau shapes (flat within $\approx15^\circ$, and zero beyond $\approx30^\circ$). A single-phase assumption leads to a temperature of $\log~T({\rm K})=6.34\pm0.03$, which is determined by the O VII/O VIII line ratio. Combining X-ray and SZ measurements, we suggest that this feature is unlikely to be the hot halo around M31 (too massive) or in the MW (too high pressure and X-ray bright). The plateau shape may be explained by a cylinder connecting the MW and M31 (the Local Hot Bridge). We constrain its length to be about 400 kpc, with a radius of 120 kpc, a density of $\approx 2\times10^{-4}-10^{-3} ~\rm cm^{-3}$, and a metallicity of $0.02-0.1~ Z_\odot$. The baryon mass is $\gtrsim10^{11}~M_\odot$, and the oxygen mass is about $\gtrsim10^8~M_\odot$, which contribute to the baryon or metal budget of the Local Group.

preprint2020arXiv

Chandra Detection of Three X-ray Bright Quasars at z>5

We report Chandra detection of three UV bright radio quiet quasars at $z\gtrsim5$. We have collected a sufficient number of photons to extract an X-ray spectrum of each quasar to measure their basic X-ray properties, such as the X-ray flux, power law photon index ($Γ$), and optical-to-X-ray spectral slope ($α_{\rm OX}$). J074749+115352 at $z=5.26$ is the X-ray brightest radio-quiet quasar at $z>5$. It may have a short timescale variation (on a timescale of $\sim3800\rm~s$ in the observer's frame, or $\sim600\rm~s$ in the rest frame) which is however largely embedded in the statistical noise. We extract phase folded spectra of this quasar. There are two distinguishable states: a "high soft" state with an average X-ray flux $\sim2.7$ times of the "low hard" state, and a significantly steeper X-ray spectral slope ($Γ=2.40_{-0.32}^{+0.33}$ vs $1.78_{-0.24}^{+0.25}$). We also compare the three quasars detected in this paper to other quasar samples. We find that J074749+115352, with a SMBH mass of $M_{\rm SMBH}\approx1.8\times10^9\rm~M_\odot$ and an Eddington ratio of $λ_{\rm Edd}\approx2.3$, is extraordinarily X-ray bright. It has an average $α_{\rm OX}=-1.46\pm0.02$ and a 2-10 keV bolometric correction factor of $L_{\rm bol}/L_{\rm2-10keV}=42.4\pm5.8$, both significantly depart from some well defined scaling relations. We compare $Γ$ of the three quasars to other samples at different redshifts, and do not find any significant redshift evolution based on the limited sample of $z>5$ quasars with reliable measurements of the X-ray spectral properties.

preprint2020arXiv

SZ Scaling Relations of Galaxy Groups and Clusters Near the North Ecliptic Pole

SZ scaling relations have been used to test the self-similar prediction for massive galaxy clusters, but little attention has been given to individual galaxy groups. We investigate the scaling relations of galaxy groups and clusters near the North Ecliptic Pole using X-ray and SZ observations. This region of the sky is where both the ROSAT and Planck satellites achieved their deepest observations, permitting the investigation of lower mass systems. Our sample consists of 62 X-ray detected groups and clusters, spanning a mass range of $10^{13.4}M_{\odot}<~M_{500}<10^{15}M_{\odot}$ and redshifts of $0.03\lesssim z \lesssim 0.82$. We extract the total SZ flux from unresolved Planck data and estimate the fraction of the SZ flux within $R_{500}$ assuming two different pressure profiles. The SZ scaling relations were derived using a Bayesian technique that accounts for censored data. We find a power law slope of $1.73^{+0.19}_{-0.18}$ for the $Y_{SZ}-M_{500}$ relation which is consistent with the self-similar prediction of 5/3. The slope of $0.89^{+0.09}_{-0.08}$ for the $Y_{SZ}-L_{X,500}$ relation is in agreement with other observational studies but not the self-similar prediction of 5/4, and the $Y_{SZ}-Y_{X}$ relation lies below the 1:1 relation when the slope is fixed to unity. The determined scaling relations are dependent on the selected pressure profile, so resolved data are needed to determine the effects of AGN feedback. In addition, we find a number of potential cluster candidates in the Planck Compton maps that were not identified in our X-ray sample.

preprint2020arXiv

The Warm Gas in the MW: A Kinematical Model

We develop a kinematical model for the Milky Way Si IV-bearing gas to determine its density distribution and kinematics. This model is constrained by a column density line shape sample extracted from the {\it HST}/COS archival data, which contains 186 AGN sight lines. We find that the Si IV ion density distribution is dominated by an extended disk along the $z$-direction (above or below the midplane), i.e., $n(z)=n_0\exp(-(z/z_0)^{0.82})$, where $z_0$ is the scale height of $6.3_{-1.5}^{+1.6}$ kpc (northern hemisphere) and $3.6_{-0.9}^{+1.0}$ kpc (southern hemisphere). The density distribution of the disk in the radial direction shows a sharp edge at $15-20$ kpc given by, $n(r_{\rm XY})=n_0\exp(-(r_{\rm XY}/r_0)^{3.36})$, where $r_0 \approx 12.5\pm0.6$ kpc. The difference of density distributions over $r_{\rm XY}$ and $z$ directions indicates that the warm gas traced by \ion{Si}{4} is mainly associated with disk processes (e.g., feedback or cycling gas) rather than accretion. We estimate the mass of the warm gas (within 50 kpc) is $\log (M(50 {\rm kpc})/M_\odot)\approx8.1$ (assuming $Z\approx0.5Z_\odot$), and a $3σ$ upper limit of $\log (M(250 {\rm kpc})/M_\odot)\approx9.1$ (excluding the Magellanic system). Kinematically, the warm gas disk is nearly co-rotating with the stellar disk at $v_{\rm rot}=215\pm3\rm~km~s^{-1}$, which lags the midplane rotation by about $8\rm~km~s^{-1}~kpc^{-1}$ (within 5 kpc). Meanwhile, we note that the warm gas in the northern hemisphere has significant accretion with $v_{\rm acc}$ of $69\pm 7\rm ~km~s^{-1}$ at 10 kpc (an accretion rate of $-0.60_{-0.13}^{+0.11}~M_\odot\rm~yr^{-1}$), while in the southern hemisphere, there is no measurable accretion, with an upper limit of $0.4~M_\odot\rm~yr^{-1}$.

preprint2016arXiv

A hot gaseous galaxy halo candidate with Mg X absorption

The hot gas in galaxy halos may account for a significant fraction of missing baryons in galaxies, and some of these gases can be traced by high ionization absorption systems in QSO UV spectra. Using high S/N ratio $ HST$/COS spectra, we discovered a high ionization state system at $z=1.1912$ in the sightline toward \object{LBQS 1435-0134}, and two-components absorption lines are matched for Mg X, Ne VIII, Ne VI, O VI, Ne V, O V, Ne IV, O IV, N IV, O III, and H I. Mg X, detected for the first time ($5.8 σ$), is a particularly direct tracer of hot galactic halos, as its peak ion fraction occurs near $10^{6.1}\rm~ K$, about the temperature of a virialized hot galaxy halo of mass $\sim 0.5 M^*$. With Mg X and NeVIII, a photoionization model cannot reproduce the observed column densities with path lengths of galaxy halos. For collisional ionization models, one or two temperature models do not produce acceptable fits, but a three temperature model or a power law model can produce the observed results. In the power law model, ${\rm d}N/{\rm d}T = 10^{4.4\pm 2.2-[Z/X]} T^{1.55\pm 0.41}$ with temperatures in the range $10^{4.39\pm0.13} {\rm~K} < T < 10^{6.04\pm0.05}~ {\rm K}$, the total hydrogen column density is $8.2 \times 10^{19} (0.3Z_{\odot}/Z) \rm~ cm^{-2}$ and the positive power law index indicates most of the mass is at the high temperature end. We suggest that this absorption system is a hot volume-filled galaxy halo rather than interaction layers between the hot halo and cool clouds. The temperature dependence of the column density is likely due to the local mixture of multiple phase gases.

preprint2016arXiv

Chandra ACIS Survey of X-ray Point Sources in Nearby Galaxies. II. X-ray Luminosity Functions and Ultraluminous X-ray Sources

Based on the recently completed {\it Chandra}/ACIS survey of X-ray point sources in nearby galaxies, we study the X-ray luminosity functions (XLFs) for X-ray point sources in different types of galaxies and the statistical properties of ultraluminous X-ray sources (ULXs). Uniform procedures are developed to compute the detection threshold, to estimate the foreground/background contamination, and to calculate the XLFs for individual galaxies and groups of galaxies, resulting in an XLF library for 343 galaxies of different types. With the large number of surveyed galaxies, we have studied the XLFs and ULX properties across different host galaxy types, and confirm with good statistics that the XLF slope flattens from lenticular ($α\sim1.50\pm0.07$) to elliptical ($\sim1.21\pm0.02$), to spirals ($\sim0.80\pm0.02$), to peculiars ($\sim0.55\pm0.30$), and to irregulars ($\sim0.26\pm0.10$). The XLF break dividing the neutron star and black hole binaries is also confirmed, albeit at quite different break luminosities for different types of galaxies. A radial dependency is found for ellipticals, with a flatter XLF slope for sources located between $D_{25}$ and 2$D_{25}$, suggesting the XLF slopes in the outer region of early-type galaxies are dominated by low-mass X-ray binaries in globular clusters. This study shows that the ULX rate in early-type galaxies is $0.24\pm0.05$ ULXs per surveyed galaxy, on a $5σ$ confidence level. The XLF for ULXs in late-type galaxies extends smoothly until it drops abruptly around $4\times10^{40}$ erg s$^{-1}$, and this break may suggest a mild boundary between the stellar black hole population possibly including 30 $M_\odot$ black holes with super-Eddington radiation and intermediate mass black holes.

preprint2016arXiv

Strategies for Detecting the Missing Hot Baryons in the Universe

About 30-50% of the baryons in the local Universe are unaccounted for and are likely in a hot phase, 10^5.5-10^8 K. A hot halo (10^6.3 K) is detected around the Milky Way through the O VII and O VIII resonance absorption and emission lines in the soft X-ray band. Current instruments are not sensitive enough to detect this gas in absorption around other galaxies and galaxy groups, the two most likely sites. We show that resonant line absorption by this hot gas can be detected with current technology, with a collecting area exceeding about 300 cm^2 and a resolution R > 2000. For a few notional X-ray telescope configurations that could be constructed as Explorer or Probe missions, we calculate the differential number of O VII and O VIII absorbers as a function of equivalent width through redshift space, dN/dz. The hot halos of individual external galaxies produce absorption that should be detectable out to about their virial radii. For the Milky Way, one can determine the radial distribution of density, temperature, and metallicity, after making optical depth corrections. Spectroscopic observations can determine the rotation of a hot gaseous halo.

preprint2016arXiv

The Circum-Galactic Medium of MASsive Spirals I: Overview and a Case Study of NGC 5908

The Circum-Galactic Medium of MASsive Spirals (CGM-MASS) is a project studying the overall content, physical and chemical properties, and spatial distributions of the multi-phase circum-galactic medium (CGM) around a small sample of the most massive isolated spiral galaxies in the local Universe. We introduce the sample and present a detailed case study of the XMM-Newton observation of the hot gas halo of NGC5908. After data calibration, point source removal, and background analysis, we find that the diffuse soft X-ray emission of NGC5908 is significantly more extended than the stellar light in the vertical direction. The 0.5-1.25keV radial intensity profile tracing hot gas emission can be detected above the background out to about $2^\prime$, or $30\rm~kpc$ from the nucleus. The unresolved soft X-ray emission can be characterized with a $β$-model with a slope of $β\approx0.68$. The unresolved 0.5-2keV luminosity within $1^\prime$ is $6.8\times10^{39}\rm~ergs~s^{-1}$, but young stellar sources may contribute significantly to it. Assuming a metallicity of 0.2solar, an upper limit (without subtracting the very uncertain young stellar contribution) to the mass of hot gas within this radius is $2.3\times10^9\rm~M_\odot$. The cooling radius is $r_{\rm cool}\approx25\rm~kpc$ or $\approx0.06r_{\rm 200}$, within which the hot gas could cool radiatively in less than 10Gyr, and the cooling of hot gas could significantly contribute in replenishing the gas consumed in star formation. The hot gas accounts for $\approx1.9\%$ of the baryon detected within the cooling radius. By comparing NGC5908 to other galaxies, we find that its X-ray luminosity per stellar mass is consistent with lower-mass non-starburst field spiral galaxies. However, a large scatter in hot gas soft X-ray emissivity is indicated for spiral galaxies with $M_*\gtrsim2\times10^{11}\rm~M_\odot$.

preprint2016arXiv

The Interaction of the Fermi Bubbles with the Milky Way's Hot Gas Halo

The Fermi bubbles are two lobes filled with non-thermal particles that emit gamma rays, extend $\approx$10 kpc vertically from the Galactic center, and formed from either nuclear star formation or accretion activity on Sgr A*. Simulations predict a range of shock strengths as the bubbles expand into the surrounding hot gas halo distribution ($T_{halo} \approx 2 \times 10^6$ K), but with significant uncertainties in the energetics, age, and thermal gas structure. The bubbles should contain thermal gas with temperatures between $10^6$ and $10^8$ K, with potential X-ray signatures. In this work, we constrain the bubbles' thermal gas structure by modeling the OVII and OVIII emission line strengths from archival XMM-Newton and Suzaku data. Our emission model includes a hot thermal volume-filled bubble component cospatial with the gamma-ray region, and a shell of compressed material. We find that a bubble/shell model with $n \approx 1 \times 10^{-3}$ cm$^{-3}$ and with log($T$) $\approx$ 6.60-6.70 is consistent with the observed line intensities. In the framework of a continuous Galactic outflow, we infer a bubble expansion rate, age, and energy injection rate of $490_{-77}^{+230}$ km s$^{-1}$, $4.3_{-1.4}^{+0.8}$ Myr, and $2.3_{-0.9}^{+5.1} \times 10^{42}$ erg s$^{-1}$. These estimates are consistent with the bubbles forming from a Sgr A* accretion event rather than from nuclear star formation.

preprint2016arXiv

The Milky Way's Hot Gas Kinematics: Signatures in Current and Future OVII Absorption Line Observations

Detections of $z \approx$ 0 oxygen absorption and emission lines indicate the Milky Way hosts a hot ($\sim 10^6$ K), low-density plasma extending $\gtrsim$50 kpc into the Mily Way's halo. Current X-ray telescopes cannot resolve the line profiles, but the variation of their strengths on the sky constrains the radial gas distribution. Interpreting the OVII K$α$ absorption line strengths has several complications, including optical depth and line of sight velocity effects. Here, we present model absorption line profiles accounting for both of these effects to show the lines can exhibit asymmetric structures and be broader than the intrinsic Doppler width. The line profiles encode the hot gas rotation curve, the net inflow or outflow of hot gas, and the hot gas angular momentum profile. We show how line of sight velocity effects impact the conversion between equivalent width and the column density, and provide modified curves of growth accounting for these effects. As an example, we analyze the LMC sight line pulsar dispersion measure and OVII equivalent width to show the average gas metallicity is $\gtrsim 0.6 Z_{\odot}$ and $b$ $\gtrsim$ 100 km s$^{-1}$. Determining these properties offers valuable insights into the dynamical state of the Milky Way's hot gas, and improves the line strength interpretation. We discuss future strategies to observe these effects with an instrument that has a spectral resolution of about 3000, a goal that is technically possible today.

preprint2016arXiv

The Rotation of the Hot Gas Around the Milky Way

The hot gaseous halos of galaxies likely contain a large amount of mass and are an integral part of galaxy formation and evolution. The Milky Way has a 2e6 K halo that is detected in emission and by absorption in the OVII resonance line against bright background AGNs, and for which the best current model is an extended spherical distribution. Using XMM-Newton RGS data, we measure the Doppler shifts of the OVII absorption-line centroids toward an ensemble of AGNs. These Doppler shifts constrain the dynamics of the hot halo, ruling out a stationary halo at about 3sigma and a corotating halo at 2sigma, and leading to a best-fit rotational velocity of 183+/-41 km/s for an extended halo model. These results suggest that the hot gas rotates and that it contains an amount of angular momentum comparable to that in the stellar disk. We examined the possibility of a model with a kinematically distinct disk and spherical halo. To be consistent with the emission-line X-ray data the disk must contribute less than 10% of the column density, implying that the Doppler shifts probe motion in the extended hot halo.

preprint2016arXiv

The Swift AGN and Cluster Survey. II. Cluster Confirmation with SDSS Data

We study 203 (of 442) Swift AGN and Cluster Survey extended X-ray sources located in the SDSS DR8 footprint to search for galaxy over-densities in three dimensional space using SDSS galaxy photometric redshifts and positions near the Swift cluster candidates. We find 104 Swift clusters with a >3sigma galaxy over-density. The remaining targets are potentially located at higher redshifts and require deeper optical follow-up observations for confirmation as galaxy clusters. We present a series of cluster properties including the redshift, BCG magnitude, BCG-to-X-ray center offset, optical richness, and X-ray luminosity. We also detect red sequences in ~85% of the 104 confirmed clusters. The X-ray luminosity and optical richness for the SDSS confirmed Swift clusters are correlated and follow previously established relations. The distribution of the separations between the X-ray centroids and the most likely BCG is also consistent with expectation. We compare the observed redshift distribution of the sample with a theoretical model, and find that our sample is complete for z <~ 0.3 and is still 80% complete up to z ~= 0.4, consistent with the SDSS survey depth. These analysis results suggest that our Swift cluster selection algorithm has yielded a statistically well-defined cluster sample for further studying cluster evolution and cosmology. We also match our SDSS confirmed Swift clusters to existing cluster catalogs, and find 42, 23 and 1 matches in optical, X-ray and SZ catalogs, respectively, so the majority of these clusters are new detections.

preprint2015arXiv

A Deep XMM-Newton Study of the Hot Gaseous Halo Around NGC 1961

We examine 11 XMM-Newton observations of the giant spiral galaxy NGC 1961, allowing us to study the hot gaseous halo of a spiral galaxy in unprecedented detail. We perform a spatial and a spectral analysis; with the former, the hot halo is detected to at least 80 kpc and with the latter its properties can be measured in detail up to 42 kpc. We find evidence for a negative gradient in the temperature profile as is common for elliptical galaxies. We measure a rough metallicity profile, which is consistent with being flat at $Z \sim 0.2 Z_{\odot}$. Converting to this metallicity, the deprojected density profile is consistent with previous parametric fits, with no evidence for a break within 42 kpc ($\sim$0.1R$_{\text{vir}}$). Extrapolating to the virial radius, we infer a hot halo mass comparable to the stellar mass of the galaxy, and a baryon fraction from the stars and hot gas of around 30%. The cooling time of the hot gas is orders of magnitude longer than the dynamical time, making the hot halo stable against cooling instabilities, and we argue that an extended stream of neutral Hydrogen seen to the NW of this galaxy is instead likely due to accretion from the intergalactic medium. The low metallicity of the hot halo suggests it too was likely accreted. We compare the hot halo of NGC 1961 to hot halos around isolated elliptical galaxies, and show that the total mass determines the hot halo properties better than the stellar mass.

preprint2015arXiv

Binary Frequencies in Globular Clusters

Binary stars are predicted to have an important role in the evolution of globular clusters, so we obtained binary fractions for 35 globular clusters that were imaged in the F606W and F814W with the ACS on the Hubble Space Telescope. When compared to the values of prior efforts (Sollima et al. 2007; Milone et al. 2012), we find significant discrepancies, despite each group correcting for contamination effects and having performed the appropriate reliability tests. The most reliable binary fractions are obtained when restricting the binary fraction to q > 0.5. Our analysis indicates that the range of the binary fractions is nearly an order of magnitude for the lowest dynamical ages, suggesting that there is a broad distribution in the binary fraction at globular cluster formation. Dynamical effects also appears to decrease the core binary fractions by a factor of two over a Hubble time, but this is a weak relationship. We confirm a correlation from previous work that the binary fraction within the core radius decreases with cluster age, indicating that younger clusters formed with higher binary fractions. The strong radial gradient in the binary fraction with cluster radius appears to be a consequence of dynamical interactions. It is often not present in dynamically young clusters but nearly always present in dynamically old clusters.

preprint2015arXiv

Constraining the Milky Way's Hot Gas Halo with OVII and OVII Emission Lines

The Milky Way hosts a hot ($\approx 2 \times 10^6$ K), diffuse, gaseous halo based on detections of z = 0 OVII and OVIII absorption lines in quasar spectra and emission lines in blank-sky spectra. Here we improve constraints on the structure of the hot gas halo by fitting a radial model to a much larger sample of OVII and OVIII emission line measurements from XMM-Newton EPIC-MOS spectra compared to previous studies ($\approx$ 650 sightlines). We assume a modified $β$-model for the halo density distribution and a constant-density Local Bubble from which we calculate emission to compare with the observations. We find an acceptable fit to the OVIII emission line observations with $χ^{2}_{red}$ (dof) = 1.08 (644) for best-fit parameters of $n_o r_c^{3β} = 1.35 \pm 0.24$ cm$^{-3}$ kpc$^{3β}$ and $β= 0.50 \pm 0.03$ for the hot gas halo and negligible Local Bubble contribution. The OVII observations yield an unacceptable $χ^{2}_{red}$ (dof) = 4.69 (645) for similar best-fit parameters, which is likely due to temperature or density variations in the Local Bubble. The OVIII fitting results imply hot gas masses of $M$(< 50 kpc) = $3.8_{-0.3}^{+0.3} \times 10^{9} M_{\odot}$ and $M$(< 250 kpc) = $4.3_{-0.8}^{+0.9} \times 10^{10} M_{\odot}$, accounting for $\lesssim$ 50% of the Milky Way's missing baryons. We also explore our results in the context of optical depth effects in the halo gas, the halo gas cooling properties, temperature and entropy gradients in the halo gas, and the gas metallicity distribution. The combination of absorption and emission line analyses implies a sub-solar gas metallicity that decreases with radius, but that also must be $\geq 0.3 Z_{\odot}$ to be consistent with the pulsar dispersion measure toward the LMC.

preprint2015arXiv

Non-Detection of X-Ray Emission From Sterile Neutrinos in Stacked Galaxy Spectra

We conduct a comprehensive search for X-ray emission lines from sterile neutrino dark matter, motivated by recent claims of unidentified emission lines in the stacked X-ray spectra of galaxy clusters and the centers of the Milky Way and M31. Since the claimed emission lines lie around 3.5 keV, we focus on galaxies and galaxy groups (masking the central regions), since these objects emit very little radiation above $\sim 2$ keV and offer a clean background against which to detect emission lines. We develop a formalism for maximizing the signal-to-noise of decaying dark matter emission lines by weighing each X-ray event according to the expected dark matter profile. In total, we examine 81 and 89 galaxies with Chandra and XMM-Newton respectively, totaling 15.0 and 14.6 Ms of integration time. We find no significant evidence of any emission lines, placing strong constraints on the mixing angle of sterile neutrinos with masses between 4.8-12.4 keV. In particular, if the 3.57 keV feature from Bulbul et al. (2014) were due to 7.1 keV sterile neutrino emission, we would have detected it at $4.4σ$ and $11.8σ$ in our two samples. The most conservative estimates of the systematic uncertainties reduce these constraints to $4.4σ$ and 7.8$σ$, or letting the line energy vary between 3.50 and 3.60 keV reduces these constraints to $2.7σ$ and $11.0σ$ respectively. Unlike previous constraints, our measurements do not depend on the model of the X-ray background or on the assumed logarithmic slope of the center of the dark matter profile.

preprint2015arXiv

The SWIFT AGN and Cluster Survey I: Number Counts of AGN and Galaxy Clusters

The Swift AGN and Cluster Survey (SACS) uses 125 deg^2 of Swift XRT serendipitous fields with variable depths surrounding gamma-ray bursts to provide a medium depth (4e-15 erg/s/cm^2) and area survey filling the gap between deep, narrow Chandra/XMM-Newton surveys and wide, shallow ROSAT surveys. Here we present a catalog of 22,563 point sources and 442 extended sources and examine the number counts of the AGN and galaxy cluster populations. SACS provides excellent constraints on the AGN number counts at the bright end with negligible uncertainties due to cosmic variance, and these constraints are consistent with previous measurements. We use Wise mid-infrared (MIR) colors to classify the sources. For AGN we can roughly separate the point sources into MIR-red and MIR-blue AGN, finding roughly equal numbers of each type in the soft X-ray band (0.5-2 keV), but fewer MIR-blue sources in the hard X-ray band (2-8 keV). The cluster number counts, with 5% uncertainties from cosmic variance, are also consistent with previous surveys but span a much larger continuous flux range. Deep optical or IR follow-up observations of this cluster sample will significantly increase the number of higher redshift (z > 0.5) X-ray-selected clusters.

preprint2014arXiv

Modeling X-ray Emission Around Galaxies

Extended X-ray emission can be studied either spatially (through its surface brightness profile) or spectrally (by analyzing the spectrum at various locations in the field). Both techniques have advantages and disadvantages, and when the emission becomes particularly faint and/or extended, the two methods can disagree. We argue that an ideal approach would be to model the events file directly, and therefore to use both the spectral and spatial information which are simultaneously available for each event. In this work we propose a first step in this direction, introducing a method for spatial analysis which can be extended to leverage spectral information simultaneously. We construct a model for the entire X-ray image in a given energy band, and generate a likelihood function to compare the model to the data. A critical goal of this modeling is disentangling vignetted and unvignetted backgrounds through their different spatial distributions. Employing either maximum likelihood or Markov Chain Monte Carlo, we can derive probability distribution functions for the source and background parameters together, or we can fit and subtract the background, leaving the description of the source non-parametric. We calibrate and demonstrate this method against a variety of simulated images, and then apply it to Chandra observations of the hot gaseous halo around the elliptical galaxy NGC 720. We are able to follow the X-ray emission below a tenth of the background, and to infer a hot gas mass within 35 kpc of 4-5x10^9 Msun, with some indication that the profile continues to at least 50 kpc and that it steepens as the radius increases. We derive much stronger constraints on the surface brightness profile than previous studies, which employed the spectral method, and we show that the density profiles inferred from these studies are in conflict with the observed surface brightness profile. (abridged)

preprint2013arXiv

Binary Frequencies in a Sample of Globular Clusters. I. Methodology and Initial Results

Binary stars are thought to be a controlling factor in globular cluster evolution, since they can heat the environmental stars by converting their binding energy to kinetic energy during dynamical interactions. Through such interaction, the binaries determine the time until core collapse. To test predictions of this model, we have determined binary fractions for 35 clusters. Here we present our methodology with a representative globular cluster NGC 4590. We use HST archival ACS data in the F606W and F814W bands and apply PSF-fitting photometry to obtain high quality color-magnitude diagrams. We formulate the star superposition effect as a Poisson probability distribution function, with parameters optimized through Monte-Carlo simulations. A model-independent binary fraction of (6.2 +- 0.3)% is obtained by counting stars that extend to the red side of the residual color distribution after accounting for the photometric errors and the star superposition effect. A model-dependent binary fraction is obtained by constructing models with a known binary fraction and an assumed binary mass-ratio distribution function. This leads to a binary fraction range of 6.8% to 10.8%, depending on the assumed shape to the binary mass ratio distribution, with the best fit occurring for a binary distribution that favors low mass ratios (and higher binary fractions). We also represent the method for radial analysis of the binary fraction in the representative case of NGC 6981, which shows a decreasing trend for the binary fraction towards the outside, consistent with theoretical predictions for the dynamical effect on the binary fraction.

preprint2013arXiv

Direct Detections of Young Stars in Nearby Elliptical Galaxies

Small amounts of star formation in elliptical galaxies are suggested by several results: surprisingly young ages from optical line indices, cooling X-ray gas, and mid-IR dust emission. Such star formation has previously been difficult to directly detect, but using UV Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) imaging, we have identified individual young stars and star clusters in four nearby ellipticals. This technique is orders of magnitude more sensitive than other methods, allowing detections of star formation to 10^(-5) Msun/yr. Ongoing star formation is detected in all galaxies, including three ellipticals that have previously exhibited potential signposts of star forming conditions (NGC 4636, NGC 4697, and NGC 4374), as well as the typical "red and dead" NGC 3379. The current star formation in our closest targets, where we are most complete, is between 1-8x10^(-5) Msun/yr. The star formation history was roughly constant from 0.5-1.5 Gyr (at 3-5x10^(-4) Msun/yr), but decreased by a factor of several in the past 0.3 Gyr. Most star clusters have a mass between 10^2 - 10^4 Msun. The specific star formation rates of ~10^(-16) yr^(-1) (at the present day) or ~10^(-14) yr^(-1) (when averaging over the past Gyr) imply that a fraction 10^(-8) of the stellar mass is younger than 100 Myr and 10^(-5) is younger than 1 Gyr, quantifying the level of frosting of recent star formation over the otherwise passive stellar population. There is no obvious correlation between either the presence or spatial distribution of postulated star formation indicators and the star formation we detect.

preprint2013arXiv

Outflow vs. Infall in Spiral Galaxies: Metal Absorption in the Halo of NGC 891

Gas accreting onto a galaxy will be of low metallicity while halo gas due to a galactic fountain will be of near-solar metallicity. We test these predictions by measuring the metal absorption line properties of halo gas 5 kpc above the plane of the edge-on galaxy NGC 891, using observations taken with HST/STIS toward a bright background quasar. Metal absorption lines of Fe II, Mg II, and Mg I in the halo of NGC 891 are clearly seen, and when combined with recent deep H I observations, we are able to place constraints on the metallicity of the halo gas for the first time. The H I line width defines the line broadening, from which we model opacity effects in these metal lines, assuming the absorbing gas is continuously distributed in the halo. The gas-phase metallicities are [Fe/H] = -1.18 +/- 0.07 and [Mg/H] = -0.23 +0.36/-0.27 (statistical errors) and this difference is probably due to differential depletion onto grains. When corrected for such depletion using Galactic gas as a guide, both elements have approximately solar or even supersolar abundances. This suggests that the gas is from the galaxy disk, probably expelled into the halo by a galactic fountain, rather than from accretion of intergalactic gas, which would have a low metallicity. The abundances would be raised by significant amounts if the absorbing gas lies in a few clouds with thermal widths smaller than the rotational velocity of the halo. If this is the case, both the abundances and [Mg/Fe] would be supersolar.

preprint2013arXiv

Puzzling accretion onto a black hole in the ultraluminous X-ray source M101 ULX-1

There are two proposed explanations for ultraluminous X-ray sources (ULXs) with luminosities in excess of $10^{39}$ erg s$^{-1}$. They could be intermediate-mass black holes (more than 100-1,000, solar masses, $M_\odot$) radiating at sub-maximal (sub-Eddington) rates, as in Galactic black-hole X-ray binaries but with larger, cooler accretion disks. Alternatively, they could be stellar-mass black holes radiating at Eddington or super-Eddington rates. On its discovery, M101 ULX-1 had a luminosity of $3\times10^{39}$ erg s$^{-1}$ and a supersoft thermal disk spectrum with an exceptionally low temperature -- uncomplicated by photons energized by a corona of hot electrons -- more consistent with the expected appearance of an accreting intermediate-mass black hole. Here we report optical spectroscopic monitoring of M101 ULX-1. We confirm the previous suggestion that the system contains a Wolf-Rayet star, and reveal that the orbital period is 8.2 days. The black hole has a minimum mass of $5M_\odot$, and more probably a mass of $20-30M_\odot$, but we argue that it is very unlikely to be an intermediate-mass black hole. Therefore its exceptionally soft spectra at high Eddinton ratios violate the expectations for accretion onto stellar-mass black holes. Accretion must occur from captured stellar wind, which has hitherto been thought to be so inefficient that it could not power an ultraluminous source.

preprint2013arXiv

The Megasecond Chandra XVP Observation of NGC 3115: Witnessing the Flow of Hot Gas within the Bondi Radius

Observational confirmation of hot accretion model predictions has been hindered by the challenge to resolve spatially the Bondi radii of black holes with X-ray telescopes. Here, we use the Megasecond Chandra X-ray Visionary Project (XVP) observation of the NGC~3115 supermassive black hole to place the first direct observational constraints on the spatially and spectroscopically resolved structures of the X-ray emitting gas inside the Bondi radius of a black hole. We measured temperature and density profiles of the hot gas from a fraction out to tens of the Bondi radius (R_B = 2.4-4.8 arcsec = 112-224 pc). The projected temperature jumps significantly from ~0.3 keV beyond 5 arcsec to ~0.7 keV within ~4-5 arcsec, but then abruptly drops back to ~0.3 keV within ~3 arcsec. This is contrary to the expectation that the temperature should rise toward the center for a radiatively inefficient accretion flow. A hotter thermal component of ~1 keV inside 3 arcsec (~150 pc) is revealed using a two component thermal model, with the cooler ~0.3 keV thermal component dominating the spectra. We argue that the softer emission comes from diffuse gas physically located within $\sim 150$~pc from the black hole. The density profile is broadly consistent with rho ~ r^{-1} within the Bondi radius for either the single temperature or the two-temperature model. The X-ray data alone with physical reasoning argue against the absence of a black hole, supporting that we are witnessing the onset of the gravitational influence of the supermassive black hole.

preprint2012arXiv

A Deep X-ray View of the Hot Halo in the Edge-on Spiral Galaxy NGC 891

NGC 891 is a nearby edge-on galaxy that is similar to the Milky Way and has a hot X-ray emitting halo that could arise from accretion, a galactic fountain, or a combination of the two. The metallicity of the gas can help distinguish between these models, and here we report on results that use 138 ks of archival Chandra data and 92 ks of new XMM-Newton data to measure the temperature and metallicity of the hot halo of the galaxy. We find good fits for a thermal model with kT ~ 0.2 keV and Z ~ 0.1 solar, and rule out solar metallicity to more than 99% confidence. This result suggests accretion from the intergalactic medium as the origin for the hot halo. However, it is also possible to fit a two-temperature thermal model with solar metallicity where kT_1 = 0.1 keV and kT_2 = 0.25 keV. A consideration of the cooling rate and scale height prefers the single-temperature model. We also find that the cooling rate in the hot gas cannot explain the massive HI halo in the steady state. In addition, a galactic fountain model cannot eject enough mass to account for the HI halo, and we speculate that the neutral halo may be gas from a prior outflow that has since cooled.

preprint2012arXiv

A New Ultraluminous X-ray Source in the Nearby Edge-on Spiral NGC 891

We report the discovery of a new candidate ultraluminous X-ray source (ULX) in the nearby edge-on spiral galaxy NGC 891. The source, which has an absorbed flux of F_X ~ 10^-12 erg/s/cm^2 (corresponding to L_X > 10^40 erg/s at 9 Mpc), must have begun its outburst in the past 5 years as it is not detected in prior X-ray observations between 1986 and 2006. We try empirical fits to the XMM-Newton spectrum, finding that the spectrum is fit very well as emission from a hot disk, a cool irradiated disk, or blurred reflection from the innermost region of the disk. The simplest physically motivated model with an excellent fit is a hot disk around a stellar-mass black hole (a super-Eddington outburst), but equally good fits are found for each model. We suggest several follow-up experiments that could falsify these models.

preprint2012arXiv

Extended Hot Halos Around Isolated Galaxies Observed in the ROSAT All-Sky Survey

We place general constraints on the luminosity and mass of hot X-ray emitting gas residing in extended "hot halos" around nearby massive galaxies. We examine stacked images of 2165 galaxies from the 2MASS Very Isolated Galaxy Catalog (2MVIG), as well as subsets of this sample based on galaxy morphology and K-band luminosity. We detect X-ray emission at high confidence (ranging up to nearly 10σ) for each subsample of galaxies. The average L_X within 50 kpc is 1.0\pm0.1 (statistical) \pm0.2 (systematic) x10^40 erg/s, although the early-type galaxies are more than twice as luminous as the late-type galaxies. Using a spatial analysis, we also find evidence for extended emission around five out of seven subsamples (the full sample, the luminous galaxies, early-type galaxies, luminous late-type galaxies, and luminous early-type galaxies) at 92.7%, 99.3%, 89.3%, 98.7%, and 92.1% confidence, respectively. Several additional lines of evidence also support this conclusion and suggest that about 1/2 of the total emission is extended, and about 1/3 of the extended emission comes from hot gas. For the sample of luminous galaxies, which has the strongest evidence for extended emission, the average hot gas mass is 4x10^9 Msun within 50 kpc and the implied accretion rate is 0.4 Msun/yr.

preprint2012arXiv

The Production Rate of SN Ia Events in Globular Clusters

In globular clusters, dynamical evolution produces luminous X-ray emitting binaries at a rate about 200 times greater than in the field. If globular clusters also produce SNe Ia at a high rate, it would account for much of the SN Ia events in early type galaxies and provide insight into their formation. Here we use archival HST images of nearby galaxies that have hosted SNe Ia to examine the rate at which globular clusters produce these events. The location of the SN Ia is registered on an HST image obtained before the event or after the supernova faded. Of the 36 nearby galaxies examined, 21 had sufficiently good data to search for globular cluster hosts. None of the 21 supernovae have a definite globular cluster counterpart, although there are some ambiguous cases. This places an upper limit to the enhancement rate of SN Ia production in globular clusters of about 42 at the 95% confidence level, which is an order of magnitude lower than the enhancement rate for luminous X-ray binaries. Even if all of the ambiguous cases are considered as having a globular cluster counterpart, the upper bound for the enhancement rate is 82 at the 95% confidence level, excluding an enhancement rate of 200. Barring unforeseen selection effects, we conclude that globular clusters are not responsible for producing a significant fraction of the SN Ia events in early-type galaxies.

preprint2012arXiv

Ultraviolet Spectra of ULX Systems

To further understand the nature of the optical counterparts associated with Ultraluminous X-Ray sources (ULXs), we obtained far ultraviolet spectra of the reported counterparts to the ULXs in NGC 1313, Holmberg II, NGC 5204, and M81. The spectral resolution of the ACS prism spectra degrades from 300 at 1300 A to 40 at 1850 A, so longer wavelength features have lower S/N. The spectra of the ULXs in NGC 1313, Ho II, and NGC 5204 are quite similar, showing the N V 1240 A line at about the same equivalent width strength. The presence of this and other emission lines confirms the presence of an accretion disk, probably diluted by the light of an early B star companion. The spectra differ strongly from high mass X-ray binaries dominated by O star winds (e.g., Cyg X-1) and are most similar to intermediate mass X-ray binary systems in the Milky Way and LMC. This indicates that the mass transfer is due to Roche Lobe overflow. The spectrum of the ULX in M81 is quite weak but suggestive of a late-type Wolf-Rayet star.

preprint2011arXiv

Detection of a Hot Gaseous Halo Around the Giant Spiral Galaxy NGC 1961

Hot gaseous halos are predicted around all large galaxies and are critically important for our understanding of galaxy formation, but they have never been detected at distances beyond a few kpc around a spiral galaxy. We used the Chandra ACIS-I instrument to search for diffuse X-ray emission around an ideal candidate galaxy: the isolated giant spiral NGC 1961. We observed four quadrants around the galaxy for 30 ks each, carefully subtracting background and point source emission, and found diffuse emission that appears to extend to 40-50 kpc. We fit $β$-models to the emission, and estimate a hot halo mass within 50 kpc of $5\times10^9 M_{\odot}$. When this profile is extrapolated to 500 kpc (the approximate virial radius), the implied hot halo mass is $1-3\times10^{11} M_{\odot}$. These mass estimates assume a gas metallicity of $Z = 0.5 Z_{\odot}$. This galaxy's hot halo is a large reservoir of gas, but falls significantly below observational upper limits set by pervious searches, and suggests that NGC 1961 is missing 75% of its baryons relative to the cosmic mean, which would tentatively place it below an extrapolation of the baryon Tully-Fisher relationship of less massive galaxies. The cooling rate of the gas is no more than 0.4 $M_{\odot}$/year, more than an order of magnitude below the gas consumption rate through star formation. We discuss the implications of this halo for galaxy formation models.

preprint2011arXiv

Dynamical mass constraints on the ultraluminous X-ray Source NGC 1313 X-2

Dynamical mass measurements hold the key to answering whether ultraluminous X-ray sources (ULXs) are intermediate mass black holes (IMBHs) or stellar mass black holes with special radiation mechanisms. NGC 1313 X-2 is so far the only ULX with HST light curves, the orbital period, and the black hole's radial velocity amplitude based on the He II $\lambda4686$Å disk emission line shift of $\sim200$ km/s. We constrain its black hole mass and other parameters by fitting observations to a binary light curve code with accommodations for X-ray heating of the accretion disk and the secondary. Given the dynamical constraints from the observed light curves and the black hole radial motion and the observed stellar environment age, the only acceptable models are those with 40-50 Myrs old intermediate mass secondaries in their helium core and hydrogen shell burning phase filling 40%-80% of their Roche lobes. The black hole can be a massive black hole of a few tens of $M_\odot$ that can be produced from stellar evolution of low metalicity stars, or an IMBH of a few hundred to above 1000$M_\odot$ if its true radial velocity $2K^\prime<40$ km/s. Further observations are required to better measure the black hole radial motion and the light curves in order to determine whether NGC 1313 X-2 is a stellar black hole or an IMBH.

preprint2011arXiv

Resolving the Bondi Accretion Flow toward the Supermassive Black Hole of NGC 3115 with Chandra

Gas undergoing Bondi accretion onto a supermassive black hole (SMBH) becomes hotter toward smaller radii. We searched for this signature with a Chandra observation of the hot gas in NGC 3115, which optical observations show has a very massive SMBH. Our analysis suggests that we are resolving, for the first time, the accretion flow within the Bondi radius of an SMBH. We show that the temperature is rising toward the galaxy center as expected in all accretion models in which the black hole is gravitationally capturing the ambient gas. There is no hard central point source that could cause such an apparent rise in temperature. The data support that the Bondi radius is at about 4 arcsec-5 arcsec (188-235 pc), suggesting an SMBH of 2 x 10^9 M_sun that is consistent with the upper end of the optical results. The density profile within the Bondi radius has a power-law index of 1.03^{+0.23}_{-0.21} which is consistent with gas in transition from the ambient medium and the accretion flow. The accretion rate at the Bondi radius is determined to be {\dot M}_B = 2.2 x 10^{-2} M_sun yr^{-1}. Thus, the accretion luminosity with 10% radiative efficiency at the Bondi radius (10^{44} erg s^{-1}) is about six orders of magnitude higher than the upper limit of the X-ray luminosity of the nucleus.

preprint2011arXiv

XMM-Newton Detects a Hot Gaseous Halo in the Fastest Rotating Spiral Galaxy UGC 12591

We present our XMM-Newton observation of the fastest rotating spiral galaxy UGC 12591. We detect hot gas halo emission out to 110 kpc from the galaxy center, and constrain the halo gas mass to be smaller than 3.5e11 solar masses. We also measure the temperature of the hot gas as T=0.64\pm0.03 keV. Combining our X-ray constraints and the near-infrared and radio measurements in the literature, we find a baryon mass fraction of 0.03-0.04 in UGC 12591, suggesting a missing baryon mass of 75% compared with the cosmological mean value. Combined with another recent measurement in NGC 1961, the result strongly argues that the majority of missing baryons in spiral galaxies does not reside in their hot halos. We also find that UGC 12591 lies significantly below the baryonic Tully-Fisher relationship. Finally, we find that the baryon fractions of massive spiral galaxies are similar to those of galaxy groups with similar masses, indicating that the baryon loss is ultimately controlled by the gravitational potential well. The cooling radius of this gas halo is small, similar to NGC 1961, which argues that the majority of stellar mass of this galaxy is not assembled as a result of cooling of this gas halo.

preprint2010arXiv

Do Hot Haloes Around Galaxies Contain the Missing Baryons?

Galaxies are missing most of their baryons, and many models predict these baryons lie in a hot halo around galaxies. We establish observationally motivated constraints on the mass and radii of these haloes using a variety of independent arguments. First, the observed dispersion measure of pulsars in the Large Magellanic Cloud allows us to constrain the hot halo around the Milky Way: if it obeys the standard NFW profile, it must contain less than 4-5% of the missing baryons from the Galaxy. This is similar to other upper limits on the Galactic hot halo, such as the soft X-ray background and the pressure around high velocity clouds. Second, we note that the X-ray surface brightness of hot haloes with NFW profiles around large isolated galaxies is high enough that such emission should be observed, unless their haloes contain less than 10-25% of their missing baryons. Third, we place constraints on the column density of hot haloes using nondetections of OVII absorption along AGN sightlines: in general they must contain less than 70% of the missing baryons or extend to no more than 40 kpc. Flattening the density profile of galactic hot haloes weakens the surface brightness constraint so that a typical L$_*$ galaxy may hold half its missing baryons in its halo, but the OVII constraint remains unchanged, and around the Milky Way a flattened profile may only hold $6-13%$ of the missing baryons from the Galaxy ($2-4 \times 10^{10} M_{\odot}$). We also show that AGN and supernovae at low to moderate redshift - the theoretical sources of winds responsible for driving out the missing baryons - do not produce the expected correlations with the baryonic Tully-Fisher relationship and so are insufficient to explain the missing baryons from galaxies. We conclude that most of missing baryons from galaxies do not lie in hot haloes around the galaxies, and that the missing baryons never fell into the potential wells of protogalaxies in the first place. They may have been expelled from the galaxies as part of the process of galaxy formation.

preprint2010arXiv

Evidence for a Stellar Disruption by an IMBH in an Extragalactic Globular Cluster

We report [O III] 5007 Angstrom and [N II] 6583 Angstrom emission from a globular cluster harboring the ultraluminous X-ray source (ULX) CXOJ033831.8-352604 in the Fornax elliptical galaxy NGC1399. No accompanying Balmer emission lines are present in the spectrum. One possibility is that the forbidden lines emanate from X-ray illuminated debris of a star that has been tidally-disrupted by an intermediate-mass black hole (IMBH), with this debris also feeding the black hole leading to the observed X-ray emission. The line strengths indicate that the minimum size of the emitting region is ~10^15 cm, and if the 70 km/s half-widths of the emission lines represent rotation around the black hole, a minimum black hole mass of 1000 solar mass is implied. The non-detection of H-alpha and H-beta emission lines suggests a white dwarf star was disrupted, although the presence of strong nitrogen emission is somewhat of a mystery.

preprint2010arXiv

Metal Production in Galaxy Clusters: The Non-Galactic Component

The metallicity in galaxy clusters is expected to originate from the stars in galaxies, with a population dominated by high mass stars likely being the most important stellar component, especially in rich clusters. We examine the relationship between the metallicity and the prominence of galaxies as measured by the star to baryon ratio, M$_*$/M$_{bary}$. Counter to expectations, we rule out a metallicity that is proportional to M$_*$/M$_{bary}$, where the best fit has the gas phase metallicity decreasing with M$_*$/M$_{bary}$, or the metallicity of the gas plus the stars being independent of M$_*$/M$_{bary}$. This implies that the population of stars responsible for the metals is largely proportional to the total baryonic mass of the cluster, not to the galaxy mass within the cluster. If generally applicable, most of the heavy elements in the universe were not produced within galaxies.

preprint2010arXiv

On the Baryon Fractions in Clusters and Groups of Galaxies

We present the baryon fractions of 2MASS groups and clusters as a function of cluster richness using total and gas masses measured from stacked ROSAT X-ray data and stellar masses estimated from the infrared galaxy catalogs. We detect X-ray emission even in the outskirts of clusters, beyond r_200 for richness classes with X-ray temperatures above 1 keV. This enables us to more accurately determine the total gas mass in these groups and clusters. We find that the optically selected groups and clusters have flatter temperature profiles and higher stellar-to-gas mass ratios than the individually studied, X-ray bright clusters. We also find that the stellar mass in poor groups with temperatures below 1 keV is comparable to the gas mass in these systems. Combining these results with individual measurements for clusters, groups, and galaxies from the literature, we find a break in the baryon fraction at ~1 keV. Above this temperature, the baryon fraction scales with temperature as f_b \propto T^0.20\pm0.03. We see significantly smaller baryon fractions below this temperature, and the baryon fraction of poor groups joins smoothly onto that of systems with still shallower potential wells such as normal and dwarf galaxies where the baryon fraction scales with the inferred velocity dispersion as f_b \propto σ^1.6. The small scatter in the baryon fraction at any given potential well depth favors a universal baryon loss mechanism and a preheating model for the baryon loss. The scatter is, however, larger for less massive systems. Finally, we note that although the broken power-law relation can be inferred from data points in the literature alone, the consistency between the baryon fractions for poor groups and massive galaxies inspires us to fit the two categories of objects (galaxies and clusters) with one relation.

preprint2009arXiv

Elemental Abundances in the X-Ray Gas of Early-Type Galaxies with XMM and Chandra Observations

The source of hot gas in elliptical galaxies is thought to be due to stellar mass loss, with contributions from supernova events and possibly from infall from a surrounding environment. This picture predicts supersolar values for the metallicity of the gas toward the inner part of the galaxy, which can be tested by measuring the gas phase abundances. We use high-quality data for 10 nearby early-type galaxy from XMM-Newton, featuring both the EPIC and the Reflection Grating Spectrometer, where the strongest emission lines are detected with little blending; some Chandra data are also used. We find excellent consistency in the elemental abundances between the different XMM instruments and good consistency with Chandra. Differences in abundances with aperture size and model complexity are examined, but large differences rarely occur. For a two-temperature thermal model plus a point source contribution, the median Fe and O abundances are 0.86 and 0.44 of the Solar value, while Si and Mg abundances are similar to that for Fe. This is similar to stellar abundances for these galaxies but supernovae were expected to enhance the gas phase abundances considerably, which is not observed.

preprint2009arXiv

Oxygen Metallicity Determinations from Optical Emission Lines in Early-type Galaxies

We measured the oxygen abundances of the warm (T$\sim 10^{4}K$) phase of gas in seven early-type galaxies through long-slit observations. A template spectra was constructed from galaxies void of warm gas and subtracted from the emission-line galaxies, allowing for a clean measurement of the nebular lines. The ratios of the emission lines are consistent with photoionization, which likely originates from the UV flux of post-asymototic giant branch (PAGB) stars. We employ H II region photoionization models to determine a mean oxygen metallicity of $1.01\pm0.50$ solar for the warm interstellar medium (ISM) in this sample. This warm ISM 0.5 to 1.5 solar metallicity is consistent with modern determinations of the metallicity in the hot (T$\sim 10^{6}-10^{7}K$) ISM and the upper range of this warm ISM metallicity is consistent with stellar population metallicity determinations. A solar metallicity of the warm ISM favors an internal origin for the warm ISM such as AGB mass loss within the galaxy.

preprint2009arXiv

Redshift Evolution in the Iron Abundance of the Intracluster Medium

Clusters of galaxies provide a closed box within which one can determine the chemical evolution of the gaseous baryons with cosmic time. We studied this metallicity evolution in the hot X-ray emitting baryons through an analysis of XMM-Newton observations of 29 galaxy clusters in the redshift range 0.3 < z < 1.3. Taken alone, this data set does not show evidence for significant evolution. However, when we also include a comparable sample of 115 clusters observed with Chandra (Maughan et al. 2008) and a lower redshift sample of 70 clusters observed with XMM at z < 0.3 (Snowden et al. 2008), there is definitive evidence for a decrease in the metallicity. This decrease is approximately a factor of two from z = 0 to z \approx 1, over which we find a least-squares best-fit line Z(z) / Z_{\odot} = (0.46 \pm 0.05) - (0.38 \pm 0.03)z. The greatest uncertainty in the evolution comes from poorly constrained metallicities in the highest redshift bin.

preprint2008arXiv

High-Velocity Clouds in the Nearby Spiral Galaxy M 83

We present deep HI 21-cm and optical observations of the face-on spiral galaxy M 83 obtained as part of a project to search for high-velocity clouds (HVCs) in nearby galaxies. Anomalous-velocity neutral gas is detected toward M 83, with 5.6x10^7 Msolar of HI contained in a disk rotating 40-50 km/s more slowly in projection than the bulk of the gas. We interpret this as a vertically extended thick disk of neutral material, containing 5.5% of the total HI within the central 8 kpc. Using an automated source detection algorithm to search for small-scale HI emission features, we find eight distinct, anomalous-velocity HI clouds with masses ranging from 7x10^5 to 1.5x10^7 Msolar and velocities differing by up to 200 km/s compared to the HI disk. Large on-disk structures are coincident with the optical spiral arms, while unresolved off-disk clouds contain no diffuse optical emission down to a limit of 27 r' mag per square arcsec. The diversity of the thick HI disk and larger clouds suggests the influence of multiple formation mechanisms, with a galactic fountain responsible for the slowly-rotating disk and on-disk discrete clouds, and tidal effects responsible for off-disk cloud production. The mass and kinetic energy of the HI clouds are consistent with the mass exchange rate predicted by the galactic fountain model. If the HVC population in M 83 is similar to that in our own Galaxy, then the Galactic HVCs must be distributed within a radius of less than 25 kpc.

preprint2008arXiv

Mass Loss From Evolved Stars in Elliptical Galaxies

Most of the X-ray emitting gas in early-type galaxies probably originates from red giant mass loss and here we model the interaction between this stellar mass loss and the hot ambient medium. Using two-dimensional hydrodynamic simulations, we adopt a temperature for the ambient medium of 3E6 K along with a range of ambient densities and stellar velocities. When the stellar velocity is supersonic relative to the ambient medium, a bow shock occurs, along with a shock driven into the stellar ejecta, which heats only a fraction of the gas. Behind the bow shock, a cool wake develops but the fast flow of the hot medium causes Kelvin-Helmholtz instabilities to grow and these fingers are shocked and heated (without radiative cooling). Along with the mixing of this wake material with the hot medium, most of the stellar ejecta is heated to approximately the temperature of the hot ambient medium within 2 pc of the star. With the addition of radiative cooling, some wake material remains cool (< 1E5 K), accounting for up to 25% of the stellar mass loss. Less cooled gas survives when the ambient density is lower or when the stellar velocity is higher than in our reference case. These results suggest that some cooled gas should be present in the inner part of early-type galaxies that have a hot ambient medium. These calculations may explain the observed distributed optical emission line gas as well as the presence of dust in early-type galaxies.

preprint2000arXiv

Velocity Gradients in the Intracluster Gas of the Perseus Cluster

We report the results of spatially resolved X-ray spectroscopy of 8 different ASCA pointings distributed symmetrically around the center of the Perseus cluster. The outer region of the intracluster gas is roughly isothermal, with temperature ~ 6-7 keV, and metal abundance ~ 0.3 Solar. Spectral analysis of the central pointing is consistent with the presence of a cooling flow and a central metal abundance gradient. A significant velocity gradient is found along an axis at a position angle of \~135 deg, which is ~ 45 deg discrepant with the major axis of the X-ray elongation. The radial velocity difference is found to be greater than 1000 km/s/Mpc at the 90% confidence level. Simultaneous fittings of GIS 2 & 3 indicate that the velocity gradient is significant at the 95% confidence level and the F-test rules out constant velocities at the 99% level. Intrinsic short and long term variations of gain are unlikely (P < 0.03) to explain the velocity discrepancies.

preprint1998arXiv

Measuring Molecular, Neutral Atomic, and Warm Ionized Galactic Gas Through X-Ray Absorption

We study the column densities of neutral atomic, molecular, and warm ionized Galactic gas through their continuous absorption of extragalactic X-ray spectra at |b| > 25 degrees. For N(H,21cm) < 5x10^20 cm^-2 there is an extremely tight relationship between N(H,21cm) and the X-ray absorption column, N(xray), with a mean ratio along 26 lines of sight of N(xray)/N(H,21cm) = 0.972 +- 0.022. This is significantly less than the anticpated ratio of 1.23, which would occur if He were half He I and half He II in the warm ionized component. We suggest that the ionized component out of the plane is highly ionized, with He being mainly He II and He III. In the limiting case that H is entirely HI, we place an upper limit on the He abundance in the ISM of He/H <= 0.103. At column densities N(xray) > 5x10^20 cm^-2, which occurs at our lower latitudes, the X-ray absorption column N(xray) is nearly double N(H,21cm). This excess column cannot be due to the warm ionized component, even if He were entirely He I, so it must be due to a molecular component. This result implies that for lines of sight out of the plane with |b| ~ 30 degrees, molecular gas is common and with a column density comprable to N(H,21cm). This work bears upon the far infrared background, since a warm ionized component, anticorrelated with N(H,21cm), might produce such a background. Not only is such an anticorrelation absent, but if the dust is destroyed in the warm ionized gas, the far infrared background may be slightly larger than that deduced by Puget et al. (1996).

preprint1995arXiv

ROSAT Observations of Compact Groups of Galaxies

We have systematically analyzed a sample of 13 new and archival ROSAT PSPC observations of compact groups of galaxies: 12 Hickson Compact Groups plus the NCG 2300 group. We find that approximately two-thirds of the groups have extended X-ray emission and, in four of these, the emission is resolved into diffuse emission from gas at a temperature of $k$T $\sim 1$ keV in the group potential. All but one of the groups with extended emission have a spiral fraction of less than 50\%. The baryon fraction of groups with diffuse emission is 5--19\%, similar to the values in clusters of galaxies. However, with a single exception (HCG 62), the gas-to-stellar-mass ratio in our groups has a median value near 5\%, somewhat greater than the values for individual early-type galaxies and two orders of magnitude lower than in clusters of galaxies. The X-ray luminosities of individual group galaxies are comparable to those of similar field galaxies, although the L$_X$-L$_B$ relation for early-type galaxies may be flatter in compact groups than in the field.