Consequences for some dark energy candidates from the type Ia supernova SN 1997ff

We examine the status of various dark energy models in light of the recently observed SN 1997ff at   z ≈1.7  . The modified data still fit a pure cosmological constant Λ or a quintessence with an equation of state similar to that of Λ. The kinematical Λ models,  Λ∼ S ^-2  and  Λ∼ H ²  , also fit the data reasonably well and require less dark energy density (hence more matter energy density) than is required by the constant Λ model. However, the model  Λ∼ S ^-2  with low energy density becomes unphysical as it cannot accommodate higher redshift objects.
We also examine an alternative explanation of the data, namely the absorption by the intervening whisker-like dust, and find that the quasi-steady state (QSS) model and the Friedmann–Robertson–Walker (FRW) model  Ω_m0=0.33  without any dark energy also fit the data reasonably well.
We notice that the addition of SN 1997ff to the old data has worsened the fit to most of the models, except a closed FRW model with a constant Λ and a closed quintessence model with   ω _φ =-0.82  , and the models have started departing from each other as we go above   z =1  . However, to make a clear discrimination possible, a few more supernovae with   z >1  are required.
We have also calculated the age of the Universe in these models and find that, in the models with a constant Λ, the expansion age is uncomfortably close to the age of the globular clusters. Quintessence models show even lower age. The kinematical Λ models are, however, interesting in this connection (especially the model  Λ∼ H ²)  , as they give a remarkably large age of the Universe.     

A reanalysis of the X-ray luminosities of clusters of galaxies with 0.3<z<0.6 in the EMSS sample

The X-ray luminosities of the Einstein Extended Medium Sensitivity Survey (EMSS) clusters of galaxies with redshifts  0.3< z <0.6  are remeasured using ROSAT Position Sensitive Proportional Counter (PSPC) data. It is found that the new luminosities are on average  1.18±0.08  times higher than previously measured, but also that this ratio depends strongly on the X-ray core radii we measure. For the clusters with small core radii, in general we confirm the EMSS luminosities, but for clusters with core radii >250 kpc (the constant value assumed in the EMSS), the new luminosities are  2.2±0.15  times the previous measurements. The X-ray luminosity function (XLF) at  0.3< z <0.6  is recalculated and is found to be consistent with the local XLF. The constraints on the updated properties of the  0.3< z <0.6  EMSS sample, including a comparison with the number of clusters predicted from local XLFs, indicate that the space density of luminous, massive clusters either has not evolved or has increased by a small factor ∼2 since   z =0.4  . The implications of this result are discussed in terms of constraints on the cosmological parameter Ω₀.     

Cluster strong lensing in the Millennium simulation: the effect of galaxies and structures along the line-of-sight

We use ray-tracing through the Millennium simulation to study how secondary matter structures along the line-of-sight and the stellar mass in galaxies affect strong cluster lensing, in particular the cross-section for giant arcs. Furthermore, we investigate the distribution of the cluster Einstein radii and the radial distribution of giant arcs. We find that additional structures along the line-of-sight increase the strong-lensing optical depth by  ∼10–25 per cent  , while strong-lensing cross-sections of individual clusters are frequently boosted by as much as  ∼50 per cent  . The enhancement is mainly due to structures that are not correlated with the lens. Cluster galaxies increase the strong-lensing optical depth by up to a factor of 2, while interloping galaxies are not significant. We conclude that these effects need to be taken into account for predictions of the giant arc abundance, but they are not large enough to fully account for the reported discrepancy between predicted and observed abundances.
Furthermore, we find that Einstein radii defined via the area enclosed by the critical curve are 10–30 per cent larger than those defined via radial surface mass density profiles. The contributions of radial and tangential arcs to the radial distribution of arcs can be clearly distinguished. The radial distribution of tangential arcs is very broad and extends out to several Einstein radii. Thus, individual arcs are not well suited for constraining Einstein radii.     

Does the Fornax dwarf spheroidal have a central cusp or core?

The dark matter dominated Fornax dwarf spheroidal has five globular clusters orbiting at ∼1 kpc from its centre. In a cuspy cold dark matter halo the globulars would sink to the centre from their current positions within a few Gyr, presenting a puzzle as to why they survive undigested at the present epoch. We show that a solution to this timing problem is to adopt a cored dark matter halo. We use numerical simulations and analytic calculations to show that, under these conditions, the sinking time becomes many Hubble times; the globulars effectively stall at the dark matter core radius. We conclude that the Fornax dwarf spheroidal has a shallow inner density profile with a core radius constrained by the observed positions of its globular clusters. If the phase space density of the core is primordial then it implies a warm dark matter particle and gives an upper limit to its mass of ∼0.5 keV, consistent with that required to significantly alleviate the substructure problem.     

Primordial non-Gaussianities in the intergalactic medium

We present results from the first high-resolution hydrodynamical simulations of non-Gaussian cosmological models. We focus on the statistical properties of the transmitted Lyman-α flux in the high-redshift intergalactic medium. Imprints of non-Gaussianity are present and are larger at high redshifts. Differences larger than 20 per cent at   z > 3  in the flux probability distribution function for high-transmissivity regions (voids) are expected for values of the non-linearity parameter   f _NL=±100  when compared to a standard Λ cold dark matter cosmology with   f _NL= 0  . We also investigate the one-dimensional flux bispectrum: at the largest scales (corresponding to tens of Mpc), we expect deviations in the flux bispectrum up to 20 per cent at   z ∼ 4  (for   f _NL=±100  ), significantly larger than deviations of ∼3 per cent in the flux power spectrum. We briefly discuss possible systematic errors that can contaminate the signal. Although challenging, a detection of non-Gaussianities in the interesting regime of scales and redshifts probed by the Lyman-α forest could be possible with future data sets.     

Cluster contribution to the X-ray background as a cosmological probe

Extensive measurements of the X-ray background (XRB) yield a reasonably reliable characterization of its basic properties. Having resolved most of the cosmic XRB into discrete sources, the levels and spectral shapes of its main components can be used to probe both the source populations and also alternative cosmological and large-scale structure models. Recent observations of clusters seem to provide evidence that clusters formed earlier and are more abundant than predicted in the standard Λ cold dark matter (ΛCDM) model. This motivates interest in alternative models that predict enhanced power on cluster scales. We calculate predicted levels and spectra of the superposed emission from groups and clusters of galaxies in ΛCDM and in two viable alternative non-Gaussian  (χ²)  and early dark energy models. The predicted levels of the contribution of clusters to the XRB in the non-Gaussian models exceed the measured level at low energies and levels of the residual XRB in the 2–8 keV band; these particular models are essentially ruled out. Our work demonstrates the diagnostic value of the integrated X-ray emission from clusters, by considering also its dependences on different metallicities, gas and temperature profiles, Galactic absorption, merger scenarios and on a non-thermal pressure component. We also show that the XRB can be used for an upper limit for the concentration parameter value.     

Analytical approach to subhalo population in dark matter haloes

In the standard model of cosmic structure formation, dark matter haloes form by gravitational instability. The process is hierarchical: smaller systems collapse earlier, and later merge to form larger haloes. The galaxy clusters, hosted by the largest dark matter haloes, are at the top of this hierarchy and representing the largest as well as the last structures formed in the Universe, while the smaller and first haloes are those Earth-sized dark subhaloes that have been both predicted by theoretical considerations and found in numerical simulations, though there do not exist any observational hints of their existence. The probability that a halo of mass m at redshift z will be part of a larger halo of mass M at the present time can be described in the frame of the extended Press & Schecter theory making use of the progenitor (conditional) mass function. Using the progenitor mass function, we calculate analytically, at redshift zero, the distribution of subhaloes in mass, formation epoch and rarity of the peak of the density field at the formation epoch. That is done for a Milky Way size system, assuming both a spherical and an ellipsoidal collapse model. Our calculation assumes that small progenitors do not lose mass due to dynamical processes after entering the parent halo, and that they do not interact with other subhaloes. For a Λ cold dark matter power spectrum, we obtain a subhalo mass function  d n /d m   proportional to   m ^−α  with a model-independent  α∼ 2  . Assuming that the dark matter is a weakly interacting massive particle, the inferred distributions are used to test the feasibility of an indirect detection in the γ-ray energy band of such a population of subhaloes with a Gamma-ray Large Area Space Telescope like satellite.     

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.     

Dark matter subhaloes in numerical simulations

We use cosmological Λ cold dark matter (CDM) numerical simulations to model the evolution of the substructure population in 16 dark matter haloes with resolutions of up to seven million particles within the virial radius. The combined substructure circular velocity distribution function (VDF) for hosts of 10¹¹ to  10¹⁴ M_⊙  at redshifts from zero to two or higher has a self-similar shape, is independent of host halo mass and redshift, and follows the relation  d n /d v = (1/8)( v _cmax/ v _cmax,host)⁻⁴  . Halo to halo variance in the VDF is a factor of roughly 2 to 4. At high redshifts, we find preliminary evidence for fewer large substructure haloes (subhaloes). Specific angular momenta are significantly lower for subhaloes nearer the host halo centre where tidal stripping is more effective. The radial distribution of subhaloes is marginally consistent with the mass profile for   r ≳ 0.3 r _vir  , where the possibility of artificial numerical disruption of subhaloes can be most reliably excluded by our convergence study, although a subhalo distribution that is shallower than the mass profile is favoured. Subhalo masses but not circular velocities decrease towards the host centre. Subhalo velocity dispersions hint at a positive velocity bias at small radii. There is a weak bias towards more circular orbits at lower redshift, especially at small radii. We additionally model a cluster in several power-law cosmologies of   P ∝ kⁿ   , and demonstrate that a steeper spectral index, n , results in significantly less substructure.     

A new independent limit on the cosmological constant/dark energy from the relativistic bending of light by Galaxies and clusters of Galaxies

We derive new limits on the value of the cosmological constant, Λ, based on the Einstein bending of light by systems where the lens is a distant galaxy or a cluster of galaxies. We use an amended lens equation in which the contribution of Λ to the Einstein deflection angle is taken into account and use observations of Einstein radii around several lens systems. We use in our calculations a Schwarzschild–de Sitter vacuole exactly matched into a Friedmann–Robertson–Walker background and show that a Λ-contribution term appears in the deflection angle within the lens equation. We find that the contribution of the Λ-term to the bending angle is larger than the second-order term for many lens systems. Using these observations of bending angles, we derive new limits on the value of Λ. These limits constitute the best observational upper bound on Λ after cosmological constraints and are only two orders of magnitude away from the value determined by those cosmological constraints.     

Are fossil groups a challenge of the cold dark matter paradigm?

We study six groups and clusters of galaxies suggested in the literature to be 'fossil' systems (i.e. to have luminous diffuse X-ray emission and a magnitude gap of at least 2 mag R between the first and the second ranked member within half of the virial radius), each having good quality X-ray data and Sloan Digital Sky Survey (SDSS) spectroscopic or photometric coverage out to the virial radius. The poor cluster AWM 4 is clearly established as a fossil system, and we confirm the fossil nature of four other systems (RX J1331.5+1108, RX J1340.6+4018, RX J1256.0+2556 and RX J1416.4+2315), while the cluster RX J1552.2+2013 is disqualified as fossil system. For all systems, we present the luminosity functions within 0.5 and 1 virial radius that are consistent, within the uncertainties, with the universal luminosity function of clusters. For the five bona fide fossil systems, having a mass range  2 × 10¹³–3 × 10¹⁴ M_⊙  , we compute accurate cumulative substructure distribution functions (CSDFs) and compare them with the CSDFs of observed and simulated groups/clusters available in the literature. We demonstrate that the CSDFs of fossil systems are consistent with those of normal observed clusters and do not lack any substructure with respect to simulated galaxy systems in the cosmological Λ cold dark matter (ΛCDM) framework. In particular, this holds for the archetype fossil group RX J1340.6+4018 as well, contrary to earlier claims.     

Voids in a ΛCDM universe

We study the formation and evolution of voids in the dark matter distribution using various simulations of the popular Λ cold dark matter cosmogony. We identify voids by requiring them to be regions of space with a mean overdensity of −0.8 or less – roughly the equivalent of using a spherical overdensity group finder for haloes. Each of the simulations contains thousands of voids. The distribution of void sizes in the different simulations shows good agreement when differences in particle and grid resolution are accounted for. Voids very clearly correspond to minima in the smoothed initial density field. Apart from a very weak dependence on the mass resolution, the rescaled mass profiles of voids in the different simulations agree remarkably well. We find a universal void mass profile of the form  ρ(< r )/ρ( r _eff) ∝ exp[( r / r _eff)^α]  , where r _eff is the effective radius of a void and  α∼ 2  . The mass function of haloes in voids is steeper than that of haloes that populate denser regions. In addition, the abundances of void haloes seem to evolve somewhat more strongly between redshifts ∼1 and 0 than the global abundances of haloes.     

Tracing the nature of dark energy with galaxy distribution

Dynamical dark energy (DE) is a viable alternative to the cosmological constant. Constructing tests to discriminate between Λ and dynamical DE models is difficult, however, because the differences are not large. In this paper we explore tests based on the galaxy mass function, the void probability function (VPF), and the number of galaxy clusters. At high z , the number density of clusters shows large differences between DE models, but geometrical factors reduce the differences substantially. We find that detecting a model dependence in the cluster redshift distribution is a significant challenge. We show that the galaxy redshift distribution is potentially a more sensitive characteristic. We do this by populating dark matter haloes in N -body simulations with galaxies using well-tested halo occupation distributions. We also estimate the VPF and find that samples with the same angular surface density of galaxies, in different models, exhibition almost model-independent VPF which therefore cannot be used as a test for DE. Once again, geometry and cosmic evolution compensate each other. By comparing VPFs for samples with fixed galaxy mass limits, we find measurable differences.     

The elliptical galaxy colour–magnitude relation as a discriminant between the monolithic and merger paradigms

The colour–magnitude relation (CMR) of cluster elliptical galaxies has been widely used to constrain their star formation histories (SFHs) and to discriminate between the monolithic collapse and merger paradigms of elliptical galaxy formation. We use a Λ cold dark matter hierarchical merger model of galaxy formation to investigate the existence and redshift evolution of the elliptical galaxy CMR in the merger paradigm. We show that the SFH of cluster ellipticals predicted by the model is quasi-monolithic , with only ∼10 per cent of the total stellar mass forming after   z ∼ 1  . The quasi-monolithic SFH results in a predicted CMR that agrees well with its observed counterpart in the redshift range  0 < z < 1.27  . We use our analysis to argue that the elliptical-only CMR can be used to constrain the SFHs of present-day cluster ellipticals only if we believe a priori in the monolithic collapse model. It is not a meaningful tool for constraining the SFH in the merger paradigm, since a progressively larger fraction of the progenitor set of present-day cluster ellipticals is contained in late-type star-forming systems at higher redshift, which cannot be ignored when deriving the SFHs. Hence, the elliptical-only CMR is not a useful discriminant between the two competing theories of elliptical galaxy evolution.     

Simulations of deep pencil-beam redshift surveys

We create mock pencil-beam redshift surveys from very large cosmological N -body simulations of two cold dark matter (CDM) cosmogonies, an Einstein–de Sitter model ( τ CDM) and a flat model with Ω₀=0.3 and a cosmological constant (ΛCDM). We use these to assess the significance of the apparent periodicity discovered by Broadhurst et al. Simulation particles are tagged as 'galaxies' so as to reproduce observed present-day correlations. They are then identified along the past light-cones of hypothetical observers to create mock catalogues with the geometry and the distance distribution of the Broadhurst et al. data. We produce 1936 (2625) quasi-independent catalogues from our τ CDM (ΛCDM) simulation. A couple of large clumps in a catalogue can produce a high peak at low wavenumbers in the corresponding one-dimensional power spectrum, without any apparent large-scale periodicity in the original redshift histogram. Although the simulated redshift histograms frequently display regularly spaced clumps, the spacing of these clumps varies between catalogues and there is no 'preferred' period over our many realizations. We find only a 0.72 (0.49) per cent chance that the highest peak in the power spectrum of a τ CDM (ΛCDM) catalogue has a peak-to-noise ratio higher than that in the Broadhurst et al. data. None of the simulated catalogues with such high peaks shows coherently spaced clumps with a significance as high as that of the real data. We conclude that in CDM universes, the regularity on a scale of ∼130  h ⁻¹ Mpc observed by Broadhurst et al. has a priori probability well below 10⁻³.     

An explanation for long flares from extragalactic globular cluster X-ray sources

We compare orbits in a thin axisymmetric disc potential in Modified Newtonian Dynamics (MOND) with those in a thin disc plus near-spherical dark matter halo predicted by a ΛCDM cosmology. Remarkably, the amount of orbital precession in MOND is nearly identical to that which occurs in a mildly oblate CDM Galactic halo (potential flattening   q = 0.9  ), consistent with recent constraints from the Sagittarius stream. Since very flattened mass distributions in MOND produce rounder potentials than in standard Newtonian mechanics, we show that it will be very difficult to use the tidal debris from streams to distinguish between a MOND galaxy and a standard CDM galaxy with a mildly oblate halo.
If a galaxy can be found with either a prolate halo or one that is more oblate than   q ∼ 0.9  this would rule out MOND as a viable theory. Improved data from the leading arm of the Sagittarius dwarf – which samples the Galactic potential at large radii – could rule out MOND if the orbital pole precession can be determined to an accuracy of the order of  ±1°  .     

Supersonic motions of galaxies in clusters

We study motions of galaxies in galaxy clusters formed in the concordance Λ cold dark matter cosmology. We use high-resolution cosmological simulations that follow the dynamics of dark matter and gas and include various physical processes critical for galaxy formation: gas cooling, heating and star formation. Analysing the motions of galaxies and the properties of intracluster gas in a sample of eight simulated clusters at z = 0, we study the velocity dispersion profiles of the dark matter, gas and galaxies. We measure the mean velocity of galaxy motions and gas sound speed as a function of radius and calculate the average Mach number of galaxy motions. The simulations show that galaxies, on average, move supersonically with the average Mach number of ≈1.4, approximately independent of the cluster-centric radius. The supersonic motions of galaxies may potentially provide an important source of heating for the intracluster gas by driving weak shocks and via dynamical friction, although these heating processes appear to be inefficient in our simulations. We also find that galaxies move slightly faster than the dark matter particles. The magnitude of the velocity bias,   b _v ≈ 1.1  , is, however, smaller than the bias estimated for subhaloes in dissipationless simulations. Interestingly, we find velocity bias in the tangential component of the velocity dispersion, but not in the radial component. Finally, we find significant random bulk motions of gas. The typical gas velocities are of order ≈20–30 per cent of the gas sound speed. These random motions provide about 10 per cent of the total pressure support in our simulated clusters. The non-thermal pressure support, if neglected, will bias measurements of the total mass in the hydrostatic analyses of the X-ray cluster observations.     

Evolution of the luminosity function and colours of galaxies in a Λ cold dark matter universe

The luminosity function of galaxies is derived from a cosmological hydrodynamic simulation of a Λ cold dark matter universe with the aid of a stellar population synthesis model. At     , the resulting B -band luminosity function has a flat faint-end slope of     with the characteristic luminosity and the normalization in fair agreement with observations, while the dark matter halo mass function is steep with a slope of     . The colour distribution of galaxies also agrees well with local observations. We also discuss the evolution of the luminosity function, and the colour distribution of galaxies from     to 5. A large evolution of the characteristic mass in the stellar mass function as a result of number evolution is compensated by luminosity evolution; the characteristic luminosity increases only by 0.8 mag from     to 2, and then declines towards higher redshift, while the B -band luminosity density continues to increase from     to 5 (but only slowly at     .     

On the importance of high-redshift intergalactic voids

We investigate the properties of 1D flux 'voids' (connected regions in the flux distribution above the mean-flux level) by comparing hydrodynamical simulations of large cosmological volumes with a set of observed high-resolution spectra at z ∼ 2. After addressing the effects of box size and resolution, we study how the void distribution changes when the most significant cosmological and astrophysical parameters are varied. We find that the void distribution in the flux is in excellent agreement with predictions of the standard Λcold dark matter (ΛCDM) cosmology, which also fits other flux statistics remarkably well. We then model the relation between flux voids and the corresponding 1D gas-density field along the line of sight and make a preliminary attempt to connect the 1D properties of the gas-density field to the 3D dark matter distribution at the same redshift. This provides a framework that allows statistical interpretations of the void population at high redshift using observed quasar spectra, and eventually it will enable linking the void properties of the high-redshift universe with those at lower redshifts, which are better known.