We constrain the assembly history of galaxies in the projected central 0.5 Mpc of the Coma cluster by performing structural decomposition on 69 massive (M_star >= 10^9 M_sun) galaxies using high-resolution F814W images from the HST Treasury Survey of Coma. Each galaxy is modeled with up to three Sersic components having a free Sersic index n. After excluding the two cDs in the projected central 0.5 Mpc of Coma, 57% of the galactic stellar mass in the projected central 0.5 Mpc of Coma resides in classical bulges/ellipticals while 43% resides in cold disk-dominated structures. Most of the stellar mass in Coma may have been assembled through major (and possibly minor) mergers. Hubble types are assigned based on the decompositions, and we find a strong morphology-density relation; the ratio of (E+S0):spirals is (91.0%):9.0%. In agreement with earlier work, the size of outer disks in Coma S0s/spirals is smaller compared with lower-density environments captured with SDSS (Data Release 2). Among similar-mass clusters from a hierarchical semi-analytic model, no single cluster can simultaneously match all the global properties of the Coma cluster. The model strongly overpredicts the mass of cold gas and underpredicts the mean fraction of stellar mass locked in hot components over a wide range of galaxy masses. We suggest that these disagreements with the model result from missing cluster physics (e.g., ram-pressure stripping), and certain bulge assembly modes (e.g., mergers of clumps). Overall, our study of Coma underscores that galaxy evolution is not solely a function of stellar mass, but also of environment.

We present a new semi-analytic model for the common growth of black holes (BHs) and galaxies within a hierarchical Universe. The model is tuned to match the mass function of BHs at z=0 and the luminosity functions of active galactic nuclei (AGNs) at z<4. We use a new observational constraint, which relates the luminosity of AGNs to the star-formation rate (SFR) of their host galaxies. We show that this new constraint is important in various aspects: a) it indicates that BH accretion events are episodic; b) it favours a scenario in which BH accretion is triggered by merger events of all mass ratios; c) it constrains the duration of both merger-induced star-bursts and BH accretion events. The model reproduces the observations once we assume that only 4 per cent of the merger events trigger BH accretion; BHs accretion is not related to secular evolution; and only a few per cent of the mass made in bursts goes into the BH. We find that AGNs with low or intermediate luminosity are mostly being triggered by minor merger events, in broad agreement with observations. Our model matches various observed properties of galaxies, such as the stellar mass function at z<4 and the clustering of galaxies at redshift zero. This allows us to use galaxies as a reliable backbone for BH growth, with reasonable estimates for the frequency of merger events. Other modes of BH accretion, such as disk-instability events, were not considered here, and should be further examined in the future.

We use an empirical approach to model the stellar mass of galaxies according to their host dark-matter haloes and subhaloes ('HASH'), where each galaxy resides in a subhalo taken from a large N-body cosmological simulation. This approach allows us to study the mass relation between subhaloes and galaxies (MR) using: the weak lensing signal (WL), the two point auto-correlation function (CF), and the stellar mass function of galaxies (SMF). Our method is based on modeling the lensing signal directly from the cosmological simulation, and should thus be more accurate than other methods based on analytic halo models. We find that the WL does not provide a strong constraint on the MR, allowing for more than a factor of 10 freedom in the subhalo mass of central galaxies, for a given stellar mass. The freedom for satellite galaxies is much larger, providing a very poor constraint on the fraction of satellite galaxies (0.05-0.8). These results are not significantly modified when using both the SMF and WL as constraints. We show that for the most massive galaxies, observational constraints based on the CF with 0.1 dex errors, are equivalent to 0.05 dex error in the WL. For intermediate and low mass galaxies the WL and CF constrain the MR in a different way. Although the WL is currently not adding much information at these masses, it has the potential of being important using future, more accurate measurements. The models found here do not match simultaneously the observed CF and WL signals, and show a limited ability to match the WL & SMF. We suspect that this is partially due to the cosmological model assumed here and we therefore adopt a mock WL signal through most of this work. In comparison to previous models in the literature, the method presented here is probably more general, as it reveals a larger range of solutions for a given set of observational constraints [abridged].

We study for the first time the environment of massive black hole (BH) seeds (~10^4-5 Msun) formed via the direct collapse of pristine gas clouds in massive haloes (>10^7 Msun) at z>6. Our model is based on the evolution of dark matter haloes within a cosmological N-body simulation, combined with prescriptions for the formation of BH along with both Pop III and Pop II stars. We calculate the spatially-varying intensity of Lyman Werner (LW) radiation from stars and identify the massive pristine haloes in which it is high enough to shut down molecular hydrogen cooling. In contrast to previous BH seeding models with a spatially constant LW background, we find that the intensity of LW radiation due to local sources, J_local, can be up to 10^6 times the spatially averaged background in the simulated volume and exceeds the critical value, J_crit, for the complete suppression of molecular cooling, in some cases by 4 orders of magnitude. Even after accounting for possible metal pollution in a halo from previous episodes of star formation, we find a steady rise in the formation rate of direct collapse (DC) BHs with decreasing redshift from 10^{-3}/Mpc^3/z at z=12 to 10^{-2}/Mpc^3/z at z=6. The onset of Pop II star formation at z~16 simultaneously marks the onset of the epoch of DCBH formation, as the increased level of LW radiation from Pop II stars is able to elevate the local levels of the LW intensity to J_local > J_crit while Pop III stars fail to do so at any time. The number density of DCBHs is sensitive to the number of LW photons and can vary by an order of magnitude at z=6 after accounting for reionisation feedback. Haloes hosting DCBHs are more clustered than similar massive counterparts that do not host DCBHs, especially at redshifts z>10. We also show that planned surveys with JWST should be able to detect the supermassive stellar precursors of DCBHs.

In order to reproduce the low mass end of the stellar mass function, most current models of galaxy evolution invoke very efficient supernova feedback. This solution seems to suffer from several shortcomings however, like predicting too little star formation in low mass galaxies at z=0. In this work, we explore modifications to the star formation (SF) law as an alternative solution to achieve a match to the stellar mass function. This is done by applying semi-analytic models based on De Lucia & Blaizot, but with varying SF laws, to the Millennium and Millennium-II simulations, within the formalism developed by Neistein & Weinmann. Our best model includes lower SF efficiencies than predicted by the Kennicutt-Schmidt law at low stellar masses, no sharp threshold of cold gas mass for SF, and a SF law that is independent of cosmic time. These simple modifications result in a model that is more successful than current standard models in reproducing various properties of galaxies less massive than 10^{10}Msun. The improvements include a good match to the observed auto-correlation function of galaxies, an evolution of the stellar mass function from z=3 to z=0 similar to observations, and a better agreement with observed specific star formation rates. However, our modifications also lead to a dramatic over-prediction of the cold mass content of galaxies. This shows that finding a successful model may require fine-tuning of both star formation and supernovae feedback, as well as improvements on gas cooling, or perhaps the inclusion of a yet unknown process which efficiently heats or expels gas at high redshifts.

In this work we develop a new method to turn a state-of-the-art hydrodynamical cosmological simulation of galaxy formation (HYD) into a simple semi-analytic model (SAM). This is achieved by summarizing the efficiencies of accretion, cooling, star formation, and feedback given by the HYD, as functions of the halo mass and redshift. Surprisingly, by turning the HYD into a SAM, we conserve the mass of individual galaxies, with deviations at the level of 0.1 dex, on an object-by-object basis. This is true for all redshifts, and for the mass of stars and gas components, although the agreement reaches 0.2 dex for satellite galaxies at low redshift. We show that the same level of accuracy is obtained even in case the SAM uses only one phase of gas within each galaxy. Moreover, we demonstrate that the formation history of one massive galaxy provides sufficient information for the SAM to reproduce the population of galaxies within the entire cosmological box. The reasons for the small scatter between the HYD and SAM galaxies are: a) The efficiencies are matched as functions of the halo mass and redshift, meaning that the evolution within merger-trees agrees on average. b) For a given galaxy, efficiencies fluctuate around the mean value on time scales of 0.2-2 Gyr. c) The various mass components of galaxies are obtained by integrating the efficiencies over time, averaging out these fluctuations. We compare the efficiencies found here to standard SAM recipes and find that they often deviate significantly. For example, here the HYD shows smooth accretion that is less effective for low mass haloes, and is always composed of hot or dilute gas; cooling is less effective at high redshift; and star formation changes only mildly with cosmic time. The method developed here can be applied in general to any HYD, and can thus serve as a common language for both HYDs and SAM (Abridged).

The observational indications for a constant specific star-formation rate (sSFR) in the redshift range z=2-7 are puzzling in the context of current galaxy-formation models. Despite the tentative nature of the data, their marked conflict with theory motivates a study of the possible implications. The plateau at sSFR ~ 2 Gyr^-1 is hard to reproduce because (a) its level is low compared to the cosmological specific accretion rate at z > 6, (b) it is higher than the latter at z ~ 2, (c) the natural correlation between SFR and stellar mass makes it difficult to manipulate their ratio, and (d) a low SFR at high z makes it hard to produce enough massive galaxies by z ~ 2. Using a flexible semi-analytic model, we explore ad-hoc modifications to the standard physical recipes trying to obey the puzzling observational constraints. Successful models involve non-trivial modifications, such as (a) a suppressed SFR at z > 4 in galaxies of all masses, by enhanced feedback or reduced SFR efficiency, following an initial active phase at z > 7, (b) a delayed gas consumption into stars, allowing the gas that was prohibited from forming stars or ejected at high z to form stars later in more massive galaxies, and (c) enhanced growth of massive galaxies, in terms of either faster assembly or more efficient starbursts in mergers, or by efficient star formation in massive haloes.

We develop a new empirical methodology to study the relation between the stellar mass of galaxies and the mass of their host subhaloes. Our approach is similar to abundance matching, and is based on assigning a stellar mass to each subhalo within a large cosmological N-body simulation. However, we significantly extend the method to include a different treatment for central and satellite galaxies, allowing the stellar mass of satellite galaxies to depend on both the host subhalo mass, and on the halo mass. Unlike in previous studies, our models are constrained by two observations: the stellar mass function of galaxies, and their auto-correlation functions (CFs). We present results for ~10^6 different successful models, showing that the uncertainty in the host subhalo mass reaches a factor of ~10 for a given stellar mass. Our results thus indicate that the host subhalo mass of central and satellite galaxies is poorly constrained by using information coming solely from the abundance and CFs of galaxies. In addition, we demonstrate that the fraction of stellar mass relative to the universal baryon fraction can vary between different models, and can reach ~0.6 for a specific population of low mass galaxies. We conclude that additional observational constraints are needed, in order to better constrain the mass relation between haloes and galaxies. These might include weak lensing, satellite kinematics, or measures of clustering other than the CFs. Moreover, CFs based on larger surveys with a better estimate of the systematic uncertainties are needed.

Recent studies emphasize that an empirical relation between the stellar mass of galaxies and the mass of their host dark matter subhaloes can predict the clustering of galaxies and its evolution with cosmic time. In this paper we study the assumptions made by this methodology using a semi-analytical model (SAM). To this end, we randomly swap between the locations of model galaxies within a narrow range of subhalo mass (M_infall). We find that shuffled samples of galaxies have different auto-correlation functions in comparison with the original model galaxies. This difference is significant even if central and satellite galaxies are allowed to follow a different relation between M_infall and stellar mass, and can reach a factor of 2 for massive galaxies at redshift zero. We analyze three features within SAMs that contribute to this effect: a) The relation between stellar mass and subhalo mass evolves with redshift for central galaxies, affecting satellite galaxies at the time of infall. b) The stellar mass of galaxies falling into groups and clusters at high redshift is different from the mass of central galaxies at the same time. c) The stellar mass growth for satellite galaxies after infall can be significant and depends on the infall redshift and the group mass. We show that the above is true for differing SAMs, and that the effect is sensitive to the treatment of dynamical friction and stripping of gas in satellite galaxies. We find that by using the FoF group mass at redshift zero in addition to M_infall, an empirical model is able to accurately reproduce the clustering properties of galaxies. On the other hand, using the infall redshift as a second parameter does not yield as good results because it is less correlated with stellar mass. Our analysis indicates that environmental processes are important for modeling the clustering and abundance of galaxies. (Abridged)

We develop a new formalism for modeling the formation and evolution of galaxies within a hierarchical universe. Similarly to standard semi-analytical models we trace galaxies inside dark-matter merger-trees. The formalism includes treatment of feedback, star-formation, cooling, smooth accretion, gas stripping in satellite galaxies, and merger-induced star bursts. However, unlike in other models, each process is assumed to have an efficiency which depends only on the host halo mass and redshift. This allows us to describe the various components of the model in a simple and transparent way. By allowing the efficiencies to have any value for a given halo mass and redshift, we can easily encompass a large range of scenarios. To demonstrate this point, we examine several different galaxy formation models, which are all consistent with the observational data. Each model is characterized by a different unique feature: cold accretion in low mass haloes, zero feedback, stars formed only in merger-induced bursts, and shutdown of star-formation after mergers. Using these models we are able to examine the degeneracy inherent in galaxy formation models, and look for observational data that will help to break this degeneracy. We show that the full distribution of star-formation rates in a given stellar mass bin is promising in constraining the models. We compare our approach in detail to the semi-analytical model of De Lucia & Blaizot. It is shown that our formalism is able to produce a very similar population of galaxies once the same median efficiencies per halo mass and redshift are being used. We provide a public version of the model galaxies on our web-page, along with a tool for running models with user-defined parameters. Our model is able to provide results for a 62.5 h^{-1} Mpc box within just a few seconds.

We study merger histories of dark-matter haloes in a suite of N-body simulations that span different cosmological models. The simulated cases include the up-to-date WMAP5 cosmology and other test cases based on the Einstein-deSitter cosmology with different power spectra. We provide a robust fitting function for the conditional mass function (CMF) of progenitor haloes of a given halo. This fit is valid for the different cosmological models and for different halo masses and redshifts, and it is a significant improvement over earlier estimates. Based on this fit, we develop a simple and accurate technique for transforming the merger history of a given simulated halo into haloes of different mass, redshift and cosmology. Other statistics such as main-progenitor history and merger rates are accurately transformed as well. This method can serve as a useful tool for studying galaxy formation. It is less sensitive to the low accuracy of the fit at small time-steps, and it can thus replace the more elaborate task of construction Monte-Carlo realizations. As an alternative approach, we confirm the earlier finding by Neistein & Dekel that the main-progenitor follows a log-normal distribution. This property of merger trees allows us to better capture their behaviour as a function of time and descendant mass, but a broader suite of simulations is required for evaluating the dependence of the log-normal parameters on the cosmological model.

The massive galaxies in the young universe, ten billion years ago, formed stars at surprising intensities. Although this is commonly attributed to violent mergers, the properties of many of these galaxies are incompatible with such events, showing gas-rich, clumpy, extended rotating disks not dominated by spheroids (Genzel et al. 2006, 2008). Cosmological simulations and clustering theory are used to explore how these galaxies acquired their gas. Here we report that they are stream-fed galaxies, formed from steady, narrow, cold gas streams that penetrate the shock-heated media of massive dark matter haloes (Dekel & Birnboim 2006; Keres et al. 2005). A comparison with the observed abundance of star-forming galaxies implies that most of the input gas must rapidly convert to stars. One-third of the stream mass is in gas clumps leading to mergers of mass ratio greater than 1:10, and the rest is in smoother flows. With a merger duy cycle of 0.1, three-quarters of the galaxies forming stars at a given rate are fed by smooth streams. The rarer, submillimetre galaxies that form stars even more intensely are largely merger-induced starbursts. Unlike destructive mergers, the streams are likely to keep the rotating disk configuration intact, although turbulent and broken into giant star-forming clumps that merge into a central spheroid (Noguchi 1999; Genzel et al. 2008, Elmegreen, Bournaud & Elmegreen 2008, Dekel, Sari & Ceverino 2009). This stream-driven scenario for the formation of disks and spheroids is an alternative to the merger picture.

We derive analytic merger rates for dark-matter haloes within the framework of the Extended Press-Schechter (EPS) formalism. These rates become self-consistent within EPS once we realize that the typical merger in the limit of a small time-step involves more than two progenitors, contrary to the assumption of binary mergers adopted in earlier studies. We present a general method for computing merger rates that span the range of solutions permitted by the EPS conditional mass function, and focus on a specific solution that attempts to match the merger rates in N-body simulations. The corrected EPS merger rates are more accurate than the earlier estimates of Lacey & Cole, by ~20% for major mergers and by up to a factor of ~3 for minor mergers of mass ratio 1:10^4. Based on the revised merger rates, we provide a new algorithm for constructing Monte-Carlo EPS merger trees, that could be useful in Semi-Analytic Modeling. We provide analytic expressions and plot numerical results for several quantities that are very useful in studies of galaxy formation. This includes (a) the rate of mergers of a given mass ratio per given final halo, (b) the fraction of mass added by mergers to a halo, and (c) the rate of mergers per given main progenitor. The creation and destruction rates of haloes serve for a self-consistency check. Our method for computing merger rates can be applied to conditional mass functions beyond EPS, such as those obtained by the ellipsoidal collapse model or extracted from $N$-body simulations.

We present a simple and efficient empirical algorithm for constructing dark-matter halo merger trees that reproduce the distribution of trees in the Millennium cosmological $N$-body simulation. The generated trees are significantly better than EPS trees. The algorithm is Markovian, and it therefore fails to reproduce the non-Markov features of trees across short time steps, except for an accurate fit to the evolution of the average main progenitor. However, it properly recovers the full main progenitor distribution and the joint distributions of all the progenitors over long-enough time steps, $\Delta \omega \simeq \Delta z>0.5$, where $\omega \simeq 1.69/D(t)$ is the self-similar time variable and $D(t)$ refers to the linear growth of density fluctuations. We find that the main progenitor distribution is log-normal in the variable $\sigma^2(M)$, the variance of linear density fluctuations in a sphere encompassing mass $M$. The secondary progenitors are successfully drawn one by one from the remaining mass using a similar distribution function. These empirical findings may be clues to the underlying physics of merger-tree statistics. As a byproduct, we provide useful, accurate analytic time-invariant approximations for the main progenitor accretion history and for halo merger rates.

We simulate the buildup of galaxies by spherical gas accretion through dark matter haloes, subject to the development of virial shocks. We find that a uniform cosmological accretion turns into a rapidly varying disc buildup rate. The generic sequence of events (Shocked-Accretion Massive Burst & Shutdown: SAMBA) consists of four distinct phases: a) continuous cold accretion while the halo is below a threshold mass Msh~10^12Mo, b) tentative quenching of gas supply for ~2Gyr, starting once the halo is ~Msh and growing a rapidly expanding shock, c) a massive burst due to the big crunch of ~10^11Mo gas in \~0.5Gyr, when the heated gas cools and joins new infalling gas, and d) a long-term shutdown, enhanced by a temporary shock instability in late SAMBAs, those that quench at z~2, burst at z~1 and end up quenched in 10^12-13Mo haloes today. These events occur at all redshifts in galaxies of baryonic mass \~10^11Mo and involve a substantial fraction of this mass. They arise from rather smooth accretion, or minor mergers, which, unlike major mergers, may leave the disc intact while being built in a rapid pace. The early bursts match observed maximum starbursting discs at z>~2, predicted to reside in <~10^13Mo haloes. The late bursts resemble discy LIRGs at z<~1. The tentative quenching gives rise to a substantial population of ~10^11Mo galaxies with a strongly suppressed star-formation rate at z~1-3. The predicted long-term shutdown leads to red & dead galaxies in groups. A complete shutdown in more massive clusters requires an additional quenching mechanism, as may be provided by clumpy accretion. Alternatively, the SAMBA bursts may trigger the AGN activity that couples to the hot gas above Msh and helps the required quenching. The SAMBA phenomenon is yet to be investigated using cosmological simulations.

Stellar-population analyses of today's galaxies show "downsizing", where the stars in more massive galaxies tend to have formed earlier and over a shorter time span. We show that this phenomenon is not necessarily "anti-hierarchical" but rather has its natural roots in the bottom-up clustering process of dark-matter haloes. While the main progenitor does indeed show an opposite effect, the integrated mass in all the progenitors down to a given minimum mass shows a robust downsizing that is qualitatively similar to what has been observed. These results are derived analytically from the standard extended Press Schechter (EPS) theory, and are confirmed by merger trees based on EPS or drawn from N-body simulations. The downsizing is valid for any minimum mass, as long as it is the same for all haloes at any given time, but the effect is weaker for smaller minimum mass. If efficient star formation is triggered by atomic cooling, then a minimum halo mass arises naturally from the minimum virial temperature for cooling, T~10^4 K, though for such a small minimum mass the effect is weaker than observed. Baryonic feedback effects, which are expected to stretch the duration of star formation in small galaxies and shut it down in massive haloes at late epochs, are likely to play a subsequent role in shaping up the final downsizing behaviour. Other appearances of downsizing, such as the decline with time of the typical mass of star-forming galaxies, may not be attributed to the gravitational clustering process but rather arise from the gas processes.

We present an I-band Tully-Fisher relation (TFR) for 18 nearby S0 galaxies using kinematics derived from long slit spectroscopy of stellar absorption lines. Our estimates of the circular velocity, V_c, at 2-3 exponential disk scale lengths account for line-of-sight projection and for the stellar random motions through an asymmetric drift correction. Uniform and accurate distance calibration for all objects is available from surface brightness fluctuation measurements of Tonry et al. (1998). Despite the care taken in estimating both V_c and M_I, the TFR shows an intrinsic scatter, $\sim 0.7$ mag in M_I, or 0.15 in log_10(V_c). This result is surprising, as S0 galaxies appear to have both the simple kinematics of disk galaxies, and the simple stellar populations of early-type galaxies. Remarkably, in this sample of overall rotation-dominated galaxies, the central stellar velocity dispersion is a better predictor of the total I-band luminosity (through the Fundamental Plane relations) than the circular speed at several exponential scale lengths. Furthermore, the TFR zeropoint, or the mean stellar I-band luminosity at a given V_c, differs by only $\sim 0.5$ mag between our sample of S0s and Mathewson et al.'s (1992) sample of late-type spirals, once both data sets are brought onto a consistent distance scale. This offset is less than expected if S0s are former spiral galaxies with prematurely truncated star-formation ($\gtorder 4$ Gyrs ago).