No Arabic abstract
Simple but flexible dynamical models are useful for many purposes, including serving as the starting point for more complex models or numerical simulations of galaxies, clusters, or dark matter haloes. We present SpheCow, a new light-weight and flexible code that allows one to easily explore the structure and dynamics of any spherical model. Assuming an isotropic or Osipkov-Merritt anisotropic orbital structure, the code can automatically calculate the dynamical properties of any model with either an analytical density profile or an analytical surface density profile as starting point. We have extensively validated SpheCow using a combination of comparisons to analytical and high-precision numerical calculations, as well as the calculation of inverse formulae. SpheCow contains readily usable implementations for many standard models, including the Plummer, Hernquist, NFW, Einasto, Sersic, and Nuker models. The code is publicly available as a set of C++ routines and as a Python module, and it is designed to be easily extendable, in the sense that new models can be added in a straightforward way. We demonstrate this by adding two new families of models in which either the density slope or the surface density slope is described by an algebraic sigmoid function. We advocate the use of the SpheCow code to investigate the full dynamical structure for models for which the distribution function cannot be expressed analytically and to explore a much wider range of models than is possible using analytical models alone.
Galaxy kinematics and gravitational lensing are two complementary ways to constrain the distribution of dark matter on galaxy scales. The typical dark matter density profiles adopted in dynamical studies cannot easily be adopted in lensing studies. Ideally, a mass model should be used that has the global characteristics of realistic dark matter distributions, and that allows for an analytical calculation of the magnifications and deflection angles. A simple model with these properties, the broken-power-law (BPL) model, has very recently been introduced. We examine the dynamical structure of the family of BPL models. We derive simple closed expressions for basic dynamical properties, and study the distribution function under the assumption of velocity isotropy. We find that none of the BPL models with realistic parameters has an isotropic distribution function that is positive over the entire phase space, implying that the BPL models cannot be supported by an isotropic velocity distribution, or models with a more radially anisotropic orbital structure. This result limits the attractiveness of the BPL family as a tool for lensing studies to some degree. More generally, we find that not all members of the general family of double power-law or Zhao models, often used to model dark matter haloes, can be supported by an isotropic or radially anisotropic distribution function. In other words, the distribution function may become negative even for spherically symmetric models with a well-behaved density profile.
Using estimates of dark halo masses from satellite kinematics, weak gravitational lensing, and halo abundance matching, combined with the Tully-Fisher and Faber-Jackson relations, we derive the mean relation between the optical, V_opt, and virial, V_200, circular velocities of early- and late-type galaxies at redshift z~0. For late-type galaxies V_opt ~ V_200 over the velocity range V_opt=90-260 km/s, and is consistent with V_opt = V_maxh (the maximum circular velocity of NFW dark matter haloes in the concordance LCDM cosmology). However, for early-type galaxies V_opt e V_200, with the exception of early-type galaxies with V_opt simeq 350 km/s. This is inconsistent with early-type galaxies being, in general, globally isothermal. For low mass (V_opt < 250 km/s) early-types V_opt > V_maxh, indicating that baryons have modified the potential well, while high mass (V_opt > 400 km/s) early-types have V_opt < V_maxh. Folding in measurements of the black hole mass - velocity dispersion relation, our results imply that the supermassive black hole - halo mass relation has a logarithmic slope which varies from ~1.4 at halo masses of ~10^{12} Msun/h to ~0.65 at halo masses of 10^{13.5} Msun/h. The values of V_opt/V_200 we infer for the Milky Way and M31 are lower than the values currently favored by direct observations and dynamical models. This offset is due to the fact that the Milky Way and M31 have higher V_opt and lower V_200 compared to typical late-type galaxies of the same stellar masses. We show that current high resolution cosmological hydrodynamical simulations are unable to form galaxies which simultaneously reproduce both the V_opt/V_200 ratio and the V_opt-M_star (Tully-Fisher/Faber-Jackson) relation.
The development of methods and algorithms to solve the $N$-body problem for classical, collisionless, non-relativistic particles has made it possible to follow the growth and evolution of cosmic dark matter structures over most of the Universes history. In the best studied case $-$ the cold dark matter or CDM model $-$ the dark matter is assumed to consist of elementary particles that had negligible thermal velocities at early times. Progress over the past three decades has led to a nearly complete description of the assembly, structure and spatial distribution of dark matter haloes, and their substructure in this model, over almost the entire mass range of astronomical objects. On scales of galaxies and above, predictions from this standard CDM model have been shown to provide a remarkably good match to a wide variety of astronomical data over a large range of epochs, from the temperature structure of the cosmic background radiation to the large-scale distribution of galaxies. The frontier in this field has shifted to the relatively unexplored subgalactic scales, the domain of the central regions of massive haloes, and that of low-mass haloes and subhaloes, where potentially fundamental questions remain. Answering them may require: (i) the effect of known but uncertain baryonic processes (involving gas and stars), and/or (ii) alternative models with new dark matter physics. Here we present a review of the field, focusing on our current understanding of dark matter structure from $N$-body simulations and on the challenges ahead.
Comparison of observed satellite galaxies of the Milky Way (hereafter MW) with dark matter subhaloes in cosmological $N$-body simulations of MW-mass haloes suggest that such subhaloes, if they exist, are occupied by satellites in a stochastic fashion. We examine how inefficient massive star formation and associated supernova feedback in high-redshift progenitors of present-day low-mass subhaloes might contribute to this stochasticity. Using a Monte Carlo approach to follow the assembly histories of present-day low-mass haloes with $10^7 lesssim M leq 10^{10}$ ${rm M}_{odot}$, we identify when cooling and star formation is likely to proceed, and observe that haloes with present-day masses $lesssim 10^9 {rm M}_{odot}$ never grow sufficiently massive to support atomic hydrogen line cooling. Noting that the star formation timescale decreases sharply with stellar mass as $t_{rm PMS} propto m_{ast}^{-2.5}$, we argue that, should the conditions for high mass star formation arise in low-mass haloes, the ensuing supernovae are likely to disrupt ongoing lower-mass star formation and unbind gas within the halo. This potentially star-forming gas is unlikely to be replenished in lower mass haloes because of, e.g. cosmological reionization, and so we expect galaxy formation to be stymied in a manner that depends on host halo assembly history and the efficiency and timing of star formation in proto-galaxies, which we illustrate using a Monte Carlo model. Based on these simple physical arguments, we assert that stochasticity of star formation and feedback is an essential but overlooked ingredient in modelling galaxy formation on the smallest scales.
Galaxy-galaxy weak lensing is a direct probe of the mean matter distribution around galaxies. The depth and sky coverage of the CFHT Legacy Survey yield statistically significant galaxy halo mass measurements over a much wider range of stellar masses ($10^{8.75}$ to $10^{11.3} M_{odot}$) and redshifts ($0.2 < z < 0.8$) than previous weak lensing studies. At redshift $z sim 0.5$, the stellar-to-halo mass ratio (SHMR) reaches a maximum of $4.0pm0.2$ percent as a function of halo mass at $sim 10^{12.25} M_{odot}$. We find, for the first time from weak lensing alone, evidence for significant evolution in the SHMR: the peak ratio falls as a function of cosmic time from $4.5 pm 0.3$ percent at $z sim 0.7$ to $3.4 pm 0.2$ percent at $z sim 0.3$, and shifts to lower stellar mass haloes. These evolutionary trends are dominated by red galaxies, and are consistent with a model in which the stellar mass above which star formation is quenched downsizes with cosmic time. In contrast, the SHMR of blue, star-forming galaxies is well-fit by a power law that does not evolve with time. This suggests that blue galaxies form stars at a rate that is balanced with their dark matter accretion in such a way that they evolve along the SHMR locus. The redshift dependence of the SHMR can be used to constrain the evolution of the galaxy population over cosmic time.