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.     

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.     

The signature of dark energy on the local Hubble flow

Using N -body simulations of flat, dark energy-dominated cosmologies, we show that galaxies around simulated binary systems resembling the Local Group (LG) have low peculiar velocities, in good agreement with observational data. We have compared results for LG-like systems selected from large, high-resolution simulations of three cosmologies: a ΛCDM model, a ΛWDM model with a 2-keV warm dark matter candidate, and a quintessence (QCDM) model with an equation-of-state parameter   w =−0.6  . The Hubble flow is significantly colder around LGs selected in a flat, Λ-dominated cosmology than around LGs in open or critical models, showing that a dark energy component manifests itself on the scales of nearby galaxies, cooling galaxy peculiar motions. Flows in the ΛWDM and QCDM models are marginally colder than in the ΛCDM one.
The results of our simulations have been compared to existing data and to a new data set of 28 nearby galaxies with robust distance measures (Cepheids and surface brightness fluctuations). The measured line-of-sight velocity dispersion is given by  σ_H= (88 ± 20  km s⁻¹) × ( R /7 Mpc)  . The best agreement with observations is found for LGs selected in the ΛCDM cosmology in environments with  −0.1 < δρ/ρ < 0.6  on scales of 7 Mpc, in agreement with existing observational estimates on the local matter density. These results provide new, independent evidence for the presence of dark energy on scales of a few megaparsecs, corroborating the evidence gathered from observations of distant objects and the early Universe.     

How galaxies lose their angular momentum

The processes are investigated by which gas loses its angular momentum during the protogalactic collapse phase, leading to disc galaxies that are too compact with respect to the observations. High-resolution N -body/SPH simulations in a cosmological context are presented including cold gas and dark matter (DM). A halo with quiet merging activity since redshift   z ∼ 3.8  and with a high-spin parameter is analysed that should be an ideal candidate for the formation of an extended galactic disc. We show that the gas and the DM have similar specific angular momenta until a merger event occurs at   z ∼ 2  with a mass ratio of 5:1. All the gas involved in the merger loses a substantial fraction of its specific angular momentum due to tidal torques and dynamical friction processes falls quickly into the centre. In contrast, gas infall through small subclumps or accretion does not lead to catastrophic angular momentum loss. In fact, a new extended disc begins to form from gas that was not involved in the 5:1 merger event and that falls in subsequently. We argue that the angular momentum problem of disc galaxy formation is a merger problem: in cold dark matter cosmology substantial mergers with mass ratios of 1:1 to 6:1 are expected to occur in almost all galaxies. We suggest that energetic feedback processes could in principle solve this problem, however only if the heating occurs at the time or shortly before the last substantial merger event. Good candidates for such a coordinated feedback would be a merger-triggered starburst or central black hole heating. If a large fraction of the low angular momentum gas would be ejected, late-type galaxies could form with a dominant extended disc component, resulting from late infall, a small bulge-to-disc ratio and a low baryon fraction, in agreement with observations.     

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.     

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.     

Satellite systems around galaxies in hydrodynamic simulations

We investigate the properties of satellite galaxies formed in N -body/SPH simulations of galaxy formation in the ΛCDM cosmology. The simulations include the main physical effects thought to be important in galaxy formation and, in several cases, produce realistic spiral discs. In total, a sample of nine galaxies of luminosity comparable to the Milky Way was obtained. At magnitudes brighter than the resolution limit,   M_V =−12  , the luminosity function of the satellite galaxies in the simulations is in excellent agreement with data for the Local Group. The radial number density profile of the model satellites, as well as their gas fractions also match observations very well. In agreement with previous N -body studies, we find that the satellites tend to be distributed in highly flattened configurations whose major axis is aligned with the major axis of the (generally triaxial) dark halo. In two out of three systems with sufficiently large satellite populations, the satellite system is nearly perpendicular to the plane of the galactic disc, a configuration analogous to that observed in the Milk Way. The discs themselves are perpendicular to the minor axis of their host haloes in the inner parts, and the correlation between the orientation of the galaxy and the shape of the halo persists even out to the virial radius. However, in one case the disc's minor axis ends up, at the virial radius, perpendicular to the minor axis of the halo. The angular momenta of the galaxies and their host halo tend to be well aligned.     

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.     

Faint Lyman-break galaxies as a crucial test for galaxy formation models

It has recently been shown that galaxy formation models within the Λ cold dark matter cosmology predict that, compared to the observed population, small galaxies (with stellar masses  <10¹¹ M_⊙  ) form too early, are too passive since   z ∼ 3  and host too old stellar populations at   z = 0  . We then expect an overproduction of small galaxies at   z ≳ 4  that should be visible as an excess of faint Lyman-break galaxies. To check whether this excess is present, we use the morgana galaxy formation model and grasil spectrophotometric  +  radiative transfer code to generate mock catalogues of deep fields observed with Hubble Space Telescope Advanced Camera for Surveys. We add observational noise and the effect of Lyman α emission, and perform colour–colour selections to identify Lyman-break galaxies. The resulting mock candidates have plausible properties that closely resemble those of observed galaxies. We are able to reproduce the evolution of the bright tail of the luminosity function of Lyman-break galaxies (with a possible underestimate of the number of the brightest i -dropouts), but uncertainties and degeneracies in dust absorption parameters do not allow to give strong constraints to the model. Besides, our model shows a clear excess with respect to observations of faint Lyman-break galaxies, especially of   z ₈₅₀∼ 27 V   -dropouts at   z ∼ 5  . We quantify the properties of these 'excess' galaxies and discuss the implications: these galaxies are hosted in dark matter haloes with circular velocities in excess of 100 km s⁻¹, and their suppression may require a deep rethinking of stellar feedback processes taking place in galaxy formation.     

Mass-to-light ratio gradients in early-type galaxy haloes

Owing to the fact that the near future should see a rapidly expanding set of probes of the halo masses of individual early-type galaxies, we introduce a convenient parameter for characterizing the halo masses from both observational and theoretical results:  ∇_ℓϒ  , the logarithmic radial gradient of the mass-to-light ratio. Using halo density profiles from Λ-cold dark matter (CDM) simulations, we derive predictions for this gradient for various galaxy luminosities and star formation efficiencies  ε_SF  . As a pilot study, we assemble the available  ∇_ℓϒ  data from kinematics in early-type galaxies – representing the first unbiased study of halo masses in a wide range of early-type galaxy luminosities – and find a correlation between luminosity and  ∇_ℓϒ  , such that the brightest galaxies appear the most dark-matter dominated. We find that the gradients in most of the brightest galaxies may fit in well with the ΛCDM predictions, but that there is also a population of fainter galaxies whose gradients are so low as to imply an unreasonably high star formation efficiency  ε_SF > 1  . This difficulty is eased if dark haloes are not assumed to have the standard ΛCDM profiles, but lower central concentrations.     

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).     

The evolution of the galaxy red sequence in simulated clusters and groups

N -body/hydrodynamical simulations of the formation and evolution of galaxy groups and clusters in a Λ cold dark matter (ΛCDM) cosmology are used in order to follow the building-up of the colour–magnitude relation in two clusters and in 12 groups. We have found that galaxies, starting from the more massive, move to the red sequence (RS) as they get aged over times and eventually set upon a 'dead sequence' (DS) once they have stopped their bulk star formation activity. Fainter galaxies keep having significant star formation out to very recent epochs and lie broader around the RS. Environment plays a role as galaxies in groups and cluster outskirts hold star formation activity longer than the central cluster regions. However, galaxies experiencing infall from the outskirts to the central parts keep star formation on until they settle on to the DS of the core galaxies. Merging contributes to mass assembly until z ∼ 1, after which major events only involve the brightest cluster galaxies.
The emerging scenario is that the evolution of the colour–magnitude properties of galaxies within the hierarchical framework is mainly driven by star formation activity during dark matter haloes assembly. Galaxies progressively quenching their star formation settle to a very sharp 'red and dead' sequence, which turns out to be universal, its slope and scatter being almost independent of the redshift (since at least z ∼ 1.5) and environment.
Differently from the DS, the operatively defined RS evolves more evidently with z , the epoch when it changes its slope being closely corresponding to that at which the passive galaxies population takes over the star-forming one: this goes from z ≃ 1 in clusters down to 0.4 in normal groups.     

A sample of 6C radio sources designed to find objects at redshift z>4– III. Imaging and the radio galaxy K–z relation

In this paper, the third and final of a series, we present complete K -band imaging and some complementary I -band imaging of the filtered 6C* sample. We find no systematic differences between the K – z relation of 6C* radio galaxies and those from complete samples, so the near-infrared properties of luminous radio galaxies are not obviously biased by the additional 6C* radio selection criteria (steep spectral index and small angular size). The 6C* K – z data significantly improve delineation of the K – z relation for radio galaxies at high redshift ( z >2) . Accounting for non-stellar contamination, and for correlations between radio luminosity and stellar mass, we find little support for previous claims that the underlying scatter in the stellar luminosity of radio galaxies increases significantly at z >2 . In a particular spatially flat universe with a cosmological constant (Ω_M=0.3 and Ω_Λ=0.7) , the most luminous radio sources appear to be associated with galaxies with a luminosity distribution with a high mean (≈5  L *), and a low dispersion ( σ ∼0.5 mag) which formed their stars at epochs corresponding to z ≳2.5 . This result is in line with recent submillimetre studies of high-redshift radio galaxies and the inferred ages of extremely red objects from faint radio samples.     

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.     

Inferring dark halo structure from observed scaling laws of late-type galaxies and LSBs

We re-examine the Fall & Efstathiou scenario for galaxy formation, including the dark halo gravitational reaction to the formation of the baryon disc, as well as continuous variations in the intrinsic halo density profile. The recently published rotation curves of low surface brightness (LSB) and dwarf galaxies together with previously known scaling relations provide sufficient information on the present-day structure of late-type disc galaxies to invert the problem. By requiring that the models reproduce all the observational restrictions we can fully constrain the initial conditions of galaxy formation, with a minimum of assumptions, in particular without the need to specify a cold dark matter (CDM) halo profile. This allows one to solve for all the initial conditions, in terms of the halo density profile, the baryon fraction and the total angular momentum. We find that a unique initial halo shape is sufficient to accurately reproduce the rotation curves of both LSB and normal late-type spiral galaxies. This unique halo profile differs substantially from that found in standard CDM models. A galactic baryon fraction of 0.065 is found. The initial value of the dimensionless angular momentum is seen to be the principal discriminator between the galaxy classes we examine. The present-day scalings between structural parameters are seen to originate in the initial conditions.     

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     .     

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.     

Survival of substructure within dark matter haloes

Using high-resolution cosmological N -body simulations, we investigate the survival of dark matter satellites falling into larger haloes. Satellites preserve their identity for some time after merging. We compute their loss of mass, energy and angular momentum as they are dissolved by dynamical friction, tidal forces and collisions with other satellites. We also analyse the evolution of their internal structure. Satellites with less than a few per cent of the mass of the main halo may survive for several billion years, whereas larger satellites rapidly sink into the centre of the main halo potential well and lose their identity. Penetrating encounters between satellites are frequent and may lead to significant mass loss and disruption. Only a minor fraction of cluster mass (10–15 per cent on average) is bound to substructure at most redshifts of interest. We discuss the application of these results to the survival and extent of dark matter haloes associated with galaxies in clusters, and to their interactions. We find that a minor fraction of galaxy-size dark matter haloes are disrupted by redshift z  = 0. The fraction of satellites undergoing close encounters is similar to the observed fraction of interacting or merging galaxies in clusters at moderate redshift.     

Understanding the halo-mass and galaxy-mass cross-correlation functions

We use the Millennium Simulation (MS) to measure the cross-correlation between halo centres and mass (or equivalently the average density profiles of dark haloes) in a Lambda cold dark matter (ΛCDM) cosmology. We present results for radii in the range  10  h ⁻¹ kpc < r < 30  h ⁻¹ Mpc  and for halo masses in the range  4 × 10¹⁰ < M ₂₀₀ < 4 × 10¹⁴  h ⁻¹ M_⊙  . Both at   z = 0  and at   z = 0.76  these cross-correlations are surprisingly well fitted if the inner region is approximated by a density profile of NFW or Einasto form, the outer region by a biased version of the linear mass autocorrelation function, and the maximum of the two is adopted where they are comparable. We use a simulation of galaxy formation within the MS to explore how these results are reflected in cross-correlations between galaxies and mass. These are directly observable through galaxy–galaxy lensing. Here also we find that simple models can represent the simulation results remarkably well, typically to ≲10 per cent. Such models can be used to extend our results to other redshifts, to cosmologies with other parameters, and to other assumptions about how galaxies populate dark haloes. Our galaxy formation simulation already reproduces current galaxy–galaxy lensing data quite well. The characteristic features predicted in the galaxy–galaxy lensing signal should provide a strong test of the ΛCDM cosmology as well as a route to understanding how galaxies form within it.