Completely dark galaxies: their existence, properties and strategies for finding them

There are a number of theoretical and observational hints that large numbers of low-mass galaxies composed entirely of dark matter exist in the field. The theoretical considerations follow from the prediction of cold dark matter theory that there exist many low-mass galaxies for every massive one. The observational considerations follow from the observed paucity of these low-mass galaxies in the field but not in dense clusters of galaxies; this suggests that the lack of small galaxies in the field is due to the inhibition of star formation in the galaxies as opposed to the fact that their small dark matter haloes do not exist. In this work we outline the likely properties of low-mass dark galaxies, and describe observational strategies for finding them, and where in the sky to search. The results are presented as a function of the global properties of dark matter, in particular the presence or absence of a substantial baryonic dark matter component. If the dark matter is purely cold and has a Navarro, Frenk & White density profile, directly detecting dark galaxies will only be feasible with present technology if the galaxy has a maximum velocity dispersion in excess of 70 km s⁻¹, in which case the dark galaxies could strongly lens background objects. This is much higher than the maximum velocity dispersions in most dwarf galaxies. If the dark matter in galaxy haloes has a baryonic component close to the cosmic ratio, the possibility of directly detecting dark galaxies is much more realistic; the optimal method of detection will depend on the nature of the dark matter. A number of more indirect methods are also discussed.     

Apocalypse soon

Based upon a simple vacuum Lagrangian, comprising cosmological and quadratic scalar field terms, a cosmological model is presented the history of which is indistinguishable from that of an innocuous low-density cold dark matter (CDM) universe, but the future of which is very much shorter. For sensible values of the deceleration parameter (0< q ₀<1), its age is greater than 85 per cent of the Hubble time, thus resolving the current version of the age crisis, which appears to be that t ₀∼1/ H ₀ while q ₀ is significantly positive.     

Determination of cosmology and evolution from K-band magnitude–redshift and number count observations

We determine cosmological and evolutionary parameters from the 3CR K -band Hubble diagram and K -band number counts, assuming that the galaxies in question undergo pure luminosity evolution. Separately the two data sets are highly degenerate with respect to choice of cosmological and evolutionary parameters, but in combination the degeneracy is resolved. Of models that either are flat or have  Ω_Λ=0  , the preferred ones are close to the canonical case  Ω_{cold  matter}=1  ,  Ω_Λ=0  , with luminosity evolution amounting to 1 mag brighter at   z =1  .     

Clustering of galaxy clusters in cold dark matter universes

We use very large cosmological N -body simulations to obtain accurate predictions for the two-point correlations and power spectra of mass-limited samples of galaxy clusters. We consider two currently popular cold dark matter (CDM) cosmogonies, a critical density model ( τ CDM) and a flat low density model with a cosmological constant (ΛCDM). Our simulations each use 10⁹ particles to follow the mass distribution within cubes of side 2  h ⁻¹ Gpc ( τ CDM) and 3  h ⁻¹ Gpc (ΛCDM) with a force resolution better than 10⁻⁴ of the cube side. We investigate how the predicted cluster correlations increase for samples of increasing mass and decreasing abundance. Very similar behaviour is found in the two cases. The correlation length increases from     for samples with mean separation     to     for samples with     The lower value here corresponds to τ CDM and the upper to ΛCDM. The power spectra of these cluster samples are accurately parallel to those of the mass over more than a decade in scale. Both correlation lengths and power spectrum biases can be predicted to better than 10 per cent using the simple model of Sheth, Mo & Tormen. This prediction requires only the linear mass power spectrum and has no adjustable parameters. We compare our predictions with published results for the automated plate measurement (APM) cluster sample. The observed variation of correlation length with richness agrees well with the models, particularly for ΛCDM. The observed power spectrum (for a cluster sample of mean separation     ) lies significantly above the predictions of both models.     

The direct detection of non-baryonic dark matter in the Galaxy?

It has been argued in a number of recent papers that dark matter is in the form of Jupiter-mass primordial black holes that betray their presence by microlensing quasars. This lensing accounts for a number of characteristic properties of quasar light curves, in both single quasars and gravitationally lensed multiple systems, that are not explained on the basis of intrinsic variation. One prediction of this idea is that Jupiter-mass bodies will be detected by the MACHO experiment as short events of about 2 d duration, although the expected frequency of detection is still very hard to estimate. However, the recent report by the MACHO group of the detection of a Jupiter-mass body in the direction of the Galactic bulge is consistent with this prediction, and is possibly the first direct detection of non-baryonic matter in the Galaxy.     

On the errors on Ω0: Monte Carlo simulations of the EMSS cluster sample

We perform Monte Carlo simulations of synthetic EMSS cluster samples, to quantify the systematic errors and the statistical uncertainties on the estimate of Ω₀ derived from fits to the cluster number density evolution and to the X-ray temperature distribution up to z =0.83 . We identify the scatter around the relation between cluster X-ray luminosity and temperature to be a source of systematic error, of the order of Δ_systΩ₀=0.09 , if not properly taken into account in the modelling. After correcting for this bias, our best Ω₀ is 0.66. The uncertainties on the shape and normalization of the power spectrum of matter fluctuations imply relatively large uncertainties on this estimate of Ω₀, of the order of Δ_statΩ₀=0.1 at the 1 σ level. On the other hand, the statistical uncertainties due to the finite size of the high-redshift sample are twice as small. Therefore, what is needed in order to improve the accuracy of Ω₀ estimates based on cluster number density evolution is a more reliable measure of the local temperature function and a better understanding of the cluster observed properties both in the local Universe and at high redshift, that is the relation between cluster mass, temperature and luminosity. This requires detailed observations of X-ray selected cluster samples, in comparison with hydrodynamic simulations including refined physics.     

An improved model for contraction of dark matter haloes in response to condensation of baryons

The cooling of gas in the centres of dark matter haloes is expected to lead to a more concentrated dark matter distribution. The response of dark matter to the condensation of baryons is usually calculated using the model of adiabatic contraction, which assumes spherical symmetry and circular orbits. Following Gnedin et al., we improve this model by modifying the assumed invariant from M ( r ) r to     , where r and     are the current and orbit-averaged particle positions. We explore the effect of the bulge in the inner regions of the halo for different values of the bulge-to-disc mass ratio. We find that the bulge makes the velocity curve rise faster in the inner regions of the halo. We present an analytical fitting curve that describes the velocity curve of the halo after dissipation. The results should be useful for dark matter detection studies.     

The formation of the Hubble sequence of disc galaxies: the effects of early viscous evolution

We investigate a model of disc galaxies whereby viscous evolution of the gaseous disc drives material inwards to form a protobulge. We start from the standard picture of disc formation through the settling of gas into a dark halo potential well, with the disc initially coming into centrifugal equilibrium with detailed conservation of angular momentum. We derive generic analytic solutions for the disc–halo system after adiabatic compression of the dark halo, with free choice of the input virialized dark halo density profile and of the specific angular momentum distribution. We derive limits on the final density profile of the halo in the central regions. Subsequent viscous evolution of the disc is modelled by a variation of the specific angular momentum distribution of the disc, providing analytic solutions to the final disc structure. The assumption that the viscous evolution time-scale and the star formation time-scale are similar leads to predictions of the properties of the stellar components. Focusing on small 'exponential' bulges, i.e., ones that may be formed through a disc instability, we investigate the relationship between the assumed initial conditions, such as halo 'formation', or assembly, redshift z _f, spin parameter λ , baryonic fraction F , and final disc properties such as global star formation time-scale, gas fraction, and bulge-to-disc ratio. We find that the present properties of discs, such as the scalelength, are compatible with a higher initial formation redshift if the redistribution by viscous evolution is included than if it is ignored. We also quantify the dependence of final disc properties on the ratio F λ , thus including the possibility that the baryonic fraction varies from galaxy to galaxy, as perhaps may be inferred from the observations.     

The coarse geometry of merger trees in Λ cold dark matter

We introduce the contour process to describe the geometrical properties of merger trees. The contour process translates a tree structure into a one-dimensional object: the contour walk. We characterize the contour walk measuring its length and action. The length is proportional to the number of progenitors in the tree, and the action is a proxy for the mean length of a branch in the tree.
We construct the contour walk for merger trees extracted from the public data base of the Millennium Run and also for merger trees constructed with a public Monte Carlo code which implements a Markovian algorithm. The trees correspond to haloes of final masses between 10¹¹ and  10¹⁴  h ⁻¹ M_⊙  . We study how the length and action scale with the mass of the final halo.
In all cases, except for the action measured from Markovian trees, we find a transitional scale around  1–6 × 10¹²  h ⁻¹ M_⊙  , which can be interpreted as the point where the halo growth starts to be dominated by mergers and not accretion. The transitional scale depends weakly on the mass resolution and time discretization of the merging trees.