Next Contents Previous

5. THE BRIGHT SIDE OF GALAXY FORMATION AND EVOLUTION

The LambdaCDM scenario of cosmic structure formation has been well tested for perturbations that are still in the linear or quasilinear phase of evolution. These tests are based, among other cosmological probes, on accurate measurements of:

Although these cosmological probes are based on observations of luminous (baryonic) objects, the physics of baryons plays a minor or indirect role in the properties of the linear mass perturbations. The situation is different at small (galaxy) scales, where perturbations went into the non-linear regime and the dissipative physics of baryons becomes relevant. The interplay of DM and baryonic processes is crucial for understanding galaxy formation and evolution. The progress in this field was mostly heuristic; the LambdaCDM scenario provides the initial and boundary conditions for modeling galaxy evolution, but the complex physics of the baryonic processes, in the absence of fundamental theories, requires a model adjustment through confrontation with the observations.

Following, I will outline some key concepts, ingredients, and results of the galaxy evolution study based on the LambdaCDM scenario. Some of the pioneer papers in this field are those of Gunn [57], White & Reese [131], Fall & Efstathiou [43], Blumental et al. [15], Davis et al. [36], Katz & Gunn [65], White & Frenk [130], Kauffmann et al. [66]. For useful lecture notes and recent reviews see e.g., Longair [76, 77], White [129], Steinmetz [113], Firmani & Avila-Reese [46].

The main methods of studying galaxy formation and evolution in the LambdaCDM context are:

5.1. Disks

The formation of galaxy disks deep inside the CDM halos is a generic process in the LambdaCDM scenario. Let us outline the (simplified) steps of disk galaxy formation in this scenario:

1. DM halo growth. The "mold" for disk formation is provided by the mass and AM distributions of the virialized halo, which grows hierarchically. A description of these aspects were presented in the previous Section.

2. Gas cooling and infall, and the maximum mass of galaxies. It is common to assume that the gas in a halo is shock-heated during collapse to the virial temperature [131]. The gas then cools radiatively and falls in a free-fall time to the center. The cooling function Lambda(n, Tk; Z) depends on the gas density, temperature, and composition 14. Since the seminal work by White & Frenk (1990) [130], the rate infall of gas available to form the galaxy is assumed to be driven either by the free-fall time, tff, if tff > tcool or by the cooling time tcool if tff < tcool. The former case applies to halos of masses smaller than approximately 5 × 1011 Modot, whilst the latter applies to more massive halos. The cooling flow from the quasistatic hot atmosphere is the process that basically limits the baryonic mass of galaxies [105], and therefore the bright end of the galaxy luminosity function; for the outer, dilute hot gas in large halos, tcool becomes larger than the Hubble time. However, detailed calculations show that even so, in massive halos too much gas cools, and the bright end of the predicted luminosity function results with a decrease slower than the observed one [12]. Below we will see some solutions proposed to this problem.

More recently it was shown that the cooling of gas trapped in filaments during the halo collapse may be so rapid that the gas flows along the filaments to the center, thus avoiding shock heating [69]. However, this process is efficient only for halos less massive than 2.5 × 1011 Modot, which in any case (even if shock-heating happens), cool their gas very rapidly [19]. Thus, for modeling the formation of disks, and for masses smaller than ~ 5 × 1011 Modot, we may assume that gas infalls in a dynamical time since the halo has virialized, or in two dynamical times since the protostructure was at its maximum expansion.

3. Disk formation, the origin of exponentially, and rotation curves. The gas, originally distributed in mass and AM as the DM, cools and collapses until it reaches centrifugal balance in a disk. Therefore, assuming detailed AM conservation, the radial mass distribution of the disk can be calculated by equating its specific AM to the AM of its final circular orbit in centrifugal equilibrium. The typical collapse factor of the gas within a DM halo is ~ 10 - 15 15, depending on the initial halo spin parameter lambda; the higher the lambda, the more extended (lower surface density) is the resulting disk. The surface density profile of the disks formed within CDM halos is nearly exponential, which provides an explanation to the long-standing question of why galaxy disks are exponential. This is a direct consequence of the AM distribution acquired by the halos by tidal torques and mergers. In more detail, however, the profiles are more concentrated in the center and with a slight excess in the periphery than the exponential law [45, 22]. The cusp in the central disk could give rise to either a photometrical bulge [120] or to a real kinematical bulge due to disk gravitational instability enhanced by the higher central surface density [2] (bulge secular formation). In a few cases (high-lambda, low-concentrated halos), purely exponential disks can be formed.

Baryons are a small mass fraction in the CDM halos, however, the disk formed in the center is very dense (recall the high collapse factors), so that the contribution of the baryonic disk to the inner gravitational potential is important or even dominant. The formed disk will drag gravitationally DM, producing an inner halo contraction that is important to calculate for obtaining the rotation curve decomposition. The method commonly used to calculate it is based on the approximation of radial adiabatic invariance, where spherical symmetry and circular orbits are assumed (e.g., [47, 82]). However, the orbits in CDM halos obtained in N-body simulations are elliptical rather than circular; by generalizing the adiabatic invariance to elliptical orbits, the halo contraction becomes less efficient [132, 52].

The rotation curve decomposition of disks within contracted LambdaCDM halos are in general consistent with observations [82, 45, 132] (nearly-flat total rotation curves; maximum disk for high-surface brightness disks; submaximum disk for the LSB disks; in more detail, the outer rotation curve shape depends on surface density, going from decreasing to increasing at the disk radius for higher to lower densities, respectively). However, there are important non-solved issues. For example, from a large sample of observed rotation curves, Persic et al. [95] inferred that the rotation curve shapes are described by an "universal" profile that (i) depends on the galaxy luminosity and (ii) implies a halo profile different from the CDM (NFW) profile. Other studies confirm only partially these claims [123, 132, 25]. Statistical studies of rotation curves are very important for testing the LambdaCDM scenario.

In general, the structure and dynamics of disks formed within LambdaCDM halos under the assumption of detailed AM conservation seem to be consistent with observations. An important result to remark is the successful prediction of the infrared Tully-Fisher relation and its scatter 16. The core problem mentioned in Section 4.2 is the most serious potential difficulty. Other potential difficulties are: (i) the predicted disk size (surface brightness) distribution implies a P(lambda) distribution narrower than that corresponding to LambdaCDM halos by almost a factor of two [74]; (ii) the internal AM distribution inferred from observations of dwarf galaxies seems not to be in agreement with the LambdaCDM halo AM distribution [122]; (iii) the inference of the halo profile from the statistical study of rotation curve shapes seems not to be agreement with CMD halos. In N-body+hydrodynamical simulations of disk galaxy formation there was common another difficulty called the 'angular momentum catastrophe': the simulated disks ended too much concentrated, apparently due to AM transference of baryons to DM during the gas collapse. The formation of highly concentrated disks also affects the shape of the rotation curve (strongly decreasing), as well as the zero-point of the Tully-Fisher relation. Recent numerical simulations are showing that the 'angular momentum catastrophe', rather than a physical problem, is a problem related to the resolution of the simulations and the correct inclusion of feedback effects.

4. Star formation and feedback. We are coming to the less understood and most complicated aspects of the models of galaxy evolution, which deserve separate notes. The star formation (SF) process is studied at two levels (each one by two separated communities!): (i) the small-scale physics, related to the complex processes by which the cold gas inside molecular clouds fragments and collapses into stars, and (ii) the large-scale physics, related to the disk global instabilities that give rise to the largest unities of SF, the molecular clouds. The SF physics incorporated to galaxy evolution models is still oversimplified, phenomenological and refers to the latter item. The large-scale SF cycle in normal galaxies is believed to be self-regulated by a balance between the energy injection due to SF (mainly SNe) and dissipation (radiative or turbulent). Two main approaches have been used to describe the SF self-regulation in models of galaxy evolution: (a) the halo cooling-feedback approach [130]), (b) the disk turbulent ISM approach [44, 124].

According to the former, the cool gas is reheated by the "galaxy" SF feedback and driven back to the intrahalo medium until it again cools radiatively and collapses into the galaxy. This approach has been used in semi-analytical models of galaxy formation where the internal structure and hydrodynamics of the disks are not treated in detail. The reheating rate is assumed to depend on the halo circular velocity Vc: dot{M}rh propto dot{M}s / Vcalpha, where dot{M}s is the SF rate (SFR) and alpha geq 2. Thus, the galaxy SFR, gas fraction and luminosity depend on Vc. In these models, the disk ISM is virtually ignored and the SN-energy injection is assumed to be as efficient as to reheat the cold gas up to the virial temperature of the halo. A drawback of the model is that it predicts hot X-ray halos around disk galaxies much more luminous than those observed.

Approach (b) is more appropriate for models where the internal processes of the disk are considered. In this approach, the SF at a given radius r is assumed to be triggered by disk gravitational instabilities (Toomre criterion) and self-regulated by a balance between energy injection (mainly by SNe) and dissipation in the turbulent ISM in the direction perpendicular to the disk plane:

Equations 15-16 (15)

(16)

where vg and Sigmag are the gas velocity dispersion and surface density, kappa is the epicyclic frequency, Qcrit is a critical value for instability, gammaSN and epsilonSN are the kinetic energy injection efficiency of the SN into the gas and the SN energy generated per gram of gas transformed into stars, respectively, dot{Sigma}* is the surface SFR, and dot{Sigma}E,accr is the kinetic energy input due to mass accretion rate (or eventually any other energy source as AGN feedback). The key parameter in the self-regulating process is the dissipation time td. The disk ISM is a turbulent, non-isothermal, multi-temperature flow. Turbulent dissipation in the ISM is typically efficient (td ~ 107 - 108yr) in such a way that self-regulation happens at the characteristic vertical scales of the disk. Thus, there is not too much room for strong feedback with the gas at heights larger than the vertical scaleheigth of normal present-day disks: self-regulation is at the level of the disk, but not at the level of the gas corona around. With this approach the predicted SFR is proportional to Sigmagn (Schmidt law), with n approx 1.4 - 2 varying along the disk, in good agreement with observational inferences. The typical SF timescales are not longer than 3 - 4 Gyr. Therefore, to keep active SFRs in the disks, gas infall is necessary, a condition perfectly fulfilled in the LambdaCDM scenario.

Given the SFR radius by radius and time by time, and assuming an IMF, the corresponding luminosities in different color bands can be calculated with stellar population synthesis models. The final result is then an evolving inside-out luminous disk with defined global and local colors.

5. Secular evolution The "quiet" evolution of galaxy disks as described above can be disturbed by minor mergers (satellite accretion) and interactions with close galaxy companions. However, as several studies have shown, the disk may suffer even intrinsic instabilities which lead to secular changes in its structure, dynamics, and SFR. The main effects of secular evolution, i.e. dynamical processes that act in a timescale longer than the disk dynamical time, are the vertical thickening and "heating" of the disk, the formation of bars, which are efficient mechanisms of radial AM and mass redistribution, and the possible formation of (pseudo)bulges (see for recent reviews [71, 33]). Models of disk galaxy evolution should include these processes, which also can affect disk properties, for example increasing the disk scale radii [117].

5.2. Spheroids

As mentioned in Section 2, the simple appearance, the dominant old stellar populations, the alpha-elements enhancement, and the dynamically hot structure of spheroids suggest that they were formed by an early (z gtapprox 4) single violent event with a strong burst of star formation, followed by passive evolution of their stellar population (monolithic mechanism). Nevertheless, both observations and theory point out to a more complex situation. There are two ways to define the formation epoch of a spheroid: when most of its stars formed or when the stellar spheroid acquired its dynamical properties in violent or secular processes. For the monolithic collapse mechanism both epochs coincide.

In the context of the LambdaCDM scenario, spheroids are expected to be formed basically as the result of major mergers of disks. However,

Besides, stellar disks may develop spheroids in their centers (bulges) by secular evolution mechanisms, both intrinsic or enhanced by minor mergers and interactions; this channel of spheroid formation should work for late-type galaxies and it is supported by a large body of observations [71]. But the picture is even more complex in the hierarchical cosmogony as galaxy morphology may be continuously changing, depending on the MAH (smooth accretion and violent mergers) and environment. An spheroid formed early should continue accreting gas so that a new, younger disk grows around. A naive expectation in the context of the LambdaCDM scenario is that massive elliptical galaxies should be assembled mainly by late major mergers of the smaller galaxies in the hierarchy. It is also expected that the disks in galaxies with small bulge-to-disk ratios should be on average redder than those in galaxies with large bulge-to-disk ratios, contrary to observations.

Although it is currently subject of debate, a more elaborate picture of spheroid formation is emerging now in the context of the LambdaCDM hierarchical scenario (see [106, 46, 39] and the references therein). The basic ideas are that massive ellipticals formed early (z gtapprox 3) and in a short timescale by the merging of gas-rich disks in rare high-peak, clustered regions of the Universe. The complex physics of the merging implies (i) an ultraluminous burst of SF obscured by dust (cool ULIRG phase) and the establishment of a spheroidal structure, (ii) gas collapse to the center, a situation that favors the growth of the preexisting massive black hole(s) through an Eddington or even super-Eddington regime (warm ULIRG phase), (iii) the switch on of the AGN activity associated to the supermassive black hole when reaching a critical mass, reverting then the gas inflow to gas outflow (QSO phase), (iv) the switch off of the AGN activity leaving a giant stellar spheroid with a supermassive black hole in the center and a hot gas corona around (passive elliptical evolution). In principle, the hot corona may cool by cooling flows and increase the mass of the galaxy, likely renewing a disk around the spheroid. However, it seems that recurrent AGN phases (less energetic than the initial QSO phase) are possible during the life of the spheroid. Therefore, the energy injected from AGN in the form of radio jets (feedback) can be responsible for avoiding the cooling flow. This way is solved the problem of disk formation around the elliptical, as well as the problem of the extended bright end in the luminosity function. It is also important to note that as soon as the halo hosting the elliptical becomes a subhalo of the group or cluster, the MAH is truncated (Section 4). According to the model just described, massive elliptical galaxies were in place at high redshifts, while less massive galaxies (collapsing from more common density peaks) assembled later. This model was called downsizing or anti-hierarchical. In spite of the name, it fits perfectly within the hierarchical LambdaCDM scenario.

5.3. Drivers of the Hubble sequence



13 Recall that linear theory relates the peculiar velocity, that is the velocity deviation from the Hubble flow, to the density contrast. It is said that the cosmological velocity field is potential; any primordial rotational motion able to give rise to a density perturbation decays as the Universe expands due to angular momentum conservation. Back.

14 The main cooling processes for the intrahalo gas are collisional excitation and ionization, recombination, and bremsstrahlung. The former is the most efficient for kinetic temperatures Tk approx 104 - 105 K and for neutral hydrogen and single ionized helium; for a meta-enriched gas, cooling is efficient at temperatures between 105 - 107 K. At higher temperatures, where the gas is completely ionized, the dominant cooling process is bremsstrahlung. At temperatures lower than 104 K (small halos) and in absence of metals, the main cooling process is by H2 and HD molecule line emission. Back.

15 It is interesting to note that in the absence of a massive halo around galaxies, the collapse factor would be larger by ~ M / Md approx 20, where M and Md are the total halo and disk masses, respectively [90]. Back.

16 In Section 4.1 we have shown that the basis of the Tully-Fisher relation is the CDM halo M - Vm relation. From the pure halo to the disk+halo system there are several intermediate processes that could distort the original M - Vm relation. However, it was shown that the way in which the CDM halo couples with the disk and the way galaxies transform their gas into stars "conspire" to keep the relation. Due to this conspiring, the Tully-Fisher relation is robust to variations in the baryon fraction fB (or mass-to-luminosity ratios) and in the spin parameter lambda [45]. Back.

Next Contents Previous