We present the results of a joint XMM-Newton and NuSTAR observation (200 ks) of the galaxy cluster Abell 523 at $z=0.104$. The peculiar morphology of the cluster radio halo and its outlier position in the radio power P(1.4 GHz) - X-ray luminosity plane make it an ideal candidate for the study of radio-X-ray correlations and for the search of inverse Compton (IC) emission. We constructed thermodynamic maps derived from the XMM observations to describe in detail the physical and dynamical state of the ICM. We performed a point-to-point comparison in terms of surface brightness between the X-ray and radio emissions, to quantify their morphological discrepancies. Making use of NuSTAR's hard X-ray focusing capability, we looked for IC emission both globally and locally, after modeling the purely thermal component with a multi-temperature description. The thermodynamic maps suggest the presence of a secondary merging process that could be responsible for the peculiar radio halo morphology. This hypothesis is supported by the comparison between the X-ray and radio surface brightnesses, which shows a broad intrinsic scatter and a series of outliers. The global NuSTAR spectrum can be explained by purely thermal gas emission, and there is no convincing evidence that an IC component is needed. The $3\sigma$ upper limit on the IC flux in the 20-80 keV band is in the range $\left[2.2 - 4.0\right] \times 10^{-13} \, \mathrm{erg} \, \mathrm{s}^{-1} \, \mathrm{cm}^{-2}$, implying a lower limit on the magnetic field strength in the range $B > [0.23 - 0.31] \, \mu G$. Locally, we looked for IC emission in the central region of the cluster radio halo finding a $3\sigma$ upper limit on the 20-80 keV non-thermal flux of $3.17 \times 10^{-14} \, \mathrm{erg} \, \mathrm{s}^{-1} \, \mathrm{cm}^{-2}$, corresponding to a lower limit on the magnetic field strength of $B \gtrsim 0.81 \, \mu G$.

We present $\sim0.10^{\prime\prime}-$resolution Atacama Large Millimeter/submillimeter Array (ALMA) CO(2$-$1) imaging of the arcsecond-scale ($r \approx 150$ pc) dusty molecular disk in the giant elliptical galaxy NGC 3258. The data provide unprecedented resolution of cold gas disk kinematics within the dynamical sphere of influence of a supermassive black hole, revealing a quasi-Keplerian central increase in projected rotation speed rising from 280 km s$^{-1}$ at the disk's outer edge to $>400$ km s$^{-1}$ near the disk center. We construct dynamical models for the rotating disk and fit beam-smeared model CO line profiles directly to the ALMA data cube. Our models incorporate both flat disks and tilted-ring disks that provide a better fit of the mildly warped structure in NGC 3258. We show that the exceptional angular resolution of the ALMA data makes it possible to infer the host galaxy's mass profile within $r=150$ pc solely from the ALMA CO kinematics, without relying on optical or near-infrared imaging data to determine the stellar mass profile. Our model therefore circumvents any uncertainty in the black hole mass that would result from the substantial dust extinction in the galaxy's central region. The best model fit yields $M_\mathrm{BH} = 2.249\times10^9$ $M_\odot$ with a statistical model-fitting uncertainty of just 0.18\%, and systematic uncertainties of 0.62\% from various aspects of the model construction and 12\% from uncertainty in the distance to NGC 3258. This observation demonstrates the full potential of ALMA for carrying out highly precise measurements of $M_\mathrm{BH}$ in early-type galaxies containing circumnuclear gas disks

Compact elliptical galaxies (CEGs) are candidates for local analogs of the high-redshift "red nuggets" thought to represent the progenitors of today's early-type galaxies (ETGs). To address whether the structure of the dark matter (DM) halo in a CEG also reflects the extremely quiescent and isolated evolution of its stars, we use a new $\approx 122$ ks Chandra observation together with a shallow $\approx 13$ ks archival observation of the CEG Mrk 1216 to perform a hydrostatic equilibrium analysis of the luminous and relaxed X-ray plasma emission extending out to a radius $0.85r_{2500}$. We examine several DM model profiles and in every case obtain a halo concentration $(c_{200})$ that is a large positive outlier in the theoretical $\Lambda$CDM $c_{200}-M_{200}$ relation; i.e., ranging from $3.4\,\sigma - 6.3\, \sigma$ above the median $\Lambda$CDM relation in terms of the intrinsic scatter. The high value of $c_{200}$ we measure implies an unusually early formation time that firmly establishes the relic nature of the DM halo in Mrk 1216. The highly concentrated DM halo leads to a higher DM fraction and smaller total mass slope at $1R_e$ compared to nearby normal ETGs. In addition, the highly concentrated total mass profile of Mrk 1216 cannot be described by MOND without adding DM, and it deviates substantially from the Radial Acceleration Relation. Mrk 1216 contains $\approx 80\%$ of the cosmic baryon fraction within $r_{200}$. The radial profile of the ratio of cooling time to free-fall time varies within a narrow range $(t_c/t_{\rm ff}\approx 14-19)$ over a large central region ($r\le 10$ kpc) suggesting "precipitation-regulated AGN feedback" for a multiphase plasma, though presently there is little evidence for cool gas in Mrk 1216. The properties of Mrk 1216 are remarkably similar to those of the nearby fossil group NGC 6482.

There is mounting evidence that compact elliptical galaxies (CEGs) are local analogs of the high-redshift "red nuggets" thought to represent progenitors of today's early-type galaxies (ETGs). We report the discovery of extended X-ray emission from a hot interstellar / intragroup medium in two CEGs, Mrk 1216 and PGC 032873, using shallow archival Chandra observations. We find that PGC 032873 has an average gas temperature $k_BT=0.67\pm 0.06$ keV within a radius of 15 kpc, and a luminosity $L_{\rm x} = (1.8\pm 0.2)\times 10^{41}$ erg s$^{-1}$ within a radius of 100kpc. For Mrk 1216, which is closer and more luminous $[L_{\rm x}(\rm <100~kpc) = (12.1\pm 1.9)\times 10^{41}$ erg s$^{-1}]$, we performed a spatially resolved spectral analysis in 7 annuli out to a radius of 73 kpc. Using an entropy-based hydrostatic equilibrium (HE) procedure, we obtain a good constraint on the $H$-band stellar mass-to-light ratio, $M_{\rm stars}/L_H=1.33\pm 0.21$ solar, in good agreement with stellar dynamical (SD) studies, which supports the HE approximation. We obtain a density slope $2.22\pm 0.08$ within $R_e$ consistent with other CEGs and normal local ETGs, while the dark matter (DM) fraction within $R_e$, $f_{\rm DM}=0.20\pm 0.07$, is similar to local ETGs. We place a constraint on the SMBH mass, $M_{\rm BH} = (5\pm 4)\times 10^{9}\, M_{\odot}$, with a 90% upper limit of $M_{\rm BH} = 1.4\times 10^{10}\, M_{\odot}$, consistent with a recent SD measurement. We obtain a halo concentration $(c_{200}=17.5\pm 6.7)$ and mass [$M_{200} = (9.6\pm 3.7)\times 10^{12}\, M_{\odot}$], where $c_{200}$ exceeds the mean $\Lambda$CDM value ($\approx 7$), consistent with a system that formed earlier than the general halo population. We suggest that these galaxies, which reside in group-scale halos, should be classified as fossil groups. (Abridged)

We present results from an Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 2 program to map CO(2-1) emission in nearby early-type galaxies (ETGs) that host circumnuclear gas disks. We obtained $\sim0.3''-$resolution Band 6 observations of seven ETGs selected on the basis of dust disks in Hubble Space Telescope images. We detect CO emission in five at high signal-to-noise ratio with the remaining two only faintly detected. All CO emission is coincident with the dust and is in dynamically cold rotation. Four ETGs show evidence of rapid central rotation; these are prime candidates for higher-resolution ALMA observations to measure the black hole masses. In this paper we focus on the molecular gas and continuum properties. Total gas masses and H$_2$ column densities for our five CO-bright galaxies are on average $\sim10^8$ $M_\odot$ and $\sim10^{22.5}$ cm$^{-2}$ over the $\sim$kpc-scale disks, and analysis suggests that these disks are stabilized against gravitational fragmentation. The continuum emission of all seven galaxies is dominated by a central, unresolved source, and in five we also detect a spatially extended component. The $\sim$230 GHz nuclear continua are modeled as power laws ranging from $S_\nu \sim \nu^{-0.4}$ to $\nu^{1.6}$ within the observed frequency band. The extended continuum profiles of the two radio-bright (and CO-faint) galaxies are roughly aligned with their radio jet and suggests resolved synchrotron jets. The extended continua of the CO-bright disks are coincident with optically thick dust absorption and have spectral slopes that are consistent with thermal dust emission.

We revisit the massive isolated elliptical galaxy / fossil group NGC 6482 for which previous X-ray studies of a modest Chandra observation obtained a very uncertain, but also possibly very high, halo concentration. We present new measurements of the hot gas surface brightness, temperature, and iron abundance using the modest Chandra observation and a previously unpublished Suzaku observation, the latter of which allows measurements of the gas properties to be extended out to ~r_2500. By constructing hydrostatic equilibrium models of the gas with separate components for the gas, BCG stellar mass, and the dark matter (DM), we measure c_200 = 32.2 +/- 7.1 and M_200 = (4.5 +/- 0.6 x 10^12 M_sun using an NFW DM profile. For a halo of this mass, c_200 exceeds the mean value (7.1) expected for relaxed LCDM halos by $3.5 \sigma$ in terms of the observational error, and by $6 \sigma$ considering the intrinsic scatter in the LCDM c-M relation, which situates NGC 6482 as the most extreme outlier known for a fossil system. We explored several variants of adiabatic contraction (AC) models and, while the AC models provide fits of the same quality as the un-contracted models, they do have the following advantages: (1) smaller c_200 that is less of an outlier in the LCDM c-M relation, and (2) baryon fractions that agree better with the mean cosmic value. While the standard AC prescriptions yield a BCG stellar mass that is uncomfortably small compared to results from stellar population synthesis (SPS) models, a weaker AC variant that artificially shuts off cooling and star formation at z=2 yields the same stellar mass as the un-contracted models. For these reasons, we believe our X-ray analysis prefers this weaker AC variant applied to either an NFW or Einasto DM halo. Finally, the BCG stellar mass strongly favors SPS models with a Chabrier or Kroupa IMF over a Salpeter IMF. (Abridged)

In this second paper on the entire virial region of the relaxed fossil cluster RXJ1159+5531, we present a hydrostatic analysis of the hot intracluster medium (ICM). For a model consisting of ICM, stellar mass from the central galaxy (BCG), and an NFW dark matter (DM) halo, we obtain good descriptions of the projected radial profiles of ICM emissivity and temperature. The BCG stellar mass is clearly detected with M_star/L_K = 0.61 +/- 0.11 solar, consistent with stellar population synthesis models for a Milky-Way IMF. We obtain a halo concentration, c_200 =8.4 +/- 1.0, and virial mass, M_200 = 7.9 +/- 0.6 x 10^13 M_sun. For its mass, the inferred concentration is larger than most relaxed halos produced in cosmological simulations with Planck parameters, consistent with RXJ1159+5531 forming earlier than the general halo population. The baryon fraction at r_200, f_b,200 = 0.134 +/- 0.007, is slightly below the Planck value (0.155) for the universe. When we account for the stellar baryons associated with non-central galaxies and the uncertain intracluster light, f_b,200 increases by ~0.015, consistent with the cosmic value. Performing our analysis in the context of MOND still requires a large DM fraction (85.0% +/- 2.5% at r=100 kpc) similar to that obtained using the standard Newtonian approach. The detection of a plausible stellar BCG mass component distinct from the NFW DM halo in the total gravitational potential suggests that ~10^14 M_sun represents the mass scale above which dissipation is unimportant in the formation of the central regions of galaxy clusters. (Abridged)

We present Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 3 observations of CO(2-1) emission from the circumnuclear disk in the E/S0 galaxy NGC 1332 at 0.044" resolution. The disk exhibits regular rotational kinematics and central high-velocity emission (+/-500 km/s) consistent with the presence of a compact central mass. We construct models for a thin, dynamically cold disk in the gravitational potential of the host galaxy and black hole, and fit the beam-smeared model line profiles directly to the ALMA data cube. Model fits successfully reproduce the disk kinematics out to r=200 pc. Fitting models just to spatial pixels within projected r=50 pc of the nucleus (two times larger than the black hole's gravitational radius of influence), we find M_BH=6.64(-0.63,+0.65)*10^8 solar masses. This observation demonstrates ALMA's powerful capability to determine the masses of supermassive black holes by resolving gas kinematics on small angular scales in galaxy nuclei.

We present first results from a program of Atacama Large Millimeter/submillimeter Array (ALMA) CO(2-1) observations of circumnuclear gas disks in early-type galaxies. The program was designed with the goal of detecting gas within the gravitational sphere of influence of the central black holes. In NGC 1332, the 0.3"-resolution ALMA data reveal CO emission from the highly inclined (i~ 83 degrees) circumnuclear disk, spatially coincident with the dust disk seen in Hubble Space Telescope images. The disk exhibits a central upturn in maximum line-of-sight velocity reaching +-500 km/s relative to the systemic velocity, consistent with the expected signature of rapid rotation around a supermassive black hole. Rotational broadening and beam smearing produce complex and asymmetric line profiles near the disk center. We constructed dynamical models for the rotating disk and fitted the modeled CO line profiles directly to the ALMA data cube. Degeneracy between rotation and turbulent velocity dispersion in the inner disk precludes the derivation of strong constraints on the black hole mass, but model fits allowing for a plausible range in the magnitude of the turbulent dispersion imply a central mass in the range ~(4-8)x10^8 Msun. We argue that gas-kinematic observations resolving the black hole's projected radius of influence along the disk's minor axis will have the capability to yield black hole mass measurements that are largely insensitive to systematic uncertainties in turbulence or in the stellar mass profile. For highly inclined disks, this is a much more stringent requirement than the usual sphere-of-influence criterion.

Abell~1142 is a low-mass galaxy cluster at low redshift containing two comparable Brightest Cluster Galaxies (BCG) resembling a scaled-down version of the Coma Cluster. Our Chandra analysis reveals an X-ray emission peak, roughly 100 kpc away from either BCG, which we identify as the cluster center. The emission center manifests itself as a second beta-model surface brightness component distinct from that of the cluster on larger scales. The center is also substantially cooler and more metal rich than the surrounding intracluster medium (ICM), which makes Abell 1142 appear to be a cool core cluster. The redshift distribution of its member galaxies indicates that Abell 1142 may contain two subclusters with each containing one BCG. The BCGs are merging at a relative velocity of ~1200 km/s. This ongoing merger may have shock-heated the ICM from ~ 2 keV to above 3 keV, which would explain the anomalous L_X--T_X scaling relation for this system. This merger may have displaced the metal-enriched "cool core" of either of the subclusters from the BCG. The southern BCG consists of three individual galaxies residing within a radius of 5 kpc in projection. These galaxies should rapidly sink into the subcluster center due to the dynamical friction of a cuspy cold dark matter halo.

Previous analysis of the fossil-group/cluster RXJ1159+5531 with X-ray observations from a central Chandra pointing and an offset-North Suzaku pointing indicate a radial intracluster medium (ICM) entropy profile at the virial radius ($R_{\rm vir}$) consistent with predictions from gravity-only cosmological simulations, in contrast to other cool-core clusters. To examine the generality of these results, we present three new Suzaku observations that, in conjunction with the North pointing, provide complete azimuthal coverage out to $R_{\rm vir}$. With two new Chandra ACIS-I observations overlapping the North Suzaku pointing, we have resolved $\gtrsim$50\% of the cosmic X-ray background there. We present radial profiles of the ICM density, temperature, entropy, and pressure obtained for each of the four directions. We measure only modest azimuthal scatter in the ICM properties at $R_{\rm 200}$ between the Suzaku pointings: 7.6\% in temperature and 8.6\% in density, while the systematic errors can be significant. The temperature scatter, in particular, is lower than that studied at $R_{\rm 200}$ for a small number of other clusters observed with Suzaku. These azimuthal measurements verify that RXJ1159+5531 is a regular, highly relaxed system. The well-behaved entropy profiles we have measured for RXJ1159+5531 disfavor the weakening of the accretion shock as an explanation of the entropy flattening found in other cool-core clusters but is consistent with other explanations such as gas clumping, electron-ion non-equilibrium, non-thermal pressure support, and cosmic ray acceleration. Finally, we mention that the large-scale galaxy density distribution of RXJ1159+5531 seems to have little impact on its gas properties near $R_{\rm vir}$.

Highly ionized, z=0 metal absorption lines detected in the X-ray spectra of background active galactic nuclei (AGNs) provide an effective method to probe the hot ($T\sim10^6$ K) gas and its metal content in and around the Milky Way. We present an all-sky survey of the $K_{\alpha}$ transition of the local O VII absorption lines obtained by Voigt-profile fitting archival XMM-Newton observations. A total of 43 AGNs were selected, among which 12 are BL Lac-type AGNs, and the rest are Seyfert 1 galaxies. At above the $3\sigma$ level the local O VII absorption lines were detected in 21 AGNs, among which 7 were newly discovered in this work. The sky covering fraction, defined as the ratio between the number of detections and the sample size, increases from at about 40% for all targets to 100% for the brightest targets, suggesting a uniform distribution of the O VII absorbers. We correlate the line equivalent width with the Galactic coordinates and do not find any strong correlations between these quantities. Some AGNs have warm absorbers that may complicate the analysis of the local X-ray absorber since the recession velocity can be compensated by the outflow velocity, especially for the nearby targets. We discuss the potential impact of the warm absorbers on our analysis. A comprehensive theoretical modelling of the X-ray absorbers will be presented in a later paper.

Weak spectral features in BL Lacertae objects (BL Lac) often provide a unique opportunity to probe the inner region of this rare type of active galactic nucleus. We present a Hubble Space Telescope/Cosmic Origins Spectrograph observation of the BL Lac H 2356-309. A weak Ly$\alpha$ emission line was detected. This is the fourth detection of a weak Ly$\alpha$ emission feature in the ultraviolet (UV) band in the so-called "high energy peaked BL Lacs", after Stocke et al. Assuming the line-emitting gas is located in the broad line region (BLR) and the ionizing source is the off-axis jet emission, we constrain the Lorentz factor ($\Gamma$) of the relativistic jet to be $\geq 8.1$ with a maximum viewing angle of 3.6$^\circ$. The derived $\Gamma$ is somewhat larger than previous measurements of $\Gamma \approx 3 - 5$, implying a covering factor of $\sim$ 3% of the line-emitting gas. Alternatively, the BLR clouds could be optically thin, in which case we constrain the BLR warm gas to be $\sim 10^{-5}\rm\ M_{\odot}$. We also detected two HI and one OVI absorption lines that are within $|\Delta v| < 150\rm\ km\ s^{-1}$ of the BL Lac object. The OVI and one of the HI absorbers likely coexist due to their nearly identical velocities. We discuss several ionization models and find a photoionization model where the ionizing photon source is the BL Lac object can fit the observed ion column densities with reasonable physical parameters. This absorber can either be located in the interstellar medium of the host galaxy, or in the BLR.

We report on the X-ray observation of a strong lensing selected group, SL2S J08544-0121, with a total mass of $2.4 \pm 0.6 \times 10^{14}$ $\rm{M_\odot}$ which revealed a separation of $124\pm20$ kpc between the X-ray emitting collisional gas and the collisionless galaxies and dark matter (DM), traced by strong lensing. This source allows to put an order of magnitude estimate to the upper limit to the interaction cross section of DM of 10 cm$^2$ g$^{-1}$. It is the lowest mass object found to date showing a DM-baryons separation and it reveals that the detection of bullet-like objects is not rare and confined to mergers of massive objects opening the possibility of a statistical detection of DM-baryons separation with future surveys.

The "missing baryons", in the form of warm-hot intergalactic medium (WHIM), are expected to reside in cosmic filamentary structures that can be traced by signposts such as large-scale galaxy superstructures. The clear detection of an X-ray absorption line in the Sculptor Wall demonstrated the success of using galaxy superstructures as a signpost to search for the WHIM. Here we present an XMM-Newton Reflection Grating Spectrometer (RGS) observation of the blazar Mkn 501, located in the Hercules Supercluster. We detected an O VII K$\alpha$ absorption line at the 98.7 % level ($2.5\sigma$) at the redshift of the foreground Hercules Supercluster. The derived properties of the absorber are consistent with theoretical expectations of the WHIM. We discuss the implication of our detection for the search for the "missing baryons". While this detection shows again that using signposts is a very effective strategy to search for the WHIM, follow-up observations are crucial both to strengthen the statistical significance of the detection and to rule out other interpretations. A local, $z\sim0$ O VII K$\alpha$ absorption line was also clearly detected at the $4\sigma$ level, and we discuss its implications for our understanding of the hot gas content of our Galaxy.

We explore the utility of narrow band X-ray surface photometry as a tool for making fully Bayesian, hydrostatic mass measurements of clusters of galaxies, groups and early-type galaxies. We demonstrate that it is sufficient to measure the surface photometry with the Chandra X-ray observatory in only three (rest frame) bands (0.5--0.9 keV, 0.9--2.0 keV and 2.0--7.0 keV) in order to constrain the temperature, density and abundance of the hot interstellar medium (ISM). Adopting parametrized models for the mass distribution and radial entropy profile and assuming spherical symmetry, we show that the constraints on the mass and thermodynamic properties of the ISM that are obtained by fitting data from all three bands simultaneously are comparable to those obtained by fitting similar models to the temperature and density profiles derived from spatially resolved spectroscopy, as is typically done. We demonstrate that the constraints can be significantly tightened when exploiting a recently derived, empirical relationship between the gas fraction and the entropy profile at large scales, eliminating arbitrary extrapolations at large radii. This "Scaled Adiabatic Model" (ScAM) is well suited to modest signal-to-noise data, and we show that accurate, precise measurements of the global system properties are inferred when employing it to fit data from even very shallow, snapshot X-ray observations. The well-defined asymptotic behaviour of the model also makes it ideally suited for use in Sunyaev-Zeldovich studies of galaxy clusters.

We present a combined X-ray, optical, and radio analysis of the galaxy group IC 1860 using the currently available Chandra and XMM data, literature multi-object spectroscopy data and GMRT data. The Chandra and XMM imaging and spectroscopy reveal two surface brightness discontinuities at 45 and 76 kpc shown to be consistent with a pair of cold fronts. These features are interpreted as due to sloshing of the central gas induced by an off-axis minor merger with a perturber. This scenario is further supported by the presence of a peculiar velocity of the central galaxy IC 1860 and the identification of a possible perturber in the optically disturbed spiral galaxy IC 1859. The identification of the perturber is consistent with the comparison with numerical simulations of sloshing. The GMRT observation at 325 MHz shows faint, extended radio emission contained within the inner cold front, as seen in some galaxy clusters hosting diffuse radio mini-halos. However, unlike mini-halos, no particle reacceleration is needed to explain the extended radio emission, which is consistent with aged radio plasma redistributed by the sloshing. There is strong analogy of the X-ray and optical phenomenology of the IC 1860 group with two other groups, NGC 5044 and NGC 5846, showing cold fronts. The evidence presented in this paper is among the strongest supporting the currently favored model of cold-front formation in relaxed objects and establishes the group scale as a chief environment to study this phenomenon.

Recent hydrostatic X-ray studies of the hot interstellar medium (ISM) in early-type galaxies underestimate the gravitating mass as compared to stellar dynamics, implying modest, but significant deviations from exact hydrostatic equilibrium. We present a method for combining X-ray measurements and stellar dynamical constraints in the context of Bayesian statistics that allows the radial distribution of the implied nonthermal pressure or bulk motions in the hot ISM to be constrained. We demonstrate the accuracy of the method with hydrodynamical simulations tailored to produce a realistic galaxy model. Applying the method to the nearby elliptical galaxy NGC4649, we find a significant but subdominant nonthermal pressure fraction (0.27+/-0.06) in the central (<5 kpc) part of the galaxy, similar to the level of deviations from hydrostatic equilibrium expected in galaxy clusters. Plausible sources of systematic error, if important, may reduce this fraction. This would imply >360 km/s random turbulence or a magnetic field B=(39+/-6)(n_e/0.1 cm^{-3})^{0.59+/-0.09} muG, whereas gas rotation alone is unlikely to explain the detailed nonthermal profile. Future observations with Astro-H will allow turbulence or gas rotation at this level to be detected.

The Elliptical Isolated X-ray (ElIXr) Galaxy Survey is a volume-limited (<110Mpc) study of optically selected, isolated, Lstar elliptical galaxies, to provide an X-ray census of galaxy-scale (virial mass, Mvir < 1e13 Msun) objects, and identify candidates for detailed hydrostatic mass modelling. In this paper, we present a Chandra and XMM study of one such candidate, NGC1521, and constrain its distribution of dark and baryonic matter. We find a morphologically relaxed hot gas halo, extending almost to R500, that is well described by hydrostatic models similar to the benchmark, baryonically closed, Milky Way-mass elliptical galaxy NGC720. We obtain good constraints on the enclosed gravitating mass (M500=3.8e12+/-1e12 Msun, slightly higher than NGC\thin 720), and baryon fraction (fb500=0.13+/-0.03). We confirm at 8.2-sigma the presence of a dark matter (DM) halo consistent with LCDM. Assuming a Navarro-Frenk-White DM profile, our self-consistent, physical model enables meaningful constraints beyond R500, revealing that most of the baryons are in the hot gas. Within the virial radius, fb is consistent with the Cosmic mean, suggesting that the predicted massive, quasi-hydrostatic gas halos may be more common than previously thought. We confirm that the DM and stars conspire to produce an approximately powerlaw total mass profile (rho \propto r^-alpha) that follows the recently discovered scaling relation between alpha and optical effective radius. Our conclusions are insensitive to modest, observationally motivated, deviations from hydrostatic equilibrium. Finally, after correcting for the enclosed gas fraction, the entropy profile is close to the self-similar prediction of gravitational structure formation simulations, as observed in massive galaxy clusters.

We present a Chandra, Suzaku and Rosat study of the hot Intra Group Medium (IGrM) of the relaxed fossil group/ poor cluster RXJ1159+5531. This group exhibits an advantageous combination of flat surface brightness profile, high luminosity and optimal distance, allowing the gas to be detected out to the virial radius (Rvir=R108=1100 kpc) in a single Suzaku pointing, while the complementary Chandra data reveal a round morphology and relaxed IGrM image down to kpc scales. We measure the IGrM entropy profile over 3 orders of magnitude in radius, including 3 data bins beyond 0.5R200 with azimuthal coverage (>30%). We find no evidence that the profile flattens at large scales (>R500), and when corrected for the enclosed gas fraction, the entropy profile is very close to the predictions from self-similar structure formation simulations, as seen in massive clusters. Within Rvir, we measure a baryon fraction of 0.17+/-0.02, consistent with the Cosmological value. These results are in sharp contrast to the gas behaviour at large scales recently reported in the Virgo and Perseus clusters, and indicate that substantial gas clumping cannot be ubiquitous near Rvir, at least in highly evolved (fossil) groups.

This is the second of two papers investigating the spherical averaging of ellipsoidal galaxy clusters in the context of X-ray and Sunyaev-Zel'dovich (SZ) observations. In the present study we quantify the orientation-average bias and scatter in observables that result from spherically averaging clusters described by ellipsoidal generalizations of the NFW profile or a nearly scale-free logarithmic potential. Although the mean biases are small and mostly <1%, the flattest cluster models generally have a significant mean bias; i.e., averaging over all orientations does not always eliminate projection biases. Substantial biases can result from different viewing orientations, where the integrated Compton-y parameter (Y_SZ) and the concentration have the largest scatter (as large as sigma ~10% for Y_SZ), and the emission-weighted temperature (T_X) has the smallest (sigma < ~0.5%). The very small scatter for T_X leads to Y_X and M_gas having virtually the same orientation biases. Substantial scatter is expected for individual clusters (up to sigma ~8%) in the correlation between Y_SZ and Y_X in comparison to the small mean bias (sigma < ~1%) applicable to a random sample of clusters of sufficient size. For ellipsoidal NFW models we show that the orientation bias for the total cluster mass attains a minimum near the radius r_2500 so that the spherically averaged mass computed at this radius is always within ~0.5% of the true value for any orientation. Finally, to facilitate the accounting for orientation bias in X-ray and SZ cluster studies, we provide cubic polynomial approximations to the mean orientation bias and 1-sigma scatter for each cluster observable as a function of axial ratio for the ellipsoidal NFW models.

This is the first of two papers investigating the deprojection and spherical averaging of ellipsoidal galaxy clusters. We specifically consider applications to hydrostatic X-ray and Sunyaev-Zel'dovich (SZ) studies, though many of the results also apply to isotropic dispersion-supported stellar dynamical systems. Here we present analytical formulas for galaxy clusters described by a gravitational potential that is a triaxial ellipsoid of constant shape and orientation. For this model type we show that the mass bias due to spherically averaging X-ray observations is independent of the temperature profile, and for the special case of a scale-free logarithmic potential, there is exactly zero mass bias for any shape, orientation, and temperature profile. The ratio of spherically averaged intracluster medium (ICM) pressures obtained from SZ and X-ray measurements depends only on the ICM intrinsic shape, projection orientation, and H_0, which provides another illustration of how cluster geometry can be recovered through a combination of X-ray and SZ measurements. We also demonstrate that Y_SZ and Y_X have different biases owing to spherical averaging, which leads to an offset in the spherically averaged Y_SZ - Y_X relation. A potentially useful application of the analytical formulas presented is to assess the error range of an observable (e.g., mass, Y_SZ) accounting for deviations from assumed spherical symmetry, without having to perform the ellipsoidal deprojection explicitly. Finally, for dedicated ellipsoidal studies, we also generalize the spherical onion peeling method to the triaxial case for a given shape and orientation.

(Abriged) Assuming that the hydrostatic equilibrium holds between the intracluster medium and the gravitational potential, we constrain the NFW profiles in a sample of 44 X-ray luminous galaxy clusters observed with XMM-Newton in the redshift range 0.1-0.3. We evaluate several systematic uncertainties that affect our reconstruction of the X-ray masses. We measure the concentration c200, the dark mass M200 and the gas mass fraction within R500 in all the objects of our sample, providing the largest dataset of mass parameters for galaxy clusters in this redshift range. We confirm that a tight correlation between c200 and M200 is present and in good agreement with the predictions from numerical simulations and previous observations. When we consider a subsample of relaxed clusters that host a Low-Entropy-Core (LEC), we measure a flatter c-M relation with a total scatter that is lower by 40 per cent. From the distribution of the estimates of c200 and M200, with associated statistical (15-25%) and systematic (5-15%) errors, we use the predicted values from semi-analytic prescriptions calibrated through N-body numerical runs and measure sigma_8*Omega_m^(0.60+-0.03)= 0.45+-0.01 (at 2 sigma level, statistical only) for the subsample of the clusters where the mass reconstruction has been obtained more robustly, and sigma_8*Omega_m^(0.56+-0.04) = 0.39+-0.02 for the subsample of the 11 more relaxed LEC objects. With the further constraint from the fgas distribution in our sample, we break the degeneracy in the sigma_8-Omega_m plane and obtain the best-fit values sigma_8~1.0+-0.2 (0.75+-0.18 when the subsample of the more relaxed objects is considered) and Omega_m = 0.26+-0.01.

We review X-ray constraints on dark matter in giant elliptical galaxies (10^{12} M_sun <~ M_vir <~ 10^{13} M_sun) obtained using the current generation of X-ray satellites, beginning with an overview of the physics of the hot interstellar medium and mass modeling methodology. Dark matter is now firmly established in many galaxies, with inferred NFW concentration parameters somewhat larger than the mean theoretical relation. X-ray observations confirm that the total mass profile (baryons+DM) is close to isothermal (M ~ r), and new evidence suggests a more general power-law relation for the slope of the total mass profile that varies with the stellar half-light radius. We also discuss constraints on the baryon fraction, super-massive black holes, and axial ratio of the dark matter halo. Finally, we review constraints on non-thermal gas motions and discuss the accuracy of the hydrostatic equilibrium approximation in elliptical galaxies.

Since the launch of the Einstein X-ray Observatory in the 1970s, a number of broad absorption features have been reported in the X-ray spectra of BL Lac objects. These features are often interpreted as arising from high velocity outflows intrinsic to the BL Lac object, therefore providing important information about the inner environment around the central engine. However, such absorption features have not been observed more recently with high-resolution X-ray telescopes such as Chandra and XMM-Newton. In this paper, we report the detection of a transient X-ray absorption feature intrinsic to the BL Lac object H 2356-309 with the Chandra X-ray Telescope. This BL Lac object was observed during XMM cycle 7, Chandra cycle 8 and 10, as part of our campaign to investigate X-ray absorption produced by the warm-hot intergalactic medium (WHIM) residing in the foreground large scale superstructure. During one of the 80 ksec, Chandra cycle 10 observations, a transient absorption feature was detected at 3.3-sigma (or 99.9% confidence level, accounting for the number of trials), which we identify as the OVIII K-alpha line produced by an absorber intrinsic to the BL Lac object. None of the other 11 observations showed this line. We constrain the ionization parameter (25 <~ Xi <~ 40) and temperature (10^5 < T < 2.5 10^7 K) of the absorber. This absorber is likely produced by an outflow with a velocity up to 1,500 km/s. There is a suggestion of possible excess emission on the long-wavelength side of the absorption line; however, the derived properties of the emission material are very different from those of the absorption material, implying it is unlikely a typical P Cygni-type profile.

We present a study of the dark and luminous matter in the isolated elliptical galaxy NGC720, based on deep X-ray observations made with Chandra and Suzaku. The gas is reliably measured to ~R2500, allowing us to place good constraints on the enclosed mass and baryon fraction (fb) within this radius (M2500=1.6e12+/-0.2e12 Msun, fb(2500)=0.10+/-0.01; systematic errors are <~20%). The data indicate that the hot gas is close to hydrostatic, which is supported by good agreement with a kinematical analysis of the dwarf satellite galaxies. We confirm a dark matter (DM) halo at ~20-sigma. Assuming an NFW DM profile, our physical model for the gas distribution enables us to obtain meaningful constraints at scales larger than R2500, revealing that most of the baryons are in the hot gas. We find that fb within Rvir is consistent with the Cosmological value, confirming theoretical predictions that a ~Milky Way-mass (Mvir=3.1e12+/-0.4e12 Msun) galaxy can sustain a massive, quasi-hydrostatic gas halo. While fb is higher than the cold baryon fraction typically measured in similar-mass spiral galaxies, both the gas fraction (fg) and fb in NGC720 are consistent with an extrapolation of the trends with mass seen in massive galaxy groups and clusters. After correcting for fg, the entropy profile is close to the self-similar prediction of gravitational structure formation simulations, as observed in galaxy clusters. Finally, we find a strong heavy metal abundance gradient in the ISM similar to those observed in massive galaxy groups.

Recent work based on a global measurement of the ICM properties find evidence for an increase of the iron abundance in galaxy clusters with temperature around 2-4 keV up to a value about 3 times larger than that typical of very hot clusters. We have started a study of the metal distribution in these objects from the sample of Baumgartner et al. (2005), aiming at resolving spatially the metal content of the ICM. We report here on a 42ks XMM observation of the first object of the sample, the cluster Abell 2028. The XMM observation reveals a complex structure of the cluster over scale of 300 kpc, showing an interaction between two sub-clusters in cometary-like configurations. At the leading edges of the two substructures cold fronts have been detected. The core of the main subcluster is likely hosting a cool corona. We show that a one-component fit for this region returns a biased high metallicity. This inverse iron bias is due to the behavior of the fitting code in shaping the Fe-L complex. In presence of a multi-temperature structure of the ICM, the best-fit metallicity is artificially higher when the projected spectrum is modeled with a single temperature component and it is not related to the presence of both Fe-L and Fe-K emission lines in the spectrum. After accounting for the bias, the overall abundance of the cluster is consistent with the one typical of hotter, more massive clusters. We caution the interpretation of high abundances inferred when fitting a single thermal component to spectra derived from relatively large apertures in 3-4 keV clusters, because the inverse iron bias can be present. Most of the inferences trying to relate high abundances in 3-4 keV clusters to fundamental physical processes will likely have to be revised.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We make use of a 500ks Chandra HRC-S/LETG spectrum of the blazar H2356-309, combined with a lower S/N spectrum of the same target, to search for the presence of warm-hot absorbing gas associated with two Large-Scale Structures (LSSs) crossed by this sightline at z=0.062 (the Pisces-Cetus Supercluster, PCS) and at z=0.128 ("Farther Sculptor Wall", FSW). No statistically significant (>=3sigma) individual absorption is detected from any of the strong He- or H-like transitions of C, O and Ne at the redshifts of the structures. However we are still able to constrain the physical and geometrical parameters of the associated putative absorbing gas, by performing joint spectral fit of marginal detections and upper limits of the strongest expected lines with our self-consistent hybrid ionization WHIM spectral model. At the redshift of the PCS we identify a warm phase with logT=5.35_-0.13^+0.07 K and log N_H =19.1+/-0.2 cm^-2 possibly coexisting with a hotter and less significant phase with logT=6.9^+0.1_-0.8 K and log N_H=20.1^+0.3_-1.7 cm^-2 (1sigma errors). For the FSW we estimate logT=6.6_-0.2^+0.1 K and log N_H=19.8_-0.8^+0.4 cm^-2. Our constraints allow us to estimate the cumulative number density per unit redshifts of OVII WHIM absorbers. We also estimate the cosmological mass density obtaining Omega_b(WHIM)=(0.021^+0.031_-0.018) (Z/Z_sun)^-1, consistent with the mass density of the intergalactic 'missing baryons' for high metallicities.

In a previous paper we reported a 3-sigma detection of an absorption line from the Warm-Hot Intergalactic Medium (WHIM) using the Chandra and XMM X-ray grating spectra of the blazar H2356-309, the sight-line of which intercepts the Sculptor Wall, a large-scale superstructure of galaxies at z ~ 0.03. To verify our initial detection, we obtained a deep (500 ks), follow-up exposure of H2356-309 as part of the Cycle-10 Chandra Large Project Program. From a joint analysis of the Cycle-10 and previous (Cycle-8) Chandra grating data we detect the redshifted OVII WHIM line at a significance level of 3.4-sigma, a substantial improvement over the 1.7-sigma level reported previously when using only the Cycle-8 data. The significance increases to 4.0-sigma when the existing XMM grating data are included in the analysis, thus confirming at higher significance the existence of the line at the redshift of the Sculptor Wall with an equivalent width of 28.5+/-10.5 mA (90% confidence). We obtain a 90% lower limit on the OVII column density of 0.8 10^16 cm^-2 and a 90% upper limit on the Doppler-b parameter of 460 km/s. Assuming the absorber is uniformly distributed throughout the ~ 15 Mpc portion of the blazar's sight-line that intercepts the Sculptor Wall, that the OVII column density is ~ 2 10^16 cm^-2 (corresponding to b > 150 km/s where the inferred column density is only weakly dependent on b), and that the oxygen abundance is 0.1 solar, we estimate a baryon over-density of ~ 30 for the WHIM, which is consistent with the peak of the WHIM mass fraction predicted by cosmological simulations. The clear detection of OVII absorption in the Sculptor Wall demonstrates the viability of using current observatories to study WHIM in the X-ray absorption spectra of blazars behind known large-scale structures.

We present a survey, using the Chandra X-ray observatory, of the central gravitating mass profiles in a sample of 10 galaxies, groups and clusters, spanning ~2 orders of magnitude in virial mass. We find the total mass distributions from ~0.2--10Re, where Re is the optical effective radius of the central galaxy, are remarkably similar to powerlaw density profiles. The negative logarithmic slope of the mass density profiles, alpha, systematically varies with Re, from alpha=2, for systems with Re~4kpc to alpha=1.2 for systems with Re>30kpc. Departures from hydrostatic equilibrium are likely to be small and cannot easily explain this trend. We show that the conspiracy between the baryonic (Sersic) and dark matter (NFW/ Einasto) components required to maintain a powerlaw total mass distribution naturally predicts an anti-correlation between alpha and Re that is very close to what is observed. The systematic variation of alpha with Re implies a dark matter fraction within Re that varies systematically with the properties of the galaxy in such a manner as to reproduce, without fine tuning, the observed tilt of the fundamental plane. We speculate that establishing a nearly powerlaw total mass distribution is therefore a fundamental feature of galaxy formation and the primary factor which determines the tilt of the fundamental plane.

We describe two-dimensional gasdynamical computations of the X-ray emitting gas in the rotating elliptical galaxy NGC 4649 that indicate an inflow of about one solar mass per year at every radius. Such a large instantaneous inflow cannot have persisted over a Hubble time. The central constant-entropy temperature peak recently observed in the innermost 150 parsecs is explained by compressive heating as gas flows toward the central massive black hole. Since the cooling time of this gas is only a few million years, NGC 4649 provides the most acutely concentrated known example of the cooling flow problem in which the time-integrated apparent mass that has flowed into the galactic core exceeds the total mass observed there. This paradox can be resolved by intermittent outflows of energy or mass driven by accretion energy released near the black hole. Inflowing gas is also required at intermediate kpc radii to explain the ellipticity of X-ray isophotes due to spin-up by mass ejected by stars that rotate with the galaxy and to explain local density and temperature profiles. We provide evidence that many luminous elliptical galaxies undergo similar inflow spin-up. A small turbulent viscosity is required in NGC 4649 to avoid forming large X-ray luminous disks that are not observed, but the turbulent pressure is small and does not interfere with mass determinations that assume hydrostatic equilibrium.

Recent work based on a global measurement of the ICM properties find evidence for an increase of the iron abundance in galaxy clusters with temperature around 2-4 keV. We have undertaken a study of the metal distribution in nearby clusters in this temperature range, aiming at resolving spatially the metal content of the ICM. The XMM observation of the first object of the sample, the cluster Abell 2028, reveals a complex structure of the cluster over scale of ~ 300 kpc, showing an interaction between two sub-clusters in a ``cometary'' configuration. We show that a naive one-component fit for the core of Abell 2028 returns a biased high metallicity. This is due to the inverse iron-bias, which is not related to the presence in the spectrum of both Fe-L and Fe-K emission lines but to the behavior of the fitting code in shaping the Fe-L complex of a one temperature component to adjust to the multi-temperature structure of the projected spectrum.

We present AGN feedback in the interesting cases of two groups: AWM 4 and NGC 5044. AWM 4 is characterized by a combination of properties which seems to defy the paradigm for AGN heating in cluster cores: a flat inner temperature profile indicative of a past, major heating episode which completely erased the cool core, as testified by the high central cooling time (> 3 Gyrs) and by the high central entropy level (~ 50 keV cm^2), and yet an active central radio galaxy with extended radio lobes out to 100 kpc, revealing recent feeding of the central massive black hole. A recent Chandra observation has revealed the presence of a compact cool corona associated with the BCG, solving the puzzle of the apparent lack of low entropy gas surrounding a bright radio source, but opening the question of its origin. NGC 5044 shows in the inner 10 kpc a pair of cavities together with a set of bright filaments. The cavities are consistent with a recent AGN outburst as also indicated by the extent of dust and H_alpha emission even though the absence of extended 1.4 GHz emission remains to be explained. The soft X-ray filaments coincident with H_alpha and dust emission are cooler than those which do not correlate with optical and infrared emission, suggesting that dust-aided cooling can contribute to the overall cooling. For the first time sloshing cold fronts at the scale of a galaxy group have been observed in this object.

We present new mass measurements for the supermassive black holes (SMBHs) in the centres of three early-type galaxies. The gas pressure in the surrounding, hot interstellar medium (ISM) is measured through spatially resolved spectroscopy with the Chandra X-ray observatory, allowing the SMBH mass (Mbh) to be inferred directly under the hydrostatic approximation. This technique does not require calibration against other SMBH measurement methods and its accuracy depends only on the ISM being close to hydrostatic, which is supported by the smooth X-ray isophotes of the galaxies. Combined with results from our recent study of the elliptical galaxy NGC4649, this brings to four the number of galaxies with SMBHs measured in this way. Of these, three already have mass determinations from the kinematics of either the stars or a central gas disc, and hence join only a handful of galaxies with Mbh measured by more than one technique. We find good agreement between the different methods, providing support for the assumptions implicit in both the hydrostatic and the dynamical models. The stellar mass-to-light ratios for each galaxy inferred by our technique are in agreement with the predictions of stellar population synthesis models assuming a Kroupa initial mass function (IMF). This concurrence implies that no more than ~10-20% of the ISM pressure is nonthermal, unless there is a conspiracy between the shape of the IMF and nonthermal pressure. Finally, we compute Bondi accretion rates, finding that the two galaxies with the highest rates exhibit little evidence of X-ray cavities, suggesting that the correlation with the AGN jet power takes time to be established.

We present XMM RGS and Chandra LETG observations of the blazar, H 2356-309, located behind the Sculptor Wall, a large-scale galaxy structure expected to harbor high-density Warm-Hot Intergalactic Medium (WHIM). Our simultaneous analysis of the RGS and LETG spectra yields a 3-sigma detection of the crucial redshifted O vii K-alpha line with a column density (>~ 10^{16} cm^{-2}) consistent with similar large-scale structures produced in cosmological simulations. This represents the first detection of non-local WHIM from X-ray absorption studies where XMM and Chandra data are analyzed simultaneously and the absorber redshift is already known, thus providing robust evidence for the expected repository of the "missing baryons".

We demonstrate that two approximations to the chi^2 statistic as popularly employed by observational astronomers for fitting Poisson-distributed data can give rise to intrinsically biased model parameter estimates, even in the high counts regime, unless care is taken over the parameterization of the problem. For a small number of problems, previous studies have shown that the fractional bias introduced by these approximations is often small when the counts are high. However, we show that for a broad class of problem, unless the number of data bins is far smaller than \sqrt{N_c}, where N_c is the total number of counts in the dataset, the bias will still likely be comparable to, or even exceed, the statistical error. Conversely, we find that fits using Cash's C-statistic give comparatively unbiased parameter estimates when the counts are high. Taking into account their well-known problems in the low count regime, we conclude that these approximate chi^2 methods should not routinely be used for fitting an arbitrary, parameterized model to Poisson-distributed data, irrespective of the number of counts per bin, and instead the C-statistic should be adopted. We discuss several practical aspects of using the C-statistic in modelling real data. We illustrate the bias for two specific problems, measuring the count-rate from a lightcurve and obtaining the temperature of a thermal plasma from its X-ray spectrum measured with the Chandra X-ray observatory. In the context of X-ray astronomy, we argue the bias could give rise to systematically mis-calibrated satellites and a ~5-10% shift in galaxy cluster scaling relations.

We present a two-dimensional analysis of the bright nearby galaxy group NGC 5044 using the currently available Chandra and XMM data. In the inner 10 kpc a pair of cavities are evident together with a set of bright X-ray filaments. If the cavities are interpreted as gas displaced by relativistic plasma inflated by an AGN, even in the absence of extended 1.4 GHz emission, this would be consistent with a recent outburst as also indicated by the extent of dust and H_alpha emission. The soft X-ray filaments coincident with H_alpha and dust emission are cooler than the ones which do not correlate with optical and infrared emission. We suggest that dust-aided cooling contributes to form warm (T =10^4 K) gas, emitting H_alpha radiation. At 31 kpc and 67 kpc a pair of cold fronts are present, indicative of sloshing due to a dynamical perturbation caused by accretion of a less massive group, also suggested by the peculiar velocity of the brightest galaxy NGC 5044 with respect to the mean group velocity.

Cosmological simulations of galaxy clusters typically find that the weight of a cluster at a given radius is not balanced entirely by the thermal gas pressure of the hot ICM, with theoretical studies emphasizing the role of random turbulent motions to provide the necessary additional pressure support. Using a set of high-resolution, hydrodynamical simulations of galaxy clusters that include radiative cooling and star formation, we find instead that in the most relaxed clusters rotational support exceeds that from random turbulent motions for radii 0.1 - 0.5 r_500, and that the observed clusters are much rounder than the simulated, relaxed clusters within ~ 0.4 r_500. Moreover, while the observed clusters display an average ellipticity profile that does not vary significantly with radius, the ellipticity of the relaxed CDM clusters declines markedly with increasing radius, suggesting that the ICM of the observed clusters rotates less rapidly than that of the relaxed CDM clusters out to ~ 0.6 r_500. We also find the ellipticity profile of a simulated cluster without radiative cooling is in much better agreement with the observations, implying that over-cooling has a substantial impact on the gas dynamics and morphology out to larger radii than previously recognized. It also suggests that the 10%-20% systematic errors from non-thermal gas pressure support reported for simulated cluster masses, obtained from fitting simulated X-ray data over large radial ranges within r_500, may need to be revised downward. These results demonstrate the utility of X-ray ellipticity profiles as a probe of ICM rotation and over-cooling which should be used to constrain future cosmological cluster simulations.

We present a Chandra survey of LMXBs in 24 early-type galaxies. Correcting for detection incompleteness, the X-ray luminosity function (XLF) of each galaxy is consistent with a powerlaw with negative logarithmic differential slope, beta~2.0. However, beta strongly correlates with incompleteness, indicating the XLF flattens at low-Lx. The composite XLF is well-fitted by a powerlaw with a break at 2.21(+0.65,-0.56)E38 erg/s and beta=1.40(+0.10,-0.13) and 2.84(+0.39,-0.30) below and above it, respectively. The break is close to the Eddington limit for a 1.4Msun neutron-star, but the XLF shape rules out its representing the division between neutron-star and black-hole systems. Although the XLFs are similar, we find evidence of some variation between galaxies. The high-Lx XLF slope does not correlate with age, but may correlate with [alpha/Fe]. Considering only LMXBs with Lx>1E37 erg/s, matching the LMXBs with globular clusters (GCs) identified in HST observations of 19 of the galaxies, we find the probability a GC hosts an LMXB is proportional to LGC^alpha ZFe^gamma} where alpha=1.01+/-0.19 and gamma=0.33+/-0.11. Correcting for GC luminosity and colour effects, and detection incompleteness, we find no evidence that the fraction of LMXBs with Lx>1e37 erg/s in GCs (40%), or the fraction of GCs hosting LMXBs (~6.5%) varies between galaxies. The spatial distribution of LMXBs resembles that of GCs, and the specific frequency of LMXBs is proportional to the GC specific luminosity, consistent with the hypothesis that all LMXBs form in GCs. If the LMXB lifetime is tau and the duty cycle is Fd, our results imply ~1.5 (tau/1E8 yr)^-1 /Fd LMXBs are formed per Gyr per GC and we place an upper limit of 1 active LMXB in the field per 3.4E9Lsun of V-band luminosity.

We present a Chandra study of the hot ISM in the giant elliptical galaxy NGC4649. In common with other group-centred ellipticals, its temperature profile rises with radius in the outer parts of the galaxy, from ~0.7keV at 2kpc to ~0.9keV by 20kpc. However, within the central ~2kpc the trend reverses and the temperature peaks at ~1.1keV within the innermost 200pc. Under the assumption of hydrostatic equilibrium, we demonstrate that the central temperature spike arises due to the gravitational influence of a quiescent central super-massive black hole. We constrain the black hole mass (MBH) to $(3.35^{+0.67}_{-0.95})\times 10^9$Msun (90% confidence), in good agreement with stellar kinematics measurements. This is the first direct measurement of MBH based on studies of hydrostatic X-ray emitting gas, which are sensitive to the most massive black holes, and is a crucial validation of both mass-determination techniques. This agreement clearly demonstrates the gas must be close to hydrostatic, even in the very centre of the galaxy, which is consistent with the lack of morphological disturbances in the X-ray image. NGC4649 is now one of only a handful of galaxies for which MBH has been measured by more than one method. At larger radii, we were able to decompose the gravitating mass profile into stellar and dark matter (DM) components. Unless one accounts for the DM, a standard Virial analysis of the stars dramatically over-estimates the stellar mass of the galaxy. We find the measured J-band stellar mass-to-light ratio, 1.37+/-0.10 Msun/Lsun, is in good agreement with simple stellar population model calculations for this object.

Although very successful in explaining the observed conspiracy between the baryonic distribution and the gravitational field in spiral galaxies without resorting to dark matter (DM), the modified Newtonian dynamics (MOND) paradigm still requires DM in X-ray bright systems. Here, to get a handle on the distribution and importance of this DM, and thus on its possible form, we deconstruct the mass profiles of 26 X-ray emitting systems in MOND, with temperatures ranging from 0.5 to 9 keV. Initially we compute the MOND dynamical mass as a function of radius, then subtract the known gas mass along with a component of galaxies which includes the cD galaxy with $M/L_K=1$. Next we test the compatibility of the required DM with ordinary massive neutrinos at the experimental limit of detection ($m_{\nu}=2$ eV), with density given by the Tremaine-Gunn limit. Even by considering that the neutrino density stays constant and maximal within the central 100 or 150 kpc (which is the absolute upper limit of a possible neutrino contribution there), we show that these neutrinos can never account for the required DM within this region. The natural corollary of this finding is that, whereas clusters (T $\ga$ 3 keV) might have most of their mass accounted for if ordinary neutrinos have a 2 eV mass, groups (T $\lsim$ 2 keV) cannot be explained by a 2 eV neutrino contribution. This means that, for instance, cluster baryonic dark matter (CBDM, Milgrom 2007) or even sterile neutrinos would present a more satisfactory solution to the problem of missing mass in MOND X-ray emitting systems.

The core of the relaxed cluster AWM 4 is characterized by a unique combination of properties which defy a popular scenario for ANG heating of cluster cores. A flat inner temperature profile is indicative of a past, major heating episode which completely erased the cool core, as testified by the high central cooling time (~ 3 Gyr) and by the high central entropy level (~ 60 keV cm^2). Yet the presence of a 1.4 GHz active central radio galaxy with extended radio lobes out to 100 kpc, reveals recent feeding of the central massive black hole. A system like AWM 4 should have no radio emission at all if only feedback from the cooling hot gas regulates the AGN activity.

We present radial mass profiles within 0.3 r_vir for 16 relaxed galaxy groups-poor clusters (kT range 1-3 keV) selected for optimal mass constraints from the Chandra and XMM data archives. After accounting for the mass of hot gas, the resulting mass profiles are described well by a two-component model consisting of dark matter (DM), represented by an NFW model, and stars from the central galaxy. The stellar component is required only for 8 systems, for which reasonable stellar mass-to-light ratios (M/L_K) are obtained, assuming a Kroupa IMF. Modifying the NFW dark matter halo by adiabatic contraction does not improve the fit and yields systematically lower M/L_K. In contrast to previous results for massive clusters, we find that the NFW concentration parameter (c_vir) for groups decreases with increasing M_vir and is inconsistent with no variation at the 3 sigma level. The normalization and slope of the c_vir-M_vir relation are consistent with the standard LambdaCDM cosmological model with sigma_8 = 0.9. The small intrinsic scatter measured about the c_vir-M_vir relation implies the groups represent preferentially relaxed, early forming systems. The mean gas fraction (f =0.05 +/- 0.01) of the groups measured within an overdensity Delta=2500 is lower than for hot, massive clusters, but the fractional scatter (sigma_f/f=0.2) for groups is larger, implying a greater impact of feedback processes on groups, as expected.

We present the concentration (c)-virial mass (M) relation of 39 galaxy systems ranging in mass from individual early-type galaxies up to the most massive galaxy clusters, (0.06-20) x 10^{14} M_sun. We selected for analysis the most relaxed systems possessing the highest quality data currently available in the Chandra and XMM public data archives. A power-law model fitted to the X-ray c-M relation requires at high significance (6.6 sigma) that c decreases with increasing M, which is a general feature of CDM models. The median and scatter of the c-M relation produced by the flat, concordance LCDM model (Omega_m=0.3, sigma_8=0.9) agrees with the X-ray data provided the sample is comprised of the most relaxed, early forming systems, which is consistent with our selection criteria. Holding the rest of the cosmological parameters fixed to those in the concordance model the c-M relation requires 0.76< sigma_8 <1.07 (99% conf.), assuming a 10% upward bias in the concentrations for early forming systems. The tilted, low-sigma_8 model suggested by a new WMAP analysis is rejected at 99.99% confidence, but a model with the same tilt and normalization can be reconciled with the X-ray data by increasing the dark energy equation of state parameter to w ~ -0.8. When imposing the additional constraint of the tight relation between sigma_8 and Omega_m from studies of cluster abundances, the X-ray c-M relation excludes (>99% conf.) both open CDM models and flat CDM models with Omega_m ~1. This result provides novel evidence for a flat, low-Omega_m universe with dark energy using observations only in the local (z << 1) universe. Possible systematic errors in the X-ray mass measurements of a magnitude ~10% suggested by CDM simulations do not change our conclusions.

We report the discovery of a galaxy cluster serendipitously detected as an extended X-ray source in an offset observation of the group NGC 5044. The cluster redshift, z=0.281, determined from the optical spectrum of the brightest cluster galaxy, agrees with that inferred from the X-ray spectrum using the Fe K alpha complex of the hot ICM (z=0.27 +/- 0.01). Based on the 50 ks XMM observation, we find that within a radius of 383 kpc the cluster has an unabsorbed X-ray flux, f_X (0.5-2 keV) = 3.34 (+0.08, -0.13) x 10^{-13} erg/cm^2/s, a bolometric X-ray luminosity, L_X = 2.21 (+0.34, -0.19) x 10^{44} erg/s, kT = 3.57 +/- 0.12 keV, and metallicity, 0.60 +/- 0.09 solar. The cluster obeys the scaling relations for L_X and T observed at intermediate redshift. The mass derived from an isothermal NFW model fit is, M_vir = 3.89 +/- 0.35 x 10^{14} solar masses, with a concentration parameter, c = 6.7 +/- 0.4, consistent with the range of values expected in the concordance cosmological model for relaxed clusters. The optical properties suggest this could be a ``fossil cluster''.

The enormous range of X-ray luminosities among elliptical galaxies with similar optical luminosities can be attributed to large scatter and systematic variations in the dark halo mass determined from X-ray observations. The mean halo mass decreases sharply with both the X-ray and K-band luminosities. Smaller halos contain less diffuse hot gas and have lower X-ray luminosities. In addition the hot gas in low mass halos is more easily depleted by supernova-driven outflows, further lowering the X-ray luminosity.

We present an X-ray analysis of the radial mass profile of the radio-quiet galaxy cluster A2589 between 0.015-0.25 r_vir using an XMM-Newton observation. Except for a ~16 kpc shift of the X-ray center of the R=45-60 kpc annulus, A2589 possesses a remarkably symmetrical X-ray image and is therefore an exceptional candidate for precision studies of its mass profile by applying hydrostatic equilibrium. The total gravitating matter profile is well described by the NFW model (fractional residuals <~10%) with c_vir=6.1 +/- 0.3 and M_vir = 3.3 +/- 0.3 x 10^{14} M_sun (r_vir = 1.74 +/- 0.05 Mpc) in excellent agreement with LCDM. When the mass of the hot ICM is subtracted from the gravitating matter profile, the NFW model fitted to the resulting dark matter (DM) profile produces essentially the same result. However, if a component accounting for the stellar mass (M_*) of the cD galaxy is included, then the NFW fit to the DM profile is substantially degraded in the central r ~50 kpc for reasonable M_*/L_V. Modifying the NFW DM halo by adiabatic contraction arising from the early condensation of stellar baryons in the cD galaxy further degrades the fit. The fit is improved substantially with a Sersic-like model recently suggested by high resolution N-body simulations but with an inverse Sersic index, alpha ~0.5, a factor of ~3 higher than predicted. We argue that neither random turbulent motions nor magnetic fields can provide sufficient non-thermal pressure support to reconcile the XMM mass profile with adiabatic contraction of a CDM halo assuming reasonable M_*/L_V. Our results support the scenario where, at least for galaxy clusters, processes during halo formation counteract adiabatic contraction so that the total gravitating mass in the core approximately follows the NFW profile.

We present a Chandra study of mass profiles in 7 elliptical galaxies, of which 3 have galaxy-scale and 4 group-scale halos demarcated at 1E13Msun. These represent the best available data for nearby objects with comparable X-ray luminosities. We measure ~flat mass-to-light (M/L) profiles within an optical half-light radius (Reff), rising by an order of magnitude at ~10Reff, which confirms the presence of dark matter (DM). The data indicate hydrostatic equilibrium, which is also supported by agreement with studies of stellar kinematics in elliptical galaxies. The data are well-fitted by a model comprising an NFW DM profile and a baryonic component following the optical light. The distribution of DM halo concentration parameters (c) versus Mvir agrees with LCDM predictions and our observations of bright groups. Concentrations are slightly higher than expected, which is most likely a selection effect. Omitting the stellar mass drastically increases c, possibly explaining large concentrations found by some past observers. The stellar M/LK agree with population synthesis models, assuming a Kroupa IMF. Allowing adiabatic compression (AC) of the DM halo by baryons made M/L more discrepant, casting some doubt on AC. Our best-fitting models imply total baryon fractions \~0.04--0.09, consistent with models of galaxy formation incorporating strong feedback. The groups exhibit positive temperature gradients, consistent with the "Universal" profiles found in other groups and clusters, whereas the galaxies have negative gradients, suggesting a change in the evolutionary history of the systems around Mvir=1E13 Msun.

Elliptical galaxies and their groups having the largest L_x/L_B lie close to the locus in the L_x,L_B diagram expected for closed systems with baryon fractions equal to the cosmic mean value, f_b = 0.16. The estimated baryon fractions for several of these galaxies/groups are also close to 0.16 when the gas density is extrapolated to the virial radius. Evidently they are the least massive baryonically closed systems. Gas retention in these groups implies that non-gravitational heating cannot exceed about 1 keV per particle, consistent with the heating required to produce the deviation of groups from the L_x - T correlation for more massive clusters. Isolated galaxies/groups with X-ray luminosities significantly lower than baryonically closed groups may have undermassive dark halos, overactive central AGNs, or higher star formation efficiencies. The virial mass and hot gas temperatures of nearly or completely closed groups correlate with the group X-ray luminosities and the optical luminosities of the group-centered elliptical galaxy, an expected consequence of their merging history. The ratio of halo mass to the mass of the central galaxy for X-ray luminous galaxy/groups is about 80.

We present the first in a series of papers studying the X-ray properties of 28 early type galaxies spanning ~3 orders of magnitude in L_x. We report emission-weighted constraints for Fe, O, Ne, Mg, Si, S and Ni. We find no evidence of the very sub-solar Z_Fe historically reported, confirming a trend in recent X-ray observations of bright galaxies and groups, nor do we find Z_Fe correlates with L_x. In general we find no evidence for multi-phase ISM, indicating multi-temperature fits found with ASCA were due to temperature gradients Chandra can resolve. We compare the stellar Z_Fe, estimated from simple stellar population model fits, to Z_Fe(ISM). Excepting one possible outlier we find no evidence that the gas is much more metal-poor than the stars and, in a few systems Z_Fe(ISM) is higher. In general the two exhibit similar Z_Fe, which is inconsistent with both galactic wind models and recent hierarchical chemical enrichment simulations. Adopting standard SNIa and SNII metal yields our abundance ratio constraints imply 66+/-11% of the Fe within the ISM was produced in SNIa, similar to the Solar neighbourhood and implying similar enrichment histories for the ISM in a spiral and in elliptical galaxies. Although these values are sensitive to systematic uncertainties in the SN yields, they are in good agreement with observations of more massive systems, indicating remarkable homology in the enrichment process operating from cluster scales to low-to-intermediate L_x galaxies. The data also uniformly exhibit the low Z_O/Z_Mg abundance ratios which have been reported in some other systems. This is inconsistent with metal-production calculations in SNII and may indicate an additional source of enrichment, such as Population III hypernovae.(abridged)

The Chandra and XMM data for a sample of 19 relaxed groups/poor clusters, covering the temperature range 1-3 keV and selected to have the best available data, reveal a remarkable similarity in their temperature profile: cool cores outside of which the temperatures reach a peak for radii less than 0.1 of the virial radius and then decline. We fitted the derived mass profiles using an NFW model, which provides a good fit to the data when accounting for the central galaxy in the inner region. The concentration parameters and virial masses are in the range c = 5-22 and M = 2 x 10^13 - 4 x 10^14 solar masses, in general agreement with the concentrations found in numerical simulations.

Cosmological simulations of galaxy formation predict a universal form for the mass profile of dark matter (DM) halos from cluster to galaxy scales. Remarkably few interesting constraints exist, however, on DM halos in early-type galaxies. Using Chandra we present the temperature, density and mass profiles of a small sample of early-type galaxies, revealing significant DM in each case. When a component is included to account for stellar mass and the DM halo is allowed to respond adiabatically to the baryonic condensation into stars, the mass profiles are well-fitted by the universal profile, with Virial masses and concentrations in agreement with simulations. However, only ~half, or less, of the mass within R_e seems attributable to the stars, implying stellar M_*/L_B ~1-5.

We have analyzed an XMM-Newton observation of the cluster Abell 2589. Apart from a low-level asymmetry in the central region, the cluster appears very relaxed and does not show presence of a central AGN. We derived constraints for the radial temperature, density and, assuming hydrostatic equilibrium, mass profiles. We find that the best fit to the dark matter profile is given by the Sersic-like profile proposed by Navarro et al. (2004). The NFW model does not provide a good fit. We also tested whether the central stellar component could affect the profile through the adiabatic contraction model but were unable to distinguish it from a simple "dark matter + stars" modeling.

(Abridged) We study predictions for galaxy cluster observables that can test the statistics of dark matter halo shapes expected in a flat LCDM universe. We present a simple analytical model for the prediction of cluster-scale X-ray observations, approximating clusters as isothermal systems in hydrostatic equilibrium, and dark matter haloes as ellipsoids with uniform axial ratios. We test the model against high-resolution, hydrodynamic cluster simulations to gauge its reliability. We find that this simple prescription does a good job of predicting the distribution of cluster X-ray ellipticities compared to the simulations as long as one focuses on cluster regions that are less sensitive to recent mergers. Based on this simple model, the distribution of cluster-size halo shapes expected in the concordance LCDM cosmology implies an X-ray ellipticity distribution with a mean of 0.32 +- 0.01 and a scatter of 0.14 +- 0.01 for the mass range (1-4)x10^{14} Msun/h. We find it important to include the mass dependence of halo shape to make comparisons to observational samples that contain many, very massive clusters. We analyse the systematics of four observational samples of cluster ellipticities and find that our results are statistically compatible with observations. In particular, we find remarkably good agreement between two recent ROSAT samples and LCDM predictions that DO NOT include gas cooling. We also test how well our analytical model can predict Sunyaev-Zel'dovich decrement maps and find that it is less successful although still useful; the model does not perform as well as a function of flux level in this case because of the changing triaxiality of dark matter haloes as a function of radial distance. Both this effect and the changing alignment of isodensity shells of dark matter haloes leave an imprint on cluster gas...

We present an analysis of a Chandra ACIS-I observation of the massive galaxy cluster A644. This cluster was previously classified as a cooling flow, but no radio emission has been detected from its cD galaxy. Outside the core (R ~75 kpc ~0.03R_vir) the hot ICM has properties consistent with a (relaxed) cool-core cluster out to the largest radii investigated (R ~415 kpc ~0.14 R_vir). Over this region the gravitating mass profile is described well by a Navarro-Frenk-White profile with concentration parameter, c = 6.1 +/- 1.2, and virial radius, R_vir = 2.9 +/- 0.4 Mpc. However, inside the core the temperature and entropy profiles reverse their inward radial decline and rise at the center; the inner temperature profile is inconsistent with a constant at the 2.3 sigma level. Although the core region does not display X-ray cavities or filamentary structures characteristic of radio-loud, cool-core clusters, the peak of the X-ray emission is offset from that of the centroid of the global X-ray halo by ~60 kpc. The position of the cD galaxy lies approximately between the X-ray peak and centroid, further testifying to a merger origin for the properties of the X-ray emission in the core. We discuss the implications of A644 and the small number of radio-quiet, cool-core clusters for the AGN feedback paradigm to suppress cooling flows in clusters.

Using Chandra archival data, we quantify the evolution of cluster morphology with redshift. Clusters form and grow through mergers with other clusters and groups, and the amount of substructure in clusters in the present epoch and how quickly it evolves with redshift depend on the underlying cosmology. Our sample includes 40 X-ray selected, luminous clusters from the Chandra archive, and we quantify cluster morphology using the power ratio method (Buote & Tsai 1995). The power ratios are constructed from the moments of the X-ray surface brightness and are related to a cluster's dynamical state. We find that, as expected qualitatively from hierarchical models of structure formation, high-redshift clusters have more substructure and are dynamically more active than low-redshift clusters. Specifically, the clusters with z>0.5 have significantly higher average third and fourth order power ratios than the lower redshift clusters. Of the power ratios, $P_3/P_0$ is the most unambiguous indicator of an asymmetric cluster structure, and the difference in $P_3/P_0$ between the two samples remains significant even when the effects of noise and other systematics are considered. After correcting for noise, we apply a linear fit to $P_3/P_0$ versus redshift and find that the slope is greater than zero at better than 99% confidence. This observation of structure evolution indicates that dynamical state may be an important systematic effect in cluster studies seeking to constrain cosmology, and when calibrated against numerical simulations, structure evolution will itself provide interesting bounds on cosmological models.

We present spectral analysis of the diffuse emission in the normal, isolated, moderate-Lx S0 NGC 1332, constraining both the temperature profile and the metal abundances in the ISM. The characteristics of the point source population and the gravitating matter are discussed in two companion papers. The diffuse emission comprises hot gas, with an ~isothermal temperature profile (~0.5 keV), and emission from unresolved point-sources. In contrast with the cool cores of many groups and clusters, we find a small central temperature peak. We obtain emission-weighted abundance contraints within 20 kpc for several key elements: Fe, O, Ne, Mg and Si. The measured iron abundance (Z_Fe=1.1 in solar units; >0.53 at 99% confidence) strongly excludes the very sub-solar values often historically reported for early-type galaxies but agrees with recent observations of brighter galaxies and groups. The abundance ratios, with respect to Fe, of the other elements were also found to be ~solar, although Z_o/Z_Fe was significantly lower (<0.4). Such a low O abundance is not predicted by simple models of ISM enrichment by Type Ia and Type II supernovae, and may indicate a significant contribution from primordial hypernovae. Revisiting Chandra observations of the moderate-Lx, isolated elliptical NGC 720, we obtain similar abundance constraints. Adopting standard SNIa and SNII metal yields, our abundance ratio constraints imply 73+/-5% and 85+/-6% of the Fe enrichment in NGC 1332 and NGC 720, respectively, arises from SNIa. Although these results are sensitive to the considerable systematic uncertainty in the SNe yields, they are in good agreement with observations of more massive systems. These two moderate-Lx early-type galaxies reveal a consistent pattern of metal enrichment from cluster scales to moderate Lx/Lb galaxies. (abridged)

The galaxy group NGC 5044 consists of a luminous giant elliptical galaxy surrounded by a cluster of approximately 160 low luminosity and dwarf galaxies, mostly of early type. The cumulative projected radial distribution of dwarf galaxies in the NGC 5044 group, unlike distributions of more luminous galaxies in rich clusters, does not follow a projected dark matter (NFW) profile. Instead, a deficiency or absence of low luminosity galaxies is apparent in NGC 5044 within about 350 kpc, about one third of the current group virial radius. Most of the dwarf galaxies identified in NGC 5044 entered the virial radius at redshifts less than about 2 or 3, and very few entered at higher redshifts. The peculiar, non-NFW shape of the projected cumulative dwarf galaxy distribution in NGC 5044 within 350 kpc resembles the characteristic shallow cumulative distribution of dark subhalos that are also known to be relatively young. Dynamical friction is unlikely to explain the apparent lack of group member galaxies at small radii in NGC 5044.

We describe a new type of dynamical model for hot gas in galaxy groups and clusters in which gas moves simultaneously in both radial directions. Circulation flows are consistent with (1) the failure to observe cooling gas in X-ray spectra, (2) multiphase gas observed near the centers of these flows and (3) the accumulation of iron in the hot gas from Type Ia supernovae in the central galaxy. Dense inflowing gas cools, producing a positive central temperature gradient, as in normal cooling flows. Bubbles of hot, buoyant gas flow outward. Circulation flows eventually cool catastrophically if the outward flowing gas transports mass but no heat; to maintain the circulation both mass and energy must be supplied to the inflowing gas over a large volume, extending to the cooling radius. The rapid radial recirculation of gas produces a flat central core in the gas iron abundance, similar to many observations. We believe the circulation flows described here are the first gasdynamic, long-term evolutionary models that are in good agreement with all essential features observed in the hot gas: little or no gas cools as required by XMM spectra, the gas temperature increases outward near the center, and the gaseous iron abundance is about solar near the center and decreases outward.

Chandra ACIS-S3 observations of the S0 galaxy NGC 1332 resolve much of the emission into 73 point-sources, of which 37 lie within the D_25 isophote. The unresolved emission is discussed in two companion papers. The source luminosity function (XLF) shows the break seen in other early-type galaxies at \~2x10^{38}ergs/s. After correcting for detection incompleteness due to source confusion and contamination from diffuse emission, the break vanishes and the data are well-described as a single power law. This casts further doubt on there being a ``universal'' XLF break in such galaxies dividing neutron-star and black-hole systems. The slope of the differential XLF (dN/dL), beta=2.7+/-0.5, is marginally (~2.5 sigma) steeper than has been found for analogous fits of other early-type galaxies but resembles what is generally seen at high luminosities. Two of the sources within D_25 are ULX, although neither have L_X>2x 10^{39} ergs/s. The absence of very luminous ULX in early-type galaxies suggests a break in the XLF slope at ~1-2x 10^{39} ergs/s, but we cannot constrain it in NGC 1332. The sources have a spatial distribution consistent with the optical light and display a range of properties expected for a LMXB population. The spectra of the individual sources, as well as the composite source spectra, agree with observations of other early-type galaxies, although a small number of highly-absorbed sources are seen. Two sources have very soft spectra, two show strong variability and one source shows evidence of an extended radial profile. We do not detect a central source in NGC 1332, but we find a faint (L_X=2+/-1 x 10^{38} ergs/s) point-source coincident with the centre of the companion dwarf galaxy NGC 1331. (Abridged)

A new XMM observation of the outer regions of the galaxy group NGC 5044 indicates hot gas iron abundances of only ~0.15 solar between r=0.2 - 0.4 r_vir. While the total baryon mass within the virial radius may be close to the cosmic mean value observed in rich clusters, the ratio of total iron mass to optical light in NGC 5044 is about 3 times lower than that in rich clusters. The remarkably low iron abundance over a large volume of the intergroup gas in the outer regions of NGC 5044 cannot be easily understood in terms of the outflow of enriched gas in a group wind during its early history or by the long term enrichment by the group member galaxies that currently occupy this region. It is possible that the stars in NGC 5044 did not produce iron with the same efficiency as in clusters, or the iron resides in non-luminous clouds or stars, or the entropy of the iron-enriched gas created in early galactic starburst winds was too high to penetrate the group gas of lower entropy.

We present an analysis of a Chandra--ACIS observation of the galaxy cluster A2589 to constrain the radial distribution of the total gravitating matter and the dark matter in the core of the cluster. A2589 is especially well-suited for this analysis because the hot gas in its core region (r < ~0.1 Rvir) is undisturbed by interactions with a central radio source. From the largest radius probed (r=0.07 Rvir) down to r ~0.02 Rvir dark matter dominates the gravitating mass. Over this region the radial profiles of the gravitating and dark matter are fitted well by the NFW and Hernquist profiles predicted by CDM. The density profiles are also described well by power laws, rho ~r^{-alpha}, where alpha=1.37 +/- 0.14 for the gravitating matter and alpha=1.35 +/- 0.21 for the dark matter. These values are consistent with profiles of CDM halos but are significantly larger than alpha ~0.5 found in LSB galaxies and expected from self-interacting dark matter models.

I review constraints on the radial density profiles and ellipticities of the dark matter obtained from recent X-ray observations with Chandra and XMM of elliptical galaxies and galaxy clusters and discuss their implications, especially for the self-interacting dark matter model.

Using Chandra archival data, we quantify the evolution of cluster morphology with redshift. To quantify cluster morphology, we use the power ratio method developed by Buote and Tsai (1995). Power ratios are constructed from moments of the two-dimensional gravitational potential and are, therefore, related to a cluster's dynamical state. Our sample will include 40 clusters from the Chandra archive with redshifts between 0.11 and 0.89. These clusters were selected from two fairly complete flux-limited X-ray surveys (the ROSAT Bright Cluster Sample and the Einstein Medium Sensitivity Survey), and additional high-redshift clusters were selected from recent ROSAT flux-limited surveys. Here we present preliminary results from the first 28 clusters in this sample. Of these, 16 have redshifts below 0.5, and 12 have redshifts above 0.5.

The failure of the XMM-Newton and Chandra X-ray telescopes to detect cooling gas in elliptical galaxies and clusters of galaxies has led many to adopt the position that the gas is not cooling at all and that heating by an active nucleus in the central E or cD galaxy is sufficient to offset radiative cooling. In this paper we explore an idealized limiting example of this point of view in which hot, buoyant bubbles formed near the center return the inflowing, radiatively cooling gas to distant regions in the flow. We show that idealized steady state, centrally heated non-cooling flows can indeed be constructed. In addition, the emission-weighted temperature profiles in these circulating flows resemble those of normal cooling flows. However, these solutions are valid only (1) for a range of bubble parameters for which there is no independent justification, (2) for a limited spatial region in the cooling flow and (3) for a limited period of time after which cooling seems inevitable. Our exploration of non-cooling flows is set in the context of galaxy/group flows.

Using new XMM and Chandra observations we present an analysis of the metal abundances of the hot gas within a radius of 100 kpc of the bright nearby galaxy group NGC 5044. Motivated by the inconsistent abundance and temperature determinations obtained by different observers for X-ray groups, we provide a detailed investigation of the systematic errors on the derived abundances considering the effects of the temperature distribution, calibration, plasma codes, bandwidth, Galactic Nh, and background rate. The iron abundance (Fe) drops from Fe ~1 solar within R ~50 kpc to Fe ~0.4 solar near R=100 kpc. This radial decline in Fe is highly significant: Fe=1.09 +/- 0.04 solar (stat) +/- 0.05 solar + 0.18 solar (sys) within R=48 kpc (5') compared to Fe=0.44 +/- 0.02 solar (stat) +/- 0.10 solar + 0.13 solar (sys) over R=48-96 kpc (5'-10'). The data rule out with high confidence a very sub-solar value for Fe within R=48 kpc confirming that previous claims of very sub-solar central Fe values in NGC 5044 were primarily the result of the Fe Bias: i.e., the incorrect assumption of spatially isothermal and single-phase gas when in fact temperature variations exist. Within R=48 kpc we obtain Si/Fe = 0.83 +/- 0.02 (stat) +/- 0.02 + 0.07 (sys) and S/Fe = 0.54 +/- 0.02 (stat) +/- 0.01 + 0.01 (sys) in solar units. These ratios are consistent with their values at larger radii and imply that SNIa have contributed ~80% of the iron mass within a 100 kpc radius of NGC 5044. This SNIa fraction is similar to the Sun and suggests an IMF similar to that of the Milky Way. At the very center (R ~2 kpc) the XMM and Chandra CCDs and the XMM RGS show that the Fe drops to ~50% of its value at immediately larger radius analogously to that seen in some galaxy clusters. (Abridged)

Using new XMM and Chandra observations we present an analysis of the temperature structure of the hot gas within a radius of 100 kpc of the bright nearby galaxy group NGC 5044. A spectral deprojection analysis of data extracted from circular annuli reveals that a two-temperature model (2T) of the hot gas is favored over single-phase or cooling flow (M_dot = 4.5 +/- 0.2 M_{sun}/yr) models within the central ~30 kpc. Alternatively, the data can be fit equally well if the temperature within each spherical shell varies continuously from ~T_h to T_c ~ T_h/2, but no lower. The high spatial resolution of the Chandra data allows us to determine that the temperature excursion T_h --> T_c required in each shell exceeds the temperature range between the boundaries of the same shell in the best-fitting single-phase model. This is strong evidence for a multi-phase gas having a limited temperature range. The cooler component of the 2T model has a temperature (T_c \~ 0.7 keV) similar to the kinetic temperature of the stars. The hot phase has a temperature (T_h ~ 1.4 keV) characteristic of the virial temperature of the \~10^{13} M_{sun} halo expected in the NGC 5044 group. However, in view of the morphological disturbances and X-ray holes visible in the Chandra image within R ~10 kpc, bubbles of gas heated to ~T_h in this region may be formed by intermittent AGN feedback. Some additional heating at larger radii may be associated with the evolution of the cold front near R ~50 kpc, as suggested by the sharp edge in the EPIC images.

Using Chandra archival data, we have begun a study to quantify the evolution of cluster morphology with redshift. To quantify cluster morphology, we use the power ratio method developed by Buote and Tsai (1995). Power ratios are constructed from moments of the two-dimensional gravitational potential and are, therefore, related to a cluster's dynamical state. Our sample will include around 50 clusters from the Chandra archive with redshifts between 0.11 and 1.1. These clusters were selected from two fairly complete flux-limited X-ray surveys (the ROSAT Bright Cluster Sample and the Einstein Medium Sensitivity Survey), and additional high-redshift clusters were selected from recent ROSAT flux-limited surveys. Here we present preliminary results from the first 15 clusters in this sample. Of these, eight have redshifts below 0.5, and seven have redshifts above 0.5.

We have used a high resolution Chandra observation to examine the core mass distribution of the unusually regular cD cluster Abell 2029. This system is especially well-suited for analysis of its mass distribution under the assumption of hydrostatic equilibrium: it exhibits an undisturbed X-ray morphology and a single-phase ICM. We estimate the total mass, and the contributions of the gas and dark matter components from <3 arcsec to ~3 arcmin (<4.4-260 kpc, 0.001-0.1 r_vir). The gas density is well fitted by either a double beta-model, or a ``cusped'' beta-model. The temperature data are increasing, and can be approximated by a power-law, with alpha_T = 0.27+/-0.01. Fitting smooth functions to density and T, we obtain a mass profile, which is remarkably well described down to 0.002 r_vir by the NFW model, for which we measure r_scale = 540+/-90 kpc (~0.2 r_vir) and c = 4.4 +/- 0.9. The mass profile is also approximated by a power-law fit with alpha_m = 1.81+/-0.04 (corresponding to a logarithmic density slope of -1.19+/-0.04). The density profile is too shallow to be fitted with a Moore et al. model. The consistency with NFW at <0.01 r_vir is incompatible with SIDM, and contrasts with the results from rotation curves of LSB and dwarf galaxies. This suggests that while CDM may adequately describe clusters, it does not currently account properly for small halos. Using the cD light, we observe a total M/L_V ~ 12 at r < 20 kpc, rising to > 100 beyond 200 kpc. The consistency with an NFW halo model and the large M/L suggest the cluster is DM-dominated down to <~ 0.005 r_vir. We observe a gas mass fraction = 0.139+/-0.004 at the limit of our observations, finding an upper limit to Omega_m <= 0.29+/-0.03 (abridged).

With a Chandra ACIS-S3 observation, we detect 42 X-ray point sources in the elliptical galaxy NGC 720, including a possible central source. Most of these sources will be low-mass X-ray binaries (LMXBs), and 12 are located within 2" of globular cluster candidates. We investigate both the hardness ratios and combined spectra of the sources. They exhibit a distribution of X-ray colors similar to those seen in other early-type galaxies. We find that there is a population of highly absorbed sources located at large distances from the center of the galaxy. The overall spatial distribution of sources is consistent with the ellipticity and position angle of the galaxy, but the sources appear to form several arcs. NGC 720 contains nine ultraluminous sources (L_x >= 10^39 ergs/s). This number is more than have previously been detected in an early-type galaxy but similar to the number seen in the Antennae merger system. The ratio L_ULX/L_B for NGC 720 is more than double the ratio for the S0 galaxy NGC 1553 and a factor of seven higher than for the elliptical galaxy NGC 4697, although uncertainties in the distance and the source spectral properties could bring these ratios into agreement. The X-ray source luminosity function is also nearly as flat as those seen in disk and merger systems. The large number of ULXs in NGC 720 and its relatively flat XLF may provide evidence against the association of all ULXs with young stars. We also see a possible high luminosity break in the luminosity function at 2x10^39 h_50^-2 ergs/s.

X-ray observations with Chandra and XMM are providing valuable new measurements of the baryonic and dark matter content of groups and clusters. Masses of cD clusters obtained from X-ray and gravitational lensing studies generally show good agreement, therefore providing important validation of both methods. Gas fractions have been obtained for several clusters that verify previous results for a low matter density (Omega_m ~0.3). Chandra has also provided measurements of the mass profiles deep down into several cluster cores and has generally found no significant deviations from CDM predictions in contrast to the flat core density profiles inferred from the rotation curves of low-surface brightness galaxies and dwarf galaxies; i.e., there is no evidence for self-interacting dark matter in cluster cores. Finally, initial studies of the iron and silicon abundances in centrally E-dominated groups show that they have pronounced gradients from 1-2 solar values within the central 30-50 kpc that fall to values of 0.3-0.5 solar at larger radii. The Si/Fe ratios are consistent with approximately 80% of the metals originating from Type Ia supernovae. Several cD clusters also display central Fe enhancements suggestive of Type Ia supernova enrichment, though some have central dips that may provide a vital clue for solving the cooling flow mystery.

We present an initial analysis of a new XMM observation of NGC 1399, the central elliptical galaxy of the Fornax group. Spectral fitting of the spatially resolved spectral data of the EPIC MOS and pn CCDs reveals that a two-temperature model (2T) of the hot gas is favored over single-phase and cooling flow models within the central ~20 kpc. The preference for the 2T model applies whether or not the data are deprojected. The cooler component has a temperature (~0.9 keV) similar to the kinetic temperature of the stars while the hotter component has a temperature (~1.5 keV) characteristic of the virial temperature of a ~10^{13} M_sun halo. The two-phase model (and other multitemperature models) removes the ``Fe Bias'' within r < ~20 kpc and gives Z_Fe/Z_sun 1.5-2. At larger radii the iron abundance decreases until Z_Fe/Z_sun \~0.5 for r ~50 kpc. The Si abundance is super-solar (1.2-1.7 solar) within the central regions while Z_Si/Z_Fe ~0.8 over the entire region studied. These Fe and Si abundances imply that ~80% of the Fe mass within r ~50 kpc originates from Type Ia supernovae (SNIa). This SNIa fraction is similar to that inferred for the Sun and therefore suggests a stellar initial mass function similar to the Milky Way.

(Abridged) We present an analysis of a Chandra ACIS-S observation of the elliptical galaxy NGC 720 to verify the existence of a dark matter (DM) halo and to measure its ellipticity. The ACIS-S3 image reveals over 60 point sources. For semi-major axes a<~150" (18.2h_{70}^{-1} kpc) the ellipticity of the diffuse emission is ex ~0.15, which is less than the values 0.2-0.3 obtained from ROSAT because the point sources contaminated the ROSAT values. The Chandra data confirm the ~20 deg position angle (PA) twist discovered by ROSAT, but the Chandra twist is more gradual also because of the point sources contaminating the ROSAT values. Overall the ex and PA values for a<~150" can be explained by the triaxial model of NGC 720 published by Romanowsky & Kochanek. Since the optical image displays no substantial isophote twisting, the X-ray PA twist requires a massive DM halo if the hot gas is in hydrostatic equilibrium. The mass-follows-light hypothesis is also inconsistent with the Chandra ellipticities at the 96% (98%) level for oblate (prolate) symmetry. Thus, both the PA twist and the ellipticities of the Chandra image imply a DM halo independent of the gas T profile -- evidence that cannot be obviated by alternative gravity theories such as MOND. The DM density model, rho ~a^{-2}, provides the best fit and gives ellipticities of 0.37 +/- 0.03 (0.36 +/- 0.02) for oblate (prolate) models. These moderate ellipticities for the DM halo are inconsistent with both the nearly spherical halos predicted if the DM is self-interacting and with the highly flattened halos predicted if the DM is cold molecular gas. These ellipticities may also be too large to be explained by warm DM, but are consistent with galaxy-sized halos formed in the Lambda-CDM paradigm.

We have obtained high spatial resolution temperature and abundance profiles for the galaxy cluster Abell 2029 with the Chandra ACIS-S instrument. Our observations reveal that the spectra are well-fit by a single-phase gas in each annulus. While the temperature of the intracluster medium drops from ~9 keV at 3 arcmin (260 kpc) to 3 keV in the central 5 arcsec (5 kpc) of the cluster, there is no evidence for gas emitting at temperatures below 3 keV. The addition of a cooling flow component does not improve the fits, despite previous claims for a massive cooling flow. There is also no evidence for excess absorption above the Galactic N_H value. We also observe a steep Fe abundance gradient, such that Z_Fe >= 2 Z_solar (assuming meteoritic solar Fe abundance) in the core, consistent with significant enrichment from Fe-rich Type Ia supernovae in the cD galaxy. The Fe abundance drops to ~0.5 Z_solar at 3 arcmin (260 kpc), in good agreement with previous BeppoSAX measurements. The high resolution image reveals neither a strong central point source, nor any filamentary structure related to a cooling flow or a merger. The absence of a strong merger signature argues against the creation of the wide-angle-tail radio source morphology in a merger event. (modified for astro-ph).

X-ray observations have played a key role in the study of substructure and merging in galaxy clusters. I review the evidence for cluster substructure and mergers obtained from X-ray observations with satellites that operated before Chandra and XMM. Different techniques to study cluster mergers via X-ray imaging and spectral data are discussed with an emphasis on the quantitative analysis of cluster morphologies. I discuss the implications of measurements of cluster morphologies for cosmology and the origin of radio halos.

Previously it has been recognized that radio halos in galaxy clusters are preferentially associated with merging systems as indicated by substructure in the X-ray images and temperature maps. Since, however, many clusters without radio halos also possess substructure, the role of mergers in the formation of radio halos has remained unclear. By using power ratios to relate gravitational potential fluctuations to substructure in X-ray images, we provide the first quantitative comparison of the dynamical states of clusters possessing radio halos. A correlation between the 1.4 GHz power (P_{1.4}) of the radio halo (or relic) and the magnitude of the dipole power ratio (P_1/P_0) is discovered such that approximately P_{1.4} ~ P_1/P_0; i.e., the strongest radio halos appear only in those clusters currently experiencing the largest departures from a virialized state. From additional consideration of a small number of highly disturbed clusters without radio halos detected at 1.4 GHz, and recalling that radio halos are more common in clusters with high X-ray luminosity (Giovannini, Tordi, & Feretti), we argue that radio halos form preferentially in massive (L_x >~ 0.5 x 10^{45} erg/s) clusters experiencing violent mergers (P_1/P_0 >~ 0.5 x 10^{-4}) that have seriously disrupted the cluster core. The association of radio halos with massive, large-P_1/P_0, core-disrupted clusters is able to account for both the vital role of mergers in accelerating the relativistic particles responsible for the radio emission as well as the rare occurrence of radio halos in cluster samples.

Recent XMM-Newton observations of the cooling flow gas in M87 indicate sharply decreasing oxygen, iron and silicon abundances within $\sim 5$ kpc of the galactic center. This result is unexpected since stellar mass loss and Type Ia supernovae are expected to produce pronounced central abundance maxima for all three elements. However, it has been suggested that many of the strong X-ray lines are optically thick and diffuse to larger radii in the cooling flow before escaping, falsifying the central abundances. We verify with radiation transfer calculations that this effect does indeed occur in the M87 cooling flow, but that it is insufficient to account for the M87 observations. We suggest that some source of continuous opacity is required to reduce the central X-ray line emission, perhaps by warm gas at $T \sim 10^5 - 10^6$ K. The radial surface brightness profiles of X-ray resonance lines are also sensitive to turbulence in cooling flows which reduces the line center optical depths considerably. Turbulence may provide sufficient energy to continuously heat the warm absorbing gas.

We summarize the discovery of oxygen absorption and warm (10^5-10^6 K) gas in cooling flows. Special attention is given to new results for M87 for which we find the strongest evidence to date for ionized oxygen absorption in these systems. We briefly discuss implications for observations of cooling flows with Chandra and XMM.

We present a re-analysis of the ROSAT PSPC data within the central 100 kpc of M87 to search for intrinsic oxygen absorption similar to that recently measured in several galaxies and groups. Using a spatial-spectral deprojection analysis we find the strongest evidence to date for intrinsic oxygen absorption in the hot gas of a galaxy, group, or cluster. Single-phase plasma models modified by intervening Galactic absorption cannot fit the 0.2-2.2 keV PSPC data as they under-predict the 0.2-0.4 keV region and over-predict the 0.5-0.8 keV region where the emission and absorption residuals are obvious upon visual inspection of the spectral fits. Since the excess emission between 0.2-0.4 keV rules out intrinsic absorption from cold gas or dust, the most reasonable model for the excess emission and absorption features is warm, collisionally ionized gas with a temperature of ~10^6 K. Simple multiphase models (cooling flow, two phases) modified by both intervening Galactic absorption and by a single oxygen edge provide good fits and yield temperatures and Fe abundances of the hot gas that agree with previous determinations by ASCA and SAX. The multiphase models of M87 inferred from the PSPC can account for the excess EUV emission observed with EUVE and the excess X-ray absorption inferred from EINSTEIN and ASCA data above 0.5 keV. Although the total mass of the warm gas implied by the oxygen absorption is consistent with the matter deposited by a cooling flow, the suppression of the mass deposition rate and the distortion of the X-ray isophotes in the region where the radio emission is most pronounced suggest some feedback effect from the AGN on the cooling gas. (Abridged)

Using spatially resolved, deprojected ROSAT PSPC spectra of 10 of the brightest cooling flow galaxies and groups with low Galactic column densities we have detected intrinsic absorption over energies ~0.4-0.8 keV in half of the sample. Since no intrinsic absorption is indicated for energies below ~0.4 keV, the most reasonable model for the absorber is collisionally ionized gas at temperatures T=10^{5-6} K with most of the absorption arising from ionized states of oxygen but with a significant contribution from carbon and nitrogen. The soft X-ray emission of this warm gas can explain the sub-Galactic column densities of cold gas inferred within the central regions of most of the systems. Attributing the absorption to ionized gas reconciles the large columns of cold H and He inferred from EINSTEIN and ASCA with the lack of such columns inferred from ROSAT. Within the central ~10-20 kpc, where the constraints are most secure, the estimated mass of the ionized absorber is consistent with most (perhaps all) of the matter deposited by a cooling flow over the lifetime of the flow. Since the warm absorber produces no significant H or He absorption the large absorber masses are consistent with the negligible atomic and molecular H inferred from HI and CO observations of cooling flows. It is also found that if T > ~2x10^5 K then the optical and UV emission implied by the warm gas does not violate published constraints. Finally, we discuss how the prediction of warm ionized gas as the product of mass drop-out in these and other cooling flows can be verified with new CHANDRA and XMM observations. (Abridged)

Previous studies of the Fe abundances in the hot gas of galaxies and groups have reported conflicting results with most studies finding very sub-solar Fe abundances that disagree with standard theory. To investigate the possible role of Fe abundance gradients on these measurements we present deprojection analysis of the ROSAT PSPC data of 10 of the brightest cooling flow galaxies and groups. The PSPC allows for spatially resolved spectral analysis on a half-arcminute scale, and interesting constraints on both the temperatures and Fe abundances are possible because the ~1 keV temperatures of these systems are well matched to the bandpass of the PSPC. In 9 out of 10 systems we find clear evidence that the Fe abundance decreases with increasing radius: Z_Fe ~1 solar - (several) solar within the central radial bin (r < ~10 kpc) which decreases to Z_Fe ~0.5 solar at the largest radii examined (r ~50-100 kpc). The Fe abundances (and temperatures) are consistent with the average values for these systems that we obtained in our previous analyses of the ASCA data using multitemperature models which confirms that previous inferences of very sub-solar Fe abundances from ASCA arise from the incorrect assumption of isothermal gas and not the presence of Fe abundance gradients. We discuss why this ``Fe Bias'' affects much more seriously the measurements of Z_Fe from ASCA data than from ROSAT data. We show that the Fe abundance profiles for these galaxies and groups are consistent with a gas-dynamical model where the gas is enriched by stellar ejecta and supernovae in the ``solar supernova proportion'', the stars formed with a Galactic IMF, and the gas is diluted by mixing with primordial gas at large radii.

The inhomogeneous cooling flow scenario predicts the existence of large quantities of gas in massive elliptical galaxies, groups, and clusters that have cooled and dropped out of the flow. Using spatially resolved, deprojected X-ray spectra from the ROSAT PSPC we have detected strong absorption over energies ~0.4-0.8 keV intrinsic to the central ~1 arcmin of the galaxy, NGC 1399, the group, NGC 5044, and the cluster, A1795. These systems have amongst the largest nearby cooling flows in their respective classes and low Galactic columns. Since no excess absorption is indicated for energies below ~0.4 keV the most reasonable model for the absorber is warm, collisionally ionized gas with T=10^{5-6} K where ionized states of oxygen provide most of the absorption. Attributing the absorption only to ionized gas reconciles the large columns of cold H and He inferred from Einstein and ASCA with the lack of such columns inferred from ROSAT, and also is consistent with the negligible atomic and molecular H inferred from HI, and CO observations of cooling flows. The prediction of warm ionized gas as the product of mass drop-out in these and other cooling flows can be verified by Chandra, XMM, and ASTRO-E.

In this paper we describe a self-contained method for performing the spectral-imaging deconvolution of X-ray data on clusters of galaxies observed by the ASCA satellite. Spatially-resolved spectral studies of data from this satellite require such a correction because its optics redistribute photons over regions which are of comparable size to the angular scales of interest in clusters. This scattering is a function not only of spatial position but also energy. To perform a correction for these effects we employ Maximum-Likelihood deconvolution of the image (within energy bands of 1 keV) to determine the spatial redistribution, followed by a Monte-Carlo energy reassignment of photon energies with position to determine the spectral redistribution. We present tests on simulated cluster data, convolved with the various instrumental characteristics and the X-ray background, which show that our methodology can successfully recover a variety of intrinsic temperature profiles in typical observational circumstances. In Paper-II we apply our spectral-imaging deconvolution procedure to a large sample of galaxy clusters to determine temperature profiles, some of which will be used in subsequent mass determinations, presented in Paper-III.

We analyze the ASCA spectra accumulated within ~100 kpc radii of 12 of the brightest groups of galaxies. Upon fitting isothermal models (1T) jointly to the ASCA SIS and GIS spectra we obtain fits for most groups that are of poor or at best marginal quality and give very sub-solar metallicities similar to previous studies, <Z> = 0.29 +/- 0.12 Z_sun. Two-temperature models (2T) provide significantly better fits for 11 out of the 12 groups and in every case have metallicities that are substantially larger than obtained for the 1T models, <Z> = 0.75 +/- 0.24 Z_sun. Although not very well constrained, for most of the groups absorption in excess of the Galactic value is indicated for the cooler temperature component of the 2T models. A simple multiphase cooling flow model gives results analogous to the 2T models including large metallicities, <Z> = 0.65 +/- 0.17 Z_sun. The nearly solar Fe abundances and also solar alpha/Fe ratios indicated by the 2T and cooling flow models are consistent with models of the chemical enrichment of ellipticals, groups, and clusters which assume ratios of Type Ia to Type II supernova and an IMF similar to those of the Milky Way. Thus, we have shown that the very sub-solar Fe abundances and Si/Fe enhancements obtained from most previous studies within r ~100 kpc of galaxy groups are an artifact of their fitting isothermal models to the X-ray spectra which also has been recently demonstrated for the brightest elliptical galaxies. Owing to the importance of these results for interpreting X-ray spectra, in an appendix we use simulated ASCA observations to examine in detail the ``Fe bias'' and ``Si bias'' associated with the spectral fitting of ellipticals, groups, and clusters of galaxies.

We examine the K shell emission lines produced by isothermal and simple multiphase models of the hot gas in elliptical galaxies and galaxy clusters to determine the most effective means for constraining the width of the differential emission measure (xi(T)) in these systems which we characterize by a dimensionless parameter, sigma_xi. Comparison of line ratios of two-temperature (sigma_xi << 1) and cooling flow (sigma_xi ~1) models is presented in detail. We find that a two-temperature model can approximate very accurately a cooling flow spectrum over 0.5-10 keV. We have re-analyzed the ASCA spectra of three of the brightest galaxy clusters to assess the evidence for multiphase gas in their cores: M87 (Virgo), the Centaurus cluster, and the Perseus cluster. K-alpha emission line blends of Si, S, Ar, Ca, and Fe are detected in each system as is significant Fe K-beta emission. The Fe K-beta/K-alpha ratios are consistent with optically thin plasma models and do not suggest resonance scattering in these systems. Consideration of both the ratios of H-like to He-like K-alpha lines and the local continuum temperatures clearly rules out isothermal gas in each case. To obtain more detailed constraints we fitted plasma models over 1.6-9 keV where the emission is dominated by these K shell lines and by continuum. In each case the ASCA spectra cannot determine whether the gas emits at only two temperatures or over a continuous range of temperatures as expected for a cooling flow. The metal abundances are near solar for all of the multiphase models. We discuss the implications of these results and examine the prospects for determining the temperature structure in these systems with upcoming X-ray missions.

We examine whether isothermal models of the hot gas can successfully describe the ASCA and ROSAT spectra of NGC 1399, NGC 4472, NGC 4636, and NGC 5044. Broad-band spectral fitting of the ASCA SIS and GIS data accumulated within a radius of ~5 arcmin for each galaxy shows that isothermal models are unable to fit the SIS data near 1 keV. In addition, these isothermal models typically fail to produce the large equivalent widths of the K-alpha line blends of the H-like and He-like ions of Si and S which are measured independently of the Fe L emission lines. Two-temperature models provide excellent broad-band fits to both the SIS and GIS data of each galaxy with the relative abundances (except for NGC 4636) fixed at their solar values. A simple multiphase cooling flow model fits nearly as well as the two-temperature model for NGC 1399, NGC 4472, and NGC 5044. The multiphase models also predict more accurately the Si and S equivalent widths and the ratios of Si XIV/XIII and S XVI/XV than the isothermal models. Using various approaches we find that the temperature gradients inferred from the ROSAT PSPC data of these galaxies, especially for NGC 1399 and NGC 5044, are inconsistent with the isothermal models derived from fitting the ASCA data within a single-aperture but are very consistent with the multiphase models. Therefore, models which assume isothermal gas within r ~5 arcmin are inconsistent with both the ASCA and ROSAT PSPC data of these elliptical galaxies. Simple two-temperature models and multiphase cooling flows provide good descriptions of these data sets with Fe abundances of ~1-2 solar and (except for NGC 4636) relative abundances fixed at their solar values. We discuss the implications of these nearly solar abundances. (Abridged)

We present analysis of the shape and radial mass distribution of the E4 galaxy NGC 3923 using archival X-ray data from the ROSAT PSPC and HRI. The X-ray isophotes are significantly elongated with ellipticity e_x=0.15 (0.09-0.21) (90% confidence) for semi-major axis a\sim 10h^{-1}_70 kpc and have position angles aligned with the optical isophotes within the estimated uncertainties. Applying the Geometric Test for dark matter, which is independent of the gas temperature profile, we find that the ellipticities of the PSPC isophotes exceed those predicted if M propto L at a marginal significance level of 85% (80%) for oblate (prolate) symmetry. Detailed hydrostatic models of an isothermal gas yield ellipticities for the gravitating matter, e_mass=0.35-0.66 (90% confidence), which exceed the intensity weighted ellipticity of the R-band optical light, <e_R> = 0.30 (e_R^max=0.39). We conclude that mass density profiles with rho\sim r^{-2} are favored over steeper profiles if the gas is essentially isothermal (which is suggested by the PSPC spectrum) and the surface brightness in the central regions (r<~15") is not modified substantially by a multi-phase cooling flow, magnetic fields, or discrete sources. We argue that these effects are unlikely to be important for NGC 3923. (The derived e_{mass} range is very insensitive to these issues.) Our spatial analysis also indicates that the allowed contribution to the ROSAT emission from a population of discrete sources with Sigma_x propto Sigma_R is significantly less than that indicated by the hard spectral component measured by ASCA.

We present spectral analysis of ASCA data of 17 elliptical and 3 lenticular galaxies most of which have high L_x/L_B. Single-temperature models (MEKAL and Raymond-Smith) give unacceptable fits (chi^2_red > 1.5) in most cases and, in agreement with previous studies, give very sub-solar abundances, <Z> = 0.19 +/- 0.12 Z_sun (MEKAL). The spectra for approximately half the sample are better fit with a cooling-flow model which in three cases gives a substantially better fit. The abundances derived from the cooling-flow model are also significantly larger, <Z> = 0.6 +/- 0.5 Z_sun. We empirically tested the reliability of the plasma codes in the Fe-L region and found no evidence for serious problems with the determined temperatures and abundances. Two-temperature models give substantially better fits that are formally acceptable (chi^2_red ~ 1.0) in all but a few cases. The highest S/N galaxies (which also have highest L_x/L_B) have fitted temperatures < 2 keV for both components consistent with each being distinct phases of hot gas. The lowest S/N galaxies (which also have lowest L_x/L_B) generally have a hot component with temperature, T_H > ~5 keV, which is consistent with emission from discrete sources. (We discuss the origin of these two components from analysis of L_x-L_B and L_x/L_B.) The abundances of these two-temperature models are approximately solar, <Z> = 0.9 +/- 0.7 Z_sun (MEKAL), consistent with a recent multi-phase model for the evolution of hot gas in ellipticals. Finally, for several galaxies we find evidence for absorption in excess of the Galactic value and discuss its implications.

We review the theory of how the shapes of the X-ray isophotes probe both the shape and radial distribution of gravitating matter in elliptical galaxies in a way that is more robust than the traditional spherical approach. We summarize and update the previous analyses of X-ray observations of NGC 720 and NGC 1332 and describe preliminary results for NGC 3923.

(Abridged) We examine the theoretical relationship between Omega_0 and substructure in galaxy clusters which are formed by the collapse of high density peaks in a gaussian random field. The radial mass distributions of the clusters are computed from the spherical accretion model using the adiabatic approximation following Ryden & Gunn. For a cluster of mass, M(r,t), we compute the quantity dM/M_bar at a cosmic time t and within a radius r, where dM is the accreted mass and M_bar is the average mass of the cluster during the previous relaxation time, which is computed individually for each cluster. For a real cluster in three dimensions we argue that dM/M_bar should be strongly correlated with the low order multipole ratios, Phi^{int}_l/Phi^{int}_0, of the potential due to matter interior to r. It is shown that the expected correlation between dM/M_bar and Phi^{int}_l/Phi^{int}_0 extends to the two-dimensional multipole ratios, Psi^{int}_m/Psi^{int}_0, which are well defined observables of the cluster density distribution. The strongest dependence of dM/M_bar on Omega_0 (lambda_0=0) occurs at z=0 where dM/M_bar propto Omega_0^{1/2} for relaxation times ~1-2 crossing times and only very weakly depends on mass and radius. The fractional accreted mass in CDM models with Omega_0+lambda_0=1 depends very weakly on Omega_0 and has a magnitude similar to the Omega_0=1 value. dM/M_bar evolves more rapidly with redshift in low-density universes and decreases significantly with radius for Omega_0=1 models for z > ~0.5. We discuss how to optimize constraints on Omega_0 and lambda_0 using cluster morphologies.

We summarize a series of papers reporting studies of the ellipticities of X-ray images of five Abell clusters and two early type galaxies to learn about the shape of the gravitational potential and the underlying matter distributions. For the clusters we find the dark matter is centrally concentrated with ellipticity $\epsilon_{DM} \sim 0.5$. For the galaxies we find independent evidence for extended, flattened dark matter halos. One case of a clear isophotal twist in the X-rays may indicate a triaxial halo.

We examine the sensitivity of the spatial morphologies of galaxy clusters to Omega_0 and P(k) using high-resolution N-body simulations with large dynamic range. Variants of the standard CDM model are considered having different spatial curvatures, SCDM (Omega_0=1), OCDM (Omega_0=0.35), LCDM (Omega_0=0.35, lambda_0=0.65), and different normalizations, sigma_8. We also explore critical density models with different spectral indices, n, of the scale-free power spectrum, P(k) propto k^n. Cluster X-ray morphologies are quantified with power ratios (PRs), where we take for the X-ray emissivity j_gas propto rho^2_DM, which we argue is a suitable approximation for analysis of PRs. We find that Omega_0 primarily influences the means of the PR distributions whereas the power spectrum (n and sigma_8) primarily affects their variances: log_10(P_3/P_0) is the cleanest probe of Omega_0 since its mean is very sensitive to Omega_0 but very insensitive to P(k). The PR means easily distinguish the SCDM and OCDM models, while the SCDM and LCDM means show a more modest, but significant, difference (3 sigma). The OCDM and LCDM models are largely indistinguishable in terms of the PRs. Finally, we compared these models to a sample of ROSAT clusters and find that the cluster morphologies strongly disfavor Omega_0=1, CDM while favoring low density, CDM models with or without a cosmological constant.

We present spectral analysis of $ASCA$ observations of the early-type galaxies NGC 720 (E4) and NGC 1332 (E7/S0) with emphasis on constraining the relative contribution to the X-ray emission from hot gas and the integrated emission from X-ray binaries. Single-temperature spectral models yield poor fits to the spectrum ($\chi^2_{red}\sim 3$) over the $\sim 0.5 - 5$ keV energy range. Two-temperature models significantly improve the spectral fits ($\chi^2_{red}\sim 1.5$) and have soft-component temperatures and sub-solar abundances consistent with previous $ROSAT$ single-temperature models ($T_{soft}\sim 0.6$ keV, abundances $\sim 0.1$) and hard-component temperatures ($T_{hard}\gtrsim 3$ keV) consistent with those expected from a discrete component. The soft component dominates the emission in both galaxies, especially in the 0.4 - 2.4 keV band used in previous $ROSAT$ studies. Combining these spectral results with $ROSAT$ data we updated constraints on the mass distributions for NGC 720 and NGC 1332. For NGC 720, which yields the more precise constraints, the ellipticity of the intrinsic shape of the mass is slightly reduced $(\Delta\epsilon_{mass}\approx 0.05)$ when the discrete component is added, $\epsilon_{mass}\sim 0.4-0.6$ $(90%)$. The estimates for the total mass increase with increasing discrete flux, and we find that models with $F_{hard}/F_{soft}=0.45$, the $2\sigma$ upper limit, have masses that exceed by $\sim 30% - 50%$ those where $F_{hard}/F_{soft}=0$.

We present spatial analysis of the deep (57ks) ROSAT HRI X-ray image of the E4 galaxy NGC 720. The orientation of the HRI surface brightness is consistent with the optical position angle $(PA)$ interior to semi-major axis $a\sim 60\arcsec$ (optical $R_e\sim 50\arcsec$). For larger $a$ the isophotes twist and eventually $(a\gtrsim 100\arcsec)$ orient along a direction consistent with the $PA$ measured with the PSPC data (Buote & Canizares 1994) -- the $\sim 30\arcdeg$ twist is significant at an estimated 99% confidence level. We argue that this twist is not the result of projected foreground and background sources, ram pressure effects, or tidal distortions. If spheroidal symmetry and a nearly isothermal hot gas are assumed, then the azimuthally averaged radial profile displays features which, when combined with the observed $PA$ twist, are inconsistent with the simple assumptions that the X-ray emission is due either entirely to hot gas or to the combined emission from hot gas and discrete sources. We discuss possible origins of the $PA$ twist and radial profile features (e.g., triaxiality).

We analyze with hydrodynamical simulations the evolution of galaxy clusters in a cosmological environment. Power ratios (Buote \& Tsai 1995) are used to quantitatively relate cluster morphologies to their dynamical states. The simulated clusters follow the same ``evolutionary track'' obeyed by a sample of low-redshift $(z<0.2)$ ROSAT PSPC clusters (Buote \& Tsai 1996) indicating that the detailed evolution of individual simulated clusters is consistent with observed clusters. However, the distribution of simulated clusters (for $\Omega=1$ standard Cold Dark Matter) along the evolutionary track at the present epoch, which indicates a measure of the present balance of cluster formation and relaxation rates, suggests that there are too many simulated clusters with significant amounts of substructure to be consistent with the observations, thus favoring a lower value of $\Omega$. Perpendicular to the evolutionary track the distributions of observed and simulated clusters are consistent which may indicate a success of the cosmological model (e.g., power spectrum). Analysis of high-redshift simulated clusters suggests that the distribution of clusters both along and perpendicular to the evolutionary track is effectively constant from $z\sim 0.6$ to the present but changes significantly for $z > 0.6$.

We quantify the morphologies and dynamical states of 59 galaxy clusters using the power-ratio technique of Buote & Tsai applied to ROSAT PSPC X-ray images. The clusters exhibit a particularly strong $P_2/P_0 - P_4/P_0$ correlation in the $1h^{-1}_{80}$ Mpc aperture which may be interpreted as an evolutionary track; the location of a cluster on the correlation line indicates the dynamical state of the cluster and the distribution of clusters along this track measures the rate of formation and evolution of clusters in our sample. The power ratios anti-correlate with the cooling-flow rate indicating a reasonable dependence of the flow rate on cluster morphology. The relationship of the power ratios to the optical Bautz-Morgan (BM) Type is more complex. This is because the power ratios are sensitive to unrelaxed regions of clusters within a specified scale, whereas BM types are sensitive to unrelaxed regions over many scales. We discuss further astrophysical applications exploiting the relationship between the power ratios and the evolutionary states of clusters.

We analyzed ROSAT PSPC images of the bright, nearby $(z<0.1)$ galaxy clusters A401, A1656 (Coma), A2029, A2199, and A2256 to constrain their intrinsic shapes and baryon fractions. Following Buote & Tsai we probed the aggregate structure of the clusters on scales $\sim 1.5h^{-1}_{80}$ Mpc to reduce effects of possible substructure on smaller scales ($\lesssim$ a few hundred kpc). The ellipticities of the X-ray isophotes are typically $\epsilon\approx 0.15 - 0.25$ and display negative radial gradients significant (except Coma). By assuming hydrostatic equilibrium and a variety of mass models we obtain ellipticities $\epsilon_{mass} \approx 0.40-0.55$ for isothermal models of the total gravitating matter of the clusters; the $\epsilon_{mass}$ constraints change by $<10%$ upon consideration of the small temperature gradients shown by ASCA to be typical for rich clusters; the observed X-ray ellipticity gradients require that mass models with constant ellipticity have steep density profiles, $\rho\propto r^{-4}$. Estimates of the gas masses are highly insensitive to the ellipticities of the X-ray isophotes. The clusters in our sample have increasing fractions of gas mass to total mass with radius and have $M_{gas}/M_{tot} = (4% - 11%)h^{-3/2}_{80}$ within a radius $1.5h^{-1}_{80}$ Mpc, in excellent agreement with the results of White & Fabian and the Baryon Catastrophe proposed by White et al. (1993). Finally, the ellipticities of the dark matter distributions are essentially identical to $\epsilon_{mass}$ and are consistent with the shapes of dark halos predicted by N-body simulations and the shapes of the galaxy isopleths in the clusters in contrast to our previous conclusions using {\it Einstein} data.

We have analyzed ROSAT PSPC X-ray data of the optically elongated S0 galaxy NGC 1332 with the purposes of constraining the intrinsic shape of its underlying mass and presenting a detailed investigation of the uncertainties resulting from the assumptions underlying this type of analysis. The X-ray isophotes are elongated with ellipticity $0.10 - 0.27$ (90% confidence) for semi-major axes $75\arcsec -90\arcsec$ and have orientations consistent with the optical isophotes (ellipticity $\sim 0.43$). The spectrum is poorly constrained by the PSPC data and cannot rule out sizeable radial temperature gradients or an emission component due to discrete sources equal in magnitude to the hot gas. Using (and clarifying) the "geometric test" for dark matter, we determined that the hypothesis that mass-traces-light is not consistent with the X-ray data at 68% confidence and marginally consistent at 90% confidence independent of the gas temperature profile. Detailed modeling gives constraints on the ellipticity of the underlying mass of \epsilon_{mass} = 0.47 - 0.72 (0.31 - 0.83) at 68% (90%) confidence for isothermal and polytropic models. The total mass of the isothermal models within a=43.6$ kpc (D = 20h^{-1}_{80} Mpc) is M_{tot} = (0.38 - 1.7) \times 10^{12}M_{\sun} (90% confidence) corresponding to total blue mass-to-light ratio \Upsilon_B = (31.9 - 143) \Upsilon_{\sun}. Similar results are obtained when the dark matter is fit directly using the known distributions of the stars and gas. When possible rotation of the gas and emission from discrete sources are included flattened mass distributions are still required, although the constraints on $\epsilon_{mass}$, but not the total mass, are substantially weakened.

We describe and test a method to quantitatively classify clusters of galaxies according to their projected morphologies. This method will be subsequently used to place constraints on cosmological parameters ($\Omega$ and the power spectrum of primordial fluctuations on scales at or slightly smaller than that of clusters) and to test theories of cluster formation. We specifically address structure that is easily discernible in projection and dynamically important to the cluster. The method is derived from the two-dimensional multipole expansion of the projected gravitational potential and yields dimensionless {\it power ratios} as morphological statistics. If the projected mass profile is used to characterize the cluster morphology, the power ratios are directly related to the cluster potential. However, since detailed mass profiles currently exist for only a few clusters, we use the X-ray--emitting gas as an alternative tracer of cluster morphology. In this case, the relation of the power ratios to the potential is qualitatively preserved. We demonstrate the feasibility of the method by analyzing simulated observations of simple models of X-ray clusters using the instrument parameters of the ROSAT PSPC. For illustrative purposes, we apply the method to ROSAT PSPC images of A85, A514, A1750, and A2029. These clusters, which differ substantially in their X-ray morphologies, are easily distinguished by their respective power ratios. We discuss the suitability of this method to address the connection between cluster morphology and cosmology and to assess whether an individual cluster is sufficiently relaxed for analysis of its intrinsic shape using hydrostatic methods. Approximately 50 X-ray observations of Abell clusters with the PSPC will be amenable to morphological analysis using the method of this paper.

Using the simulation of Katz & White we have tested the viability of X-ray analysis for constraining the intrinsic shapes of clusters of galaxies considering the effects of both substructure and steep temperature gradients. We restrict analysis to aggregate cluster shapes on scales of $r\sim 1-2$ Mpc to reduce our sensitivity to core subclustering. For low redshifts $(z\lesssim 0.25)$ the X-ray method accurately measures the true ellipticity of the three-dimensional cluster dark matter; assuming the gas is isothermal adds only small errors. At higher redshifts the X-ray method yields unreliable results since the gas does not trace the cluster gravitational potential. We proffer necessary conditions for the reliability of X-ray methods characterized by the amount of substructure in the X-ray surface brightness images and the isophotal shapes. We conclude that measurements of aggregate cluster shapes on scales $r\sim 1-2$ Mpc are insensitive to core substructure on scales of a few hundred kpc. Therefore our results suggest that the X-ray measurements of aggregate cluster shapes by Fabricant, Rybicki, & Gorenstein and Buote & Canizares are valid provided that they do not suffer from serious projection effects. A substantial number of Abell clusters observed with the ROSAT PSPC will be amenable to X-ray shape analysis.

(shortened for babbage) We describe (1) a new test for dark matter and alternate theories of gravitation based on the relative geometries of the X-ray and optical surface brightness distributions and an assumed form for the potential of the optical light, (2) a technique to measure the shapes of the total gravitating matter and dark matter of an ellipsoidal system which is insensitive to the precise value of the temperature of the gas and to modest temperature gradients, and (3) a method to determine the ratio of dark mass to stellar mass that is dependent on the functional forms for the visible star, gas and dark matter distributions, but independent of the distance to the galaxy or the gas temperature. We apply these techniques to X-ray data from the ROSAT Position Sensitive Proportional Counter (PSPC) of the optically-flattened elliptical galaxy NGC 720. The X-ray isophotes are significantly elongated, $\epsilon = 0.20 - 0.30$ (90% confidence) for semi-major axis $a\sim 100\arcsec$. The major axes of the optical and X-ray isophotes are misaligned by $\sim 30\deg \pm 15\deg$ (90% confidence). We conclude that matter distributed like the optical light cannot produce the observed ellipticities of the X-ray isophotes independent of the gas pressure, the gas temperature, and the value of the stellar mass.