Dark matter halo response to the disc growth

We consider the sensitivity of the circular-orbit adiabatic contraction approximation to the baryon condensation rate and the orbital structure of dark matter haloes in the Λ cold dark matter (ΛCDM) paradigm. Using one-dimensional hydrodynamic simulations including the dark matter halo mass accretion history and gas cooling, we demonstrate that the adiabatic approximation is approximately valid even though haloes and discs may assemble simultaneously. We further demonstrate the validity of the simple approximation for ΛCDM haloes with isotropic velocity distributions using three-dimensional N -body simulations. This result is easily understood: an isotropic velocity distribution in a cuspy halo requires more circular orbits than radial orbits. Conversely, the approximation is poor in the extreme case of a radial orbit halo. It overestimates the response of a core dark matter halo, where radial orbit fraction is larger. Because no astronomically relevant models are dominated by low angular momentum orbits in the vicinity of the disc and the growth time-scale is never shorter than a dynamical time, we conclude that the adiabatic contraction approximation is useful in modelling the response of dark matter haloes to the growth of a disc.     

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.     

Dwarf galaxy rotation curves and the core problem of dark matter haloes

The standard cold dark matter (CDM) model has recently been challenged by the claim that dwarf galaxies have dark matter haloes with constant-density cores, whereas CDM predicts haloes with steeply cusped density distributions. Consequently, numerous alternative dark matter candidates have recently been proposed. In this paper we scrutinize the observational evidence for the incongruity between dwarf galaxies and the CDM model. To this end, we analyse the rotation curves of 20 late-type dwarf galaxies studied by Swaters. Taking the effects of beam smearing and adiabatic contraction into account, we fit mass models to these rotation curves with dark matter haloes with different cusp slopes, ranging from constant-density cores to r ⁻² cusps. Even though the effects of beam smearing are small for these data, the uncertainties in the stellar mass-to-light ratio and the limited spatial sampling of the halo's density distribution hamper a unique mass decomposition. Consequently, the rotation curves in our sample cannot be used to discriminate between dark haloes with constant-density cores and r ⁻¹ cusps. We show that the dwarf galaxies analysed here are consistent with CDM haloes in a ΛCDM cosmology, and that there is thus no need to abandon the idea that dark matter is cold and collisionless. However, the data are also consistent with any alternative dark matter model that produces dark matter haloes with central cusps less steep than r ^−1.5. In fact, we argue that based on existing H  i rotation curves alone, at best weak limits can be obtained on cosmological parameters and/or the nature of the dark matter. In order to make progress, rotation curves with higher spatial resolution and independent measurements of the mass-to-light ratio of the disc are required.     

The formation and survival of discs in a ΛCDM universe

We study the formation of galaxies in a Λ cold dark matter (ΛCDM) universe using high-resolution hydrodynamical simulations with a multiphase treatment of gas, cooling and feedback, focusing on the formation of discs. Our simulations follow eight isolated haloes similar in mass to the Milky Way and extracted from a large cosmological simulation without restriction on spin parameter or merger history. This allows us to investigate how the final properties of the simulated galaxies correlate with the formation histories of their haloes. We find that, at   z = 0  , none of our galaxies contains a disc with more than 20 per cent of its total stellar mass. Four of the eight galaxies nevertheless have well-formed disc components, three have dominant spheroids and very small discs, and one is a spheroidal galaxy with no disc at all. The   z = 0  spheroids are made of old stars, while discs are younger and formed from the inside-out. Neither the existence of a disc at   z = 0  nor the final disc-to-total mass ratio seems to depend on the spin parameter of the halo. Discs are formed in haloes with spin parameters as low as 0.01 and as high as 0.05; galaxies with little or no disc component span the same range in spin parameter. Except for one of the simulated galaxies, all have significant discs at   z ≳ 2  , regardless of their   z = 0  morphologies. Major mergers and instabilities which arise when accreting cold gas is misaligned with the stellar disc trigger a transfer of mass from the discs to the spheroids. In some cases, discs are destroyed, while in others, they survive or reform. This suggests that the survival probability of discs depends on the particular formation history of each galaxy. A realistic ΛCDM model will clearly require weaker star formation at high redshift and later disc assembly than occurs in our models.     

Star formation and chemical evolution of damped Lyman α systems

In this paper, we investigate the star formation and chemical evolution of damped Lyman α systems (DLAs) based on the disc galaxy formation model developed by Mo, Mao & White. We propose that the DLAs are the central galaxies of less-massive dark haloes present at redshifts z ∼3, and they should inhabit haloes of moderately low circular velocity. The empirical Schmidt law of star formation rates, and closed box model of chemical evolution that an approximation known as instantaneous recycling is assumed, are adopted. In our models, when the predicted distribution of metallicity for DLAs is calculated, two cases are considered. One is that, using the closed-box model, empirical Schmidt law and star formation time, the distribution of metallicity can be directly calculated. The other is that, when the simple gravitational instability of a thin isothermal gas disc as first discussed by Toomre is considered, the star formation occurs only in the region where the surface density of gas satisfies the critical value, not everywhere of a gas disc. In this case, we first obtain the region where the star formation can occur by assuming that the disc has a flat rotation curve and rotational velocity is equal to the circular velocity of the surrounding dark matter halo, and then calculate the metallicity distribution as in case one. We assume that star formation in each DLA lasts for a period of 1 Gyr from redshifts z =3. There is only one output parameter in our models, i.e. the stellar yield, which relates to the time of star formation history and is obtained by normalizing the predicted distribution of metallicity to the mean value of 1/13 Z_⊙ as presented by Pettini et al.. The predicted metallicity distribution is consistent with the current (rather limited) observational data. A random distribution of galactic discs is taken into account.     

The cosmological dependence of galactic specific angular momenta

Hydrodynamical simulations of galaxy formation in spatially flat cold dark matter (CDM) cosmologies with and without a cosmological constant (Λ) are described. A simple star formation algorithm is employed and radiative cooling is allowed only after redshift z =1 so that enough hot gas is available to form large, rapidly rotating stellar discs if angular momentum is approximately conserved during collapse. The specific angular momenta of the final galaxies are found to be sensitive to the assumed background cosmology. This dependence arises from the different angular momenta contained in the haloes at the epoch when the gas begins to collapse and the inhomogeneity of the subsequent halo evolution. In the Λ-dominated cosmology, the ratio of stellar specific angular momentum to that of the dark matter halo (measured at the virial radius) has a median value of ∼0.24 at z =0. The corresponding quantity for the Λ=0 cosmology is over three times lower. It is concluded that the observed frequency and angular momenta of disc galaxies pose significant problems for spatially flat CDM models with Λ=0 but may be consistent with a Λ-dominated CDM universe.     

Tidal tails in cold dark matter cosmologies

We study the formation of tidal tails in pairs of merging disc galaxies with structural properties motivated by current theories of cold dark matter (CDM) cosmologies. In a recent study, Dubinski, Mihos & Hernquist showed that the formation of prominent tidal tails can be strongly suppressed by massive and extended dark haloes. For the large halo-to-disc mass ratio expected in CDM cosmologies their sequence of models failed to produce strong tails like those observed in many well-known pairs of interacting galaxies. In order to test whether this effect can constrain the viability of CDM cosmologies, we construct N ‐body models of disc galaxies with structural properties derived in analogy to the recent analytical work of Mo, Mao & White. With a series of self-consistent collisionless simulations of galaxy–galaxy mergers we demonstrate that even the discs of very massive dark haloes have no problems developing long tidal tails, provided the halo spin parameter is large enough. For our class of models, the halo-to-disc mass ratio is not a good indicator of the ability to produce tails. Instead, the relative size of disc and halo or, alternatively, the ratio of circular velocity to local escape speed at the half mass radius of the disc is a more useful criterion. This result holds in all CDM models. While tidal tails can provide useful information on the structure of galaxies, it thus appears unlikely that they are able to constrain the values of the cosmological parameters within these models.     

The properties of Lyα emitting galaxies in hierarchical galaxy formation models

We present detailed predictions for the properties of Lyα-emitting galaxies in the framework of the Λ cold dark matter cosmology, calculated using the semi-analytical galaxy formation model galform . We explore a model that assumes a top-heavy initial mass function in starbursts and that has previously been shown to explain the sub-millimetre number counts and the luminosity function of Lyman-break galaxies at high redshift. We show that this model, with the simple assumption that a fixed fraction of Lyα photons escape from each galaxy, is remarkably successful at explaining the observed luminosity function of Lyα emitters (LAEs) over the redshift range  3 < z < 6.6  . We also examine the distribution of Lyα equivalent widths and the broad-band continuum magnitudes of emitters, which are in good agreement with the available observations. We look more deeply into the nature of LAEs, presenting predictions for fundamental properties such as the stellar mass and radius of the emitting galaxy and the mass of the host dark matter halo. The model predicts that the clustering of LAEs at high redshifts should be strongly biased relative to the dark matter, in agreement with observational estimates. We also present predictions for the luminosity function of LAEs at   z > 7  , a redshift range that is starting to be be probed by near-infrared surveys and using new instruments such as the Dark Ages Z Lyman Explorer (DAzLE).     

Galaxy formation and evolution: low-surface-brightness galaxies

We investigate in detail the hypothesis that low-surface-brightness galaxies (LSBs) differ from ordinary galaxies simply because they form in haloes with large spin parameters. We compute star formation rates using the Schmidt law, assuming the same gas infall dependence on surface density as used in models of the Milky Way. We build stellar population models, predicting colours, spectra and chemical abundances. We compare our predictions with observed values of metallicity and colours for LSBs, and find excellent agreement with all observables. In particular, integrated colours, colour gradients, surface brightness and metallicity match very well to the observed values of LSBs for models with ages larger than 7 Gyr and high values (λ > 0.05) for the spin parameter of the haloes. We also compute the global star formation rate (SFR) in the Universe due to LSBs, and show that it has a flatter evolution with redshift than the corresponding SFR for normal discs. We furthermore compare the evolution in redshift of [ Zn / H ] for our models to those observed in damped Lyman α systems by Pettini et al. and show that damped Lyman α system abundances are consistent with the predicted abundances at different radii for LSBs. Finally, we show how the required late redshift of collapse of the halo may constrain the power spectrum of fluctuations.     

The association between gas and galaxies – I. CFHT spectroscopy and pair analysis

We investigate the relative distribution of the gaseous contents of the Universe (as traced by a sample of Lyα absorbers), and the luminous baryonic matter (as traced by a redshift survey of galaxies in the same volume searched for Lyα absorbers), along 16 lines of sight (LOS) between redshifts 0 and 1. Our galaxy redshift survey was made with the multi-object spectrograph on the Canada–France–Hawaii Telescope and, when combined with galaxies from the literature in the same LOS, gives us a galaxy sample of 636 objects. By combining this with an absorption-line sample of 406 absorbing systems drawn from published works, we are able to study the relationship between gas and galaxies over the latter half of the age of the Universe. A correlation between absorbers and galaxies is detected out to separation of 1.5 Mpc. This correlation is weaker than the galaxy–galaxy correlation. There is also some evidence that the absorbing systems seen in C  iv are more closely related to galaxies, although this correlation could be with column density rather than metallicity. The above results are all consistent with the absorbing gas and the galaxies coexisting in dark matter filaments and knots as predicted by current models where the column density of the absorbing gas is correlated with the underlying matter density.     

On the nature of superoutbursts in dwarf novae

We determine a crucial feature of the dark halo density distribution from the fact that the luminous matter dominates the gravitational potential at about one disc scalelength R _d, but at the optical edge     the dark matter has already become the main component of the galaxy density. From the kinematics of 137 spirals we find that the dark matter halo density profiles are self-similar at least out to R _opt and show core radii much larger than the corresponding disc scalelengths. The luminous regions of spirals consist of stellar discs embedded in dark haloes with roughly constant density. This invariant dark matter profile is very difficult to reconcile with the fundamental properties of the density distribution of cold dark matter haloes. With respect to previous work, the present evidence is obtained by means of a robust method and for a large and complete sample of normal spirals.     

Damped Lyα absorbers from dwarf galaxy ejecta

It is argued that the formation of a dwarf galaxy causes a massive burst of star formation, resulting in the ejection of most of the available gas from the galaxy as a weakly collimated wind. The ejected gas can give rise to a damped Lyα absorber (DLA). Weakly collimated outflows naturally explain the asymmetric profiles seen in low-ionization absorption lines caused by heavy elements associated with DLAs, where absorption is strongest at one edge of the absorption feature. The shape of the distribution of column densities in the model agrees reasonably well with observations. In particular, the break in slope is caused by external photoionization of the wind. A semi-analytical model for galaxy formation is used to show that, for currently acceptable cosmological parameters, dwarf galaxy outflows can account for the majority of DLA systems and their distribution with redshift. This model also predicts a correlation between velocity structure and metallicity of DLA systems, in qualitative agreement with observations. DLAs do not require many large, rapidly rotating disc galaxies to have formed early on, as in other models for their origin.     

Galaxy cluster mass profiles

Accurate measurements of the mass distribution in galaxy and cluster haloes are essential to test the cold dark matter (CDM) paradigm. The cosmological model predicts a universal shape for the density profile in all haloes, independent of halo mass. Its profile has a 'cuspy' centre, with no evidence for the constant density core. In this paper, we carry out a careful analysis of 12 galaxy clusters, using Chandra data to compute the mass distribution in each system under the assumption of hydrostatic equilibrium. Due to their low concentration, clusters provide ideal objects for studying the central cusps in dark matter haloes. The majority of the systems are consistent with the CDM model, but four objects exhibit flat inner density profiles. We suggest that the flat inner profile found for these clusters is due to an underestimation of the mass in the cluster centre (rather than any problem with the CDM model), since these objects also have a centrally peaked gas mass fraction. We discuss possible causes for erroneously low-mass measurements in the cores of some systems.     

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.     

Metallicity in damped Lyman-α systems: evolution or bias?

Assuming that damped Lyman- α (DLA) systems are galactic discs, we calculate the corresponding evolution of metal abundances. We use detailed multizone models of galactic chemical evolution (reproducing successfully the observed properties of disc galaxies) and appropriate statistics (including geometrical probability factors) to calculate the average metallicity as a function of redshift. The results are compatible with available observations, provided that observational biases are taken into account , as suggested by Boissé et al. In particular, high column density and high-metallicity systems are not detected because the light of backround quasars is severely extinguished, while low column density and low-metallicity systems are not detectable through their absorption lines by current surveys. We show that these observational constraints lead to a 'no-evolution' picture for the DLA metallicity, which does not allow us to draw strong conclusions about the nature of those systems or about their role in 'cosmic chemical evolution'.     

The origin of the density distribution of disc galaxies: a new problem for the standard model of disc formation

We present new models for the formation of disc galaxies that improve upon previous models by following the detailed accretion and cooling of the baryonic mass, and by using realistic distributions of specific angular momentum. Under the assumption of detailed angular momentum conservation, the discs that form have density distributions that are more centrally concentrated than an exponential. We examine the influence of star formation, bulge formation, and feedback on the outcome of the surface brightness distributions of the stars. Low angular momentum haloes yield disc galaxies with a significant bulge component and with a stellar disc that is close to exponential, in good agreement with observations. High angular momentum haloes, on the other hand, produce stellar discs that are much more concentrated than an exponential, in clear conflict with observations. At large radii, the models reveal distinct truncation radii in both the stars and the cold gas. The stellar truncation radii result from our implementation of star formation threshold densities, and are in excellent agreement with observations. The truncation radii in the density distribution of the cold gas reflect the maximum specific angular momentum of the gas that has cooled. We find that these truncation radii occur at H  i surface densities of roughly 1 M_⊙ pc⁻², in conflict with observations. We examine various modifications to our models, including feedback, viscosity, and dark matter haloes with constant-density cores, but show that the models consistently fail to produce bulge less discs with exponential surface brightness profiles. This signals a new problem for the standard model of disc formation: if the baryonic component of the protogalaxies out of which disc galaxies form has the same angular momentum distribution as the dark matter, discs are too compact.     

The extended rotation curve and the dark matter halo of M33

We present the 21-cm rotation curve of the nearby galaxy M33 out to a galactocentric distance of 16 kpc (13 disc scalelengths). The rotation curve keeps rising out to the last measured point and implies a dark halo mass ≳5×10¹⁰ M_⊙. The stellar and gaseous discs provide virtually equal contributions to the galaxy gravitational potential at large galactocentric radii, but no obvious correlation is found between the radial distribution of dark matter and the distribution of stars or gas.
Results of the best fit to the mass distribution in M33 picture a dark halo which controls the gravitational potential from 3 kpc outward, with a matter density which decreases radially as R ^−1.3. The density profile is consistent with the theoretical predictions for structure formation in hierarchical clustering cold dark matter (CDM) models, and favours lower mass concentrations than those expected in the standard cosmogony.     

Lyman-break galaxies: high mass or low?

We reassess the hypothesis that Lyman-break galaxies (LBGs) at redshifts   z ∼ 3  mark the centres of the most massive dark matter haloes at that epoch. First we reanalyse the kinematic measurements of Pettini et al. and Erb et al. of the rest-frame optical emission lines of LBGs. We compare the distribution of the ratio of the rotation velocity to the central line width, against the expected distribution for galaxies with random inclination angles, modelled as singular isothermal spheres. The model fits the data well. On this basis we argue that the central line width provides a predictor of the circular velocity at a radius of several kpc. Assembling a larger sample of LBGs with measured line widths, we compare these results against the theoretical ΛCDM rotation curves of Mo, Mao & White, under the hypothesis that LBGs mark the centres of the most massive dark matter haloes. We find that the circular velocities are overpredicted by a substantial factor, which we estimate conservatively as  1.8 ± 0.4  . This indicates that the model is probably incorrect. The model of LBGs as relatively low-mass starburst systems, of Somerville, Primack & Faber, provides a good fit to the data.     

Vertical dynamics of disc galaxies in modified Newtonian dynamics

We investigate the possibility of discriminating between modified Newtonian dynamics (MOND) and Newtonian gravity with dark matter, by studying the vertical dynamics of disc galaxies. We consider models with the same circular velocity in the equatorial plane (purely baryonic discs in MOND and the same discs in Newtonian gravity embedded in spherical dark matter haloes), and we construct their intrinsic and projected kinematical fields by solving the Jeans equations under the assumption of a two-integral distribution function. We find that the vertical velocity dispersion of deep MOND discs can be much larger than in the equivalent spherical Newtonian models. However, in the more realistic case of high surface density discs, this effect is significantly reduced, casting doubt on the possibility of discriminating between MOND and Newtonian gravity with dark matter by using current observations.     

Spatial fluctuations in the spectral shape of the ultraviolet background at 2 < z < 3 and the reionization of helium

The low-density hydrogen and helium in the intergalactic medium (IGM) probed by quasi-stellar object (QSO) absorption lines is sensitive to the amplitude and spectral shape of the metagalactic ultraviolet (UV) background. We use realistic H  i and He  ii Lyα forest spectra, constructed from state-of-the-art hydrodynamical simulations of a Λ cold dark matter (ΛCDM) universe to confirm the reliability of using line profile fitting techniques to infer the ratio of the metagalactic H  i and He  ii ionization rates. We further show that the large spatial variations and the anticorrelation with H  i absorber density observed in the ratio of the measured He  ii to H  i column densities can be explained in a model where the H  i ionization rate is dominated by the combined UV emission from young star-forming galaxies and QSOs and the He  ii ionization rate is dominated by emission from QSOs only. In such a model the large fluctuations in the column density ratio are due to the small number of QSOs expected to contribute at any given point to the He  ii ionization rate. A significant contribution to UV emission at the He  ii photoelectric edge from hot gas in galaxies and galaxy groups would decrease the expected fluctuations in the column density ratio. Consequently, this model appears difficult to reconcile with the large increase in He  ii opacity fluctuations towards higher redshift. Our results further strengthen previous suggestions that observed He  ii Lyα forest spectra at z ∼ 2–3.5 probe the tail end of the reionization of He  ii by QSOs.