The Compression of Dark Matter Halos by Baryonic Infall

The initial radial density profiles of dark matter halos are laid down by gravitational collapse in hierarchical structure formation scenarios and are subject to further compression as baryons cool and settle to the halo centers. We here describe an explicit implementation of the algorithm, originally developed by Young, to calculate changes to the density profile as the result of adiabatic infall in a spherical halo model. Halos with random motion are more resistant to compression than are those in which random motions are neglected, which is a key weakness of the simple method widely employed. Young's algorithm results in density profiles in excellent agreement with those from N-body simulations. We show how the algorithm may be applied to determine the original uncompressed halos of real galaxies, a step which must be computed with care in order to enable a confrontation with theoretical predictions from theories such as LCDM.Comment: Revised version for ApJ. 8 pages, 8 figures, latex uses emulateap

The Evolution of X-ray Clusters and the Entropy of the Intra Cluster Medium

The thermodynamics of the diffuse, X-ray emitting gas in clusters of galaxies is determined by gravitational processes associated with shock heating, adiabatic compression, and non-gravitational processes such as heating by SNe, stellar winds, activity in the central galactic nucleus, and radiative cooling. The effect of gravitational processes on the thermodynamics of the Intra Cluster Medium (ICM) can be expressed in terms of the ICM entropy S ~ ln(T/\rho^{2/3}). We use a generalized spherical model to compute the X-ray properties of groups and clusters for a range of initial entropy levels in the ICM and for a range of mass scales, cosmic epochs and background cosmologies. We find that the statistical properties of the X-ray clusters strongly depend on the value of the initial excess entropy. Assuming a constant, uniform value for the excess entropy, the present-day X-ray data are well fitted for the following range of values K_* = kT/\mu m_p \rho^{2/3} = (0.4\pm 0.1) \times 10^{34} erg cm^2 g^{-5/3} for clusters with average temperatures kT>2 keV; K_* = (0.2\pm 0.1) \times 10^{34} erg cm^2 g^{-5/3} for groups and clusters with average temperatures kT<2 keV. These values correspond to different excess energy per particle of kT \geq 0.1 (K_*/0.4\times 10^{34}) keV. The dependence of K_* on the mass scale can be well reproduced by an epoch dependent external entropy: the relation K_* = 0.8(1+z)^{-1}\times 10^{34} erg cm^2 g^{-5/3} fits the data over the whole temperature range. Observations of both local and distant clusters can be used to trace the distribution and the evolution of the entropy in the cosmic baryons, and ultimately to unveil the typical epoch and the source of the heating processes.Comment: 53 pages, LateX, 19 figures, ApJ in press, relevant comments and references adde

How Do Galaxies Get Their Gas?

We examine the temperature history of gas accreted by forming galaxies in SPH simulations. About half the gas shock heats to roughly the virial temperature of the galaxy potential well before cooling, condensing, and forming stars, but the other half radiates its acquired gravitational energy at much lower temperatures, typically T<10^5 K, and the histogram of maximum gas temperatures is clearly bimodal. The "cold mode" of gas accretion dominates for low mass galaxies (M_baryon < 10^{10.3}Msun or M_halo < 10^{11.4}Msun), while the conventional "hot mode" dominates the growth of high mass systems. Cold accretion is often directed along filaments, allowing galaxies to efficiently draw gas from large distances, while hot accretion is quasi-spherical. The galaxy and halo mass dependence leads to redshift and environment dependence of cold and hot accretion rates, with cold mode dominating at high redshift and in low density regions today, and hot mode dominating in group and cluster environments at low redshift. Star formation rates closely track accretion rates, and we discuss the physics behind the observed environment and redshift dependence of galactic scale star formation. If we allowed hot accretion to be suppressed by conduction or AGN feedback, then the simulation predictions would change in interesting ways, perhaps resolving conflicts with the colors of ellipticals and the cutoff of the galaxy luminosity function. The transition between cold and hot accretion at M_h ~ 10^{11.4}Msun is similar to that found by Birnboim & Dekel (2003) using 1-d simulations and analytic arguments. The corresponding baryonic mass is tantalizingly close to the scale at which Kauffmann et al. (2003) find a marked shift in galaxy properties. We speculate on connections between these theoretical and observational transitions.Comment: 1 figure added, Appendix discussing SAMs added, some text changes. Matches the version accepted by MNRAS. 31 pages (MNRAS style), 21 figures,For high resolution version of the paper (highly recommended) follow http://www.astro.umass.edu/~keres/paper/ms2.ps.g

A simple analytical model for dark matter halo structure and adiabatic contraction

A simple analytical model for describing inner parts of dark matter halo is considered. It is assumed that dark matter density is power-law. The model deals with dark matter distribution function in phase space of adiabatic invariants (radial action and angular momentum). Two variants are considered for the angular part of the distribution function: narrow and broad distribution. The model allows to describe explicitly the process of adiabatic contraction of halo due to change of gravitational potential caused by condensation of baryonic matter in the centre. The modification of dark matter density in the centre is calculated, and is it shown that the standard algorithm of adiabatic contraction calculation overestimates the compressed halo density, especially in the case of strong radial anisotropy.Comment: 5 pages, 3 figures. v3 - major improvements, another halo model introduced, discussion extende

Detection of the Entropy of the Intergalactic Medium: Accretion Shocks in Clusters, Adiabatic Cores in Groups

The thermodynamics of the diffuse, X-ray emitting gas in clusters of galaxies is linked to the entropy level of the intra cluster medium. In particular, models that successfully reproduce the properties of local X-ray clusters and groups require the presence of a minimum value for the entropy in the center of X-ray halos. Such a minimum entropy is most likely generated by non-gravitational processes, in order to produce the observed break in self-similarity of the scaling relations of X-ray halos. At present there is no consensus on the level, the source or the time evolution of this excess entropy. In this paper we describe a strategy to investigate the physics of the heating processes acting in groups and clusters. We show that the best way to extract information from the local data is the observation of the entropy profile at large radii in nearby X-ray halos (z~0.1), both at the upper and lower extremes of the cluster mass scale. The spatially and spectrally resolved observation of such X-ray halos provides information on the mechanism of the heating. We demonstrate how measurements of the size of constant entropy (adiabatic) cores in clusters and groups can directly constrain heating models, and the minimum entropy value. We also consider two specific experiments: the detection of the shock fronts expected at the virial boundary of rich clusters, and the detection of the isentropic, low surface-brightness emission extending to radii larger than the virial ones in low mass clusters and groups. Such observations will be a crucial probe of both the physics of clusters and the relationship of non-gravitational processes to the thermodynamics of the intergalactic medium.Comment: ApJ accepted, 31 pages including 8 figures. Important material added; references update

Density profiles of dark matter haloes on Galactic and Cluster scales

In the present paper, we improve the "Extended Secondary Infall Model" (ESIM) of Williams et al. (2004) to obtain further insights on the cusp/core problem. The model takes into account the effect of ordered and random angular momentum, dynamical friction and baryon adiabatic contraction in order to obtain a secondary infall model more close to the collapse reality. The model is applied to structures on galactic scales (normal and dwarf spiral galaxies) and on cluster of galaxies scales. The results obtained suggest that angular momentum and dynamical friction are able, on galactic scales, to overcome the competing effect of adiabatic contraction eliminating the cusp. The NFW profile can be reobtained, in our model only if the system is constituted just by dark matter and the magnitude of angular momentum and dynamical friction are reduced with respect to the values predicted by the model itself. The rotation curves of four LSB galaxies from de Blok & Bosma (2002) are compared to the rotation curves obtained by the model in the present paper obtaining a good fit to the observational data. On scales smaller than $\simeq 10^{11} h^{-1} M_{\odot}$ the slope $\alpha \simeq 0$ and on cluster scales we observe a similar evolution of the dark matter density profile but in this case the density profile slope flattens to $\alpha \simeq 0.6$ for a cluster of $\simeq 10^{14} h^{-1} M_{\odot}$. The total mass profile, differently from that of dark matter, shows a central cusp well fitted by a NFW model.Comment: 26 pages; 4 figures A&A Accepte

The Planes of Satellite Galaxies Problem, Suggested Solutions, and Open Questions

Satellite galaxies of the Milky Way and of the Andromeda galaxy have been found to preferentially align in significantly flattened planes of satellite galaxies, and available velocity measurements are indicative of a preference of satellites in those structures to co-orbit. There is increasing evidence that such kinematically correlated satellite planes are also present around more distant hosts. Detailed comparisons show that similarly anisotropic phase-space distributions of sub-halos are exceedingly rare in cosmological simulations based on the $\Lambda$CDM paradigm. Analogs to the observed systems have frequencies of $\leq 0.5$ per cent in such simulations. In contrast to other small-scale problems, the satellite planes issue is not strongly affected by baryonic processes because the distribution of sub-halos on scales of hundreds of kpc is dominated by gravitational effects. This makes the satellite planes one of the most serious small-scale problem for $\Lambda$CDM. This review summarizes the observational evidence for planes of satellite galaxies in the Local Group and beyond, and provides an overview of how they compare to cosmological simulations. It also discusses scenarios which aim at explaining the coherence of satellite positions and orbits, and why they all are currently unable to satisfactorily resolve the issue.Comment: Invited review for MPLA, accepted for publication. 26 pages, 3 figure

Globular Cluster Formation from Colliding Substructure

We investigate a scenario where the formation of Globular Clusters (GCs) is triggered by high-speed collisions between infalling atomic-cooling subhalos during the assembly of the main galaxy host, a special dynamical mode of star formation that operates at high gas pressures and is intimately tied to LCDM hierarchical galaxy assembly. The proposed mechanism would give origin to "naked" globulars, as colliding dark matter subhalos and their stars will simply pass through one another while the warm gas within them clashes at highly supersonic speed and decouples from the collisionless component, in a process reminiscent of the Bullet galaxy cluster. We find that the resulting shock-compressed layer cools on a timescale that is typically shorter than the crossing time, first by atomic line emission and then via fine-structure metal-line emission, and is subject to gravitational instability and fragmentation. Through a combination of kinetic theory approximation and high-resolution $N$-body simulations, we show that this model may produce: (a) a GC number-halo mass relation that is linear down to dwarf galaxy scales and agrees with the trend observed over five orders of magnitude in galaxy mass; (b) a population of old globulars with a median age of 12 Gyr and an age spread similar to that observed; (c) a spatial distribution that is biased relative to the overall mass profile of the host; and (d) a bimodal metallicity distribution with a spread similar to that observed in massive galaxies.Comment: 15 pages, 5 figures, accepted for publication by the Astrophysical Journa