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.     

Forming a large disc galaxy from a z < 1 major merger

Using high-resolution SPH simulations in a fully cosmological Λ cold dark matter context, we study the formation of a bright disc-dominated galaxy that originates from a 'wet' major merger at   z = 0.8  . The progenitors of the disc galaxy are themselves disc galaxies that formed from early major mergers between galaxies with blue colours. A substantial thin stellar disc grows rapidly following the last major merger and the present-day properties of the final remnant are typical of early-type spiral galaxies, with an i -band bulge-to-disc ratio ∼0.65, a disc scalelength of 7.2 kpc,   g − r = 0.5 mag  , an H  i linewidth ( W ₂₀/2) of 238 km s⁻¹ and total magnitude   i =−22.4  . The key ingredients for the formation of a dominant stellar disc component after a major merger are (i) substantial and rapid accretion of gas through cold flows followed at late times by cooling of gas from the hot phase, (ii) supernova feedback that is able to partially suppress star formation during mergers and (iii) relative fading of the spheroidal component. The gas fraction of the progenitors' discs does not exceed 25 per cent at   z < 3  , emphasizing that the continuous supply of gas from the local environment plays a major role in the regrowth of discs and in keeping the galaxies blue. The results of this simulation alleviate the problem posed for the existence of disc galaxies by the high likelihood of interactions and mergers for galaxy-sized haloes at relatively low z .     

The effects of gas on morphological transformation in mergers: implications for bulge and disc demographics

Transformation of discs into spheroids via mergers is a well-accepted element of galaxy formation models. However, recent simulations have shown that the bulge formation is suppressed in increasingly gas-rich mergers. We investigate the global implications of these results in a cosmological framework, using independent approaches: empirical halo-occupation models (where galaxies are populated in haloes according to observations) and semi-analytic models. In both, ignoring the effects of gas in mergers leads to the overproduction of spheroids: low- and intermediate-mass galaxies are predicted to be bulge-dominated (   B / T ∼ 0.5  at  <10¹⁰ M_⊙  , with almost no 'bulgeless' systems), even if they have avoided major mergers. Including the different physical behaviour of gas in mergers immediately leads to a dramatic change: bulge formation is suppressed in low-mass galaxies, observed to be gas-rich (giving   B / T ∼ 0.1  at  <10¹⁰ M_⊙  , with a number of bulgeless galaxies in good agreement with observations). Simulations and analytic models which neglect the similarity-breaking behaviour of gas have difficulty reproducing the strong observed morphology–mass relation. However, the observed dependence of gas fractions on mass, combined with suppression of bulge formation in gas-rich mergers, naturally leads to the observed trends. Discrepancies between observations and models that ignore the role of gas increase with redshift; in models that treat gas properly, galaxies are predicted to be less bulge-dominated at high redshifts, in agreement with the observations. We discuss implications for the global bulge mass density and future observational tests.     

The effect of galaxy mass ratio on merger-driven starbursts

We employ numerical simulations of galaxy mergers to explore the effect of galaxy mass ratio on merger-driven starbursts. Our numerical simulations include radiative cooling of gas, star formation, and stellar feedback to follow the interaction and merger of four disc galaxies. The galaxy models span a factor of 23 in total mass and are designed to be representative of typical galaxies in the local universe. We find that the merger-driven star formation is a strong function of merger mass ratio, with very little, if any, induced star formation for large mass ratio mergers. We define a burst efficiency that is useful to characterize the merger-driven star formation and test that it is insensitive to uncertainties in the feedback parametrization. In accord with previous work we find that the burst efficiency depends on the structure of the primary galaxy. In particular, the presence of a massive stellar bulge stabilizes the disc and suppresses merger-driven star formation for large mass ratio mergers. Direct, coplanar merging orbits produce the largest tidal disturbance and yield the most intense burst of star formation. Contrary to naive expectations, a more compact distribution of gas or an increased gas fraction both decrease the burst efficiency. Owing to the efficient feedback model and the newer version of smoothed particle hydrodynamics employed here, the burst efficiencies of the mergers presented here are smaller than in previous studies.     

Testing model predictions of the cold dark matter cosmology for the sizes, colours, morphologies and luminosities of galaxies with the SDSS

The huge size and uniformity of the Sloan Digital Sky Survey (SDSS) make possible an exacting test of current models of galaxy formation. We compare the predictions of the galform semi-analytical galaxy formation model for the luminosities, morphologies, colours and scalelengths of local galaxies. galform models the luminosity and size of the disc and bulge components of a galaxy, and so we can compute quantities which can be compared directly with SDSS observations, such as the Petrosian magnitude and the Sérsic index. We test the predictions of two published models set in the cold dark matter cosmology: the Baugh et al. model, which assumes a top-heavy initial mass function (IMF) in starbursts and superwind feedback, and the Bower et al. model, which uses active galactic nucleus feedback and a standard IMF. The Bower et al. model better reproduces the overall shape of the luminosity function, the morphology–luminosity relation and the colour bimodality observed in the SDSS data, but gives a poor match to the size–luminosity relation. The Baugh et al. model successfully predicts the size–luminosity relation for late-type galaxies. Both models fail to reproduce the sizes of bright early-type galaxies. These problems highlight the need to understand better both the role of feedback processes in determining galaxy sizes, in particular the treatment of the angular momentum of gas reheated by supernovae, and the sizes of the stellar spheroids formed by galaxy mergers and disc instabilities.     

Stormy weather and cluster radio galaxies

New adaptive mesh refinement N-body+hydrodynamics numerical simulations are used to illustrate the complex and changing cluster environments in which many radio galaxies live and evolve. Groups and clusters of galaxies form at the intersections of filaments where they continue to accrete gas and dark matter to the present day. The accretion process produces shocks, turbulence, and transonic bulk flows forming a kind of stormy weather within the intracluster medium (ICM). Radio sources embedded within the stormy ICM form distorted, complex morphologies as observed in recent VLA cluster surveys. We show that the bending of wide-angle tailed radio sources can be understood as the result of recent cluster–subcluster mergers. We use new MHD simulations to illustrate how cluster radio halos can be formed by the shocks and turbulence produced during cluster mergers. Finally, we discuss new observations of distant Abell clusters that reveal a class of weak radio sources, probably starbursts, likely produced during the formation of the clusters as they accrete material from the supercluster environment.     

Cold gas accretion in galaxies

Evidence for the accretion of cold gas in galaxies has been rapidly accumulating in the past years. HI observations of galaxies and their environment have brought to light new facts and phenomena which are evidence of ongoing or recent accretion: (1) A large number of galaxies are accompanied by gas-rich dwarfs or are surrounded by HI cloud complexes, tails and filaments. This suggests ongoing minor mergers and recent arrival of external gas. It may be regarded, therefore, as direct evidence of cold gas accretion in the local universe. It is probably the same kind of phenomenon of material infall as the stellar streams observed in the halos of our galaxy and M 31. (2) Considerable amounts of extra-planar HI have been found in nearby spiral galaxies. While a large fraction of this gas is undoubtedly produced by galactic fountains, it is likely that a part of it is of extragalactic origin. Also the Milky Way has extra-planar gas complexes: the Intermediate- and High-Velocity Clouds (IVCs and HVCs). (3) Spirals are known to have extended and warped outer layers of HI. It is not clear how these have formed, and how and for how long the warps can be sustained. Gas infall has been proposed as the origin. (4) The majority of galactic disks are lopsided in their morphology as well as in their kinematics. Also here recent accretion has been advocated as a possible cause. In our view, accretion takes place both through the arrival and merging of gas-rich satellites and through gas infall from the intergalactic medium (IGM). The new gas could be added to the halo or be deposited in the outer parts of galaxies and form reservoirs for replenishing the inner parts and feeding star formation. The infall may have observable effects on the disk such as bursts of star formation and lopsidedness. We infer a mean “visible” accretion rate of cold gas in galaxies of at least . In order to reach the accretion rates needed to sustain the observed star formation (), additional infall of large amounts of gas from the IGM seems to be required.     

Damped Lyman-α Systems

Damped Lyman-α systems observed in the spectra of high-redshift quasars are considered as the progenitors of present-day galaxies.Indeed, the large neutral hydrogen column densities observed and the presence of metals imply that the gas is somehow closely associated with regions of star formation. The nature of the absorbing object is unclear however.The discussion of whether high-redshift damped Lyman-α systems are produced throughlarge, fast-rotating, thick disks or through building-blocks of galaxies is important since it is related to how present-day galaxies form, either through initial formation of large disks and subsequent accretion of gas or as a result of merging of pregalactic clumps.We argue here that damped Lyman-α systems are progenitors of any kind of galaxy (dwarfs and large spirals); they trace peaks in the spatial density distribution of the gas. This revised version was published online in July 2006 with corrections to the Cover Date.     

Dynamical friction and galaxy merging time-scales

The time-scale for galaxies within merging dark matter haloes to merge with each other is an important ingredient in galaxy formation models. Accurate estimates of merging time-scales are required for predictions of astrophysical quantities such as black hole binary merger rates, the build-up of stellar mass in central galaxies and the statistical properties of satellite galaxies within dark matter haloes. In this paper, we study the merging time-scales of extended dark matter haloes using N -body simulations. We compare these results to standard estimates based on the Chandrasekhar theory of dynamical friction. We find that these standard predictions for merging time-scales, which are often used in semi-analytic galaxy formation models, are systematically shorter than those found in simulations. The discrepancy is approximately a factor of 1.7 for M _sat/ M _host≈ 0.1 and becomes larger for more disparate satellite-to-host mass ratios, reaching a factor of ∼3.3 for M _sat/ M _host≈ 0.01. Based on our simulations, we propose a new, easily implementable fitting formula that accurately predicts the time-scale for an extended satellite to sink from the virial radius of a host halo down to the halo's centre for a wide range of M _sat/ M _host and orbits. Including a central bulge in each galaxy changes the merging time-scale by ≲10 per cent. To highlight one concrete application of our results, we show that merging time-scales often used in the literature overestimate the growth of stellar mass by satellite accretion by ≈40 per cent, with the extra mass gained in low mass ratio mergers.     

HST/WFPC2 imaging of the QDOT ultraluminous infrared galaxy sample

We present HST WFPC2 V -band imaging for 23 ultraluminous infrared galaxies (ULIRGs) taken from the QDOT redshift survey. The fraction of sources observed to be interacting is 87 per cent. Most of the merging systems show a number of compact 'knots', whose colour and brightness differ substantially from their immediate surroundings. Colour maps for nine of the objects show a non-uniform colour structure. Features include blue regions located towards the centres of merging systems which are likely to be areas of enhanced star formation, and compact red regions which are likely to be dust shrouded starbursts or active galactic nuclei. The host galaxies of the quasi-stellar objects (QSOs) in the sample were found to be either interacting systems or ellipticals. Our data show no evidence that ULIRGs are a simple transition stage between galaxy mergers and QSOs. We propose an alternative model for ULIRGs based on the morphologies in our sample and previous N -body simulations. Under this model ULIRGs as a class are much more diverse than a simple transition between galaxy merger and QSO. The evolution of IR power source and merger morphology in ULIRGs is driven solely by the local environment and the morphologies of the merger progenitors.     

Globular Clusters in Blue Compact Galaxies as Tracers of the Starburst History

Representing single stellar populations, globular clusters (GCs) arerelatively easy to model, thus providing powerful tools for studyingthe evolution of galaxies. This has been demonstrated for the bluecompact galaxy ESO338-IG04. GC systems in galaxies may be fossils ofstarbursts and mergers. Thus studies of GCs in the local universe mayadd to our understanding of the formation and evolution of galaxies andthe distant universe. This revised version was published online in July 2006 with corrections to the Cover Date.     

Understanding dwarf galaxies as galactic building blocks

This is a summary of a general discussion held during the third EuroConference on galaxy evolution. Various observational features of the stellar populations in present-day dwarf galaxies were presented to introduce the discussion on the possibility that these systems be the main building blocks of spiral and elliptical galaxies. Many people in the audience turned out to think that the inconsistencies among the observed properties of large and dwarf galaxies are too many to believe that the former are built up only by means of successive accretions of the latter. However, theorists of hierarchical galaxy formation suggested that present-day dwarfs are not representative of the galactic building blocks, which may be completely invisible nowadays. Some of them suggested that, contrary to what is usually assumed in hierarchical modelling, the actual building blocks were still fully gaseous systems when their major mergers occurred. If this is the case, then most of the inconsistencies can be overcome, and the scenario of hierarchical galaxy formation becomes not too different from that of a slow gas accretion. This revised version was published online in August 2006 with corrections to the Cover Date.     

Does Deuterium Enable the Formation of Primordial Brown Dwarfs?

The time evolution of barred structures is examined under the influence of the external forces exerted by a spherical halo and by prolate halos. In particular, galaxy disks are placed in the plane including the major axis of prolate halos, whose configuration is often found in cosmological simulations. N-body disks in fixed external halo fields are simulated, so that bars are formed via dynamical instability. In the subsequent evolution, the bars in prolate halos dissolve gradually with time, while the bar pattern in a spherical halo remains almost unchanged to the end of the simulation. The decay times of the bars suggest that they can be destroyed in a time smaller than a Hubble time. Our results indicate that this dissolution process could occur in real barred galaxies, if they are surrounded by massive dark prolate halos, and the configuration adopted here is not unusual from the viewpoint of galaxy formation. For a prolate halo model, an additional simulation that is restricted to two-dimensional in-plane motions has also ended up with similar bar dissolution. This means that the vertical motions of disk stars do not play an essential role in the bar dissolution demonstrated here.     

Star formation sustained by gas accretion

Numerical simulations predict that metal-poor gas accretion from the cosmic web fuels the formation of disk galaxies. This paper discusses how cosmic gas accretion controls star formation, and summarizes the physical properties expected for the cosmic gas accreted by galaxies. The paper also collects observational evidence for gas accretion sustaining star formation. It reviews evidence inferred from neutral and ionized hydrogen, as well as from stars. A number of properties characterizing large samples of star-forming galaxies can be explained by metal-poor gas accretion, in particular, the relationship among stellar mass, metallicity, and star-formation rate (the so-called fundamental metallicity relationship). They are put forward and analyzed. Theory predicts gas accretion to be particularly important at high redshift, so indications based on distant objects are reviewed, including the global star-formation history of the universe, and the gas around galaxies as inferred from absorption features in the spectra of background sources.     

The effects of photoionization on galaxy formation – I. Model and results at z=0

We develop a coupled model for the evolution of the global properties of the intergalactic medium (IGM) and the formation of galaxies, in the presence of a photoionizing background due to stars and quasars. We use this model to predict the thermodynamic history of the IGM when photoionized by galaxies forming in a cold dark matter (CDM) universe. The evolution of the galaxies is calculated using a semi-analytical model, including a detailed treatment of the effects of tidal stripping and dynamical friction on satellite galaxies orbiting inside larger dark matter haloes. We include in the model the negative feedback on galaxy formation from the photoionizing background. Photoionization inhibits galaxy formation in low-mass dark matter haloes in two ways: (i) heating of the IGM and inhibition of the collapse of gas into dark haloes by the IGM pressure, and (ii) reduction in the rate of radiative cooling of gas within haloes. The result of our method is a self-consistent model of galaxy formation and the IGM. The IGM is reheated twice (during reionization of H  i and He  ii ), and we find that the star formation rate per unit volume is slightly suppressed after each episode of reheating. We find that galaxies brighter than L _★ are mostly unaffected by reionization, while the abundance of faint galaxies is significantly reduced, leading to present-day galaxy luminosity functions with shallow faint-end slopes, in good agreement with recent observational data. Reionization also affects other properties of these faint galaxies, in a readily understandable way.     

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.     

Gas dynamic stripping and X-ray emission of cluster elliptical galaxies

Detailed three-dimensional numerical simulations of an elliptical galaxy orbiting in a gas-rich cluster of galaxies indicate that gas dynamic stripping is less efficient than the results from previous, simpler calculations by Takeda et al. and Gaetz et al. implied. This result is consistent with X-ray data for cluster elliptical galaxies. Hydrodynamic torques and direct accretion of orbital angular momentum can result in the formation of a cold gaseous disc, even in a non-rotating galaxy. The gas lost by cluster galaxies via the process of gas dynamic stripping tends to produce a colder, chemically enriched cluster gas core. A comparison of the models with the available X-ray data of cluster galaxies shows that the X-ray luminosity distribution of cluster galaxies may reflect hydrodynamic stripping, but also that a purely hydrodynamic treatment is inadequate for the cooler interstellar medium near the centre of the galaxy.     

Milky Way type galaxies in a ΛCDM cosmology

We analyse a sample of 52 000 Milky Way (MW) type galaxies drawn from the publicly available galaxy catalogue of the Millennium Simulation with the aim of studying statistically the differences and similarities of their properties in comparison to our Galaxy. Model galaxies are chosen to lie in haloes with maximum circular velocities in the range 200–250 km s⁻¹ and to have bulge-to-disc ratios similar to that of the MW. We find that model MW galaxies formed 'quietly' through the accretion of cold gas and small satellite systems. Only ≈12 per cent of our model galaxies experienced a major merger during their lifetime. Most of the stars formed ' in situ ', with only about 15 per cent of the final mass gathered through accretion. Supernovae (SNe) and active galactic nuclei (AGN) feedback play an important role in the evolution of these systems. At high redshifts, when the potential wells of the MW progenitors are shallower, winds driven by SNe explosions blow out a large fraction of the gas and metals. As the systems grow in mass, SNe feedback effects decrease and AGN feedback takes over, playing a more important role in the regulation of the star formation activity at lower redshifts. Although model MW galaxies have been selected to lie in a narrow range of maximum circular velocities, they nevertheless exhibit a significant dispersion in the final stellar masses and metallicities. Our analysis suggests that this dispersion results from the different accretion histories of the parent dark matter haloes. Statistically, we also find evidences to support the MW as a typical Sb/Sc galaxy in the same mass range, providing a suitable benchmark to constrain numerical models of galaxy formation.