Laboratory experiments, large-scale computer simulations and observational cosmology have begun to make progress in the campaign to identify the particle responsible for gravitationally-inferred dark matter. In this contribution we discuss the dark matter density profiles in the cores of nearby galaxy clusters and estimate the gamma-ray flux expected for MSSM dark matter over a range of neutralino masses.

We present a spectroscopic deprojection analysis of a sample of ten relaxed galaxy clusters. We use an empirical F-test derived from a set of Markov Chain Monte Carlo simulations to determine if the core plasma in each cluster could contain multiple phases. We derive non-parametric baryon density and temperature profiles, and use these to construct total gravitating mass profiles. We compare these profiles with the standard halo parameterizations. We find central density slopes roughly consistent with the predictions of LCDM: $-1 \lesssim d\log(\rho)/d\log(r) \lesssim -2$. We constrain the core size of each cluster and, using the results of cosmological simulations as a calibrator, place an upper limit of ~0.1 cm^2/g = 0.2 b(GeV/c^2)^{-1} (99% confidence) on the dark matter particle self-interaction cross section.

Knowledge of the structure of galaxy clusters is essential for an understanding of large scale structure in the universe, and may provide important clues to the nature of dark matter. Moreover, the shape of the dark matter distribution in the cluster core may offer insight into the structure formation process. Unfortunately, cluster cores also tend to be the site of complicated astrophysics. X-ray imaging spectroscopy of relaxed clusters, a standard technique for mapping their dark matter distributions, is often complicated by the presence of cool components in cluster cores, and the dark matter profile one derives for a cluster is sensitive to assumptions made about the distribution of this component. In addition, fluctuations in the temperature measurements resulting from normal statistical variance can produce results which are unphysical. We present here a procedure for extracting the dark matter profile of a spherically symmetric, relaxed galaxy cluster which deals with both of these complications. We apply this technique to a sample of galaxy clusters observed with the Chandra X-ray Observatory, and comment on the resulting mass profiles. For some of the clusters we compare their masses with those derived from weak and strong gravitational measurements.

We examine the X-ray spectra and variability of the sample of X-ray sources with L_X = 10^{31}-10^{33} erg s^{-1} identified within the inner 9' of the Galaxy. Very few of the sources exhibit intra-day or inter-month variations. We find that the spectra of the point sources near the Galactic center are very hard between 2--8 keV, even after accounting for absorption. When modeled as power laws the median photon index is Gamma=0.7, while when modeled as thermal plasma we can only obtain lower limits to the temperature of kT>8 keV. The combined spectra of the point sources is similarly hard, with a photon index of Gamma=0.8. Strong line emission is observed from low-ionization, He-like, and H-like Fe, both in the average spectra and in the brightest individual sources. The line ratios of the highly-ionized Fe in the average spectra are consistent with emission from a plasma in thermal equilibrium. This line emission is observed whether average spectra are examined as a function of the count rate from the source, or as a function of the hardness ratios of individual sources. This suggests that the hardness of the spectra may in fact to due local absorption that partially-covers the X-ray emitting regions in the Galactic center systems. We suggest that most of these sources are intermediate polars, which (1) often exhibit hard spectra with prominent Fe lines, (2) rarely exhibit either flares on short time scales or changes in their mean X-ray flux on long time scales, and (3) are the most numerous hard X-ray sources with comparable luminosities in the Galaxy.

Observations of extragalactic objects need to be corrected for Galactic absorption and this is often accomplished by using the measured 21 cm HI column. However, within the beam of the radio telescope there are variations in the HI column that can have important effects in interpreting absorption line studies and X-ray spectra at the softest energies. We examine the HI and DIRBE/IRAS data for lines of sight out of the Galaxy, which show evidence for HI variations in of up to a factor of three in 1 degree fields. Column density enhancements would preferentially absorb soft X-rays in spatially extended objects and we find evidence for this effect in the ROSAT PSPC observations of two bright clusters of galaxies, Abell 119 and Abell 2142. For clusters of galaxies, the failure to include column density fluctuations will lead to systematically incorrect fits to the X-ray data in the sense that there will appear to be a very soft X-ray excess. This may be one cause of the soft X-ray excess in clusters, since the magnitude of the effect is comparable to the observed values.

We present a spectral analysis of the central X-ray emission for a sample of galaxy clusters observed with Chandra. We constrain the quantity of a second cospatial temperature component using Markov Chain Monte Carlo sampling and discuss the implications for our understanding of cooling flows.

We report Chandra and XMM-Newton observations of the transient neutron star low-mass X-ray binary GRS 1741.9-2853. Chandra detected the source in outburst on 2000 October 26 at an X-ray luminosity of ~10^{36} erg/s (2--8 keV; 8 kpc), and in quiescence on 2001 July 18 at ~10^{32} erg/s. The latter observation is the first detection of GRS 1741.9-2853 in quiescence. We obtain an accurate position for the source of 17h 45m 2.33s, -28o 54' 49.7" (J2000), with an uncertainty of 0.7". GRS 1741.9-2853 was not detected significantly in three other Chandra observations, nor in three XMM-Newton observations, indicating that the luminosity of the source in quiescence varies by at least a factor of 5 between (< 0.9 - 5.0) \times 10^{32} erg/s (2--8 keV). A weak X-ray burst with a peak luminosity of 5 \times 10^{36} erg/s above the persistent level was observed with Chandra during the outburst on 2000 October 26. The energy of this burst, 10^{38} erg, is unexpectedly low, and may suggest that the accreted material is confined to the polar caps of the neutron star. A search of the literature reveals that GRS 1741.9-2853 was observed in outburst with ASCA in Fall 1996 as well, when the BeppoSAX WFC detected the three previous X-ray bursts from this source. The lack of X-ray bursts from GRS 1741.9-2853 at other epochs suggests that it produces bursts only during transient outbursts when the accretion rate onto the surface of the neutron star is about 10^{-10} M_sun/yr. A similar situation may hold for other low-luminosity bursters recently identified from WFC data.

Galaxy cluster mass distributions offer an important test of the cold dark matter picture of structure formation, and may even contain clues about the nature of dark matter. X-ray imaging spectroscopy of relaxed systems can map cluster dark matter distributions, but are usually complicated by the presence of central cool components in the intracluster medium. Here we describe a statistically correct approach to distinguishing amongst simple alternative models of the cool component, and apply it to one cluster. We also present mass profiles and central density slopes for five clusters derived from Chandra data, and illustrate how assumptions about the cool component affect the resulting mass profiles. For four of these objects, we find that the central density profile (r < 200 h_50^-1 kpc) rho(r) = r^a with -2 < a < -1, for either of two models of the central cool component. These results are consistent with standard CDM predictions.

Determining the structure of galaxy clusters is essential for an understanding of large scale structure in the universe, and may hold important clues to the identity and nature of dark matter particles. Moreover, the core dark matter distribution may offer insight into the structure formation process. Unfortunately, cluster cores also tend to be the site of complicated astrophysics. X-ray imaging spectroscopy of relaxed clusters, a standard technique for mapping their dark matter distributions, is often complicated by the presence of their putative ``cooling flow'' gas, and the dark matter profile one derives for a cluster is sensitive to assumptions made about the distribution of this gas. Here we present a statistical analysis of these assumptions and their effect on our understanding of dark matter in galaxy clusters.

We present Chandra observations of EMSS 1358+6245, a relaxed cooling flow cluster of galaxies at z = 0.328. We employ a new deprojection technique to construct temperature, gas, and dark matter profiles. We confirm the presence of cool gas in the cluster core, and our deprojected temperature profile for the hot component is isothermal over 30 kpc < r < 0.8 Mpc. Fitting the mass profile to an NFW model yields r_s = 153 [+161,-83] kpc and c = 8.4 [+3.4,-2.3]. We find good agreement between our dark matter profile and weak gravitational lensing measurements. We place an upper limit of 42 kpc (90% confidence limit) on the size of any constant density core. We compare this result to recent simulations and place a conservative upper limit on the dark matter particle scattering cross section of 0.1 cm^2/g. This limit implies that the cross-section must be velocity dependent if the relatively shallow core mass profiles of dwarf galaxies are a direct result of dark matter self-interaction.

Chandra observations of lensing galaxy clusters have now provided accurate dark matter profiles for several objects in which the intracluster medium is likely to be in hydrostatic equilibrium. We discuss Chandra observations of the mass profile of one such cluster, EMSS1358+6245. We find no evidence for flattening of the mass density profile at radii greater than 50/h_50 kpc . This result, and similar findings from Chandra observations of other clusters, appear to rule out models in of dark matter self-interaction proposed to explain the flat cores of low-surface brightness galaxies.

We have simulated the evolution of tidal debris in the Galactic halo in order to guide our ongoing survey to determine the fraction of halo mass accreted via satellite infall. Contrary to naive expectations that the satellite debris will produce a single narrow velocity peak on a smooth distribution, there are many different signatures of substructure, including multiple peaks and broad but asymmetrical velocity distributions. Observations of the simulations show that there is a high probability of detecting the presence of tidal debris with a pencil beam survey of 100 square degrees. In the limiting case of a single 10^7 Msun satellite contributing 1% of the luminous halo mass the detection probability is a few percent using just the velocities of 100 halo stars in a single 1 square degree field. The detection probabilities scale with the accreted fraction of the halo and the number of fields surveyed. There is also surprisingly little dependence of the detection probabilities on the time since the satellite became tidally disrupted, or on the initial orbit of the satellite, except for the time spent in the survey volume.

We have developed a technique to map the three-dimensional structure of the local interstellar medium using a maximum entropy reconstruction technique. A set of column densities N to stars of known distance can in principle be used to recover a three-dimensional density field n, since the two quantities are related by simple geometry through the equation N = C n, where C is a matrix characterizing the stellar spatial distribution. In practice, however, there is an infinte number of solutions to this equation. We use a maximum entropy reconstruction algorithm to find the density field containing the least information which is consistent with the observations. The solution obtained with this technique is, in some sense, the model containing the minimum structure. We apply the algorithm to several simulated data sets to demonstrate its feasibility and success at recovering ``real'' density contrasts. This technique can be applied to any set of column densities whose end points are specified. In a subsequent paper we shall describe the application of this method to a set of stellar color excesses to derive a map of the dust distribution, and to soft X-ray absorption columns to hot stars to derive a map of the total density of the interstellar medium.

A study of the cores of galaxy clusters with the Einstein SSS indicated the presence of absorbing material corresponding to 1E+12 Msun of cold cluster gas, possibly resulting from cooling flows. Since this amount of cold gas is not confirmed by observations at other wavelengths, we examined whether this excess absorption is present in the ROSAT PSPC observations of 20 bright galaxy clusters. For 3/4 of the clusters, successful spectral fits were obtained with absorption due only to the Galaxy, and therefore no extra absorption is needed within the clusters, in disagreement with the results from the Einstein SSS data for some of the same clusters. For 1/4 of the clusters, none of our spectral fits was acceptable, suggesting a more complicated cluster medium than the two-temperature and cooling flow models considered here. However, even for these clusters, substantial excess absorption is not indicated.

An extremely soft X-ray excess throughout galaxy clusters has been claimed as a new feature of these sytems, with important physical implications. We have reexamined this feature in the five clusters for which it has been discussed, using the most recent X-ray absorption cross sections, X-ray data processing techniques, and a consistent set of HI data. For the Virgocluster, we find that the spectrum can be fit with a single-temperature thermal plasma and with an X-ray absorption column that is not significantly different than the Galactic HI column. The result for Abell 1367, Abell 1656 (Coma), Abell 1795, and Abell 2199 is similar in that the difference betweenthe X-ray absorption column and the Galactic HI column is less than 3-sigma for He/H = 0.09, and for He/H = 0.10 only one cluster location leads to a Galactic HI column more than 3-sigma above the X-ray absorption column (Coma, with one location with a 3.6-sigma difference). We conclude that there is no strong evidence for the extremely soft X-ray excess in galaxy clusters.

This paper, and its companion, investigate the evolution of dense stellar systems due to the influence of two-body gravitational encounters, physical collisions and stellar evolution. Our goal is the simulation of the densest centers of galaxies, like M32, which reach stellar densities near a million solar masses per cubic parsec and which may harbor black holes. These systems have a different Safronov number (the dimensionless ratio of stellar binding energy to mean stellar kinetic energy) than globular clusters, substantially increasing the importance of physical collisions relative to gravitational encounters. In this paper, we focus only on the gravitational encounters. We demonstrate, first, that our simulations with small N with a Hernquist tree code yield results basically in accord with years of effort studying globular clusters. Second, we investigate (crudely) core collapse in rotating systems with mass segregation, to separate out the effects purely due to two-body encounters from those seen in the more complex second paper. Consistent with previous studies, we find that systems whose constituent particles follow the Salpeter initial mass function rapidly undergo core collapse through Spitzer's mass segregation instability. All rotationally flattened systems show a decrease in flattening with time, consistent with Fokker-Planck calculations. An interesting new result concerns the well-established inability of simulators to identify a static center in their simulations of collisional systems. We find that the lump of high mass stars which condenses at the center wanders about the core in Brownian motion.

We report the results of simulations of dense rotating stellar systems whose members suffer collisions and undergo stellar evolution processes. The initial configuration for each experiment is an isotropic Kuzmin-Kutuzov model. The dynamical evolution is simulated with the N-body tree code of Hernquist, modified to incorporate physical stellar collisions, stellar evolution, stellar mass loss and compact remnant formation, and star formation. In some simulations we have added a large accreting central black hole. In all systems the velocity dispersion in the halo evolves toward a radially biased state. In systems containing a central black hole, the dispersion becomes tangentially biased in the core, whereas it remains isotropic in systems with no black hole. Collisions tend to produce a single dominant stellar merger product, as opposed to a swarm of intermediate-mass stars. In cases where we have suppressed all processes except relaxation and physical collisions, objects with greater flattening produce larger stars through mergers. In systems where stellar ejecta are allowed to escape the system, mass loss from the heavy core stars temporarily reduces the core density and collision rate. Most of the simulations performed reproduced the ratio of the central collision time scale to the central relaxation time scale found in the dwarf elliptical galaxy M32. The rapid central evolution of these systems due to collisions and relaxation, combined with scaling the results in N, suggests that we are either viewing M32 at a peculiar moment in its history, or that its dynamically-inferred central density is at least in part due to the present of a massive dark object, presumably a black hole.

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