Collapse of primordial clouds — II. The role of dark matter

In this article we extend the study performed in our previous article of the collapse of primordial objects. We here analyse the behaviour of the physical parameters for clouds ranging from 10⁷ to 10¹⁵ M_⊙. We study the dynamical evolution of these clouds in two ways: as purely baryonic clouds and as clouds with non-baryonic dark matter included. We start the calculations at the beginning of the recombination era, following the evolution of the structure until the collapse (which we defined as the time when the density contrast of the baryonic matter is greater than 10⁴). We analyse the behaviour of several physical parameters of the clouds (e.g. the density contrast and the velocities of the baryonic matter and the dark matter) as a function of time and radial position in the cloud. In this study all physical processes that are relevant to the dynamical evolution of the primordial clouds, such as for example photon drag (due to the cosmic background radiation) and hydrogen molecular production, besides the expansion of the Universe, are included in the calculations. In particular we find that the clouds with dark matter collapse at higher redshift when we compare the results with the purely baryonic models. As a general result we find that the distribution of the non-baryonic dark matter is more concentrated than the baryonic one. It is important to stress that we do not take into account the putative virialization of the non-baryonic dark matter; we just follow the time and spatial evolution of the cloud, solving its hydrodynamical equations. We also studied the role of cooling–heating processes in the purely baryonic clouds.     

A dark matter disc in three cosmological simulations of Milky Way mass galaxies

Making robust predictions for the phase-space distribution of dark matter at the solar neighbourhood is vital for dark matter direct-detection experiments. To date, almost all such predictions have been based on simulations that model the dark matter alone. Here, we use three cosmological hydrodynamic simulations of bright, disc-dominated galaxies to include the effects of baryonic matter self-consistently for the first time. We find that the addition of baryonic physics drastically alters the dark matter profile in the vicinity of the solar neighbourhood. A stellar/gas disc, already in place at high redshift, causes merging satellites to be dragged preferentially towards the disc plane where they are torn apart by tides. This results in an accreted dark matter disc that contributes ∼0.25–1.5 times the non-rotating halo density at the solar position. The dark disc, unlike dark matter streams, is an equilibrium structure that must exist in disc galaxies that form in a hierarchical cosmology. Its low rotation lag with respect to the Earth significantly boosts Weakly Interacting Massive Particle (WIMP) capture in the Earth and Sun, boosts the annual modulation signal and leads to distinct variations in the flux as a function of recoil energy that allow the WIMP mass to be determined.     

Analytic solutions for coupled linear perturbations

Analytic solutions for the evolution of cosmological linear density perturbations in the baryonic gas and collisionless dark matter are derived. The solutions are expressed in a closed form in terms of elementary functions, for arbitrary baryonic mass fraction. They are obtained assuming =1 and a time-independent comoving Jeans wavenumber, k _J. By working with a time variable ln( t ^2/3), the evolution of the perturbations is described by linear differential equations with constant coefficients. The new equations are then solved by means of Laplace transformation, assuming that the gas and dark matter trace the same density field before a sudden heating epoch. In a dark matter-dominated Universe, the ratio of baryonic to dark matter density perturbation decays with time roughly as exp(5 /4) t ^5/6 to the limiting value 1/[1+( k k _J)²]. For wavenumbers the decay is accompanied by oscillations with a period in . In comparison, as increases in a baryonic matter-dominated Universe, the ratio approaches 1( k k _J)² for k k _J, and zero otherwise.     

重子物质对暗物质晕形状和角动量的影响

利用高精度大样本的冷暗物质($\Lambda$ cold dark matter, $\Lambda$CDM)宇宙学数值模拟的数据, 对重子物质如何影响暗物质晕的形状和角动量进行了研究.使用了3种数值模拟数据, 纯暗物质模拟、含辐射冷却、恒星形成和动力学超新星反馈的模拟, 包含活动星系核反馈效应的恒星形成模拟. 对这3种模拟, 还进行了不同红移处的比较. 主要结果如下.重子物理过程会改变暗物质晕的质量分布, 特别是有活动星系核反馈机制的情况下.比如, 活动星系核反馈会减少大质量暗物质晕的形成.随着宇宙的演化, 暗物质晕的空间形态逐渐由扁变圆. 重子物质的存在会加速暗物质晕形状的变化过程, 而且会使暗物质晕形状变得更圆. 但是活动星系核的反馈会对这一加速效应产生抑制.重子物质对暗物质晕的影响与暗物质晕的质量和半径都存在一定的依赖性.暗物质晕的质量越大, 它会呈现更扁的形态. 同时, 重子物质对任意质量的暗物质晕或暗物质晕在任意半径处的变圆均有一定的促进作用,尽管活动星系核反馈会抑制这一促进作用.特别是对于暗物质晕在0.2--0.6倍维里半径处的形状, 重子物质的影响尤为明显.此外, 重子物质的存在会对暗物质晕的角动量产生显著影响, 它会增大暗物质的角动量. 暗物质晕的自旋参数与质量无相关性, 但是与暗物质晕的半径存在一定的相关性.     

The effect of cosmological background dynamics on the spherical collapse in MOND

The effect of background dynamics of the universe on formation of large scale structures in the framework of Modified Newtonian Dynamics (MOND) is investigated. A spherical collapse model is used for modeling the formation of the structures. This study is done in two extreme cases: (i) assuming a universe with a low-density baryonic matter without any cold dark matter and dark energy; (ii) a dark energy dominated universe with baryonic matter, without cold dark matter. We show that for the case (ii) the structures virialize at lower redshifts with larger radii compared to the low-density background universe. The dark energy slow downs the collapse of the structures. We show that our results are compatible with recent simulations of the structure formation in MOND.     

Re-ionizing the universe without stars

Recent observations show that the measured rates of star formation in the early universe are insufficient to produce re-ionization, and therefore, another source of ionizing photons is required. In this Letter, we examine the possibility that these can be supplied by the fast accretion shocks formed around the cores of the most massive haloes (10.5<log M/M _⊙<12) on spatial scales of order 1 kpc. We model the detailed physics of these fast accretion shocks, and apply these to a simple 1-D spherical hydrodynamic accretion model for baryonic infall in dark matter halos with an Einasto density distribution. The escape of UV photons from these halos is delayed by the time taken to reach the critical accretion shock velocity for escape of UV photons; 220 km s⁻¹, and by the time it takes for these photons to ionize the surrounding baryonic matter in the accretion flow. Assuming that in the universe at large the baryonic matter tracks the dark matter, we can estimate the epoch of re-ionization in the case that accretion shocks act alone as the source of UV photons. We find that 50% of the volume (and 5-8% of the mass) of the universe can be ionized by z∼7–8. The UV production rate has an uncertainty of a factor of about 5 due to uncertainties in the cosmological parameters controlling the development of large scale structure. Because our mechanism is a steeply rising function of decreasing redshift, this uncertainty translates to a re-ionization redshift uncertainty of less than ±0.5. We also find that, even without including the UV photon production of stars, re-ionization is essentially complete by z∼5.8. Thus, fast accretion shocks can provide an important additional source of ionizing photons in the early universe.     

Cuspy dark matter haloes and the Galaxy

The microlensing optical depth to Baade's Window constrains the minimum total mass in baryonic matter within the Solar circle to be greater than ∼     , assuming the inner Galaxy is barred with viewing angle ∼20°. From the kinematics of solar neighbourhood stars, the local surface density of dark matter is ∼     . We construct cuspy haloes normalized to the local dark matter density and calculate the circular-speed curve of the halo in the inner Galaxy. This is added in quadrature to the rotation curve provided by the stellar and ISM discs, together with a bar sufficiently massive so that the baryonic matter in the inner Galaxy reproduces the microlensing optical depth. Such models violate the observational constraint provided by the tangent-velocity data in the inner Galaxy (typically at radii     . The high baryonic contribution required by the microlensing is consistent with implications from hydrodynamical modelling and the pattern speed of the Galactic bar. We conclude that the cuspy haloes favoured by the cold dark matter cosmology (and its variants) are inconsistent with the observational data on the Galaxy.     

Do the Unidentified Egret Sources Trace Annihilating Dark Matter in the Local Group?

In a cold dark matter (CDM) framework of structure formation, the dark matter haloes around galaxies assemble through successive mergers with smaller haloes. This merging process is not completely efficient, and hundreds of surviving halo cores, or subhaloes, are expected to remain in orbit within the halo of a galaxy like the Milky Way. While the dozen visible satellites of the Milky Way may trace some of these subhaloes, the majority are currently undetected. A large number of high-velocity clouds (HVCs) of neutral hydrogen are observed around the Milky Way, and it is plausible that some of the HVCs may trace subhaloes undetected in the optical. Confirming the existence of concentrations of dark matter associated with even a few of the HVCs would represent a dramatic step forward in our attempts to understand the nature of dark matter. Supersymmetric (SUSY) extensions of the Standard Model of particle physics currently suggest neutralinos as a natural well-motivated candidate for the non-baryonic dark matter of the universe. If this is indeed the case, then it may be possible to detect dark matter indirectly as it annihilates into neutrinos, photons or positrons. In particular, the centres of subhaloes might show up as point sources in gamma-ray observations. In this work, we consider the possibility that some of the unidentified EGRET γ-ray sources trace annihilating neutralino dark matter in the dark substructure of the Local Group. We compare the observed positions and fluxes of both the unidentified EGRET sources and the HVCs with the positions and fluxes predicted by a model of halo substructure, to determine up to what extent any of these three populations could be associated.     

Double starbursts triggered by mergers in hierarchical clustering scenarios

We use cosmological smooth particle hydrodynamical (SPH) simulations to study the effects of mergers in the star formation history of galactic objects in hierarchical clustering scenarios. We find that during some merger events, gaseous discs can experience two starbursts: the first one during the orbital decay phase, owing to gas inflows driven as the satellite approaches, and the second one when the two baryonic clumps collide. A trend for these first induced starbursts to be more efficient at transforming the gas into stars is also found. We detect that systems that do not experience early gas inflows have well-formed stellar bulges and more concentrated potential wells, which seem to be responsible for preventing further gas inward transport triggered by tidal forces. The potential wells concentrate owing to the accumulation of baryons in the central regions and of dark matter as the result of the pulling in by baryons. The coupled evolution of the dark matter and baryons would lead to an evolutionary sequence during which systems with shallower total potential wells suffer early gas inflows during the orbital decay phase that help to feed their central mass concentration, pulling in dark matter and contributing to build up more stable systems. Within this scenario, starbursts triggered by early gas inflows are more likely to occur at early stages of evolution of the systems and to be an important contributor to the formation of stellar bulges. Our results constitute the first proof that bulges can form as the product of collapse, collisions and secular evolution in a cosmological framework, and they are consistent with a rejuvenation of the stellar population in bulges at intermediate z with, at least, 50 per cent of the stars (in SCDM) being formed at high z .     

Maximum feedback and dark matter profiles of dwarf galaxies

The observed rotation curves of dark matter-dominated dwarf galaxies indicate low-density cores, contrary to the predictions of CDM models. A possible solution of this problem involves stellar feedback. A strong baryonic wind driven by vigorous star formation can remove a large fraction of the gas, causing the dark matter to expand. Using both numerical and analytical techniques, we explore the maximum effect of the feedback with an instantaneous removal of the gaseous disc. The energy input depends on the compactness of the disc, hence the specific angular momentum of the disc. For the plausible cosmological parameters and a wide range of the disc angular momenta, the feedback is insufficient to destroy the central halo cusp, while the inner density is lowered only by a modest factor of 2 to 6. Any realistic modelling of the feedback would have even lesser impact on dark matter. We find that no star formation effect can resolve the problems of CDM cusps.     

Populating a cluster of galaxies – I. Results at

We simulate the assembly of a massive rich cluster and the formation of its constituent galaxies in a flat, low-density universe. Our most accurate model follows the collapse, the star formation history and the orbital motion of all galaxies more luminous than the Fornax dwarf spheroidal, while dark halo structure is tracked consistently throughout the cluster for all galaxies more luminous than the SMC. Within its virial radius this model contains about     dark matter particles and almost 5000 distinct dynamically resolved galaxies. Simulations of this same cluster at a variety of resolutions allow us to check explicitly for numerical convergence both of the dark matter structures produced by our new parallel N -body and substructure identification codes, and of the galaxy populations produced by the phenomenological models we use to follow cooling, star formation, feedback and stellar aging. This baryonic modelling is tuned so that our simulations reproduce the observed properties of isolated spirals outside clusters. Without further parameter adjustment our simulations then produce a luminosity function, a mass-to-light ratio, luminosity, number and velocity dispersion profiles, and a morphology–radius relation which are similar to those observed in real clusters. In particular, since our simulations follow galaxy merging explicitly, we can demonstrate that it accounts quantitatively for the observed cluster population of bulges and elliptical galaxies.     

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.     

The post-collapse equilibrium structure of cosmological haloes in a low-density universe

An analytical model is presented for the post-collapse equilibrium structure of virialized objects that condense out of a low-density cosmological background universe, either matter-dominated or flat with a cosmological constant. This generalizes the model we derived previously for an Einstein–de Sitter (EdS) universe. The model is based upon the assumption that cosmological haloes form from the collapse and virialization of 'top-hat' density perturbations, and are spherical, isotropic and isothermal. This leads to the prediction of a unique, non-singular, truncated isothermal sphere (TIS), a particular solution of the Lane–Emden equation (suitably modified when Λ≠0) . The size and virial temperature are unique functions of the mass and redshift of formation of the object for a given background universe. The central density is roughly proportional to the critical density of the universe at the epoch of collapse. This TIS model is in good agreement with observations of the internal structure of dark-matter-dominated haloes on scales ranging from dwarf galaxies to X-ray clusters. It also reproduces many of the average properties of haloes in simulations of the cold dark matter (CDM) model to good accuracy, suggesting that it is a useful analytical approximation for haloes that form from realistic initial conditions. Our TIS model matches the density profiles of haloes in CDM N -body simulations outside the innermost region, while avoiding the steep central cusp of the latter which is in apparent conflict with observations. The TIS model may also be relevant to non-standard CDM models, such as that for self-interacting dark matter, recently proposed to resolve this conflict.     

The angular momentum content of dwarf galaxies: new challenges for the theory of galaxy formation

We compute the specific angular momentum distributions for a sample of low-mass disc galaxies observed by Swaters. We compare these distributions to those of dark matter haloes obtained by Bullock et al. from high-resolution N -body simulations of structure formation in a ΛCDM universe. We find that although the disc mass fractions are significantly smaller than the universal baryon fraction, the total specific angular momenta of the discs are in good agreement with those of dark matter haloes. This suggests that discs form out of only a small fraction of the available baryons, but yet manage to draw most of the available angular momentum. In addition we find that the angular momentum distributions of discs are clearly distinct from those of the dark matter; discs lack predominantly both low and high specific angular momenta. Understanding these findings in terms of a coherent picture for disc formation is challenging. Cooling, feedback and stripping, which are the main mechanisms to explain the small disc mass fractions found, seem unable to simultaneously explain the angular momentum distributions of the discs. In fact, it seems that the baryons that make up the discs must have been born out of angular momentum distributions that are clearly distinct from those of ΛCDM haloes. However, the dark and baryonic mass components experience the same tidal forces, and it is therefore expected that they should have similar angular momentum distributions. Therefore, understanding the angular momentum content of disc galaxies remains an important challenge for our picture of galaxy formation.     

Constraining the cosmological baryon density with X-ray clusters

We study the properties of X-ray galaxy clusters in four cold dark matter models with different baryon fractions Ω_BM, ranging from 5 to 20 per cent. By using an original three-dimensional hydrodynamic code based on the piecewise parabolic method, we run simulations on a box with a size of 64  h ⁻¹ Mpc and we identify the clusters by selecting the peaks in the X-ray luminosity field. We analyse these mock catalogues by computing the mass function, the luminosity function, the temperature distribution and the luminosity–temperature relation. By comparing the predictions of the different models to a series of recent observational results, we find that only the models with low baryonic content agree with the data, while models with larger baryon fraction are well outside the 1σ error bars. In particular, the analysis of the luminosity functions, both bolometric and in the energy band [0.5–2] keV, requires Ω_BM ≲ 0.05 when we fix the values h  = 0.5 and n  = 0.8 for the Hubble parameter and the primordial spectral index, respectively. Moreover we find that, independently of the cosmological scenario, all the considered quantities have very little redshift evolution, particularly between z  = 0.5 and 0.     

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.     

The statistics of wide-separation lensed quasars

The absence of any wide-separation gravitational lenses in the Large Bright Quasar Survey is used to place limits on the population of cluster-sized haloes in the universe, and hence constrain a number of cosmological parameters. The results agree with previous investigations in strongly ruling out the standard cold dark matter model but they are consistent with low-density universes in which the primordial fluctuation spectrum matches both cluster abundances and cosmic microwave background measurements. These conclusions are essentially independent of the cosmological constant, which is in stark contrast to the statistics of galaxy lenses. The constraints presented here are nullified if clusters have core radii of ≳10 kpc, but are free of a number of potential systematic errors, owing to the homogeneity of the data.     

The Bullet Cluster 1E0657-558 evidence shows modified gravity in the absence of dark matter

A detailed analysis of the 2006 November 15 data release X-ray surface density Σ-map and the strong and weak gravitational lensing convergence κ-map for the Bullet Cluster 1E0657-558 is performed and the results are compared with the predictions of a modified gravity (MOG) and dark matter. Our surface density Σ-model is computed using a King β-model density, and a mass profile of the main cluster and an isothermal temperature profile are determined by the MOG. We find that the main cluster thermal profile is nearly isothermal. The MOG prediction of the isothermal temperature of the main cluster is   T = 15.5 ± 3.9 keV  , in good agreement with the experimental value   T = 14.8^+2.0_−1.7 keV  . Excellent fits to the 2D convergence κ-map data are obtained without non-baryonic dark matter, accounting for the 8σ spatial offset between the Σ-map and the κ-map reported in Clowe et al. The MOG prediction for the κ-map results in two baryonic components distributed across the Bullet Cluster 1E0657-558 with averaged mass fraction of 83 per cent intracluster medium (ICM) gas and 17 per cent galaxies. Conversely, the Newtonian dark matter κ-model has on average 76 per cent dark matter (neglecting the indeterminant contribution due to the galaxies) and 24 per cent ICM gas for a baryon to dark matter mass fraction of 0.32, a statistically significant result when compared to the predicted Λ-cold dark matter cosmological baryon mass fraction of 0.176^+0.019_−0.012.     

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.     

The orientation of galaxy dark matter haloes around cosmic voids

Using the Millennium N -body Simulation we explore how the shape and angular momentum of galaxy dark matter haloes surrounding the largest cosmological voids are oriented. We find that the major and intermediate axes of the haloes tend to lie parallel to the surface of the voids, whereas the minor axis points preferentially in the radial direction. We have quantified the strength of these alignments at different radial distances from the void centres. The effect of these orientations is still detected at distances as large as 2.2 R _void from the void centre. Taking a subsample of haloes expected to contain disc-dominated galaxies at their centres we detect, at the 99.9 per cent confidence level, a signal that the angular momentum of those haloes tends to lie parallel to the surface of the voids. Contrary to the alignments of the inertia axes, this signal is only detected in shells at the void surface  (1 < R < 1.07  R _void)  and disappears at larger distances. This signal, together with the similar alignment observed using real spiral galaxies, strongly supports the prediction of the Tidal Torque theory that both dark matter haloes and baryonic matter have acquired, conjointly, their angular momentum before the moment of turnaround.