We describe an algorithm for predicting the concentrations of dark matter halos via a random walk in energy space. Given a full merger tree for a halo, the total internal energy of each halo in that tree is determined by summing the internal and orbital energies of progenitor halos. For halos described by single-parameter density profiles (such as the NFW profile) the energy can be directly mapped to a scale radius, and so to a concentration. We show that this model can accurately reproduce the mean of the concentration mass relation measured in N-body simulations, and reproduces more of the scatter in that relation than previous models. However, our model underpredicts the kurtosis of the distribution of N-body concentrations. We test this model by examining both the autocorrelation of scale radii across time, and the correlations between halo concentration and spin, and comparing to results measured from cosmological N-body simulations. In both cases we find that our model closely matches the N-body results. Our model is implemented within the open source Galacticus toolkit.

We extend the random-walk model of Vitvitska et al. for predicting the spins of dark matter halos from their merger histories. Using updated merger rates, orbital parameter distributions, and N-body constraints we show that this model can accurately reproduce the distribution of spin parameters measured in N-body simulations when we include a weak correlation between the spins of halos and the angular momenta of infalling subhalos. We further show that this model is in approximate agreement with the correlation of the spin magnitude over time as determined from N-body simulations, while it slightly underpredicts the correlation in the direction of the spin vector measured from the same simulations. This model is useful for predicting spins from merger histories derived from non-N-body sources, thereby circumventing the need for very high resolution simulations to permit accurate measurements of spins. It may be particularly relevant to modeling systems which accumulate angular momentum from halos over time (such as galactic discs) - we show that this model makes small but significant changes in the distribution of galactic disc sizes computed using the Galacticus semi-analytic galaxy formation model.

The Tri-Agency Cosmological Simulations (TACS) Task Force was formed when Program Managers from the Department of Energy (DOE), the National Aeronautics and Space Administration (NASA), and the National Science Foundation (NSF) expressed an interest in receiving input into the cosmological simulations landscape related to the upcoming DOE/NSF Vera Rubin Observatory (Rubin), NASA/ESA's Euclid, and NASA's Wide Field Infrared Survey Telescope (WFIRST). The Co-Chairs of TACS, Katrin Heitmann and Alina Kiessling, invited community scientists from the USA and Europe who are each subject matter experts and are also members of one or more of the surveys to contribute. The following report represents the input from TACS that was delivered to the Agencies in December 2018.

We present a semi-analytic model of satellite galaxies, SatGen, which can generate large samples of satellite populations for a host halo of desired mass, redshift, and assembly history. The model combines dark-matter halo merger trees, empirical relations for the galaxy-halo connection, and analytic prescriptions for tidal effects, dynamical friction, and ram pressure stripping. SatGen emulates cosmological zoom-in hydro-simulations in certain aspects. Satellites can reside in cored or cuspy DM subhaloes, depending on the halo response to baryonic physics that can be formulated from hydro-simulations and physical modeling. The subhalo profile and the stellar mass and size of a satellite evolves depending on its tidal mass loss and initial structure. The host galaxy can include a baryonic disc and a stellar bulge, each described by a density profile that allows analytic orbit integration. SatGen complements simulations by propagating the effect of halo response found in simulated field galaxies to satellites (not properly resolved in simulations) and outperforms simulations by sampling the halo-to-halo variance of satellite statistics and overcoming artificial disruption due to insufficient resolution. As a first application, we use the model to study satellites of Milky Way sized hosts, making it emulate simulations of bursty star formation and of smooth star formation, respectively, and to experiment with a disc potential in the host halo. Our model reproduces the observed satellite statistics reasonably well. Different physical recipes make a difference in satellite abundance and spatial distribution at the 25% level, not large enough to be distinguished by current observations given the halo-to-halo variance. The MW disc depletes satellites by 20% and has a subtle effect of diversifying the internal structure of satellites, important for alleviating certain small-scale problems.

We use three semi-analytic models (SAMs) of galaxy formation and evolution, run on the same 1$h^{-1}$Gpc MultiDark Planck2 cosmological simulation, to investigate the properties of [OII] emission line galaxies in the redshift range $0.6<z<1.2$. We compare model predictions with different observational data sets, including DEEP2--Firefly galaxies with absolute magnitudes. We estimate the [OII] luminosity, L[OII], using simple relations derived both from the models and observations and also using a public code. This code ideally uses as input instantaneous star formation rates (SFRs), which are only provided by one of the SAMs under consideration. We use this SAM to study the feasibility of inferring galaxies' L[OII] for models that only provide average SFRs. We find that the post-processing computation of L[OII] from average SFRs is accurate for model galaxies with dust attenuated L[OII]$\lesssim10^{42.2}$erg s$^{-1}$ ($<5\%$ discrepancy). We also explore how to derive the [OII] luminosity from simple relations using global properties usually output by SAMs. Besides the SFR, the model L[OII] is best correlated with the observed-frame $u$ and $g$ broad-band magnitudes. These correlations have coefficients (r-values) above 0.64 and a dispersion that varies with L[OII]. We use these correlations and an observational one based on SFR and metallicity to derive L[OII]. These relations result in [OII] luminosity functions and halo occupation distributions with shapes that vary depending on both the model and the method used. Nevertheless, for all the considered models, the amplitude of the clustering at scales above 1$h^{-1}$Mpc remains unchanged independently of the method used to derive L[OII].

Understanding the non-linear dynamics of satellite halos (a.k.a. ``sub-halos'') is important for predicting the abundance and distribution of dark matter substructures and satellite galaxies, and for distinguishing among microphysical dark matter models using observations. Typically, modeling these dynamics requires large N-body simulations with high resolution. Semi-analytic models can provide a more efficient way to describe the key physical processes such as dynamical friction, tidal mass loss, and tidal heating, with only a few free parameters. In this work, we present a fast Monte Carlo Markov Chain fitting approach to explore the parameter space of such a sub-halo non-linear evolution model. We use the dynamical models described in an earlier work and calibrate the models to two sets of high-resolution cold dark matter N-body simulations, ELVIS and Caterpillar. Compared to previous calibrations that used manual parameter tuning, our approach provides a more robust way to determine the best-fit parameters and their posterior probabilities. We find that jointly fitting for the sub-halo mass and maximum velocity functions can break the degeneracy between tidal stripping and tidal heating parameters, as well as providing better constraints on the strength of dynamical friction. We show that our semi-analytic simulation can accurately reproduce N-body simulations statistics, and that the calibration results for the two sets of N-body simulations agree at 95% confidence level. Dynamical models calibrated in this work will be important for future dark matter substructure studies.

The ability to accurately predict the evolution of tidally stripped haloes is important for understanding galaxy formation and testing the properties of dark matter. Most studies of substructure evolution make predictions based on empirical models of tidal mass loss that are calibrated using numerical simulations. This approach can be accurate in the cases considered, but lacks generality and does not provide a physical understanding of the processes involved. Recently, we demonstrated that truncating NFW distribution functions sharply in energy results in density profiles that resemble those of tidally stripped systems, offering a path to constructing physically motivated models of tidal mass loss. In this work, we review calculations of mass loss based on energy truncation alone, and then consider what secondary effects may modulate mass loss beyond this. We find that a combination of dependence on additional orbital parameters and variations in individual particle energies over an orbit results in a less abrupt truncation in energy space as a subhalo loses mass. Combining the energy truncation approach with a simple prediction for the mass-loss rate, we construct a full model of mass loss that can accurately predict the evolution of a subhalo in terms of a single parameter $\eta_{\rm eff}$. This parameter can be fully determined from the initial orbital and halo properties, and does not require calibration with numerical simulations.

We introduce LATIS, the Ly$\alpha$ Tomography IMACS Survey, a spectroscopic survey at Magellan designed to map the z=2.2-2.8 intergalactic medium (IGM) in three dimensions by observing the Ly$\alpha$ forest in the spectra of galaxies and QSOs. Within an area of 1.7 deg${}^2$, we will observe approximately half of $\gtrsim L^*$ galaxies at z=2.2-3.2 for typically 12 hours, providing a dense network of sightlines piercing the IGM with an average transverse separation of 2.5 $h^{-1}$ comoving Mpc (1 physical Mpc). At these scales, the opacity of the IGM is expected to be closely related to the dark matter density, and LATIS will therefore map the density field in the $z \sim 2.5$ universe at $\sim$Mpc resolution over the largest volume to date. Ultimately LATIS will produce approximately 3800 spectra of z=2.2-3.2 galaxies that probe the IGM within a volume of $4 \times 10^6 h^{-3}$ Mpc${}^3$, large enough to contain a representative sample of structures from protoclusters to large voids. Observations are already complete over one-third of the survey area. In this paper, we describe the survey design and execution. We present the largest IGM tomographic maps at comparable resolution yet made. We show that the recovered matter overdensities are broadly consistent with cosmological expectations based on realistic mock surveys, that they correspond to galaxy overdensities, and that we can recover structures identified using other tracers. LATIS is conducted in Canada-France-Hawaii Telescope Legacy Survey fields, including COSMOS. Coupling the LATIS tomographic maps with the rich data sets collected in these fields will enable novel studies of environment-dependent galaxy evolution and the galaxy-IGM connection at cosmic noon.

Simulations of cold dark matter make robust predictions about the slope and normalization of the dark matter halo and subhalo mass functions on small scales. Recent observational advances utilizing strong gravitational lensing have demonstrated the ability of this technique to place constraints on these quantities on subgalactic scales corresponding to dark matter halo masses of $10^6$--$10^9\mathrm{M}_\odot$. On these scales the physics of baryons, which make up around 17% of the matter content of the Universe but which are not included in pure dark matter N-body simulations, are expected to affect the growth of structure and the collapse of dark matter halos. In this work we develop a semi-analytic model to predict the amplitude and slope of the dark matter halo and subhalo mass functions on subgalactic scales in the presence of baryons. We find that the halo mass function is suppressed by up to 25%, and the slope is modified, ranging from $-1.916$ to $-1.868$ in this mass range. These results are consistent with current measurements, but differ sufficiently from the expectations for a dark matter only universe that it may be testable in the near future.

A key obstacle to developing a satisfying theory of galaxy evolution is the difficulty in extending analytic descriptions of early structure formation into full nonlinearity, the regime in which galaxy growth occurs. Extant techniques, though powerful, are based on approximate numerical methods whose Monte Carlo-like nature hinders intuition building. Here, we develop a new solution to this problem and its empirical validation. We first derive closed-form analytic expectations for the evolution of fixed percentiles in the real-space cosmic density distribution, {\it averaged over representative volumes observers can track cross-sectionally\}. Using the Lagrangian forms of the fluid equations, we show that percentiles in $\delta$---the density relative to the median---should grow as $\delta(t)\propto\delta_{0}^{\alpha}\,t^{\beta}$, where $\alpha\equiv2$ and $\beta\equiv2$ for Newtonian gravity at epochs after the overdensities transitioned to nonlinear growth. We then use 9.5 sq. deg. of Carnegie-Spitzer-IMACS Redshift Survey data to map {\it galaxy\} environmental densities over $0.2<z<1.5$ ($\sim$7 Gyr) and infer $\alpha=1.98\pm0.04$ and $\beta=2.01\pm0.11$---consistent with our analytic prediction. These findings---enabled by swapping the Eulerian domain of most work on density growth for a Lagrangian approach to real-space volumetric averages---provide some of the strongest evidence that a lognormal distribution of early density fluctuations indeed decoupled from cosmic expansion to grow through gravitational accretion. They also comprise the first exact, analytic description of the nonlinear growth of structure extensible to (arbitrarily) low redshift. We hope these results open the door to new modeling of, and insight-building into, the galaxy growth and its diversity in cosmological contexts.

The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work.

The free-streaming length of dark matter depends on fundamental dark matter physics, and determines the abundance and concentration of dark matter halos on sub-galactic scales. Using the image positions and flux ratios from eight quadruply-imaged quasars, we constrain the free-streaming length of dark matter and the amplitude of the subhalo mass function (SHMF). We model both main deflector subhalos and halos along the line of sight, and account for warm dark matter (WDM) free-streaming effects on the mass function and mass-concentration relation. By calibrating the scaling of the SHMF with host halo mass and redshift using a suite of simulated halos, we infer a global normalization for the SHMF. We account for finite-size background sources, and marginalize over the mass profile of the main deflector. Parameterizing dark matter free-streaming through the half-mode mass $m_{\rm{hm}}$, we constrain the thermal relic particle mass $m_{\rm{DM}}$ corresponding to $m_{\rm{hm}}$. At $95 \%$ CI: $m_{\rm{hm}} < 10^{7.8} M_{\odot}$ ($m_{\rm{DM}} > 5.2 \ \rm{keV}$). We disfavor $m_{\rm{DM}} = 4.0 \rm{keV}$ and $ m_{\rm{DM}} = 3.0 \rm{keV}$ with likelihood ratios of 7:1 and 30:1, respectively, relative to the peak of the posterior distribution. Assuming cold dark matter, we constrain the projected mass in substructure between $10^6 - 10^{9} M_{\odot}$ near lensed images. At $68 \%$ CI, we infer $2.0 - 6.1 \times 10^{7} M_{\odot} \rm{kpc^{-2}}$, corresponding to mean projected mass fraction $\bar{f}_{\rm{sub}} = 0.035_{-0.017}^{+0.021}$. At $95 \%$ CI, we obtain a lower bound on the projected mass of $0.6 \times 10^{7} M_{\odot} \rm{kpc^{-2}}$, corresponding to $\bar{f}_{\rm{sub}} > 0.005$. These results agree with the predictions of cold dark matter.

We perform a simulation with Galacticus, a semi-analytical galaxy formation model, to predict the number counts of H$\alpha$ and [OIII] emitting galaxies. With a state-of-the-art N-body simulation, UNIT, we first calibrate Galacticus with the current observation of H$\alpha$ luminosity function. The resulting model coupled with a dust attenuation model, can reproduce the current observations, including the H$\alpha$ luminosity function from HiZELS and number density from WISP. We extrapolate the model prediction to higher redshift and the result is found to be consistent with previous investigations. We then use the same galaxy formation model to predict the number counts for [OIII] emitting galaxies. The result provides further validation of our galaxy formation model and dust model. We present number counts of H$\alpha$ and [OIII] emission line galaxies for three different line flux limits: $5\times10^{-17}$erg/s/cm$^{2}$, $1\times10^{-16}$erg/s/cm$^{2}$ (6.5$\sigma$ nominal depth for WFIRST GRS), and $2\times10^{-16}$erg/s/cm$^{2}$ (3.5$\sigma$ depth of Euclid GRS). At redshift $2<z<3$, our model predicts that WFIRST can observe hundreds of [OIII] emission line galaxies per square degree with a line flux limit of $1\times10^{-16}$erg/s/cm$^{2}$. This will provide accurate measurement of large scale structure to probe dark energy over a huge cosmic volume to an unprecedented high redshift. Finally, we compare the flux ratio of H$\alpha$/[OIII] within the redshift range of $0<z<3$. Our results show the known trend of increasing H$\alpha$/[O III] flux ratio with H$\alpha$ flux at low redshift, which becomes a weaker trend at higher redshifts.

As part of the effort to meet the needs of the Large Synoptic Survey Telescope Dark Energy Science Collaboration (LSST DESC) for accurate, realistically complex mock galaxy catalogs, we have developed GalSampler, an open-source python package that assists in generating large volumes of synthetic cosmological data. The key idea behind GalSampler is to recast hydrodynamical simulations and semi-analytic models as physically-motivated galaxy libraries. GalSampler populates a new, larger-volume halo catalog with galaxies drawn from the baseline library; by using weighted sampling guided by empirical modeling techniques, GalSampler inherits statistical accuracy from the empirical model and physically-motivated complexity from the baseline library. We have recently used GalSampler to produce the cosmoDC2 extragalactic catalog made for the LSST DESC Data Challenge 2. Using cosmoDC2 as a guiding example, we outline how GalSampler can continue to support ongoing and near-future galaxy surveys such as the Dark Energy Survey (DES), the Dark Energy Spectroscopic Instrument (DESI), WFIRST, and Euclid.

The mass-concentration relation of dark matter halos reflects the assembly history of objects in hierarchical structure formation scenarios, and depends on fundamental quantities in cosmology such as the slope of the primordial matter power-spectrum. This relation is unconstrained by observations on sub-galactic scales. We derive the first measurement of the mass-concentration relation using the image positions and flux ratios from eleven quadruple-image strong gravitational lenses (quads) in the mass range $10^{6} - 10^{10} M_{\odot}$, assuming cold dark matter. Our analysis framework includes both subhalos and line of sight halos, marginalizes over nuisance parameters describing the lens macromodel, accounts for finite source effects on lensing observables, and simultaneously constrains the normalization and logarithmic slope of the mass-concentration relation, and the normalization of the subhalo mass function. At $z=0$, we constrain the concentration of $10^{8} M_{\odot}$ halos $c=12_{-5}^{+6}$ at $68 \%$ CI, and $c=12_{-9}^{+15}$ at $95 \%$ CI. For a $10^{7} M_{\odot}$ halo, we obtain $68 \%$ ($95 \%$) constraints $c=15_{-8}^{+9}$ ($c=15_{-11}^{+18}$), while for $10^{9} M_{\odot}$ halos $c=10_{-4}^{+7}$ ($c=10_{-7}^{+14}$). These results are consistent with the theoretical predictions from mass-concentration relations in the literature, and establish strong lensing by galaxies as a powerful probe of halo concentrations on sub-galactic scales across cosmological distance.

In the cold dark matter (CDM) picture of structure formation, galaxy mass distributions are predicted to have a considerable amount of structure on small scales. Strong gravitational lensing has proven to be a useful tool for studying this small-scale structure. Much of the attention has been given to detecting individual dark matter subhalos through lens modeling, but recent work has suggested that the full population of subhalos could be probed using a power spectrum analysis. In this paper we quantify the power spectrum of small-scale structure in simulated galaxies, with the goal of understanding theoretical predictions and setting the stage for using measurements of the power spectrum to test dark matter models. We use a sample of simulated galaxies generated from the \texttt{Galacticus} semi-analytic model to determine the power spectrum distribution first in the CDM paradigm and then in a warm dark matter scenario. We find that a measurement of the slope and amplitude of the power spectrum on galaxy strong lensing scales ($k\sim 1$ kpc$^{-1}$) could be used to distinguish between CDM and alternate dark matter models, especially if the most massive subhalos can be directly detected via gravitational imaging.

This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data.

Strong lensing provides a powerful means of investigating the nature of dark matter as it probes dark matter structure on sub-galactic scales. We present an extension of a forward modeling framework that uses flux ratios from quadruply imaged quasars (quads) to measure the shape and amplitude of the halo mass function, including line of sight (LOS) halos and main deflector subhalos. We apply this machinery to 50 mock lenses --- roughly the number of known quads --- with warm dark matter (WDM) mass functions exhibiting free-streaming cutoffs parameterized by the half-mode mass $m_{\rm{hm}}$. Assuming cold dark matter (CDM), we forecast bounds on $m_{\rm{hm}}$ and the corresponding thermal relic particle masses over a range of tidal destruction severity, assuming a particular WDM mass function and mass-concentration relation. With significant tidal destruction, at $2 \sigma$ we constrain $m_{\rm{hm}}<10^{7.9} \left(10^{8.4}\right) M_{\odot}$, or a 4.4 (3.1) keV thermal relic, with image flux uncertainties from measurements and lens modeling of $2\% \left(6\%\right)$. With less severe tidal destruction we constrain $m_{\rm{hm}}<10^{7} \left(10^{7.4}\right) M_{\odot}$, or an 8.2 (6.2) keV thermal relic. If dark matter is warm, with $m_{\rm{hm}} = 10^{7.7} M_{\odot}$ (5.1 keV), we would favor WDM with $m_{\rm{hm}} > 10^{7.7} M_{\odot}$ over CDM with relative likelihoods of 22:1 and 8:1 with flux uncertainties of $2\%$ and $6\%$, respectively. These bounds improve over those obtained by modeling only main deflector subhalos because LOS objects produce additional flux perturbations, especially for high redshift systems. These results indicate that $\sim 50$ quads can conclusively differentiate between warm and cold dark matter.

We present recalibrations of the GALFORM semi-analytical model of galaxy formation in a new N-body simulation with the Planck cosmology. The Planck Millennium simulation uses more than 128 billion particles to resolve the matter distribution in a cube of $800$ Mpc on a side, which contains more than 77 million dark matter haloes with mass greater than $2.12 \times 10^{9} h^{-1} {\rm M_{\odot}}$ at the present day. Only minor changes to a very small number of model parameters are required in the recalibration. We present predictions for the atomic hydrogen content (HI) of dark matter halos, which is a key input into the calculation of the HI intensity mapping signal expected from the large-scale structure of the Universe. We find that the HI mass $-$ halo mass relation displays a clear break at the halo mass above which AGN heating suppresses gas cooling, $\approx 3 \times 10^{11} h^{-1} M_{\rm \odot}$. Below this halo mass, the HI content of haloes is dominated by the central galaxy; above this mass it is the combined HI content of satellites that prevails. We find that the HI mass - halo mass relation changes little with redshift up to $z=3$. The bias of HI sources shows a scale dependence that gets more pronounced with increasing redshift.

We forecast the linear bias for H${\rm \alpha}$-emitting galaxies at high redshift. To simulate a Euclid-like and a WFIRST-like survey, we place galaxies into a large-volume dark matter halo lightcone by sampling a library of luminosity-dependent halo occupation distributions (HODs), which is constructed using a physically motivated galaxy formation model. We calibrate the dust attenuation in the lightcones such that they are able to reproduce the H{\alpha} luminosity function or the H{\alpha} cumulative number counts. The angle-averaged galaxy correlation function is computed for each survey in redshift slices of width $\Delta z = 0.2$. In each redshift bin the linear bias can be fitted with a single, scale-independent value that increases with increasing redshift. Fitting for the evolution of linear bias with redshift, we find that our Euclid-like and WFIRST-like surveys are both consistent within error with the relation $b(z) = 0.7z + 0.7$. Our bias forecasts are consistent with bias measurements from the HiZELS survey. We find that the Euclid-like and WFIRST-like surveys yield linear biases that are broadly consistent within error, most likely due to the HOD for the WFIRST-like survey having a steeper power-law slope towards larger halo masses.

We study halo mass functions with the Press-Schechter formalism for interacting dark matter models, where matter power spectra are damped due to dark acoustic oscillations in the early universe. After adopting a smooth window function, we calibrate the analytical model with numerical simulations from the "effective theory of structure formation" (ETHOS) project and fix the model parameters in the high mass regime, $M_{\rm h}\gtrsim3\times10^{10}\;{\rm M}_{\odot}$. We also perform high-resolution cosmological simulations with halo masses down to $M_{\rm h}\sim10^8\;{\rm M}_{\odot}$ to cover a wide mass range for comparison. Although the model is calibrated with ETHOS1 and CDM simulations for high halo masses at redshift $z=0$, it successfully reproduces simulations for two other ETHOS models in the low mass regime at low and high redshifts. As an application, we compare the cumulative number density of haloes to that of observed galaxies at $z=6$, and find the interacting dark matter models with a kinetic decoupling temperature below $0.5\ \rm{keV}$ is disfavored. We also perform the abundance-matching analysis and derive the stellar-halo mass relation for these models at $z=4$. Suppression in halo abundance leads to less massive haloes that host observed galaxies in the stellar mass range $M_*\simeq 10^5-10^7\ {\rm M}_{\odot}$.

We study the shape of the gas-phase mass-metallicity relation (MZR) of a combined sample of present-day dwarf and high-mass star-forming galaxies using IZI, a Bayesian formalism for measuring chemical abundances presented in Blanc et al. 2015. We observe a characteristic stellar mass scale at $M_* \simeq 10^{9.5}$M$_{\odot}$, above which the ISM undergoes a sharp increase in its level of chemical enrichment. In the $10^{6}-10^{9.5}$M$_{\odot}$ range the MZR follows a shallow power-law ($Z\propto M^{\alpha}_*$) with slope $\alpha=0.14\pm0.08$. At approaching $M_* \simeq 10^{9.5}$M$_{\odot}$ the MZR steepens significantly, showing a slope of $\alpha=0.37\pm0.08$ in the $10^{9.5}-10^{10.5}$M$_{\odot}$ range, and a flattening towards a constant metallicity at higher stellar masses. This behavior is qualitatively different from results in the literature that show a single power-law MZR towards the low mass end. We thoroughly explore systematic uncertainties in our measurement, and show that the shape of the MZR is not induced by sample selection, aperture effects, a changing N/O abundance, the adopted methodology used to construct the MZR, secondary dependencies on star formation activity, nor diffuse ionized gas (DIG) contamination, but rather on differences in the method used to measure abundances. High resolution hydrodynamical simulations can qualitatively reproduce our result, and suggest a transition in the ability of galaxies to retain their metals for stellar masses above this threshold. The MZR characteristic mass scale also coincides with a transition in the scale height and clumpiness of cold gas disks, and a typical gas fraction below which the efficiency of star formation feedback for driving outflows is expected to decrease sharply.

The most massive and luminous galaxies in the Universe serve as powerful probes to study the formation of structure, the assembly of mass, and cosmology. However, their detailed formation and evolution is still barely understood. Here we extract a sample of massive mock galaxies from the semi-analytical model of galaxy formation (SAM) GALACTICUS from the MultiDark-Galaxies, by replicating the CMASS photometric selection from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). The comparison of the GALACTICUS CMASS-mock with BOSS-CMASS data allows us to explore different aspects of the massive galaxy population at $0.5<z<0.6$, including the galaxy-halo connection and the galaxy clustering. We find good agreement between our modelled galaxies and observations regarding the galaxy-halo connection, but our CMASS-mock over-estimates the clustering amplitude of the 2-point correlation function, due to a smaller number density compared to BOSS, a lack of blue objects, and a small intrinsic scatter in stellar mass at fixed halo mass of $<0.1$ dex. To alleviate this problem, we construct an alternative mock catalogue mimicking the CMASS colour-magnitude distribution by randomly down-sampling the SAM catalogue. This CMASS-mock reproduces the clustering of CMASS galaxies within 1$\sigma$ and shows some environmental dependency of star formation properties that could be connected to the quenching of star formation and the assembly bias.

We apply the model relating halo concentration to formation history proposed by Ludlow et al. to merger trees generated using an algorithm based on excursion set theory. We find that while the model correctly predicts the median relation between halo concentration and mass, it underpredicts the scatter in concentration at fixed mass. Since the same model applied to N-body merger trees predicts the correct scatter, we postulate that the missing scatter is due to the lack of any environmental dependence in merger trees derived from excursion set theory. We show that a simple modification to the merger tree construction algorithm, which makes merger rates dependent on environment, can increase the scatter by the required amount, and simultaneously provide a qualitatively correct correlation between environment and formation epoch in the excursion set merger trees.

We calculate the net extinction of galactic light as a function of wavelength, inclination, central optical depth, and morphology for simple galactic geometries using the Hyperion radiative transfer code. Compared to previous, similar works we tabulate extinction over a much broader range of galactic properties, and using a much finer grid in the model parameters. We expect these results to be useful for constructing dust-extinguished spectra and luminosities of model galaxies and, therefore, for synthetic survey building. Results are made available as an HDF5 file at https://doi.org/10.5281/zenodo.1442826

We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic fractions and gas scaling relations at redshift z = 0. However, there are still some (model-dependent) differences, such as central galaxies being too massive, and galaxy colours (g - r) being bluer (about 0.2 dex lower at the peak position) than in observations. The agreement in gas scaling relations down to 10^{13} h^{-1} M_{\odot} between the simulations indicates that particulars of the sub-grid modelling of the baryonic physics only has a weak influence on these relations. We also include - where appropriate - a comparison to three semi-analytical galaxy formation models as applied to the same underlying dark-matter-only simulation. All simulations and derived data products are publicly available.

We compute covariance matrices for many observed estimates of the stellar mass function of galaxies from $z=0$ to $z\approx 4$, and for one estimate of the projected correlation function of galaxies split by stellar mass at $z\lesssim 0.5$. All covariance matrices include contributions due to large scale structure, the preference for galaxies to be found in groups and clusters, and for shot noise. These covariance matrices are made available for use in constraining models of galaxy formation and the galaxy-halo connection.

We present a comparison of the observed evolving galaxy stellar mass functions with the predictions of eight semi-analytic models and one halo occupation distribution model. While most models are able to fit the data at low redshift, some of them struggle to simultaneously fit observations at high redshift. We separate the galaxies into 'passive' and 'star-forming' classes and find that several of the models produce too many low-mass star-forming galaxies at high redshift compared to observations, in some cases by nearly a factor of 10 in the redshift range $2.5 < z < 3.0$. We also find important differences in the implied mass of the dark matter haloes the galaxies inhabit, by comparing with halo masses inferred from observations. Galaxies at high redshift in the models are in lower mass haloes than suggested by observations, and the star formation efficiency in low-mass haloes is higher than observed. We conclude that many of the models require a physical prescription that acts to dissociate the growth of low-mass galaxies from the growth of their dark matter haloes at high redshift.

Cosmic acceleration is the most surprising cosmological discovery in many decades. Testing and distinguishing among possible explanations requires cosmological measurements of extremely high precision probing the full history of cosmic expansion and structure growth and, ideally, compare and contrast matter and relativistic tracers of the gravity potential. This program is one of the defining objectives of the Wide-Field Infrared Survey Telescope (WFIRST), as set forth in the New Worlds, New Horizons report (NWNH) in 2010. The WFIRST mission has the ability to improve these measurements by 1-2 orders of magnitude compared to the current state of the art, while simultaneously extending their redshift grasp, greatly improving control of systematic effects, and taking a unified approach to multiple probes that provide complementary physical information and cross-checks of cosmological results. We describe in this annual report the activities of the Science Investigation Team (SIT) "Cosmology with the High Latitude Survey (HLS)" during the year 2017. This team was selected by NASA in December 2015 in order to address the stringent challenges of the WFIRST dark energy (DE) program through the Project's formulation phase. This SIT has elected to jointly address Galaxy Redshift Survey, Weak Lensing and Cluster Growth and thus fully embrace the fact that the imaging and spectroscopic elements of the HLS will be realized as an integrated observing program, and they jointly impose requirements on performance and operations. WFIRST is designed to be able to deliver a definitive result on the origin of cosmic acceleration. It is not optimized for Figure of Merit sensitivity but for control of systematic uncertainties and for having multiple techniques each with multiple cross-checks. Our SIT work focuses on understanding the potential systematics in the WFIRST DE measurements.

The use of high-quality simulated sky catalogs is essential for the success of cosmological surveys. The catalogs have diverse applications, such as investigating signatures of fundamental physics in cosmological observables, understanding the effect of systematic uncertainties on measured signals and testing mitigation strategies for reducing these uncertainties, aiding analysis pipeline development and testing, and survey strategy optimization. The list of applications is growing with improvements in the quality of the catalogs and the details that they can provide. Given the importance of simulated catalogs, it is critical to provide rigorous validation protocols that enable both catalog providers and users to assess the quality of the catalogs in a straightforward and comprehensive way. For this purpose, we have developed the DESCQA framework for the Large Synoptic Survey Telescope Dark Energy Science Collaboration as well as for the broader community. The goal of DESCQA is to enable the inspection, validation, and comparison of an inhomogeneous set of synthetic catalogs via the provision of a common interface within an automated framework. In this paper, we present the design concept and first implementation of DESCQA. In order to establish and demonstrate its full functionality we use a set of interim catalogs and validation tests. We highlight several important aspects, both technical and scientific, that require thoughtful consideration when designing a validation framework, including validation metrics and how these metrics impose requirements on the synthetic sky catalogs.

Knowledge of the number density of H$\alpha$ emitting galaxies is vital for assessing the scientific impact of the Euclid and WFIRST missions. In this work we present predictions from a galaxy formation model, Galacticus, for the cumulative number counts of H$\alpha$-emitting galaxies. We couple Galacticus to three different dust attenuation methods and examine the counts using each method. A $\chi^2$ minimisation approach is used to compare the model predictions to observed galaxy counts and calibrate the dust parameters. We find that weak dust attenuation is required for the Galacticus counts to be broadly consistent with the observations, though the optimum dust parameters return large values for $\chi^2$, suggesting that further calibration of Galacticus is necessary. The model predictions are also consistent with observed estimates for the optical depth and the H$\alpha$ luminosity function. Finally we present forecasts for the redshift distributions and number counts for two Euclid-like and one WFIRST-like survey. For a Euclid-like survey with redshift range $0.9\leqslant z\leqslant 1.8$ and H$\alpha+{\rm [NII]}$ blended flux limit of $2\times 10^{-16}{\rm erg}\,{\rm s}^{-1}\,{\rm cm}^{-2}$ we predict a number density between 3900--4800 galaxies per square degree. For a WFIRST-like survey with redshift range $1\leqslant z\leqslant 2$ and blended flux limit of $1\times 10^{-16}{\rm erg}\,{\rm s}^{-1}\,{\rm cm}^{-2}$ we predict a number density between 10400--15200 galaxies per square degree.

We present a comparison of nine galaxy formation models, eight semi-analytical and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of co-moving width 125$h^{-1}$ Mpc, with a dark-matter particle mass of $1.24\times 10^9 h^{-1}$ Msun) and the same merger trees. While their free parameters have been calibrated to the same observational data sets using two approaches, they nevertheless retain some 'memory' of any previous calibration that served as the starting point (especially for the manually-tuned models). For the first calibration, models reproduce the observed z = 0 galaxy stellar mass function (SMF) within 3-{\sigma}. The second calibration extended the observational data to include the z = 2 SMF alongside the z~0 star formation rate function, cold gas mass and the black hole-bulge mass relation. Encapsulating the observed evolution of the SMF from z = 2 to z = 0 is found to be very hard within the context of the physics currently included in the models. We finally use our calibrated models to study the evolution of the stellar-to-halo mass (SHM) ratio. For all models we find that the peak value of the SHM relation decreases with redshift. However, the trends seen for the evolution of the peak position as well as the mean scatter in the SHM relation are rather weak and strongly model dependent. Both the calibration data sets and model results are publicly available.

We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities - noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite mass resolution of the N-body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models - in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by on-going and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available.

Dwarf galaxies are known to have remarkably low star formation efficiency due to strong feedback. Adopting the dwarf galaxies of the Milky Way as a laboratory, we explore a flexible semi-analytic galaxy formation model to understand how the feedback processes shape the satellite galaxies of the Milky Way. Using Markov-Chain Monte-Carlo, we exhaustively search a large parameter space of the model and rigorously show that the general wisdom of strong outflows as the primary feedback mechanism cannot simultaneously explain the stellar mass function and the mass--metallicity relation of the Milky Way satellites. An extended model that assumes that a fraction of baryons is prevented from collapsing into low-mass halos in the first place can be accurately constrained to simultaneously reproduce those observations. The inference suggests that two different physical mechanisms are needed to explain the two different data sets. In particular, moderate outflows with weak halo mass dependence are needed to explain the mass--metallicity relation, and prevention of baryons falling into shallow gravitational potentials of low-mass halos (e.g. "pre-heating") is needed to explain the low stellar mass fraction for a given subhalo mass.

Using a cosmological $N$-body simulation, we investigate the origin and distribution of stars in the intracluster light (ICL) of a Fornax-like cluster. In a dark matter only simulation we identify a halo which, at $z=0$, has $M_200 \simeq 4.1 \times 10^{13}M_{sun}$ and $r_{200} = 700kpc$, and replace infalling subhalos with models that include spheroid and disc components. As they fall into the cluster, the stars in some of these galaxies are stripped from their hosts, and form the ICL. We consider the separate contributions to the ICL from stars which originate in the haloes and the discs of the galaxies. We find that disc ICL stars are more centrally concentrated than halo ICL stars. The majority of the disc ICL stars are associated with one initially disc-dominated galaxy that falls to the centre of the cluster and is heavily disrupted, producing part of the cD galaxy. At radial distances greater than 200kpc, well beyond the stellar envelope of the cD galaxy, stars formerly from the stellar haloes of galaxies dominate the ICL. Therefore at large distances, the ICL population is dominated by older stars.

Any measurement made using an N-body simulation is subject to noise due to the finite number of particles used to sample the dark matter distribution function, and the lack of structure below the simulation resolution. This noise can be particularly significant when attempting to measure intrinsically small quantities, such as halo spin. In this work we develop a model to describe the effects of particle noise on halo spin parameters. This model is calibrated using N-body simulations in which the particle noise can be treated as a Poisson process on the underlying dark matter distribution function, and we demonstrate that this calibrated model reproduces measurements of halo spin parameter error distributions previously measured in N-body convergence studies. Utilizing this model, along with previous measurements of the distribution of halo spin parameters in N-body simulations, we place constraints on the noise-free distribution of halo spins. We find that the noise-free median spin is 3% lower than that measured directly from the N-body simulation, corresponding to a shift of approximately 40 times the statistical uncertainty in this measurement arising purely from halo counting statistics. We also show that measurement of the spin of an individual halo to 10% precision requires at least $4\times 10^4$ particles in the halo - for halos containing 200 particles the fractional error on spins measured for individual halos is of order unity.

We present a clustering comparison of 12 galaxy formation models (including Semi-Analytic Models (SAMs) and Halo Occupation Distribution (HOD) models) all run on halo catalogues and merger trees extracted from a single {\Lambda}CDM N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the 2-Point Correlation Functions (2PCF). We also study the implications of the different treatments of orphan (galaxies not assigned to any dark matter subhalo) and non-orphan galaxies in these measurements. Our main result is that the galaxy formation models generally agree in their clustering predictions but they disagree significantly between HOD and SAMs for the orphan satellites. Although there is a very good agreement between the models on the 2PCF of central galaxies, the scatter between the models when orphan satellites are included can be larger than a factor of 2 for scales smaller than 1 Mpc/h. We also show that galaxy formation models that do not include orphan satellite galaxies have a significantly lower 2PCF on small scales, consistent with previous studies. Finally, we show that the 2PCF of orphan satellites is remarkably different between SAMs and HOD models. Orphan satellites in SAMs present a higher clustering than in HOD models because they tend to occupy more massive haloes. We conclude that orphan satellites have an important role on galaxy clustering and they are the main cause of the differences in the clustering between HOD models and SAMs.

We propose a new method for generating equilibrium models of spherical systems of collisionless particles that are finite in extent, but whose central regions resemble dark matter halos from cosmological simulations. This method involves iteratively removing unbound particles from a Navarro-Frenk-White profile truncated sharply at some radius. The resulting models are extremely stable, and thus provide a good starting point for N-body simulations of isolated halos. We provide a code to generate such models for NFW and a variety of other common density profiles. We then develop an analytic approximation to this truncated distribution function. Our method proceeds by analogy with the King model, truncating and shifting the original distribution function of an infinitely extended Navarro-Frenk-White profile in energy space. We show that the density profiles of our models closely resemble the tidally truncated density profiles seen previously in studies of satellite evolution. Pursuing this analogy further with a series of simulations of tidal mass loss, we find that our models provide a good approximation to the full distribution function of tidally stripped systems, thus allowing theoretically motivated phase-space calculations for such systems.

The study of cosmology, galaxy formation and exoplanets has now advanced to a stage where a cosmic inventory of terrestrial planets may be attempted. By coupling semi-analytic models of galaxy formation to a recipe that relates the occurrence of planets to the mass and metallicity of their host stars, we trace the population of terrestrial planets around both solar-mass (FGK type) and lower-mass (M dwarf) stars throughout all of cosmic history. We find that the mean age of terrestrial planets in the local Universe is $7\pm{}1$ Gyr for FGK hosts and $8\pm{}1$ Gyr for M dwarfs. We estimate that hot Jupiters have depleted the population of terrestrial planets around FGK stars by no more than $\approx 10\%$, and that only $\approx 10\%$ of the terrestrial planets at the current epoch are orbiting stars in a metallicity range for which such planets have yet to be confirmed. The typical terrestrial planet in the local Universe is located in a spheroid-dominated galaxy with a total stellar mass comparable to that of the Milky Way. When looking at the inventory of planets throughout the whole observable Universe, we argue for a total of $\approx 1\times 10^{19}$ and $\approx 5\times 10^{20}$ terrestrial planets around FGK and M stars, respectively. Due to light travel time effects, the terrestrial planets on our past light cone exhibit a mean age of just $1.7\pm 0.2$ Gyr. These results are discussed in the context of cosmic habitability, the Copernican principle and searches for extraterrestrial intelligence at cosmological distances.

Here we explore the evolution of galaxy ensembles at early times by writing the in situ stellar mass growth of galaxies purely as a stationary stochastic (e.g., quasi-steady state) process. By combining the mathematics of such processes with Newtonian gravity and a mean local star formation efficiency, we show that the stellar mass evolution of galaxy ensembles is directly related to the average acceleration of baryons onto dark matter halos at the onset of star formation, with explicit dependencies on initial local matter densities and halo mass. The density term specifically implies more rapid average rates of growth in higher density regions of the universe compared to low density regions, i.e., assembly bias. With this framework, using standard cosmological parameters, a mean star formation efficiency derived by other authors, and knowledge of the shape of the cosmological matter power spectrum at small scales, we analytically derive (1) the characteristic stellar masses of galaxies (M*), (2) the power-law low-mass slope (alpha) and normalization (phi*) of the stellar mass function, and (3) the evolution of the stellar mass function in time over 12.5 > z > 2. Correspondingly, the rise in the cosmic star formation rate density over these epochs, while the universe can sustain unabated fueling of star formation, also emerges naturally. All of our findings are consistent with the deepest available data, including the expectation of alpha~-7/5; i.e., a stellar mass function low-mass slope that is notably shallower than that of the halo mass function, and with no systematic deviations from a mean star formation efficiency with density or mass, nor any explicit, additional feedback mechanisms. These derivations yield a compelling richness and complexity but also show that very few astrophysical details are required to understand the evolution of cosmic ensemble of galaxies at early times.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

A common approach in semi-analytic modeling of galaxy formation is to construct Monte Carlo realizations of merger histories of dark matter halos whose masses are sampled from a halo mass function. Both the mass function itself, and the merger rates used to construct merging histories are calibrated to N-body simulations. Typically, "backsplash" halos (those which were once subhalos within a larger halo, but which have since moved outside of the halo) are counted in both the halo mass function, and in the merger rates (or, equivalently, progenitor mass functions). This leads to a double-counting of mass in Monte Carlo merger histories which will bias results relative to N-body results. We measure halo mass functions and merger rates with this double-counting removed in a large, cosmological N-body simulation with cosmological parameters consistent with current constraints. Furthermore, we account for the inherently noisy nature of N-body halo mass estimates when fitting functions to N-body data, and show that ignoring these errors leads to a significant systematic bias given the precision statistics available from state-of-the-art N-body cosmological simulations.

Accurate modeling of galaxy formation in a hierarchical, cold dark matter universe requires the use of sufficiently high-resolution merger trees to obtain convergence in the predicted properties of galaxies. When semi-analytic galaxy formation models are applied to cosmological N-body simulation merger trees, it is often the case that those trees have insufficient resolution to give converged galaxy properties. We demonstrate a method to augment the resolution of N-body merger trees by grafting in branches of Monte Carlo merger trees with higher resolution, but which are consistent with the pre-existing branches in the N-body tree. We show that this approach leads to converged galaxy properties.

We present a new version of the GALFORM semi-analytical model of galaxy formation. This brings together several previous developments of GALFORM into a single unified model, including a different initial mass function (IMF) in quiescent star formation and in starbursts, feedback from active galactic nuclei supressing gas cooling in massive halos, and a new empirical star formation law in galaxy disks based on their molecular gas content. In addition, we have updated the cosmology, introduced a more accurate treatment of dynamical friction acting on satellite galaxies, and updated the stellar population model. The new model is able to simultaneously explain both the observed evolution of the K-band luminosity function and stellar mass function, and the number counts and redshift distribution of sub-mm galaxies selected at 850 mu. This was not previously achieved by a single physical model within the LambdaCDM framework, but requires having an IMF in starbursts that is somewhat top-heavy. The new model is tested against a wide variety of observational data covering wavelengths from the far-UV to sub-mm, and redshifts from z=0 to z=6, and is found to be generally successful. These observations include the optical and near-IR luminosity functions, HI mass function, fraction of early type galaxies, Tully-Fisher, metallicity-luminosity and size-luminosity relations at z=0, as well as far-IR number counts, and far-UV luminosity functions at z ~ 3-6. [abridged]

Many properties of the Milky Way's dark matter halo, including its mass assembly history, concentration, and subhalo population, remain poorly constrained. We explore the connection between these properties of the Milky Way and its satellite galaxy population, especially the implication of the presence of the Magellanic Clouds for the properties of the Milky Way halo. Using a suite of high-resolution $N$-body simulations of Milky Way-mass halos with a fixed final Mvir ~ 10^{12.1}Msun, we find that the presence of Magellanic Cloud-like satellites strongly correlates with the assembly history, concentration, and subhalo population of the host halo, such that Milky Way-mass systems with Magellanic Clouds have lower concentration, more rapid recent accretion, and more massive subhalos than typical halos of the same mass. Using a flexible semi-analytic galaxy formation model that is tuned to reproduce the stellar mass function of the classical dwarf galaxies of the Milky Way with Markov-Chain Monte-Carlo, we show that adopting host halos with different mass-assembly histories and concentrations can lead to different best-fit models for galaxy-formation physics, especially for the strength of feedback. These biases arise because the presence of the Magellanic Clouds boosts the overall population of high-mass subhalos, thus requiring a different stellar-mass-to-halo-mass ratio to match the data. These biases also lead to significant differences in the mass--metallicity relation, the kinematics of low-mass satellites, the number counts of small satellites associated with the Magellanic Clouds, and the stellar mass of Milky Way itself. Observations of these galaxy properties can thus provide useful constraints on the properties of the Milky Way halo.

We propose a common terminology for use in describing both temporal merger trees and spatial structure trees for dark-matter halos. We specify a unified data format in HDF5 and provide example I/O routines in C, FORTRAN and PYTHON.

We present the results of a Keck-ESI study of dwarf galaxies across a range of environment: the Perseus Cluster, the Virgo Cluster, the NGC 1407 group, and the NGC 1023 group. Eighteen dEs are targeted for spectroscopy, three for the first time. We confirm cluster membership for one Virgo dE, and group membership for one dE in the NGC 1023 group, and one dE in the NGC 1407 group for the first time. Regardless of environment, the dEs follow the same size-magnitude and $\sigma$-luminosity relation. Two of the Virgo dwarfs, VCC 1199 and VCC 1627, have among the highest central velocity dispersions ($\sigma_{0}$ = 58.4 km s$^{-1}$ and 49.2 km s$^{-1}$) measured for dwarfs of their luminosity ($M_{R}\approx -17$). Given their small sizes ($R_{e} < 300$ pc) and large central velocity dispersions, we classify these two dwarfs as compact ellipticals rather than dEs. Group dEs typically have higher mean dynamical-to-stellar mass ratios than the cluster dEs, with $M_{dyn}/M_{\star} = 5.1\pm0.6$ for the group dwarfs, vs. $M_{dyn}/M_{\star} = 2.2\pm0.5$ for the cluster sample, which includes two cEs. We also search for trends in $M_{dyn}/M_{\star}$ vs. distance from M87 for the Virgo Cluster population, and find no preference for dwarfs with high values of $M_{dyn}/M_{\star}$ to reside in the cluster outskirts vs. centre.

We develop a simple analytical model that tracks galactic metallicities governed by star formation and feedback to gain insight from the observed galaxy stellar mass-metallicity relations over a large range of stellar masses and redshifts. The model reveals the following implications of star formation and feedback processes in galaxy formation. First, the observed metallicity relations provide a stringent upper limit for the averaged outflow mass-loading factors of local galaxies, which is ~20 for M_*~10^9Msun galaxies and monotonically decreases to ~1 for M_*~10^{11}Msun galaxies. Second, the inferred upper-limit for the outflow mass-loading factor sensitively depends on whether the outflow is metal-enriched with respect to the ISM metallicity. If half of the metals ejected from SNe leave the galaxy in metal-enriched winds, the outflow mass-loading factor for galaxies at any mass can barely be higher than ~10, which puts strong constraints on galaxy formation models. Third, the relatively lower stellar-phase to gas-phase metallicity ratio for lower-mass galaxies indicate that low-mass galaxies are still rapidly enriching their metallicities in recent times, while high-mass galaxies are more settled, which seems to show a downsizing effect in the metallicity evolution of galaxies. The analysis presented in the paper demonstrates the importance of accurate measurements of galaxy metallicities and the cold gas fraction of galaxies at different redshifts for constraining star formation and feedback processes, and demonstrates the power of these relations in constraining the physics of galaxy formation.

We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.

Understanding the relationship between galaxies hosting active galactic nuclei (AGN) and the dark matter halos in which they reside is key to constraining how black-hole fueling is triggered and regulated. Previous efforts have relied on simple halo mass estimates inferred from clustering, weak gravitational lensing, or halo occupation distribution modeling. In practice, these approaches remain uncertain because AGN, no matter how they are identified, potentially live a wide range of halo masses with an occupation function whose general shape and normalization are poorly known. In this work, we show that better constraints can be achieved through a rigorous comparison of the clustering, lensing, and cross-correlation signals of AGN hosts to a fiducial stellar-to-halo mass relation (SHMR) derived for all galaxies. Our technique exploits the fact that the global SHMR can be measured with much higher accuracy than any statistic derived from AGN samples alone. Using 382 moderate luminosity X-ray AGN at z<1 from the COSMOS field, we report the first measurements of weak gravitational lensing from an X-ray selected sample. Comparing this signal to predictions from the global SHMR, we find that, contrary to previous results, most X-ray AGN do not live in medium size groups ---nearly half reside in relatively low mass halos with Mh~10^12.5 Msun. The AGN occupation function is well described by the same form derived for all galaxies but with a lower normalization---the fraction of halos with AGN in our sample is a few percent. By highlighting the relatively "normal" way in which moderate luminosity X-ray AGN hosts occupy halos, our results suggest that the environmental signature of distinct fueling modes for luminous QSOs compared to moderate luminosity X-ray AGN is less obvious than previously claimed.

We model the dynamics of dwarf early-type galaxies in the Virgo cluster when subject to a variety of environmental processes. We focus on how these processes imprint trends in dynamical state (rotational vs. pressure support as measured by the $\lambda^*_{\rm Re/2}$ statistic) with projected distance from the cluster center, and compare these results to observational estimates. We find a large scatter in the gradient of $\lambda^*_{\rm Re/2}$ with projected radius. A statistical analysis shows that models with no environmental effects produce gradients as steep as those observed in none of the 100 cluster realizations we consider, while in a model incorporating tidal stirring by the cluster potential 34% of realizations produce gradients as steep as that observed. Our results suggest that tidal stirring may be the cause of the observed radial dependence of dwarf early-type dynamics in galaxy clusters.

We analyze the stellar kinematics of 39 dwarf early-type galaxies (dEs) in the Virgo cluster. Based on the specific stellar angular momentum lambda_e and the ellipticity, we find 11 slow rotators and 28 fast rotators. The fast rotators in the outer parts of the Virgo cluster rotate significantly faster than fast rotators in the inner parts of the cluster. Moreover, 10 out of the 11 slow rotators are located in the inner 3 degrees (D < 1 Mpc) of the cluster. The fast rotators contain subtle disky structures that are visible in high-pass filtered optical images, while the slow rotators do not exhibit these structures. In addition, two of the dEs have kinematically decoupled cores and four more have emission partially filling in the Balmer absorption lines. These properties suggest that Virgo cluster dEs may have originated from late-type star-forming galaxies that were transformed by the environment after their infall into the cluster. The correlation between lambda_e and the clustercentric distance can be explained by a scenario where low luminosity star-forming galaxies fall into the cluster, their gas is rapidly removed by ram pressure stripping, although some of it can be retained in their core, their star-formation is quenched but their stellar kinematics are preserved. After a long time in the cluster and several passes through its center, the galaxies are heated up and transformed into slow rotating dEs.

We describe the methodology to include nonlinear evolution, including tidal effects, in the computation of subhalo distribution properties in both cold (CDM) and warm (WDM) dark matter universes. Using semi-analytic modeling, we include effects from dynamical friction, tidal stripping, and tidal heating, allowing us to dynamically evolve the subhalo distribution. We calibrate our nonlinear evolution scheme to the CDM subhalo mass function in the Aquarius N-body simulation, producing a subhalo mass function within the range of simulations. We find tidal effects to be the dominant mechanism of nonlinear evolution in the subhalo population. Finally, we compute the subhalo mass function for $m_\chi=1.5$ keV WDM including the effects of nonlinear evolution, and compare radial number densities and mass density profiles of subhalos in CDM and WDM models. We show that all three signatures differ between the two dark matter models, suggesting that probes of substructure may be able to differentiate between them.

A robust prediction of $\Lambda{\rm CDM}$ cosmology is the halo circular velocity function (CVF), a dynamical cousin of the halo mass function. The correspondence between theoretical and observed CVFs is uncertain, however: cluster galaxies are reported to exhibit a power-law CVF consistent with $N$-body simulations, but that of the field is distinctly Schechter-like, flattened relative to $\Lambda{\rm CDM}$ expectations at circular velocities $v_c \lesssim 200\, {\rm km\, s^{-1}}$. Groups offer a powerful probe of the role environment plays in this discrepancy as they bridge the field and clusters. Here, we construct the CVF for a large, mass- and multiplicity-complete sample of group galaxies from the Sloan Digital Sky Survey. Using independent photometric $v_c$ estimators, we find no transition from a field- to $\Lambda{\rm CDM}$-shaped CVF above $v_c = 50\, {\rm km\, s^{-1}}$ as a function of group halo mass. All groups with $12.4 \lesssim \log M_{\rm halo} / M_{\odot} \lesssim 15.1$ (Local Group analogs to rich clusters) display similar Schechter-like CVFs marginally suppressed at low-$v_c$ compared to that of the field. Conversely, some agreement with $N$-body results emerges for samples saturated with late-type galaxies, with isolated late-types displaying a CVF similar in shape to $\Lambda{\rm CDM}$ predictions. We conclude that the flattening of the low-$v_c$ slope in groups is due to their depressed late-type fractions -- environment affecting the CVF only to the extent that it correlates with this quantity -- and that previous cluster analyses may suffer from interloper contamination. These results serve as useful benchmarks for cosmological simulations of galaxy formation.

We present a simulation setup for studying the dynamical and chemical evolution of the intracluster medium (ICM) and analyze a sample of 12 galaxy clusters that are diverse both kinetically (pre-merger, merging, virialized) and in total mass (M vir = 1.17 x 10^14 - 1.06 x 10^15 M). We analyzed the metal mass fraction in the ICM as a function of redshift and discuss radial trends as well as projected 2D metallicity maps. The setup combines high mass resolution N-body simulations with the semi-analytical galaxy formation model Galacticus for consistent treatment of the subgrid physics (such as galactic winds and ram-pressure stripping) in the cosmological hydrodynamical simulations. The interface between Galacticus and the hydro simulation of the ICM with FLASH is discussed with respect to observations of star formation rate histories, radial star formation trends in galaxy clusters, and the metallicity at different redshifts. As a test for the robustness of the wind model, we compare three prescriptions from different approaches. For the wind model directly taken from Galacticus, we find mean ICM metallicities between 0.2 - 0.8Z within the inner 1Mpc at z = 0. The main contribution to the metal mass fraction comes from galactic winds. The outflows are efficiently mixed in the ICM, leading to a steady homogenization of metallicities until ram-pressure stripping becomes effective at low redshifts. We find a very peculiar and yet common drop in metal mass fractions within the inner ~200kpc of the cool cores, which is due to a combination of wind suppression by outer pressure within our model and a lack of mixing after the formation of these dense regions.

As galaxy formation and evolution over long cosmic time-scales depends to a large degree on the structure of the universe, the assembly history of galaxies is potentially a powerful approach for learning about the universe itself. In this paper we examine the merger history of dark matter halos based on the Extended Press-Schechter formalism as a function of cosmological parameters, redshift and halo mass. We calculate how major halo mergers are influenced by changes in the cosmological values of $\Omega_{\rm m}$, $\Omega_{\Lambda}$, $\sigma_{8}$, the dark matter particle temperature (warm vs. cold dark matter), and the value of a constant and evolving equation of state parameter $w(z)$. We find that the merger fraction at a given halo mass varies by up to a factor of three for halos forming under the assumption of Cold Dark Matter, within different underling cosmological parameters. We find that the current measurements of the merger history, as measured through observed galaxy pairs as well as through structure, are in agreement with the concordance cosmology with the current best fit giving $1 - \Omega_{\rm m} = \Omega_{\rm \Lambda} = 0.84^{+0.16}_{-0.17}$. To obtain a more accurate constraint competitive with recently measured cosmological parameters from Planck and WMAP requires a measured merger accuracy of $\delta f_{\rm m} \sim 0.01$, implying surveys with an accurately measured merger history over 2 - 20 deg$^{2}$, which will be feasible with the next generation of imaging and spectroscopic surveys such as Euclid and LSST.

We constrain a highly simplified semi-analytic model of galaxy formation using the $z\approx 0$ stellar mass function of galaxies. Particular attention is paid to assessing the role of random and systematic errors in the determination of stellar masses, to systematic uncertainties in the model, and to correlations between bins in the measured and modeled stellar mass functions, in order to construct a realistic likelihood function. We derive constraints on model parameters and explore which aspects of the observational data constrain particular parameter combinations. We find that our model, once constrained, provides a remarkable match to the measured evolution of the stellar mass function to $z=1$, although fails dramatically to match the local galaxy HI mass function. Several "nuisance parameters" contribute significantly to uncertainties in model predictions. In particular, systematic errors in stellar mass estimate are the dominant source of uncertainty in model predictions at $z\approx 1$, with additional, non-negligble contributions arising from systematic uncertainties in halo mass functions and the residual uncertainties in cosmological parameters. Ignoring any of these sources of uncertainties could lead to viable models being erroneously ruled out. Additionally, we demonstrate that ignoring the significant covariance between bins the observed stellar mass function leads to significant biases in the constraints derived on model parameters. Careful treatment of systematic and random errors in the constraining data, and in the model being constrained, are crucial if this methodology is to be used to test hypotheses relating to the physics of galaxy formation.

We compare predictions of cooled masses and cooling rates from three stripped-down Semi-Analytic Models (SAMs) of galaxy formation with the results of N-body+SPH simulations with gas particle mass of 3.9x10^6 Msun, where radiative cooling of a gas of primordial composition is implemented. We also run a simulation where cooling is switched on at redshift ~2, in order to test cooling models in a regime in which their approximations are expected to be valid. We confirm that cooling models implemented in SAMs are able to predict the amount of cooled mass at z=0 to within ~20 per cent. However, some relevant discrepancies are found. (i) When the contribution from poorly resolved halos is subtracted out, SAMs tend to under-predict by ~30 per cent the mass that cools in the infall-dominated regime. (ii) At large halo masses SAMs tend to over-predict cooling rates, though the numerical result may be affected by the use of SPH. (iii) As found in our previous work, cooling rates are found to be significantly affected by model details: simulations disfavour models with large cores and with quenching of cooling at major mergers. (iv) When cooling is switched on at z~2, cold gas accumulates very quickly in the simulated halos. This accumulation is reproduced by SAMs with varying degrees of accuracy.

We study the oscillations and stability of self-gravitating cylindrically symmetric fluid systems and collisionless systems. This is done by studying small perturbations to the equilibrium system and finding the normal modes, using methods similar to those used in astroseismology. We find that there is a single sequence of purely radial modes that become unstable if the adiabatic exponent is less than 1. Nonradial modes can be divided into p modes, which are stable and pressure-driven, and g modes, which are are gravity driven. The g modes become unstable if the adiabatic exponent is greater than the polytrope index. These modes are analogous to the modes of a spherical star, but their behavior is somewhat different because a cylindrical geometry has less symmetry than a spherical geometry. This implies that perturbations are classified by a radial quantum number, an azimuthal quantum number, and wavelength in the z direction, which can become arbitrarily large. We find that decreasing this wavelength increases the frequency of stable modes and increases the growth rate of unstable modes. We use use variational arguments to demonstrate that filaments of collisionless matter with ergodic distribution functions are stable to purely radial perturbations, and that filaments with ergodic power-law distribution functions are stable to all perturbations.

Well-motivated particle physics theories predict the existence of particles (such as sterile neutrinos) which acquire non-negligible thermal velocities in the early universe. These particles could behave as warm dark matter (WDM) and generate a small-scale cutoff in the linear density power spectrum which scales approximately inversely with the particle mass. If this mass is of order a keV, the cutoff occurs on the scale of dwarf galaxies. Thus, in WDM models the abundance of small galaxies, such as the satellites that orbit in the halo of the Milky Way, depends on the mass of the warm particle. The abundance also scales with the mass of the host galactic halo. We use the \galform semi-analytic model of galaxy formation to calculate the properties of galaxies in universes in which the dark matter is warm. Using this method, we can compare the predicted satellite luminosity functions to the observed data for the Milky Way dwarf spheroidals, and determine a lower bound on the thermally produced WDM particle mass. This depends strongly on the value of the Milky Way halo mass and, to some extent, on the baryonic physics assumed; we examine both of these dependencies. For our fiducial model we find that for a particle mass of 3.3 keV (the 2$\sigma$ lower limit found by Viel et al. from a recent analysis of the Lyman-$\alpha$ forest) the Milky Way halo mass is required to be $> 1.4 \times 10^{12}$ \msun. For this same fiducial model, we also find that all WDM particle masses are ruled out (at 95% confidence) if the halo of the Milky Way has a mass smaller than $1.1 \times 10^{12}$ \msun, while if the mass of the Galactic halo is greater than 1.8 $\times 10^{12}$ \msun, only WDM particle masses larger than 2 keV are allowed.

A number of recent challenges to the standard Lambda-CDM paradigm relate to discrepancies that arise in comparing the abundance and kinematics of local dwarf galaxies with the predictions of numerical simulations. Such arguments rely heavily on the assumption that the local dwarf and satellite galaxies form a representative distribution in terms of their stellar-to-halo mass ratios. To address this question, we present new, deep spectroscopy using DEIMOS on Keck for 82 low mass (10^7-10^9 solar masses) star-forming galaxies at intermediate redshift (z=0.2-1). For 50 percent of these we are able to determine resolved rotation curves using nebular emission lines and thereby construct the stellar mass Tully-Fisher relation to masses as low as 10^7 solar masses. Using scaling relations determined from weak lensing data, we convert this to a stellar-to-halo mass (SHM) relation for comparison with abundance matching predictions. We find a discrepancy between the propagated predictions from simulations compared to our observations, and suggest possible reasons for this as well as future tests that will be more effective.

The escape of ionizing radiation from galaxies plays a critical role in the evolution of gas in galaxies, and the heating and ionization history of the intergalactic medium. We present semi-analytic calculations of the escape fraction of ionizing radiation for both hydrogen and helium from galaxies ranging from primordial systems to disk-type galaxies that are not heavily dust-obscured. We consider variations in the galaxy density profile, source type, location, and spectrum, and gas overdensity/distribution factors. For sufficiently hard first-light sources, the helium ionization fronts closely track or advance beyond that of hydrogen. Key new results in this work include calculations of the escape fractions for He I and He II ionizing radiation, and the impact of partial ionization from X-rays from early AGN or stellar clusters on the escape fractions from galaxy halos. When factoring in frequency-dependent effects, we find that X-rays play an important role in boosting the escape fractions for both hydrogen and helium, but especially for He II. We briefly discuss the implications of these results for recent observations of the He II reionization epoch at low redshifts, as well as the UV data and emission-line signatures from early galaxies anticipated from future satellite missions.

We present a new method to compute the first crossing distribution in excursion set theory for the case of correlated random walks. We use a combination of the path integral formalism of Maggiore & Riotto, and the integral equation solution of Zhang & Hui, and Benson et al. to find a numerically robust and convenient algorithm to derive the first crossing distribution in terms of a perturbative expansion around the limit of an uncorrelated random walk. We apply this methodology to the specific case of a Gaussian random density field filtered with a Gaussian smoothing function. By comparing our solutions to results from Monte Carlo calculations of the first crossing distribution we demonstrate that our method accurate for power spectra $P(k)\propto k^n$ for $n=1$, becoming less accurate for smaller values of $n$. It is therefore complementary to the method of Musso & Sheth, which will therefore be more useful for standard $\Lambda$CDM power spectra. Our approach is quite general, and can be adapted to other smoothing functions, and also to non-Gaussian density fields.

We present a clustering analysis of X-ray selected AGN by compiling X-ray samples from the literature and re-estimating the dark matter (DM) halo masses of AGN in a uniform manner. We find that moderate luminosity AGN (Lx(2-10 keV)=10^42-10^44 erg/sec) in the z=0-1.3 Universe are typically found in DM haloes with masses of ~10^13 Msun. We then compare our findings to the theoretical predictions of the coupled galaxy and black hole formation model GALFORM. We find good agreement when our calculation includes the hot-halo mode of accretion onto the central black hole. This type of accretion, which is additional to the common cold accretion during disk instabilities and galaxy mergers, is tightly coupled to the AGN feedback in the model. The hot-halo mode becomes prominent in DM haloes with masses greater than ~10^12.5 Msun, where AGN feedback typically operates, giving rise to a distinct class of moderate luminosity AGN that inhabit rich clusters and superclusters. Cold gas fuelling of the black hole cannot produce the observationally inferred DM halo masses of X-ray AGN. Switching off AGN feedback in the model results in a large population of luminous quasars (Lx(2-10 keV) > 10^44 erg/sec) in DM haloes with masses up to ~10^14 Msun, which is inconsistent with the observed clustering of quasars. The abundance of hot-halo AGN decreases significantly in the z~3-4 universe. At such high redshifts, the cold accretion mode is solely responsible for shaping the environment of moderate luminosity AGN. Our analysis supports two accretion modes (cold and hot) for the fuelling of supermassive black holes and strongly underlines the importance of AGN feedback in cosmological models both of galaxy formation and black hole growth.

Large-scale structures (LSS) out to z $< 3.0$ are measured in the Cosmic Evolution Survey (COSMOS) using extremely accurate photometric redshifts (photoz). The Ks-band selected sample (from Ultra-Vista) is comprised of 155,954 galaxies. Two techniques -- adaptive smoothing and Voronoi tessellation -- are used to estimate the environmental densities within 127 redshift slices. Approximately 250 statistically significant overdense structures are identified out to z $= 3.0$ with shapes varying from elongated filamentary structures to more circularly symmetric concentrations. We also compare the densities derived for COSMOS with those based on semi-analytic predictions for a $\Lambda$CDM simulation and find excellent overall agreement between the mean densities as a function of redshift and the range of densities. The galaxy properties (stellar mass, spectral energy distributions (SEDs) and star formation rates (SFRs)) are strongly correlated with environmental density and redshift, particularly at z $< 1.0 - 1.2$. Classifying the spectral type of each galaxy using the rest-frame b-i color (from the photoz SED fitting), we find a strong correlation of early type galaxies (E-Sa) with high density environments, while the degree of environmental segregation varies systematically with redshift out to z $\sim 1.3$. In the highest density regions, 80% of the galaxies are early types at z=0.2 compared to only 20% at z = 1.5. The SFRs and the star formation timescales exhibit clear environmental correlations. At z $> 0.8$, the star formation rate density (SFRD) is uniformly distributed over all environmental density percentiles, while at lower redshifts the dominant contribution is shifted to galaxies in lower density environments.

We introduce a method for constructing end-to-end mock galaxy catalogues using a semi-analytical model of galaxy formation, applied to the halo merger trees extracted from a cosmological N-body simulation. The mocks that we construct are lightcone catalogues, in which a galaxy is placed according to the epoch at which it first enters the past lightcone of the observer, and incorporate the evolution of galaxy properties with cosmic time. We determine the position between the snapshot outputs at which a galaxy enters the observer's lightcone by interpolation. As an application, we consider the effectiveness of the BzK colour selection technique, which was designed to isolate galaxies in the redshift interval 1.4<z<2.5. The mock catalogue is in reasonable agreement with the observed number counts of all BzK galaxies, as well as with the observed counts of the subsample of BzKs that are star-forming galaxies. We predict that over 75 per cent of the model galaxies with K_{AB}<=23, and 1.4<z<2.5, are selected by the BzK technique. Interloper galaxies, outside the intended redshift range, are predicted to dominate bright samples of BzK galaxies (i.e. with K_{AB}<=21). Fainter K-band cuts are necessary to reduce the predicted interloper fraction. We also show that shallow B-band photometry can lead to confusion in classifying BzK galaxies as being star-forming or passively evolving. Overall, we conclude that the BzK colour selection technique is capable of providing a sample of galaxies that is representative of the 1.4<z<2.5 galaxy population.

We study the impact of star formation and stellar feedback prescriptions on galaxy properties predicted by means of "stripped-down" versions of independently developed semi-analytic models (SAMs). These include cooling, star formation, feedback from supernovae (SNe) and simplified prescriptions for galaxy merging, but no chemical evolution, disc instabilities or AGN feedback. We run these versions on identical samples of dark matter (DM) haloes extracted from high-resolution N-body simulations in order to perform both statistical analysis and object-by-object comparisons. We compare our results with previous work based on stripped-down versions of the same SAMs including only gas cooling, and show that all feedback models provide coherent modifications in the distribution of baryons between the various gas phases. In particular, we find that the predicted hot gas fractions are considerably increased by up to a factor of three, while the corresponding cold gas fractions are correspondingly decreased, and a significant amount of mass is ejected from the DM halo. Nonetheless, we also find relevant differences in the predicted properties of model galaxies among the three SAMs: these deviations are more relevant at mass scales comparable to that of our own Galaxy, and are reduced at larger masses, confirming the varying impact of stellar feedback at different mass scales. We also check the effect of enhanced star formation events (i.e. starbursts modes), defined in connection with galaxy mergers. We find that, in general, these episodes have a limited impact in the overall star formation histories of model galaxies, even in massive DM halos where merger-driven star formation has often been considered very important.

We describe a methodology to accurately compute halo mass functions, progenitor mass functions, merger rates and merger trees in non-cold dark matter universes using a self-consistent treatment of the generalized extended Press-Schechter formalism. Our approach permits rapid exploration of the subhalo population of galactic halos in dark matter models with a variety of different particle properties or universes with rolling, truncated, or more complicated power spectra. We make detailed comparisons of analytically derived mass functions and merger histories with recent warm dark matter cosmological N-body simulations, and find excellent agreement. We show that, once the accretion of smoothly distributed matter is accounted for, coarse-grained statistics such as the mass accretion history of halos can be almost indistinguishable between cold and warm dark matter cases. However, the halo mass function and progenitor mass functions differ significantly, with the warm dark matter cases being strongly suppressed below the free-streaming scale of the dark matter. We demonstrate the importance of using the correct solution for the excursion set barrier first-crossing distribution in warm dark matter - if the solution for a flat barrier is used instead the truncation of the halo mass function is much slower, leading to an overestimate of the number of low mass halos.

We present a spectro-photometric survey of 2522 extragalactic globular clusters (GCs) around twelve early-type galaxies, nine of which have not been published previously. Combining space-based and multi-colour wide field ground-based imaging, with spectra from the Keck DEIMOS instrument, we obtain an average of 160 GC radial velocities per galaxy, with a high velocity precision of 15 km/s per GC. After studying the photometric properties of the GC systems, such as their spatial and colour distributions, we focus on the kinematics of metal-poor (blue) and metal-rich (red) GC subpopulations to an average distance of ~8 effective radii from the galaxy centre. Our results show that for some systems the bimodality in GC colour is also present in GC kinematics. The kinematics of the red GC subpopulations are strongly coupled with the host galaxy stellar kinematics. The blue GC subpopulations are more dominated by random motions, especially in the outer regions, and decoupled from the red GCs. Peculiar GC kinematic profiles are seen in some galaxies: the blue GCs in NGC 821 rotate along the galaxy minor axis, whereas the GC system of the lenticular galaxy NGC 7457 appears to be strongly rotation supported in the outer region. We supplement our galaxy sample with data from the literature and carry out a number of tests to study the kinematic differences between the two GC subpopulations. We confirm that the GC kinematics are coupled with the host galaxy properties and find that the velocity kurtosis and the slope of their velocity dispersion profiles is different between the two GC subpopulations in more massive galaxies.

We report on the discovery of strong tidal features around a dwarf spheroidal galaxy in the Hydra I galaxy cluster, indicating its ongoing tidal disruption. This very low surface brightness object, HCC-087, was originally classified as an early-type dwarf in the Hydra Cluster Catalogue (HCC), but our re-analysis of the ESO-VLT/FORS images of the HCC unearthed a clear indication of an S-shaped morphology and a large spatial extent. Its shape, luminosity (M_V=-11.6 mag), and physical size (at a half-light radius of 3.1 kpc and a full length of ~5.9 kpc) are comparable to the recently discovered NGC 4449B and the Sagittarius dwarf spheroidal, all of which are undergoing clear tidal disruption. Aided by N-body simulations we argue that HCC-087 is currently at its first apocenter, at 150 kpc, around the cluster center and that it is being tidally disrupted by the galaxy cluster's potential itself. An interaction with the near-by (50 kpc) S0 cluster galaxy HCC-005, at M* ~ 3 x 10^10 M_sun is rather unlikely, as this constellation requires a significant amount of dynamical friction and thus low relative velocities. The S-shaped morphology and large spatial extent of the satellite would, however, also appear if HCC-087 would orbit the cluster center. These features appear to be characteristic properties of satellites that are seen in the process of being tidally disrupted, independent of the environment of the destruction. An important finding of our simulations is an orientation of the tidal tails perpendicular to the orbit.

NGC 4449 is a nearby Magellanic irregular starburst galaxy with a B-band absolute magnitude of -18 and a prominent, massive, intermediate-age nucleus at a distance from Earth of 3.8 megaparsecs. It is wreathed in an extraordinary neutral hydrogen (H I) complex, which includes rings, shells and a counter-rotating core, spanning 90 kiloparsecs. NGC 4449 is relatively isolated, although an interaction with its nearest known companion-the galaxy DDO 125, some 40 kpc to the south-has been proposed as being responsible for the complexity of its HI structure. Here we report the presence of a dwarf galaxy companion to NGC 4449, namely NGC 4449B. This companion has a V-band absolute magnitude of -13.4 and a half-light radius of 2.7 kpc, with a full extent of around 8 kpc. It is in a transient stage of tidal disruption, similar to that of the Sagittarius dwarf near the Milky Way. NGC 4449B exhibits a striking S-shaped morphology that has been predicted for disrupting galaxies but has hitherto been seen only in a dissolving globular cluster. We also detect an additional arc or disk ripple embedded in a two-component stellar halo, including a component extending twice as far as previously known, to about 20 kpc from the galaxy's centre.

The observed stellar mass function (SMF) is very different to the halo mass function predicted by Lambda-CDM, and it is widely accepted that this is due to energy feedback from supernovae and black holes. However, the strength and form of this feedback is not understood. In this paper, we use the phenomenological model GALFORM to explore how galaxy formation depends on the strength and halo mass dependence of feedback. We focus on 'expulsion' models in which the wind mass loading, beta, is proportional to 1/\vdisk^n, with n=0,1,2 and contrast these models with the successful Bower et al.\ 2008 model (B8W7). A crucial development is that our code explicitly accounts for the recapture of expelled gas as the system's halo mass (and thus gravitational potential) increases. We find that a model with modest wind speed but high mass loading matches the flat portion of the SMF. When combined with AGN feedback, the model provides a good description of the observed SMF above 10^9 h^-1 Msol. However, in the expulsion models, the brightest galaxies are assembled more recently than in B8W7, and the specific star formation rates of galaxies decrease strongly with decreasing stellar mass. The expulsion models also tend to have a cosmic star formation density that is dominated by lower mass galaxies at z=1-3, and dominated high mass galaxies at low redshift. These trends are in conflict with observational data, but the comparison highlights some deficiencies of the B8W7 model also. The experiments in this paper give us important physical insight to the impact of the feedback process on the formation histories of galaxies, but the strong mass dependence of feedback adopted in B8W7 still appears to provide the most promising description of the observed universe.

We study the evolution of the cold gas content of galaxies by splitting the interstellar medium into its atomic and molecular hydrogen components, using the galaxy formation model GALFORM in the LCDM framework. We calculate the molecular-to-atomic hydrogen mass ratio, H2/HI, in each galaxy using two different approaches; the pressure-based empirical relation of Blitz & Rosolowsky and the theoretical model of Krumholz, McKeee & Tumlinson, and apply them to consistently calculate the star formation rates of galaxies. We find that the model based on the Blitz & Rosolowsky law predicts an HI mass function, CO(1-0) luminosity function, correlations between the H2/HI ratio and stellar and cold gas mass, and infrared-CO luminosity relation in good agreement with local and high redshift observations. The HI mass function evolves weakly with redshift, with the number density of high mass galaxies decreasing with increasing redshift. In the case of the H2 mass function, the number density of massive galaxies increases strongly from z=0 to z=2, followed by weak evolution up to z=4. We also find that the H2/HI ratio of galaxies is strongly dependent on stellar and cold gas mass, and also on redshift. The slopes of the correlations between H2/HI and stellar and cold gas mass hardly evolve, but the normalisation increases by up to two orders of magnitude from z=0-8. The strong evolution in the H2 mass function and the H2/HI ratio is primarily due to the evolution in the sizes of galaxies and secondarily, in the gas fractions. The predicted cosmic density evolution of HI agrees with the observed evolution inferred from DLAs, and is dominated by low/intermediate mass halos. We find that previous theoretical studies have largely overestimated the redshift evolution of the global H2/HI ratio due to limited resolution. We predict a maximum of rho_H2/rho_HI~1.2 at z~3.5.

Dark matter halo merger trees are now routinely extracted from cosmological simulations of structure formation. These trees are frequently used as inputs to semi-analytic models of galaxy formation to provide the backbone within which galaxy formation takes place. By necessity, these merger trees are constructed from a finite set of discrete "snapshots" of the N-body simulation and so have a limited temporal resolution. To date, there has been little consideration of how this temporal resolution affects the properties of galaxies formed within these trees. In particular, the question of how many snapshots are needed to achieve convergence in galaxy properties has not be answered. Therefore, we study the convergence in the stellar and total baryonic masses of galaxies, distribution of merger times, stellar mass functions and star formation rates in the Galacticus model of galaxy formation as a function of the number of "snapshot" times used to represent dark matter halo merger trees. When utilizing snapshots between z=20 and z=0, we find that at least 128 snapshots are required to achieve convergence to within 5% for galaxy masses. This convergence is obtained for mean quantities averaged over large samples of galaxies - significant variance for individual galaxies remains even when using very large numbers of snapshots. We find only weak dependence of the rate of convergence on the distribution of snapshots in time - snapshots spaced uniformly in the expansion factor, uniformly in the logarithm of expansion factor or uniformly in the logarithm of critical overdensity for collapse work equally well in almost all cases. We provide input parameters to Galacticus which allow this type of convergence study to be tuned to other simulations and to be carried out for other galaxy properties.

We present a model for the satellites of the Milky Way in which galaxy formation is followed using semi-analytic techniques applied to the six high-resolution N-body simulations of galactic halos of the Aquarius project. The model, calculated using the Galform code, incorporates improved treatments of the relevant physics in the LambdaCDM cosmogony, particularly a self-consistent calculation of reionization by UV photons emitted by the forming galaxy population, including the progenitors of the central galaxy. Along the merger tree of each halo, the model calculates gas cooling (by Compton scattering off cosmic microwave background photons, molecular hydrogen and atomic processes), gas heating (from hydrogen photoionization and supernova energy), star formation and evolution. The evolution of the intergalactic medium is followed simultaneously with that of the galaxies. Star formation in the more massive progenitor subhalos is suppressed primarily by supernova feedback, while for smaller subhalos it is suppressed primarily by photoionization due to external and internal sources. The model is constrained to match a wide range of properties of the present day galaxy population as a whole, but at high redshift it requires an escape fraction of UV photons near unity in order completely to reionize the universe by redshift z ~ 8. In the most successful model the local sources photoionize the pre-galactic region completely by z ~ 10. In addition to the luminosity function of Milky Way satellites, the model matches their observed luminosity-metallicity relation, their radial distribution and the inferred values of the mass within 300 pc, which in the models increase slowly but significantly with luminosity. There is a large variation in satellite properties from halo to halo, with the luminosity function, for example, varying by a factor of ~ 2 among the six simulations.

The first stars and quasars are known sources of hard ionizing radiation in the first billion years of the Universe. We examine the joint effects of X-rays and hard UV radiation from such first-light sources on the hydrogen and helium reionization of the intergalactic medium (IGM) at early times, and the associated heating. We study the growth and evolution of individual HII, HeII and HeIII regions around early galaxies with first stars and/or QSO populations. We find that in the presence of helium-ionizing radiation, X-rays may not dominate the ionization and thermal history of the IGM at redshifts, z, of 10-20, contributing relatively modest increases to IGM ionization, and heating up to about 10^3--10^5 K in IGM temperatures. We also calculate the 21 cm signal expected from a number of scenarios with metal-free starbursts and quasars at these redshifts. The peak values for the spin temperature reach about 10^4 to 10^5 K in such cases. The maximum values for the 21 cm brightness temperature are around 30-40 mK in emission, while the net values of the 21 cm absorption signal range from about a few to 60 mK on scales of 0.01-1 Mpc. We find that the 21 cm signature of X-ray versus UV ionization could be distinct, with the emission signal expected from X-rays alone occurring at smaller scales than that from UV radiation, resulting from the inherently different spatial scales at which X-ray and UV ionization/heating manifest. This difference is time-dependent, and becomes harder to distinguish with an increasing X-ray contribution to the total ionizing photon production. Such differing scale-dependent contributions from X-ray and UV photons may therefore "blur" the 21 cm signature of the percolation of ionized bubbles around early halos (depending on whether a cosmic X-ray or UV background built up first), and affect the interpretation of 21 cm data constraints on reionization.

We present a global analysis of kinematics and metallicity in the nearest S0 galaxy, NGC 3115, along with implications for its assembly history. The data include high-quality wide-field imaging from Suprime-Cam on the Subaru telescope, and multi-slit spectra of the field stars and globular clusters (GCs) obtained using Keck-DEIMOS/LRIS and Magellan-IMACS. Within two effective radii, the bulge (as traced by the stars and metal-rich GCs) is flattened and rotates rapidly (v/sigma > 1.5). At larger radii, the rotation declines dramatically to v/sigma ~ 0.7, but remains well-aligned with the inner regions. The radial decrease in characteristic metallicity of both the metal-rich and metal-poor GC subpopulations produces strong gradients with power law slopes of -0.17 +/- 0.04 and -0.38 +/- 0.06 dex per dex, respectively. We argue that this pattern is not naturally explained by a binary major merger, but instead by a two-phase assembly process where the inner regions have formed in an early violent, dissipative phase, followed by the protracted growth of the outer parts via minor mergers with typical mass ratios of ~ 15-20:1.

We investigate the consequences of applying different star formation laws in the galaxy formation model GALFORM. Three broad star formation laws are implemented: the empirical relations of Kennicutt and Schmidt and Blitz & Rosolowsky and the theoretical model of Krumholz, McKee & Tumlinson. These laws have no free parameters once calibrated against observations of the star formation rate (SFR) and gas surface density in nearby galaxies. We start from published models, and investigate which observables are sensitive to a change in the star formation law, without altering any other model parameters. We show that changing the star formation law (i) does not significantly affect either the star formation history of the universe or the galaxy luminosity functions in the optical and near-IR, due to an effective balance between the quiescent and burst star formation modes; (ii) greatly affects the cold gas contents of galaxies; (iii) changes the location of galaxies in the SFR versus stellar mass plane, so that a second sequence of "passive" galaxies arises, in addition to the known "active" sequence. We show that this plane can be used to discriminate between the star formation laws.

We predict the formation histories, properties and descendants of Lyman-break galaxies (LBGs) in the Lambda-CDM cosmology. In our model, which incorporates a top-heavy IMF in starbursts, we find that most LBGs are starbursts triggered by minor mergers of galaxies. The duration of the LBG phase is predicted to be quite short, ~20-60 Myr. We investigate the distributions of stellar and halo masses and morphologies for bright (L_UV > L*_UV) and faint (L_UV > 0.1 L*_UV) LBGs at z=3, z=6 and z=10 (where we classify LBGs according to their rest-frame UV luminosities relative the observed characteristic luminosity L*_UV at z \approx 3). Bright LBGs at z=3 are predicted to have median stellar masses ~ 1x10^9 Msun/h and host halo masses ~ 3x10^{11} Msun/h, and to be typically mildly disk-dominated in stellar mass. On the other hand, faint LBGs at z=10 are predicted to have median stellar masses of only ~ 1x10^7 Msun/h and host halo masses 2x10^{10} Msun/h, and to be generally bulge-dominated. Bright LBGs at z=3 evolve to present-day galaxies with median stellar mass ~ 5x^{10} Msun/h (comparable to the Milky Way), consisting of roughly equal numbers of disk- and bulge-dominated systems, and hosted by halos with median mass ~2x10^{13} Msun/h (corresponding to medium-size galaxy groups). The model predicts that 40% of Milky Way mass galaxies at the present-day have a bright LBG progenitor in the redshift range 3<z<4, while 95% have a faint LBG progenitor in the same redshift range, and 7% have a faint LBG progenitor at 10<z<11. With our multiwavelength model, we also investigate the overlap between the LBG population and that of submillimetre selected galaxies (SMGs); at z=3, only ~1% of bright LBGs are also predicted to also be bright SMGs (with an 850 mum flux in excess of 5 mJy).

Using data from the COSMOS survey, we perform the first joint analysis of galaxy-galaxy weak lensing, galaxy spatial clustering, and galaxy number densities. Carefully accounting for sample variance and for scatter between stellar and halo mass, we model all three observables simultaneously using a novel and self-consistent theoretical framework. Our results provide strong constraints on the shape and redshift evolution of the stellar-to-halo mass relation (SHMR) from z=0.2 to z=1. At low stellar mass, we find that halo mass scales as Mh M*^0.46 and that this scaling does not evolve significantly with redshift to z=1. We show that the dark-to-stellar ratio, Mh/M*, varies from low to high masses, reaching a minimum of Mh/M*~27 at M*=4.5x10^10 Msun and Mh=1.2x10^12 Msun. This minimum is important for models of galaxy formation because it marks the mass at which the accumulated stellar growth of the central galaxy has been the most efficient. We describe the SHMR at this minimum in terms of the "pivot stellar mass", M*piv, the "pivot halo mass", Mhpiv, and the "pivot ratio", (Mh/M*)piv. Thanks to a homogeneous analysis of a single data set, we report the first detection of mass downsizing trends for both Mhpiv and M*piv. The pivot stellar mass decreases from M*piv=5.75+-0.13x10^10 Msun at z=0.88 to M*piv=3.55+-0.17x10^10 Msun at z=0.37. Intriguingly, however, the corresponding evolution of Mhpiv leaves the pivot ratio constant with redshift at (Mh/M*)piv~27. We use simple arguments to show how this result raises the possibility that star formation quenching may ultimately depend on Mh/M* and not simply Mh, as is commonly assumed. We show that simple models with such a dependence naturally lead to downsizing in the sites of star formation. Finally, we discuss the implications of our results in the context of popular quenching models, including disk instabilities and AGN feedback.

While the dark matter content within the most massive giant and smallest dwarf galaxies has been probed -- spanning a range of over one million in mass -- an important observational gap remains for galaxies of intermediate mass. This gap covers K band magnitudes of approximately -16 > M_K > -18 (for which dwarf galaxies have B--K ~ 2). On the high mass side of the gap are dwarf elliptical (dE) galaxies, that are dominated by stars in their inner regions. While the low mass side includes dwarf spheroidal (dSph) galaxies that are dark matter-dominated and ultra compact dwarf (UCD) objects that are star-dominated. Evolutionary pathways across the gap have been suggested but remain largely untested because the `gap' galaxies are faint, making dynamical measurements very challenging. With long exposures on the Keck telescope using the ESI instrument we have succeeded in bridging this gap by measuring the dynamical mass for five dwarf galaxies with M_K ~ -17.5 (M_B ~ --15.5). With the exception of our brightest dwarf galaxy, they possess relatively flat velocity dispersion profiles of around 20 km/s. By examining their 2D scaling relations and 3D fundamental manifold, we found that the sizes and velocity dispersions of these gap galaxies reveal continuous trends from dE to dSph galaxies. We conclude that low-luminosity dwarf elliptical galaxies are dominated by stars, not by dark matter, within their half light radii. This finding can be understood if internal feedback processes are operating most efficiently in gap galaxies, gravitationally heating the centrally-located dark matter to larger radii. Whereas external environmental processes, which can strip away stars, have a greater influence on dSph galaxies resulting in their higher dark matter fractions. Abridged.

Recent observations of galaxy clusters have shown that environmental effects apparently associated with the cluster begin to lower the star formation rates of galaxies at distances as great as three times the cluster virial radius. These observations may indicate preprocessing of cluster galaxies in groups or in the cluster core for galaxies on highly elliptical orbits, but may also imply that the environmental effects due to the cluster are directly affecting galaxies on their first infall. To explore these issues, we investigate different models of ram pressure stripping as it acts on satellite galaxies in clusters, and compare to observations of the radial star formation gradient in clusters. We calculate the location of the accretion shock around model clusters, and use this as the radius of onset of ram pressure stripping in the GALFORM semi-analytic model of galaxy formation. Comparison of the results of our model, and previously considered, simpler ram pressure models, with recent observations indicates that current data is unable to strongly discriminate between models of ram pressure stripping due to the complex interplay of preprocessing effects at work. However, future observations of a larger sample of clusters will likely be able to place stronger constraints on the process of ram pressure stripping and its role in shaping radial trends in and around clusters.

We have analyzed high resolution N-body simulations of dark matter halos, focusing specifically on the evolution of angular momentum. We find that not only is individual particle angular momentum not conserved, but the angular momentum of radial shells also varies over the age of the Universe by up to factors of a few. We find that torques from external structure are the most likely cause for this distribution shift. Since the model of adiabatic contraction that is often applied to model the effects of galaxy evolution on the dark-matter density profile in a halo assumes angular momentum conservation, this variation implies that there is a fundamental limit on the possible accuracy of the adiabatic contraction model in modeling the response of DM halos to the growth of galaxies.

We use a coupled model of the formation and evolution of galaxies and black holes (BH) to study the evolution of active galactic nuclei (AGN) in a cold dark matter universe. The model predicts the BH mass, spin and mass accretion history. BH mass grows via accretion triggered by discs becoming dynamically unstable or galaxy mergers (called the starburst mode) and accretion from quasi-hydrostatic hot gas haloes (called the hot-halo mode). By taking into account AGN obscuration, we obtain a very good fit to the observed luminosity functions (LF) of AGN (optical, soft and hard X-ray, and bolometric) for a wide range of redshifts (0<z<6). The model predicts a hierarchical build up of BH mass, with the typical mass of actively growing BHs increasing with decreasing redshift. Remarkably, despite this, we find downsizing in the AGN population, in terms of the differential growth with redshift of the space density of faint and bright AGN. This arises naturally from the interplay between the starburst and hot-halo accretion modes. The faint end of the LF is dominated by massive BHs experiencing quiescent accretion via a thick disc, primarily during the hot-halo mode. The bright end of the LF, on the other hand, is dominated by AGN which host BHs accreting close to or in excess of the Eddington limit during the starburst mode. The model predicts that the comoving space density of AGN peaks at z~3, similar to the star formation history. However, when taking into account obscuration, the space density of faint AGN peaks at lower redshift (z<2) than that of bright AGN (z~2-3). This implies that the cosmic evolution of AGN is shaped in part by obscuration.

We explore a dark-matter model in which there are two dark-matter species nearly degenerate in mass, with epsilon = Delta M/M << 1. The heavier particle undergoes two-body decay with a half-life tau, to the lighter dark-matter particle and a noninteracting massless particle. Unlike previous work on decaying dark matter, we explore the regime tau > 100 Myr and non-relativistic kick speeds vk / c = epsilon. Using a set of N-body simulations of isolated dark-matter halos, we show how halos change as a function of tau and vk. We find that tau < 40 Gyr is ruled out for vk > 20 km s^{-1} (epsilon > 10^{-4}) when we compare the simulations to observations of dwarf-galaxy- to cluster-mass dark matter halos. We highlight which set of observations should provide better future constraints for decays and other types of dark-matter physics.

We make a detailed investigation of the properties of Lyman-break galaxies (LBGs) in the LambdaCDM model. We present predictions for two published variants of the GALFORM semi-analytical model: the Baugh et al. (2005) model, which has star formation at high redshifts dominated by merger-driven starbursts with a top-heavy IMF, and the Bower et al. (2006) model, which has AGN feedback and a standard Solar neighbourhood IMF throughout. We show predictions for the evolution of the rest-frame far-UV luminosity function in the redshift range z=3-20, and compare with the observed luminosity functions of LBGs at z=3-10. We find that the Baugh et al. model is in excellent agreement with these observations, while the Bower et al. model predicts too many high-luminosity LBGs. Dust extinction, which is predicted self-consistently based on galaxy gas contents, metallicities and sizes, is found to have a large effect on LBG luminosities. We compare predictions for the size evolution of LBGs at different luminosities with observational data for 2<z<7, and find the Baugh et al. model to be in good agreement. We present predictions for stellar, halo and gas masses, star formation rates, circular velocities, bulge-to-disk ratios, gas and stellar metallicities and clustering bias, as functions of far-UV luminosity and redshift. We find broad consistency with current observational constraints. We then present predictions for the abundance and angular sizes of LBGs out to very high redshift (z<20), finding that planned deep surveys with JWST should detect objects out to z<15. The typical UV luminosities of galaxies are predicted to be very low at high redshifts, which has implications for detecting the galaxies responsible for reionizing the IGM; for example, at z=10, 50% of the ionizing photons are expected to be produced by galaxies fainter than M_AB(1500A)-5logh ~ -15.

We consider constraints on a phenomenological dark-matter model consisting of two nearly degenerate particle species using observed properties of the Milky Way satellite galaxy population. The two parameters of this model, assuming the particle masses are >~ GeV, are v_k, the recoil speed of the daughter particle, and tau, the lifetime of the parent particle. The satellite constraint that spans the widest range of v_k is the number of satellites that have a mass within 300 pc M300 > 5 x 10^6 solar masses, although constraints based on M300 in the classical dwarfs and the overall velocity function are competitive for v_k >~ 50 km/s. In general, we find that tau <~ 30 Gyr is ruled out for 20 km/s <~ v_k <~ 200 km/s, although we find that the limits on tau for fixed v_k can change constraints by a factor of ~3 depending on the star-formation histories of the satellites. We advocate using the distribution of M300 in Milky Way satellites determined by next-generation all-sky surveys and follow-up spectroscopy as a probe of dark-matter properties.

We track the co-evolution of supermassive black holes (SMBHs) and their host galaxies. The calculation is embedded in the GALFORM semi-analytical model which simulates the formation and evolution of galaxies in a cold dark matter (CDM) universe. During the evolution of the host galaxy, hot and cold gas are added to the SMBH by flows triggered by halo gas cooling, disc instabilities and galaxy mergers. This builds up the mass and spin of the BH, and the resulting accretion power regulates the gas cooling and subsequent star formation. The accretion flow is assumed to form a geometrically thin cool disc when the accretion rate exceeds 0.01\dot{M}_Edd, and a geometrically thick, radiatively inefficient hot flow when the accretion rate falls below this value. The resulting quasar optical luminosity function matches observations very well, and the mass of the SMBH correlates with the mass of the galaxy bulge as observed. The BH spin distribution depends strongly on whether the gas in any given accretion episode remains in the same plane (prolonged accretion) or whether, due to self-gravity, it fragments into multiple, randomly aligned accretion episodes (chaotic accretion). In the chaotic accretion model there is a clear correlation of spin with SMBH mass. Massive BHs (M>5\times10^8\Msun) are hosted by giant elliptical galaxies and are rapidly spinning, while lower mass BHs are hosted in spiral galaxies and have much lower spin. Using the Blandford-Znajek mechanism for jet production to calculate the jet power, our model is able to reproduce the radio loudness of radio galaxies, LINERS and Seyferts. This is the first confirmation that a CDM galaxy formation model can reproduce the observed phenomenology of AGN.

We describe a new, free and open source semi-analytic model of galaxy formation, Galacticus. The Galacticus model was designed to be highly modular to facilitate expansion and the exploration of alternative descriptions of key physical ingredients. We detail the Galacticus engine for evolving galaxies through a merging hierarchy of dark matter halos and give details of the specific implementations of physics currently available in Galacticus. Finally, we show results from an example model that is in reasonably good agreement with several observational datasets. We use this model to explore numerical convergence and to demonstrate the types of information which can be extracted from Galacticus.

We review the current theory of how galaxies form within the cosmological framework provided by the cold dark matter paradigm for structure formation. Beginning with the pre-galactic evolution of baryonic material we describe the analytical and numerical understanding of how baryons condense into galaxies, what determines the structure of those galaxies and how internal and external processes (including star formation, merging, active galactic nuclei etc.) determine their gross properties and evolution. Throughout, we highlight successes and failings of current galaxy formation theory. We include a review of computational implementations of galaxy formation theory and assess their ability to provide reliable modelling of this complex phenomenon. We finish with a discussion of several "hot topics" in contemporary galaxy formation theory and assess future directions for this field.

Semi-analytic models are a powerful tool for studying the formation of galaxies. However, these models inevitably involve a significant number of poorly constrained parameters that must be adjusted to provide an acceptable match to the observed universe. In this paper, we set out to quantify the degree to which observational data-sets can constrain the model parameters. By revealing degeneracies in the parameter space we can hope to better understand the key physical processes probed by the data. We use novel mathematical techniques to explore the parameter space of the GALFORM semi-analytic model. We base our investigation on the Bower et al. 2006 version of GALFORM, adopting the same methodology of selecting model parameters based on an acceptable match to the local bJ and K luminosity functions. The model contains 16 parameters that are poorly constrained, and we investigate this parameter space using the Model Emulator technique, constructing a Bayesian approximation to the GALFORM model that can be rapidly evaluated at any point in parameter space. By combining successive waves of emulation, we show that only 0.26% of the initial volume is of interest for further exploration. However, within this region we show that the Bower et al. 2006 model is only one choice from an extended sub-space of model parameters that can provide equally acceptable fits. We explore the geometry of this region and begin to explore the physical connections between parameters that are exposed by this analysis. We also consider the impact of adding additional observational data to further constrain the parameter space.

(Abridged) We use N-body simulations to study the effects that a divergent (i.e. "cuspy") dark matter (DM) profile introduces on the tidal evolution of dwarf spheroidal galaxies (dSphs). Our models assume cosmologically-motivated initial conditions where dSphs are DM-dominated systems on eccentric orbits about a host galaxy composed of a dark halo and a baryonic disc. We find that the resilience of dSphs to tidal stripping is extremely sensitive to the halo cuspiness; whereas dwarfs with a cored profile can be easily destroyed by the host disc, those with cusps always retain a bound remnant. For a given halo profile the evolution of the structural parameters as driven by tides is controlled solely by the total amount of mass lost. This information is used to construct a semi-analytic code that simulates the hierarchical build-up of spiral galaxies assuming different halo profiles and disc masses. We find that tidal encounters with discs tend to decrease the average mass of satellites at all galactocentric radii. Interestingly, satellites accreted before re-ionization (z>6), which may be singled out by anomalous metallicity patterns, survive only if haloes are cuspy. We show that the size-mass relation established from Milky Way (MW) dwarfs strongly supports the presence of cusps in the majority of these systems, as cored models systematically underestimate the masses of the known Ultra-Faint dSphs. Our models also indicate that a massive M31 disc may explain why many of its dSphs fall below the size-mass relationship derived from MW dSphs. We use our models to constrain the mass threshold below which star formation is suppressed in DM haloes, finding that luminous satellites must be accreted with masses above 10^8--10^9 M_sol in order to explain the size-mass relation observed in MW dwarfs.

We present six simulations of galactic stellar haloes formed by the tidal disruption of accreted dwarf galaxies in a fully cosmological setting. Our model is based on the Aquarius project, a suite of high resolution N-body simulations of individual dark matter haloes. We tag subsets of particles in these simulations with stellar populations predicted by the Galform semi-analytic model. Our method self-consistently tracks the dynamical evolution and disruption of satellites from high redshift. The luminosity function and structural properties of surviving satellites, which agree well with observations, suggest that this technique is appropriate. We find that accreted stellar haloes are assembled between 1<z<7 from less than 5 significant progenitors. These progenitors are old, metal-rich satellites with stellar masses similar to the brightest Milky Way dwarf spheroidals. In contrast to previous stellar halo simulations, we find that several of these major contributors survive as self-bound systems to the present day. Both the number of these significant progenitors and their infall times are inherently stochastic. This results in great diversity among our stellar haloes, which amplifies small differences between the formation histories of their dark halo hosts. The masses and density/surface-brightness profiles of the stellar haloes are consistent with expectations from the Milky Way and M31. Each halo has a complex structure, consisting of well-mixed components, tidal streams, shells and other subcomponents. This structure is not adequately described by smooth models. We find one example of an accreted thick disk. Contrasts in age and metallicity between halo stars and those in surviving satellites are in broad agreement with recent observations. [Abridged]

We use stripped-down versions of three semi-analytic galaxy formation models to study the influence of different assumptions about gas cooling and galaxy mergers. By running the three models on identical sets of merger trees extracted from high-resolution cosmological N-body simulations, we are able to perform both statistical analyses and halo-by-halo comparisons. Our study demonstrates that there is a good statistical agreement between the three models used here, when operating on the same merger trees, reflecting a general agreement in the underlying framework for semi-analytic models. We also show, however, that various assumptions that are commonly adopted to treat gas cooling and galaxy mergers can lead to significantly different results, at least in some regimes. In particular, we find that the different models adopted for gas cooling lead to similar results for mass scales comparable to that of our own Galaxy. Significant differences, however, arise at larger mass scales. These are largely (but not entirely) due to different treatments of the `rapid cooling' regime, and different assumptions about the hot gas distribution. At this mass regime, the predicted cooling rates can differ up to about one order of magnitude, with important implications on the relative weight that these models give to AGN feedback in order to counter-act excessive gas condensation in relatively massive haloes at low redshift. Different assumptions in the modelling of galaxy mergers can also result in significant differences in the timings of mergers, with important consequences for the formation and evolution of massive galaxies.

A generic expectation for gas accreted by high mass haloes is that it is shock heated to the virial temperature of the halo. In low mass haloes, or at high redshift, however, the gas cooling rate is sufficiently rapid that an accretion shock is unlikely to form. Instead, gas can accrete directly into the centre of the halo in a `cold mode' of accretion. Although semi-analytic models have always made a clear distinction between hydrostatic and rapid cooling they have not made a distinction between whether or not an accretion shock forms. Starting from the well-established Galform code, we investigate the effect of explicitly accounting for cold mode accretion using the shock stability model of Birnboim & Dekel. When we modify the code so that there is no effective feedback from galaxy formation, we find that cold mode accretion is the dominant channel for feeding gas into the galaxies at high redshifts. However, this does not translate into a significant difference in the star formation history of the universe compared to the previous code. When effective feedback is included in the model, we find that the the cold mode is much less apparent because of the presence of gas ejected from the galaxy. Thus the inclusion of the additional cold mode physics makes little difference to basic results from earlier semi-analytic models which used a simpler treatment of gas accretion. For more sophisticated predictions of its consequences, we require a better understanding of how the cold mode delivers angular momentum to galaxies and how it interacts with outflows.

The distribution of cold gas in dark matter haloes is driven by key processes in galaxy formation: gas cooling, galaxy mergers, star formation and reheating of gas by supernovae. We compare the predictions of four different galaxy formation models for the spatial distribution of cold gas. We find that satellite galaxies make little contribution to the abundance or clustering strength of cold gas selected samples, and are far less important than they are in optically selected samples. The halo occupation distribution function of present-day central galaxies with cold gas mass > 10^9 h^-1 Msun is peaked around a halo mass of ~ 10^11 h^-1 Msun, a scale that is set by the AGN suppression of gas cooling. The model predictions for the projected correlation function are in good agreement with measurements from the HI Parkes All-Sky Survey. We compare the effective volume of possible surveys with the Square Kilometre Array with those expected for a redshift survey in the near-infrared. Future redshift surveys using neutral hydrogen emission will be competitive with the most ambitious spectroscopic surveys planned in the near-infrared.

Over the past several decades, galaxy formation theory has met with significant successes. In order to test current theories thoroughly we require predictions for as yet unprobed regimes. To this end, we describe a new implementation of the Galform semi-analytic model of galaxy formation. Our motivation is the success of the model described by Bower et al. in explaining many aspects of galaxy formation. Despite this success, the Bower et al. model fails to match some observational constraints and certain aspects of its physical implementation are not as realistic as we would like. The model described in this work includes substantially updated physics, taking into account developments in our understanding over the past decade, and removes certain limiting assumptions made by this (and most other) semi-analytic models. This allows it to be exploited reliably in high-redshift and low mass regimes. Furthermore, we have performed an exhaustive search of model parameter space to find a particular set of model parameters which produce results in good agreement with a wide range of observational data (luminosity functions, galaxy sizes and dynamics, clustering, colours, metal content) over a wide range of redshifts. This model represents a solid basis on which to perform calculations of galaxy formation in as yet unprobed regimes.

We investigate the assembly of groups and clusters of galaxies using the Millennium dark matter simulation and the associated gas simulations and semi-analytic catalogues of galaxies. In particular, in order to find an observable quantity that could be used to identify early-formed groups, we study the development of the difference in magnitude between their brightest galaxies to assess the use of magnitude gaps as possible indicators. We select galaxy groups and clusters at redshift z=1 with dark matter halo mass M(R200) > 1E13/h Msun, and trace their properties until the present time (z=0). We consider only the systems with X-ray luminosity L_X> 0.25E42/h^2 erg/s at z=0. While it is true that a large magnitude gap between the two brightest galaxies of a particular group often indicates that a large fraction of its mass was assembled at an early epoch, it is not a necessary condition. More than 90% of fossil groups defined on the basis of their magnitude gaps (at any epoch between 0<z<1) cease to be fossils within 4 Gyr, mostly because other massive galaxies are assembled within their cores, even though most of the mass in their haloes might have been assembled at early times. We show that, compared to the conventional definition of fossil galaxy groups based on the magnitude gap Delta m(12)> 2 (in the R-band, within 0.5R200 of the centre of the group), an alternative criterion Delta m(14)>2.5 (within the same radius) finds 50% more early-formed systems, and those that on average retain their fossil phase longer. However, the conventional criterion performs marginally better at finding early-formed groups at the high-mass end of groups. Nevertheless, both criteria fail to identify a majority of the early-formed systems.

We use a model for the evolution of galaxies in the far-IR based on the LambdaCDM cosmology to make detailed predictions for upcoming cosmological surveys with the Herschel Space Observatory. We use the combined GALFORM semi-analytical galaxy formation model and GRASIL spectrophotometric code to compute galaxy SEDs including the reprocessing of radiation by dust. The model, which is the same as that in Baugh et al. (2005), assumes two different IMFs: a normal solar neighbourhood IMF for quiescent star formation in disks, and a very top-heavy IMF in starbursts triggered by galaxy mergers. We have shown previously that the top-heavy IMF appears necessary to explain the number counts and redshifts of faint sub-mm galaxies. In this paper, we present predictions for galaxy luminosity functions, number counts and redshift distributions in the Herschel imaging bands. We find that source confusion will be a serious problem in the deepest planned surveys. We also show predictions for physical properties such as star formation rates and stellar, gas and halo masses, together with fluxes at other wavelengths (from the far-UV to the radio) relevant for multi-wavelength follow-up observations. We investigate what fraction of the total IR emission from dust and of the high-mass star formation over the history of the Universe should be resolved by planned surveys with Herschel, and find a fraction ~30-50%, depending on confusion. Finally, we show that galaxies in Herschel surveys should be significantly clustered.

Recent focus on the importance of cold, unshocked gas accretion in galaxy formation -- not explicitly included in semi-analytic studies -- motivates the following detailed comparison between two inherently different modelling techniques: direct hydrodynamical simulation and semi-analytic modelling. By analysing the physical assumptions built into the Gasoline simulation, formulae for the emergent behaviour are derived which allow immediate and accurate translation of these assumptions to the Galform semi-analytic model. The simulated halo merger history is then extracted and evolved using these equivalent equations, predicting a strikingly similar galactic system. This exercise demonstrates that it is the initial conditions and physical assumptions which are responsible for the predicted evolution, not the choice of modelling technique. On this level playing field, a previously published Galform model is applied (including additional physics such as chemical enrichment and feedback from active galactic nuclei) which leads to starkly different predictions.

The Galform semi-analytic model of galaxy formation is used to explore the mechanisms primarily responsible for the three types of galaxies seen in the local universe: bulge, bulge+disk and disk, identified with the visual morphological types E, S0/a-Sbc, and Sc-Scd, respectively. With a suitable choice of parameters the Galform model can accurately reproduce the observed local K_s-band luminosity function (LF) for galaxies split by visual morphological type. The successful set of model parameters is used to populate the Millennium Simulation with 9.4 million galaxies and their dark matter halos. The resulting catalogue is then used to explore the evolution of galaxies through cosmic history. The model predictions concur with recent observational results including the galaxy merger rate, the star formation rate and the seemingly anti-hierarchical evolution of ellipticals. However, the model also predicts significant evolution of the elliptical galaxy LF that is not observed. The discrepancy raises the possibility that samples of z~1 galaxies which have been selected using colour and morphological criteria may be contaminated with galaxies that are not actually ellipticals.

Unbounded outflows in the form of highly collimated jets and broad winds appear to be a ubiquitous feature of accreting black hole systems. The most powerful jets are thought to derive a significant fraction, if not the majority, of their power from the rotational energy of the black hole. Whatever the precise mechanism that causes them, these jets must therefore exert a braking torque on the black hole. We calculate the spin-up function for an accreting black hole, accounting for this braking torque. We find that the predicted black hole spin-up function depends only on the black hole spin and dimensionless parameters describing the accretion flow. Using recent relativistic magnetohydrodynamical numerical simulation results to calibrate the efficiency of angular momentum transfer in the flow, we find that an ADAF flow will spin a black hole up (or down) to an equilibrium value of about 96% of the maximal spin value in the absence of jets. Combining our ADAF system with a simple model for jet power, we demonstrate that an equilibrium is reached at approximately 93% of the maximal spin value, as found in the numerical simulation studies of the spin-up of accreting black holes, at which point the spin-up of the hole by accreted material is balanced by the braking torque arising from jet production. Our model also yields a relationship between jet efficiency and black hole spin that is in surprisingly good agreement with that seen in the simulation studies, indicating that our simple model is a useful and convenient description of ADAF inflow - jet outflow about a spinning black hole for incorporation in models of the formation and evolution of galaxies, groups and clusters of galaxies.

The huge size and uniformity of the Sloan Digital Sky Survey makes 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 scale-lengths of local galaxies. GALFORM models the luminosity and size of the disk 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 Sersic index. We test the predictions of two published models set in the cold dark matter cosmology: the Baugh et al. (2005) model, which assumes a top-heavy initial mass function (IMF) in starbursts and superwind feedback, and the Bower et al. (2006) model, which uses AGN 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 disk instabilities.

The clustering amplitude of galaxies depends on their intrinsic luminosity. We compare the properties of publicly available galaxy formation models with clustering measurements from the two-degree field galaxy redshift survey. The model predictions show the same qualitative behaviour as the data but fail to match the observations at the level of accuracy at which current measurements can be made. We demonstrate that this is due to the model producing too many satellite galaxies in massive haloes. We implement simple models to describe two new processes, satellite-satellite mergers and the tidal dissolution of satellites to investigate their impact on the predicted clustering. We find that both processes need to be included in order to produce a model which matches the observations.

Models of disk galaxy formation commonly predict the existence of an extended reservoir of accreted hot gas surrounding massive spirals at low redshift. As a test of these models, we use X-ray and H-alpha data of the two massive, quiescent edge-on spirals NGC 5746 and NGC 5170 to investigate the amount and origin of any hot gas in their halos. Contrary to our earlier claim, the Chandra analysis of NGC 5746, employing more recent calibration data, does not reveal any significant evidence for diffuse X-ray emission outside the optical disk, with a 3-sigma upper limit to the halo X-ray luminosity of 4e39 erg/s. An identical study of the less massive NGC 5170 also fails to detect any extraplanar X-ray emission. By extracting hot halo properties of disk galaxies formed in cosmological hydrodynamical simulations, we compare these results to expectations for cosmological accretion of hot gas by spirals. For Milky Way-sized galaxies, these high-resolution simulations predict hot halo X-ray luminosities which are lower by a factor of ~2 compared to our earlier results reported by Toft et al. (2002). We find the new simulation predictions to be consistent with our observational constraints for both NGC 5746 and NGC 5170, while also confirming that the hot gas detected so far around more actively star-forming spirals is in general probably associated with stellar activity in the disk. Observational results on quiescent disk galaxies at the high-mass end are nevertheless providing powerful constraints on theoretical predictions, and hence on the assumed input physics in numerical studies of disk galaxy formation and evolution.

There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data.

We present new measurements of the evolution in the Lyman break galaxy (LBG) population between z~4 and z~6. By utilizing the extensive multiwavelength datasets available in the GOODS fields, we identify 2443 B, 506 V, and 137 i'-band dropout galaxies likely to be at z~4, 5, and 6. With the goal of understanding the duration of typical star formation episodes in galaxies at z>4, we examine the distribution of stellar masses and ages as a function of cosmic time. We find that at a fixed rest-UV luminosity, the average stellar masses and ages of galaxies do not increase significantly between z~6 and 4. In order to maintain this near equilibrium in the average properties of high redshift LBGs, we argue that there must be a steady flux of young, newly-luminous objects at each successive redshift. When considered along with the short duty cycles inferred from clustering measurements, these results may suggest that galaxies are undergoing star formation episodes lasting only several hundred million years. In contrast to the unchanging relationship between the average stellar mass and rest-UV luminosity, we find that the number density of massive galaxies increases considerably with time over 4<z<6. Given this rapid increase of UV luminous massive galaxies, we explore the possibility that a significant fraction of massive (M*>1e11 Msun) z~2-3 distant red galaxies (DRGs) were in part assembled in an LBG phase at earlier times. Integrating the growth in the stellar mass function of actively forming LBGs over 4<z<6 down to z~2, we find that z>3 LBGs could have contributed significantly to the quiescent DRG population, indicating that the intense star-forming systems probed by current sub-millimeter observations are not the only route toward the assembly of DRGs at z~2.

We present a method for extracting the expected cosmological 21-cm signal from the epoch of reionization, taking into account contaminating radiations and random instrumental noise. The method is based on the maximum a-posteriori probability (MAP) formalism and employs the coherence of the contaminating radiation along the line-of-sight and the three-dimensional correlations of the cosmological signal. We test the method using a detailed and comprehensive modeling of the cosmological 21-cm signal and the contaminating radiation. The signal is obtained using a high resolution N-body simulation where the gas is assumed to trace the dark matter and is reionized by stellar radiation computed from semi-analytic galaxy formation recipes. We model contaminations to the cosmological signal from synchrotron and free-free galactic foregrounds and extragalactic sources including active galactic nuclei, radio haloes and relics, synchrotron and free-free emission from star forming galaxies, and free-free emission from dark matter haloes and the intergalactic medium. We provide tests of the reconstruction method for several rms values of instrumental noise from $\sigma_{N}=1$ to 250 mK. For low instrumental noise, the recovered signal, along individual lines-of-sight, fits the true cosmological signal with a mean rms difference of $d_{rms}\approx 1.7\pm 0.6$ for $\sigma_{N}=1$ mK, and $d_{rms}\approx 4.2\pm 0.4$ for $\sigma_{N}=5$ mK. The one-dimensional power spectrum is nicely reconstructed for all values of $\sigma_{N}$ considered here, while the reconstruction of the two-dimensional power spectrum and the Minkowski functionals is good only for noise levels of the order of few mK.

Only ~10% of baryons in the universe are in the form of stars, yet most models of luminous structure formation have concentrated on the properties of the luminous stellar matter. In this paper we focus on the "flip side" of galaxy formation and investigate the properties of the material that is not presently locked up in galaxies. This "by-product" of galaxy formation can be observed as an X-ray emitting plasma (the intracluster medium, hereafter ICM) in groups and clusters, and we present a version of the Durham semi-analytic galaxy formation model GALFORM that allows us to investigate the properties of the ICM. As we would expect on the basis of gravitational scaling arguments, the previous model (presented in Bower et al. 2006) fails to reproduce even the most basic observed properties of the ICM; however, we present a simple modification to the model to allow for heat input into the ICM from the AGN "radio mode" feedback. This heating acts to expel gas from the X-ray luminous central regions of the host halo. With this modification, the model reproduces the observed gas mass fractions and luminosity-temperature relation of groups and clusters. Introducing the heating process into the model requires changes to a number of model parameters in order to retain a good match to the observed galaxy properties. With the revised parameters, the best fitting luminosity function is comparable to that presented in Bower et al. (2006). The new model makes a fundamental step forward, providing a unified model of galaxy and cluster ICM formation. However, the detailed comparison with the data is not completely satisfactory, and we highlight key areas for improvement.

Current models of galaxy formation predict satellite galaxies in groups and clusters that are redder than observed. We investigate the effect on the colours of satellite galaxies produced by the ram pressure stripping of their hot gaseous atmospheres as the satellites orbit within their parent halo. We incorporate a model of the stripping process based on detailed hydrodynamic simulations within the Durham semi-analytic model of galaxy formation. The simulations show that the environment in groups and clusters is less aggressive than previously assumed. The main uncertainty in the model is the treatment of gas expelled by supernovae. With reasonable assumptions for the stripping of this material, we find that satellite galaxies are able to retain a significant fraction of their hot gas for several Gigayears, thereby replenishing their reservoirs of cold, star forming gas and remaining blue for a relatively long period of time. A bimodal distribution of galaxy colours, similar to that observed in SDSS data, is established and the colours of the satellite galaxies are in good agreement with the data. In addition, our model naturally accounts for the observed dependence of satellite colours on environment, from small groups to high mass clusters.

Only by incorporating various forms of feedback can theories of galaxy formation reproduce the present-day luminosity function of galaxies. It has also been argued that such feedback processes might explain the counter-intuitive behaviour of "downsizing" witnessed since redshifts z=1-2. To examine this question, observations spanning 0.4 < z < 1.4 from the DEEP2/Palomar survey are compared with a suite of equivalent mock observations derived from the Millennium Simulation, populated with galaxies using the Galform code. Although the model successfully reproduces the observed total mass function and the general trend of downsizing, it fails to accurately reproduce the colour distribution and type-dependent mass functions at all redshifts probed. This failure is shared by other semi-analytical models which collectively appear to "over-quench" star formation in intermediate-mass systems. These mock lightcones are also a valuable tool for investigating the reliability of the observational results in terms of cosmic variance. Using variance estimates derived from the lightcones we confirm the significance of the decline since z=1 in the observed number density of massive blue galaxies which, we argue, provides the bulk of the associated growth in the red sequence. We also assess the limitations arising from cosmic variance in terms of our ability to observe mass-dependent growth since z=1.

We have traced the past 7 Gyr of red galaxy stellar mass growth within dark matter halos. We have determined the halo occupation distribution, which describes how galaxies reside within dark matter halos, using the observed luminosity function and clustering of 40,696 0.2<z<1.0 red galaxies in Bootes. Half of 10^{11.9} Msun/h halos host a red central galaxy, and this fraction increases with increasing halo mass. We do not observe any evolution of the relationship between red galaxy stellar mass and host halo mass, although we expect both galaxy stellar masses and halo masses to evolve over cosmic time. We find that the stellar mass contained within the red population has doubled since z=1, with the stellar mass within red satellite galaxies tripling over this redshift range. In cluster mass halos most of the stellar mass resides within satellite galaxies and the intra-cluster light, with a minority of the stellar mass residing within central galaxies. The stellar masses of the most luminous red central galaxies are proportional to halo mass to the power of a third. We thus conclude that halo mergers do not always lead to rapid growth of central galaxies. While very massive halos often double in mass over the past 7 Gyr, the stellar masses of their central galaxies typically grow by only 30%.

Luminous red galaxies (LRGs) are much rarer and more massive than L* galaxies. Coupled with their extreme colours, LRGs therefore provide a demanding testing ground for the physics of massive galaxy formation. We present the first self-consistent predictions for the abundance and properties of LRGs in hierarchical structure formation models. We test two published models which use quite different mechanisms to suppress the formation of massive galaxies: the Bower et al. (2006) model, which invokes ``AGN-feedback'' to prevent gas from cooling in massive haloes, and the Baugh et al. (2005) model which relies upon a ``superwind'' to eject gas before it is turned into stars. Without adjusting any parameters, the Bower et al. model gives an excellent match to the observed luminosity function of LRGs in the SDSS (with a median redshift of z=0.24) and to their clustering; the Baugh et al. model is less successful in these respects. Both models fail to match the observed abundance of LRGs at z=0.5 to better than a factor of ~2. In the models, LRGs are typically bulge dominated systems with M* of ~2x10^11 h^{-1} M_sun and velocity dispersions of ~250 km s^{-1}. Around half of the stellar mass in the model LRGs is already formed by z~2.2 and is assembled into one main progenitor by z~1.5; on average, only 25% of the mass of the main progenitor is added after z~1. LRGs are predicted to be found in a wide range of halo masses, a conclusion which relies on properly taking into account the scatter in the formation histories of haloes. Remarkably, we find that the correlation function of LRGs is predicted to be a power law down to small pair separations, in excellent agreement with observational estimates. Neither the Bower et al. nor the Baugh et al. model is able to reproduce the observed radii of LRGs.

We place additional constraints on the three parameters of the dark matter halo merger rate function recently proposed by Parkinson, Cole & Helly by utilizing Smoluchowski's coagulation equation, which must be obeyed by any binary merging process which conserves mass. We find that the constraints from Smoluchowski's equation are degenerate, limiting to a thin plane in the three dimensional parameter space. This constraint is consistent with those obtained from fitting to N-body measures of progenitor mass functions, and provides a better match to the evolution of the overall dark matter halo mass function, particularly for the most massive halos. We demonstrate that the proposed merger rate function does not permit an exact solution of Smoluchowski's equation and, therefore, the choice of parameters must reflect a compromise between fitting various parts of the mass function. The techniques described herein are applicable to more general merger rate functions, which may permit a more accurate solution of Smoluchowski's equation. The current merger rate solutions are most probably sufficiently accurate for the vast majority of applications.

A follow-up survey using the Submillimetre High-Angular Resolution Camera (SHARC-II) at 350 microns has been carried out to map the regions around several 850 micron-selected sources from the Submillimetre HAlf Degree Extragalactic Survey (SHADES). These observations probe the infrared luminosities and hence star-formation rates in the largest existing, most robust sample of submillimetre galaxies (SMGs). We measure 350 micron flux densities for 24 850 micron sources, seven of which are detected at >2.5-sigma within a 10 arcsec search radius of the 850 micron positions. When results from the literature are included the total number of 350 micron flux density constraints of SHADES SMGs is 31, with 15 detections. We fit a modified blackbody to the far-infrared (FIR) photometry of each SMG, and confirm that typical SMGs are dust-rich (Mdust~9x10^8 Msun), luminous (Lfir~2x10^12 Lsun), star-forming galaxies with intrinsic dust temperatures of ~35 K and star-formation rates of ~400 Msun/yr. We have measured the temperature distribution of SMGs and find that the underlying distribution is slightly broader than implied by the error bars, and that most SMGs are at 28 K with a few hotter. We also place new constraints on the 350 micron source counts, N350(>25mJy)~200-500 deg^-2.

Inspired by recent work on feedback in disk galaxies (Efstathiou 2000, Silk 2003) and on the angular momentum distribution in simulated gas halos (Sharma and Steinmetz 2005), a fully dynamic model of disk galaxy formation and evolution has been developed. This is used to demonstrate how observed galactic systems could have formed from halos similar to those found in simulations and applies physically motivated models of star formation and feedback to explore whether the true nature of these processes would be manifest from local and cosmological observables. This is made possible by computational integration with the galaxy formation model developed originally by the group at Durham University (Cole et al. 2000).

The formation of galactic discs and the efficiency of star formation within them are issues central to our understanding of galaxy formation. We have developed a detailed and versatile model of disc formation which combines the strengths of previous studies of isolated discs with those of hierarchical galaxy formation models. Disc structure is inferred from the distribution of angular momentum in hot halo gas and the hierarchical build-up of dark matter, leading to theoretically generated systems where the evolution of surface density, rotation, velocity dispersion, stability and metallicity is predicted for annular regions of width 20-100 pc. The model will be used to establish whether the accepted theory of large-scale structure formation in the universe is consistent with observed trends in the properties of disc galaxies. This first paper explicitly examines the importance of embedding such calculations within a merging hierarchy of dark matter haloes, finding that this leads to dramatically different formation histories compared to models in which discs grow in isolation. Different models of star formation are explored, and are found to have only a secondary influence on the properties of the resulting galaxy discs, the main governing factor being the infalling gas supply from the hot halo.

We present the redshift distribution of the SHADES galaxy population based on the rest-frame radio-mm-FIR colours of 120 robustly detected 850um sources in the Lockman Hole East (LH) and Subaru XMM-Newton Deep Field (SXDF). The redshift distribution derived from the full SED information is shown to be narrower than that determined from the radio-submm spectral index, as more photometric bands contribute to a higher redshift accuracy. The redshift distribution of sources derived from at least two photometric bands peaks at z ~ 2.4 and has a near-Gaussian distribution, with 50 per cent (interquartile range) of sources at z=1.8-3.1. We find a statistically-significant difference between the measured redshift distributions in the two fields; the SXDF peaking at a slightly lower redshift (median z ~ 2.2) than the LH (median z ~ 2.7), which we attribute to the noise-properties of the radio observations. We demonstrate however that there could also be field-to-field variations that are consistent with the measured differences in the redshift distributions, and hence, that the incomplete area observed by SHADES with SCUBA, despite being the largest sub-mm survey to date, may still be too small to fully characterize the bright sub-mm galaxy population. Finally we present a brief comparison with the predicted, or assumed, redshift distributions of sub-mm galaxy formation and evolution models, and we derive the contribution of these SHADES sources and the general sub-mm galaxy population to the star formation-rate density at different epochs.

Using the Galactica code of Benson et al., we obtain quantitative measurements of spheroid-to-disk ratios for a sample of 8839 galaxies observed in the Sloan Digital Sky Survey. We carry out extensive tests of this code and of Gim2D, finding that they perform similarly in all respects. From the spheroid and disk luminosities, we construct luminosity and stellar mass functions for each component and estimate the relative luminosity and stellar mass densities of disks and spheroids in the local Universe. Assuming a simple one-to-one mapping between between spheroid mass and the mass of a central supermassive black hole, we provide the most accurate determination so far of the black hole mass function in the local universe. From this, we infer a cosmological mass density of black holes of rho_BH=3.77 +/- 0.97 x 10^5 h M_sun Mpc^-3. We compare our results to predictions from current hierarchical models of galaxy formation and these are found to fare well in predicting the qualitative trends observed. We find that stars in disks contribute 35-51% of the local stellar mass density.

We present the first identification of large-scale structures (LSS) at z $< 1.1$ in the Cosmic Evolution Survey (COSMOS). The structures are identified from adaptive smoothing of galaxy counts in the pseudo-3d space ($\alpha,\delta$,z) using the COSMOS photometric redshift catalog. The technique is tested on a simulation including galaxies distributed in model clusters and a field galaxy population -- recovering structures on all scales from 1 to 20\arcmin without {\it a priori} assumptions for the structure size or density profile. Our procedure makes {\bf no} {\it a priori} selection on galaxy spectral energy distribution (SED, for example the Red Sequence), enabling an unbiased investigation of environmental effects on galaxy evolution. The COSMOS photometric redshift catalog yields a sample of $1.5\times10^5$ galaxies with redshift accuracy, $\Delta z_{FWHM}/(1+z) \leq 0.1$ at z $< 1.1$ down to I$_{AB} \leq 25$ mag. Using this sample of galaxies, we identify 42 large-scale structures and clusters. abstract truncated for astroph 25 line limit -- see preprint

The Cosmic Evolution Survey (COSMOS) was initiated with an extensive allocation (590 orbits in Cycles 12-13) using the Hubble Space Telescope (HST) for high resolution imaging. Here we review the characteristics of the HST imaging with the Advanced Camera for Surveys (ACS) and parallel observations with NICMOS and WFPC2. A square field (1.8$\sq$\deg) has been imaged with single-orbit ACS I-F814W exposures with 50% completeness for sources 0.5\arcsec in diameter at I$_{AB} $ = 26.0 mag. The ACS imaging is a key part of the COSMOS survey, providing very high sensitivity and high resolution (0.09\arcsec FWHM, 0.05\arcsec pixels) imaging and detecting 1.2 million objects to a limiting magnitude of 26.5 (AB). These images yield resolved morphologies for several hundred thousand galaxies. The small HST PSF also provides greatly enhanced sensitivity for weak lensing investigations of the dark matter distribution.

We examine the problem tidally-induced mass loss from collisionless systems such as dark matter haloes. We develop a model for tidal mass loss, based upon the phase space distribution of particles, which accounts for how both tidal and Coriolis torques perturb the angular momentum of each particle in the system. This allows us to study how both the density profile and velocity anisotropy affect the degree of mass loss--we present basic results from such a study. Our model predicts that mass loss is a continuous process even in a static tidal field, a consequence of the fact that mass loss weakens the potential of the system making it easier for further mass loss to occur. We compare the predictions of our model with N-body simulations of idealized systems in order to check its validity. We find reasonable agreement with the N-body simulations except for in the case of very strong tidal fields, where our results suggest that a higher-order perturbation analysis may be required. The continuous tidally-induced mass loss predicted by our model can lead to substantial reduction in satellite mass in cases where the traditional treatment predicts no mass loss. As such, our results may have important consequences for the orbits and survival of low mass satellites in dark matter haloes.

Models of disk galaxy formation commonly predict the existence of an extended reservoir of hot gas surrounding massive spirals at low redshift. As a test of these models, we have obtained X-ray and optical data of the two massive edge-on spirals NGC 5746 and NGC 5170, in order to investigate the amount and origin of hot gas in their disks and halos. Chandra observations of NGC 5746 reveal evidence for diffuse X-ray emission with a total luminosity of ~7 x 10^39 erg/s surrounding this galaxy out to at least ~20 kpc from the disk, whereas an identical study of the less massive NGC 5170 fails to detect any extraplanar X-ray emission. Unlike the case for other disk galaxies with detected X-ray halos, the halo emission around NGC 5746 is not accompanied by extraplanar H-alpha or radio emission, and there is no evidence for significant nuclear or starburst activity in the disk. In contrast to these other cases, the emission around NGC 5746 therefore appears to arise from the cooling of externally accreted material rather than from disk outflows. To verify this idea, we present results of cosmological simulations of galaxy formation and evolution, showing our observations to be in good agreement with expectations for cosmological accretion, while also confirming that the X-ray halos of other spirals do not fit well into an accretion scenario. We find that the estimated cooling rate of hot halo gas around NGC 5746 would provide sufficient material for star formation in the disk to proceed at its present rate. This lends support to the idea that a supply of hot ambient gas is potentially available as fuel for star formation in massive, nearby spirals, and suggests that accretion of hot gas could be important for maintaining the stellar disks of such galaxies.

We trace the assembly history of red galaxies since z=1, by measuring their evolving space density with the B-band luminosity function. Our sample of 39599 red galaxies, selected from 6.96 square degrees of imaging from the NOAO Deep Wide-Field and Spitzer IRAC Shallow surveys, is an order of magnitude larger, in size and volume, than comparable samples in the literature. We measure a higher space density of z=0.9 red galaxies than some of the recent literature, in part because we account for the faint yet significant galaxy flux which falls outside of our photometric aperture. The B-band luminosity density of red galaxies, which effectively measures the evolution of ~L* galaxies, increases by only 36 percent from z=0 to z=1. If red galaxy stellar populations have faded by 1.24 B-band magnitudes since z=1, the stellar mass contained within the red galaxy population has roughly doubled over the past 8 Gyr. This is consistent with star-forming galaxies being transformed into ~L* red galaxies after a decline in their star formation rates. In contrast, the evolution of 4L* red galaxies differs only slightly from a model with negligible star formation and no galaxy mergers since z=1. If this model approximates the luminosity evolution of red galaxy stellar populations, then 80 percent of the stellar mass contained within today's 4L* red galaxies was already in place at z=0.7. While red galaxy mergers have been observed, such mergers do not produce rapid growth of 4L* red galaxy stellar masses between z=1 and the present day.

The merger of two supermassive black holes is expected to produce a gravitational-wave signal detectable by the satellite LISA. The rate of supermassive-black-hole mergers is intimately connected to the halo merger rate, and the extended Press-Schechter formalism is often employed when calculating the rate at which these events will be observed by LISA. This merger theory is flawed and provides two rates for the merging of the same pair of haloes. We show that the two predictions for the LISA supermassive-black-hole-merger event rate from extended Press-Schechter merger theory are nearly equal because mergers between haloes of similar masses dominate the event rate. An alternative merger rate may be obtained by inverting the Smoluchowski coagulation equation to find the merger rate that preserves the Press-Schechter halo abundance, but these rates are only available for power-law power spectra. We compare the LISA event rates derived from the extended Press-Schechter merger formalism to those derived from the merger rates obtained from the coagulation equation and find that the extended Press-Schechter LISA event rates are thirty percent higher for a power spectrum spectral index that approximates the full Lambda-CDM result of the extended Press-Schechter theory.

Recent observations of the distant Universe suggest that much of the stellar mass of bright galaxies was already in place at $z>1$. This presents a challenge for models of galaxy formation because massive halos are assembled late in hierarchical cosmologies such as cold dark matter (CDM). In this paper, we discuss a new implementation of the Durham semi-analytic model in which feedback due to active galactic nuclei (AGN) is assumed to quench cooling flows in massive halos. This mechanism naturally creates a break in the local galaxy luminosity function at bright magnitudes. The model is implemented within the Millennium N-body simulation; the accurate dark matter merger trees and large number of realizations of the galaxy formation process that the simulation provides results in highly accurate statistics. After adjusting the values of the physical parameters in the model by reference to the properties of local galaxies, we use it to investigate the evolution of the K-band luminosity and galaxy stellar mass functions. We also calculate the volume averaged star formation rate density of the Universe as a function of redshift and the way in which this is apportioned amongst galaxies of different mass. The model robustly predicts a substantial population of massive galaxies out to redshift $z\sim 5$ and a star formation rate density which rises with increasing redshift in objects of all masses. Although observational data on these properties have been cited as evidence for ``anti-hierarchical'' galaxy formation, we find that when AGN feedback is taken into account, the fundamentally hierachical CDM model provides a very good match to these observations.

We present a hierarchical galaxy formation model which can account for the number counts of sources detected through their emission at sub-millimetre wavelengths. The first stage in our approach is an ab initio calculation of the star formation histories for a representative sample of galaxies, which is carried out using the semi-analytical galaxy formation model GALFORM. These star formation histories are then input into the spectro-photometric code GRASIL, to produce a spectral energy distribution for each galaxy. Dust extinction and emission are treated self consistently in our model, without having to resort to ad-hoc assumptions about the amount of attenuation by dust or the temperature at which the dust radiates. We argue that it is necessary to modify the form of the stellar initial mass function in starbursts in order to match the observed number of sub-mm sources, if we are to retain the previous good matches enjoyed between observations and model predictions in the local universe. We also list some other observational tests that have been passed by our model.

We have modelled the process of reionization of the IGM by photoionization by galaxies, in order to predict the epoch of reionization. We use a sophisticated semi-analytic model to track the formation of galaxies. Our study represents a much more complete and physically consistent modelling of reionization than has been conducted in the past, containing significant improvements in the modelling of the collapse of baryons into dark matter halos, and in the model for the cooling and condensation of gas within halos (including photoheating from a self-consistently computed ionizing background and cooling due to H2). We find that reionization can be achieved by z~10-20 in a LCDM cosmological model with sigma8~0.9. However, a cosmological model with a running spectral index is only able to achieve reionization before z~9 if extreme assumptions are made about the physics of feedback at high redshifts. We also consider the galaxy formation model discussed by Baugh et al., and find that it is able to reionize the Universe by z~12. The previous results assume that all of the ionizing photons produced by stars in galaxies are able to escape and ionize the IGM. If this is not the case, then the redshift of reionization could be substantially reduced. We find that extended periods of partial reionization and double reionizations can occur in models in which the first stars formed via cooling by H2 molecules, are very massive, and in which the escape fraction of ionizing photons ~10-30%. Such models do not fully reionize until z~6-7, but predict an electron scattering optical depth as large as 0.15. Models with lower sigma8=0.7-0.8 as suggested by the recent WMAP three year data have reduced redshifts of reionization, but can be consistent with the lower optical depth also suggested by the WMAP three year data.

Hot gaseous haloes surrounding galaxies and extending well beyond the distribution of stars are a ubiquitous prediction of galaxy formation scenarios. The haloes are believed to consist of gravitationally trapped gas with a temperature of millions of Kelvin. The existence of such hot haloes around massive elliptical galaxies has been established through their X-ray emission. While gas out-flowing from starburst spiral galaxies has been detected, searches for hot haloes around normal, quiescent spiral galaxies have so far failed, casting doubts on the fundamental physics in galaxy formation models. Here we present the first detection of a hot, large-scale gaseous halo surrounding a normal, quiescent spiral galaxy, NGC 5746, alleviating a long-standing problem for galaxy formation models. In contrast to starburst galaxies, where the X-ray halo can be powered by the supernova energy, there is no such power source in NGC 5746. The only compelling explanation is that we are here witnessing a galaxy forming from gradually in-flowing hot and dilute halo gas.

We explore whether stellar tidal streams can provide information on the secular, cosmological evolution of the Milky Way's gravitational potential and on the presence of subhalos. We carry out long-term (~t_hubble) N-body simulations of disrupting satellite galaxies in a semi-analytic Galaxy potential where the dark matter halo and the subhalos evolve according to a LCDM cosmogony. All simulations are constrained to end up with the same position and velocity at present. Our simulations account for: (i) the secular evolution of the host halo's mass, size and shape, (ii) the presence of subhalos and (iii) dynamical friction. We find that tidal stream particles respond adiabatically to the Galaxy growth so that, at present, the energy and angular momentum distribution is exclusively determined by the present Galaxy potential. In other words, all present-day observables can only constrain the present mass distribution of the Galaxy independent of its past evolution. We also show that, if the full phase-space distribution of a tidal stream is available, we can accurately determine (i) the present Galaxy's shape and (ii) the amount of mass loss from the stream's progenitor, even if this evolution spanned a cosmologically significant epoch.

The halo occupation distribution (HOD) describes the relation between galaxies and dark matter at the level of individual dark matter halos. The properties of galaxies residing at the centers of halos differ from those of satellite galaxies because of differences in their formation histories. Using a smoothed particle hydrodynamics (SPH) simulation and a semi-analytic (SA) galaxy formation model, we examine the separate contributions of central and satellite galaxies to the HOD, more specifically to the probability P(N|M) that a halo of virial mass M contains N galaxies of a particular class. In agreement with earlier results for dark matter subhalos, we find that the mean occupation function <N> for galaxies above a baryonic mass threshold can be approximated by a step function for central galaxies plus a power law for satellites, and that the distribution of satellite numbers is close to Poisson at fixed halo mass. For galaxy samples defined by different baryonic mass thresholds, there is a nearly linear relation between the minimum halo mass Mmin required to host a central galaxy and the mass M1 at which an average halo hosts one satellite, with M1 ~ 14 Mmin (SPH) or M1 ~ 18 Mmin (SA). The mean occupation number of young galaxies exhibits a local minimum at M ~ 10 Mmin where halos are too massive to host a young central galaxy but not massive enough to host satellites. We show that the conditional galaxy mass function at fixed halo mass cannot be described by a Schechter function because central galaxies produce a "bump" at high masses. We suggest parameterizations for the HOD and the conditional luminosity function that can be used to model observed galaxy clustering. Many of our predictions are in good agreement with recent results inferred from clustering in the Sloan Digital Sky Survey.

We present a Monte-Carlo model of He II reionization by QSOs and its effect on the thermal state of the clumpy intergalactic medium (IGM). The model assumes that patchy reionization develops as a result of the discrete distribution of QSOs. It includes various recipes for the propagation of the ionizing photons, and treats photo-heating self-consistently. The model provides the fraction of He III, the mean temperature in the IGM, and the He II mean optical depth -- all as a function of redshift. It also predicts the evolution of the local temperature versus density relation during reionization. Our findings are as follows: The fraction of He III increases gradually until it becomes close to unity at $z\sim 2.8-3.0$. The He II mean optical depth decreases from $\tau\sim 10$ at $z\geq 3.5$ to $\tau\leq 0.5$ at $z\leq 2.5$. The mean temperature rises gradually between $z\sim 4$ and $z\sim 3$ and declines slowly at lower redshifts. The model predicts a flattening of the temperature-density relation with significant increase in the scatter during reionization at $z\sim 3$. Towards the end of reionization the scatter is reduced and a tight relation is re-established. This scatter should be incorporated in the analysis of the Ly$\alpha$ forest at $z\leq 3$. Comparison with observational results of the optical depth and the mean temperature at moderate redshifts constrains several key physical parameters.

We present predictions for the abundance of sub-mm galaxies (SMGs) and Lyman-break galaxies (LBGs) in the $\Lambda$CDM cosmology. A key feature of our model is the self-consistent calculation of the absorption and emission of radiation by dust. The new model successfully matches the LBG luminosity function, as well reproducing the properties of the local galaxy population in the optical and IR. The model can also explain the observed galaxy number counts at $850\mum$, but only if we assume a top-heavy IMF for the stars formed in bursts. The predicted redshift distribution of SMGs depends relatively little on their flux over the range 1-$10\mjy$, with a median value of $z\approx 2.0$ at a flux of $5\mjy$, in very good agreement with the recent measurement by Chapman et al The counts of SMGs are predicted to be dominated by ongoing starbursts. However, in the model these bursts are responsible for making only a few per cent of the stellar mass locked up in massive ellipticals at the present day.

We develop a semi-analytical model that determines the evolution of the mass, position and internal structure of dark matter substructures orbiting in dark matter haloes. We apply this model to the case of the Milky Way. We focus in particular on the effects of mass loss, dynamical friction and substructure--substructure interactions, the last of which has previously been ignored in analytic models of substructure evolution. Our semi-analytical treatment reproduces both the spatial distribution of substructures and their mass function as obtained from the most recent N-body cosmological calculations of Gao et al. (2004). We find that, if mass loss is taken into account, the present distribution of substructures is practically insensitive to dynamical friction and scatterings from other substructures. Implementing these phenomena leads to a slight increase (~5%) in the number of substructures at r<0.35 r_vir, whereas their effects on the mass function are negligible. We find that mass loss processes lead to the disruption of substructures before dynamical friction and gravitational scattering can significantly alter their orbits. Our results suggest that the present substructure distribution at r>0.35 r_vir reflects the orbital properties at infall and is, therefore, purely determined by the dark matter environment around the host halo and has not been altered by dynamical evolution.

We present distributions of the orbital parameters of dark matter substructures at the time of merging into their host halo. Accurate knowledge of the orbits of dark matter substructures is a crucial input to studies which aim to assess the effects of the cluster environment on galaxies, the heating of galaxy disks and many other topics. Orbits are measured for satellites in a large number of N-body simulations. We focus on the distribution of radial and tangential velocities, but consider also distributions of orbital eccentricity and semi-major axis. We show that the distribution of radial and tangential velocities has a simple form and provide a fitting formula for this distribution. We also search for possible correlations between the infall directions of pairs of satellites, finding evidence for positive correlation at small angular separations as expected if some infall occurs along filaments. We also find (weak) evidence for correlations between the direction of the infall and infall velocity and the spin of the host halo.

The rate of merging of dark-matter halos is an absolutely essential ingredient for studies of both structure and galaxy formation. Remarkably, however, our quantitative understanding of the halo merger rate is still quite limited, and current analytic descriptions based upon the extended Press-Schechter formalism are fundamentally flawed. We show that a mathematically self-consistent merger rate must be consistent with the evolution of the halo abundance in the following sense: The merger rate must, when inserted into the Smoluchowski coagulation equation, yield the correct evolution of the halo abundance. We then describe a numerical technique to find merger rates that are consistent with this evolution. We present results from a preliminary study in which we find merger rates that reproduce the evolution of the halo abundance according to Press-Schechter for power-law power spectra. We discuss the limitations of the current approach and outline the questions that must still be answered before we have a fully consistent and correct theory of halo merger rates.

We develop an analytic model to calculate the rate at which galaxy disks are heated by dark matter substructures orbiting in their halos. The model takes into account the internal structure, mass function and accretion rate of satellites expected in the LambdaCDM cosmology, as well as the growth of the disk by accretion and mergers, but it ignores resonant heating of the disk and the dynamical effects of spiral arms and bars. We calibrate this model against N-body simulations and demonstrate that it is able to reproduce the N-body heating rates to within a factor of 3 in the majority of cases. Our model gives the distribution of disk scale-heights for galaxies of different luminosities. For L* spiral galaxies, it predicts a median disk thickness of only 5% of the radial scale-length if substructure is the only source of heating. The median disk thickness increases to nearly 20% of the radial scale-length when heating due to gravitational scattering of stars by molecular clouds is also included. The latter value is close to the thickness estimated observationally for the disk of the Milky Way galaxy. The distribution of disk thickness predicted by the model is also consistent with a recent observational determination for sub-L* galaxies by Bizyaev & Mitronova. Thus, the observed thickness of the stellar disks of spiral galaxies seems to be entirely compatible with the abundance of substructure in dark matter halos predicted by the standard Lambda-dominated cold dark matter model of structure formation. In an Omega_0=1 universe, our best model of galaxy formation produces similar scale-heights, a consequence of the fact that similar amounts of substructure are accreted by halos during the lifetime of the disk in Omega_0=1 and Omega_0=0.3, Lambda_0=0.7 cold dark matter cosmologies.

We investigate the physical mechanisms that shape the luminosity function. Beginning with the mass function of dark matter halos, we show how gas cooling, photoionization, feedback, merging and thermal conduction affect the shape of the luminosity function. We consider three processes whereby supernovae can affect the forming galaxy: (1) reheating of disk gas to the halo temperature; (2) expansion of the diffuse halo gas; (3) expulsion of cold disk gas from the halo. While feedback of form (1) is able to flatten the faint end of the luminosity function, this alone does not produce the sharp cut-off observed at large luminosities. Feedback of form (2) is also unable to solve this problem. The relative paucity of very bright galaxies can only be explained if cooling in massive halos is strongly suppressed. Conduction is a promising mechanism, but an uncomfortably high efficiency is required to suppress cooling to the desired level. If, instead, superwinds are responsible for the lack of bright galaxies, then the total energy budget required to obtain a good match to the galaxy luminosity function greatly exceeds the energy available from supernova explosions. The mechanism is only viable if the formation of central supermassive black holes play a crucial role in limiting the amount of stars that form in the host galaxy. (abridged)

Winds from pregalactic starbursts and 'miniquasars' may pollute the IGM with metals and raise its temperature to a high adiabat, and so inhibit the formation of early galaxies. We compute the thermal history of the IGM when it experiences a period of rapid, homogeneous "preheating" at high redshifts. Measurements of the temperature of the Lyamn-alpha forest at z~3 constrain the redshift and energy of preheating, and rule out models that preheat too late or to too high a temperature. We predict galaxy luminosity functions in preheated universes. The results depend crucially on whether the baryonic smoothing scale in the IGM is computed globally, or in a local, density-dependent fashion. Using a globally averaged smoothing scale, we find that models with excessive preheating produce too few L_* and fainter galaxies, and are therefore inconsistent with observational data. A density-dependent smoothing scale requires more energetic preheating to achieve the same degree of suppression in the faint-end slope. All models, however, appear unable to explain the sharp cut-off in the luminosity function at bright magnitudes. Supernova-driven preheating scenarios tend to raise the mean metallicity of the universe well above the minimum levels observed in the Lyman-alpha clouds. We find that ionizing photon escape fractions must be significantly higher than 10% in order to explain the low inferred HI fraction at z~6. While early preheating causes strong suppression of dwarf galaxy formation we show that it is not able to reproduce the observed abundance of satellite galaxies in the Local Group in detail.

The halo occupation distribution (HOD) describes the bias between galaxies and dark matter by specifying (a) the probability P(N|M) that a halo of virial mass M contains N galaxies of a particular class and (b) the relative distributions of galaxies and dark matter within halos. We calculate predicted HODs for a Lambda-CDM cosmological model using an SPH hydrodynamic simulation and a semi-analytic (SA) galaxy formation model. Although the two methods predict different galaxy mass functions, their HOD predictions agree remarkably well. For mass-selected samples, the mean occupation <N(M)> exhibits a sharp cutoff at low halo masses, a slowly rising plateau for <N>~1-2, and a more steeply rising high occupancy regime. At low <N>, the mean pair and triple counts are well below Poisson expectations, with important consequences for small scale behavior of 2- and 3-point correlation functions. The HOD depends strongly on galaxy age, with high mass halos populated mainly by old galaxies and low mass halos by young galaxies. The SPH simulation supports several simplifying assumptions about HOD bias: the most massive galaxy in a halo usually lies close to the center and moves near the halo's mean velocity; satellite galaxies have the same radial profile and velocity dispersion as the dark matter; and the mean occupation at fixed halo mass is independent of the halo's larger scale environment. By applying the SPH and SA HODs to a large volume N-body simulation, we show that both methods predict slight, observable departures from a power-law galaxy correlation function. The predicted HODs are closely tied to the underlying galaxy formation physics, they offer useful guidance to theoretical models of galaxy clustering, and they will be tested empirically by ongoing analyses of galaxy redshift surveys. (Shortened)

Both preheating of the intergalactic medium and radiative cooling of low entropy gas have been proposed to explain the deviation from self-similarity in the cluster L_x-T_x relation and the observed entropy floor in these systems. However, severe overcooling of gas in groups is necessary for radiative cooling alone to explain the observations. Non-gravitational entropy injection must therefore still be important in these systems. We point out that on scales of groups and below, gas heated to the required entropy floor cannot cool in a Hubble time, regardless of its subsequent adiabatic compression. Preheating therefore shuts off the gas supply to galaxies, and should be an important global feedback mechanism for galaxy formation. Constraints on global gas cooling can be placed from the joint evolution of the comoving star formation rate and neutral gas density. Preheating at high redshift can be ruled out; however the data does not rule out passive gas consumption without inflow since z~2. Since for preheated gas t_cool > t_dyn, we speculate that preheating could play a role in determining the Hubble sequence: at a given mass scale, high sigma peaks in the density field collapse early to form ellipticals, while low sigma peaks collapse late and quiescently accrete preheated gas to form spirals. The entropy produced by large scale shock-heating of the intergalatic medium is significant only at late times, z<1, and cannot produce these effects.

Using semi-analytic modeling techniques, we calculate the luminosity function of galaxy populations residing in cold dark matter halos of different mass. We pay particular attention to the influence of the reionization of the Universe on the number of faint galaxies and to the effects of dynamical friction and tidal limitation of satellites on the number of bright galaxies. We find substantial differences in the shapes of the galaxy luminosity functions in halos of different mass which reflect generic features of the cold dark matter model of galaxy formation and thus offer the opportunity to test it. We then consider how the individual halo luminosity functions combine together to produce the global luminosity function. Surprisingly, the global function ends up having a shallower faint end slope than those of the constituent halo luminosity functions. We compare our model predictions with the limited datasets compiled by Trentham & Hodgkin. We find good agreement with the luminosity functions measured in the Virgo and Coma clusters but significant disagreement with the luminosity functions measured in the Local Group and Ursa Minor cluster. We speculate on possible inadequacies in our modeling and in the existing observational samples. The luminosity functions of galaxies in groups and clusters identified in the 2dF and SDSS galaxy redshift surveys offer the prospect of testing galaxy formation models in detail.

We develop a hybrid galaxy formation model which uses outputs from an N-body simulation to follow the merger histories (or ``merger trees'') of dark matter halos and treats baryonic processes, such as the cooling of gas within halos and subsequent star formation, using the semi-analytic model of Cole et al. We compare this hybrid model to an otherwise identical model which utilises merger tree realisations generated by a Monte-Carlo algorithm and find that, apart from the limited mass resolution imposed by the N-body particle mass, the only significant differences between the models are due to the known discrepancy between the distribution of halo progenitor masses predicted by extended Press Schechter theory and that found in N-body simulations. We investigate the effect of limited mass resolution on the hybrid model by comparing to a purely semi-analytic model with greatly improved mass resolution. We find that the mass resolution of the simulation we use, which has a particle mass of 1.4x10^10 Msolar/h, is insufficient to produce a reasonable luminosity function for galaxies with magnitudes in the b_J band fainter than -17.

We compare the results of two techniques used to calculate the evolution of cooling gas during galaxy formation: Smooth Particle Hydrodynamics (SPH) simulations and semi-analytic modelling. We improve upon the earlier statistical comparison of Benson et al. by taking halo merger histories from the dark matter component of the SPH simulation, which allows us to compare the evolution of galaxies on an object-by-object basis in the two treatments. We use a ``stripped-down'' version of the semi-analytic model described by Helly et al. which includes only shock heating and radiative cooling of gas and which is adjusted to mimic the resolution and other parameters of a comparison SPH simulation as closely as possible. We compare the total mass of gas that cools in halos of different mass as a function of redshift as well as the masses and spatial distribution of individual ``galaxies.'' At redshift z=0, the cooled gas mass in well-resolved halos agrees remarkably well (to better than ~20%) in the SPH simulation and stripped-down semi-analytic model. At high redshift, resolution effects in the simulation become increasingly important and, as a result, more gas tends to cool in low mass halos in the SPH simulation than in the semi-analytic model. The cold gas mass function of individual galaxies in the two treatments at z=0 also agrees very well and, when the effects of mergers are accounted for, the masses of individual galaxies and their 2-point correlation functions are also in excellent agreement in the two treatments. Thus, our comparison confirms and extends the earlier conclusion of Benson et al. that SPH simulations and semi-analytic models give consistent results for the evolution of cooling galactic gas.

We present predictions for numerous statistics related to the presence of voids in the distribution of galaxies in a cold dark matter model of structure formation using a semi-analytic model of galaxy formation. Our study is able to probe galaxies with masses as low as 10^9Msun/h corresponding to absolute magnitudes of M_bJ-5log(h)=-18.1 and M_r-5log(h)=-18.7. We quantify the void and underdense probability functions, distributions of nearest neighbour distances and void sizes and compute the density profiles of voids. These results are contrasted with the expectations for dark matter (and the difference examined in terms of the galaxy/dark matter biasing relation) and are compared to analytic predictions and observational data where available. The predicted void probability functions are consistent with those measured from the Center for Astrophysics redshift surveys given the rather large uncertainties in this relatively small (for studies of voids) observational sample. The size of the observational sample is too small to probe the bias between galaxies and dark matter that we predict. We also examine the predicted properties of galaxies living within voids and contrast these with the general galaxy population. Our predictions are aimed at forthcoming large galaxy redshift surveys which should for the first time provide statistically accurate measures of the void population.

We present predictions for two statistical measures of the hydrogen reionization process at high redshift. The first statistic is the number of neutral segments identified in spectra of high redshift QSOs as a function of their length. The second is the cross-correlation of neutral regions with possible sources of ionizing radiation. These independent probes are sensitive to the topology of the ionized regions. If reionization proceeded from high to low density regions then the cross-correlation will be negative, while if voids were ionized first then we expect a positive correlation and a relatively small number of long neutral segments. We test the sensitivity of these statistics for reionization by stars in high redshift galaxies. The flux of ionizing radiation emitted from stars is estimated by identifying galaxies in an N-body simulation using a semi-analytic galaxy formation model. The spatial distribution of ionized gas is traced in various models for the propagation of the ionization fronts. A model with ionization proceeding from high to low density regions is consistent with the observations of Becker et al. (2001), while models in which ionization begins in the lowest density regions appear to be inconsistent with the present data.

We present a method to obtain quantitative measures of galaxy morphology and apply it to a spectroscopic sample of field galaxies in order to determine the luminosity and stellar mass functions of galactic disks and spheroids. We estimate, for each galaxy, the bulge-to-disk luminosity ratio in the I-band using a two-dimensional image fitting procedure. Monte Carlo simulations indicate that reliable determinations are only possible for galaxies approximately two magnitudes brighter than the photometric completeness limit, leaving a sample of 90 galaxies with well determined bulge-to-total light ratios. We construct the luminosity functions of disks and spheroids and, using a stellar population synthesis model, we estimate the stellar mass functions of each of these components. The disk and spheroid luminosity functions are remarkably similar. We do, however, find evidence in the bi-variate luminosity function that spheroid-dominated galaxies occur only among the brightest spheroids, while disk-dominated galaxies span a much wider range of disk luminosities. Remarkably, the total stellar mass residing in disks and spheroids is approximately the same. For our sample, we find the ratio of stellar masses in disks and spheroids to be 1.3+/-0.2. Ongoing large photometric and redshift surveys will lead to a large increase in the number of galaxies to which our techniques can be applied and thus to an improvement in the current estimates. (abridged)

Advances in extragalactic astronomy have prompted the development of increasingly realistic models which aim to describe the formation and evolution of galaxies. We review the philosophy behind one such technique, called semi-analytic modelling, and explain the relation between this approach and direct simulations of gas dynamics. Finally, we present model predictions for the evolution of the stellar mass of galaxies in a universe in which structure formation is hierarchical.

We re-examine the assembly history of field spheroidals as a potentially powerful discriminant of galaxy formation models. Whereas monolithic collapse and hierarchical, merger-driven, models suggest radically different histories for these galaxies, neither the theoretical predictions nor the observational data for field galaxies have been sufficiently reliable for precise conclusions to be drawn. A major difficulty in interpreting the observations, reviewed here, concerns the taxonomic definition of spheroidals in merger-based models. Using quantitative measures of recent star formation activity drawn from the internal properties of a sample of distant field galaxies in the Hubble Deep Fields, we undertake a new analysis to assess the continuous formation of spheroidal galaxies. Whereas abundances and redshift distributions of modelled spheroidals are fairly insensitive to their formation path, we demonstrate that the distribution and amount of blue light arising from recent mergers provides a more sensitive approach. With the limited resolved data currently available, the rate of mass assembly implied by the observed colour inhomogeneities is compared to that expected in popular Lambda-dominated cold dark matter models of structure formation. These models produce as many highly inhomogeneous spheroidals as observed, but underpredict the proportion of homogeneous, passive objects. We conclude that colour inhomogeneities, particularly when combined with spectroscopic diagnostics for large, representative samples of field spheroidals, will be a more valuable test of their physical assembly history than basic source counts and redshift distributions. Securing such data should be a high priority for the Advanced Camera for Surveys on Hubble Space Telescope.

We show that the distribution of the sizes and temperatures of clusters can be used to constrain cosmological models. The size-temperature (ST) distribution predicted in a flat Gaussian cluster-abundance-normalized Omega_0=0.3 model agrees well with the fairly tight ST relation observed. A larger power-spectrum amplitude sigma_8 would give rise to a larger scatter about the ST relation as would a larger value of Omega_0 and/or long non-Gaussian high-density tails in the probability density function. For Gaussian initial conditions, the ST distribution suggests a constraint sigma_8 Omega_0^{0.26} \simeq 0.76. The ST relation is expected to get tighter at high redshifts. In the process, we derive a simple formula for the halo formation-redshift distribution for non-Gaussian models. We also suggest that the discrepancy between the naive zero-redshift ST relation and that observed may be due, at least in part, to the fact that lower-mass clusters form over a wider range of redshifts. An Appendix derives an equation for the formation-redshift distribution of halos.

We describe a detailed analytic model for predicting statistical quantities (such as number counts, redshift distributions and sizes) of clusters detected in blank-field, thermal Sunyaev-Zel'dovich effect experiments. We include in this model the possibility of non-Gaussian density perturbations in the early Universe and also describe a simple model for the effects of preheating on cluster Sunyaev-Zel'dovich effect fluxes. We use this model to explore the current state of the theoretical uncertainties present in this type of analytic modelling, highlighting where further improvement will be necessary to fully exploit forthcoming surveys. We then go on to explore the constraints on cosmological parameters, the presence of any non-Gaussianity and the degree of cluster preheating which may be obtained from both the BOLOCAM and PLANCK experiments. We find that, providing redshifts can be measured for all detected clusters, the BOLOCAM experiment may provide detections of non-Gaussianity or preheating and could give approximate measurements of these effects if prior knowledge of the various cosmological parameters is taken into account. The PLANCK experiment Sunyaev-Zel'dovich effect cluster survey is predicted to provide highly accurate (~5%) measurements of the degree of non-Gaussianity and preheating while also providing measurements of several cosmological parameters to accuracies of a few percent independent from those constraints that will be derived from its detections of primordial cosmic microwave background anisotropies.

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 halos. We include in the model the negative feedback on galaxy formation from the photoionizing background. Photoionization inhibits galaxy formation in low-mass dark matter halos in two ways: (i) heating of the IGM and inhibition of the collapse of gas into dark halos by the IGM pressure, and (ii) reduction in the radiative cooling of gas within halos. The result of our method is a self-consistent model of galaxy formation and the IGM. The IGM is reheated twice (during reionization of HI and HeII), 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.

We use a self-consistent model of galaxy formation and the evolution of the intergalactic medium to study the effects of the reionization of the universe at high redshift on the properties of satellite galaxies like those seen around the Milky Way. Photoionization suppresses the formation of small galaxies, so that surviving satellites are preferentially those that formed before the universe reionized. As a result, the number of satellites expected today is about an order of magnitude smaller than the number inferred by identifying satellites with subhalos in high-resolution simulations of the dark matter. The resulting satellite population has an abundance and a distribution of circular velocities similar to those observed in the Local Group. We explore many other properties of satellite galaxies, including their gas content, metallicity and star formation rate, and find generally good agreement with available data. Our model predicts the existence of many as yet undetected satellites in the Local Group. We quantify their observability in terms of their apparent magnitude and surface brightness and also in terms of their constituent stars. A near-complete census of the Milky Way's satellites would require imaging to V~20 and to a surface brightness fainter than 26 V-band magnitudes per square arcsecond. Satellites with integrated luminosity V=15 should contain of order 100 stars brighter than B=26, with central stellar densities of a few tens per square arcminute. Discovery of a large population of faint satellites would provide a strong test of current models of galaxy formation.

We study the signal in the Cosmic Microwave Background (CMB) polarization anisotropy resulting from patchy reionization. It is well known that the primordial polarization of the CMB is very sensitive to the details of reionization. Combining a semi-analytic model of galaxy formation, in which the optical depth to the reionization epoch is in the range 0.014 to 0.048, with a high resolution N-body simulation we find that reionization generates a peak with amplitude 0.05~0.15 \mu K at large angular scales. The position of this peak reveals the size of the horizon at reionization, whilst its amplitude is a measure of the optical depth to reionization. On small scales, ionized patches prior to full reionization create a second order polarization signal due to the coupling of the free electron density fluctuation with the quadrupole moment of the temperature anisotropy. Careful study reveals that the coupling generates the same power spectra for electric and magnetic modes, whose amplitudes of polarization anisotropies from this process are predicted to be ~10nK. The amplitude depends strongly on the total baryon density $\Omega_b$ and on the spatial correlations of the free electron density in the ionized regions, and weakly on the fraction of ionizing photons able to escape their source galaxy, $f_{esc}$. The first- and second-order signals are therefore sensitive to the details of how the reionization occurred. Detection of these signals will place important constraints on the reionization history of the Universe.

There is now compelling evidence in favour of the hierarchical structure formation paradigm. Semi-analytic modelling is a powerful tool which allows the formation and evolution of galaxies to be followed in a hierarchical framework. We review some of the latest developments in this area before discussing how such models can help us to interpret observations of the high redshift Universe.

We follow the evolution of the galaxy population in a Lambda-CDM cosmology by means of high-resolution N-body simulations in which the formation of galaxies and their observable properties are calculated using a semi-analytic model. We display images of the spatial distribution of galaxies in the simulations that illustrate its evolution and provide a qualitative understanding of the processes responsible for various biases that develop. We consider three specific statistical measures of clustering at z=1 and z=0: the correlation length (in real- and redshift-space) of galaxies of different luminosity, the morphology-density relation and the genus curve of the topology of galaxy isodensity surfaces. For galaxies with luminosity below L*, the z=0 correlation length depends very little on the luminosity of the sample, but for brighter galaxies it increases very rapidly, reaching values in excess of 10Mpc/h. The "accelerated" dynamical evolution experienced by galaxies in rich clusters, which is partly responsible for this effect, also results in a strong morphology-density relation. Remarkably, this relation is already well-established at z=1. The genus curves of the galaxies are significantly different from the genus curves of the dark matter but this is not due to genuine topological differences but rather to the sparse sampling of the density field provided by galaxies. The predictions of our model at z=0 will be tested by forthcoming data from the 2dF and Sloan galaxy surveys, and those at z=1 by the DEEP and VIRMOS surveys.

We present a new model for the X-ray properties of the intracluster medium that explicitly includes heating of the gas by the energy released during the evolution of cluster galaxies. We calculate the evolution of clusters by combining the semi-analytic model of galaxy formation of Cole et al. with a simple model for the radial profile of the intracluster gas. We focus on the cluster X-ray luminosity function and on the relation between X-ray temperature and luminosity. We show that cooling alone is not enough to account for the flatness of the observed T-L relation or for the lack of strong redshift evolution in the observed X-ray luminosity function. Gas heating, on the other hand, can solve these two problems: in the \Lambda-CDM cosmology, our model reproduces fairly well the T-L relation and the X-ray luminosity function. Furthermore, it predicts only weak evolution in these two properties out to z=0.5, in agreement with recent observational data. A successful model requires an energy input of 1--2 x 10^49 ergs per solar mass of stars formed. This is comparable to the total energy released by the supernovae associated with the formation of the cluster galaxies. Thus, unless the transfer of supernovae energy to the intracluster gas is very (perhaps unrealistically) efficient, additional sources of energy, such as mechanical energy from AGN winds are required. However, the amplification of an initial energy input by the response of the intracluster medium to protocluster mergers might ease the energy requirements. Our model makes definite predictions for the X-ray properties of groups and clusters at high redshift which may soon be tested with data from the Chandra and Newton satellites.

We develop a method to measure the probability, P(N;M), of finding N galaxies in a dark-matter halo of mass M from the theoretically determined clustering properties of dark-matter halos and the observationally measured clustering properties of galaxies. Knowledge of this function and the distribution of the dark matter completely specifies all clustering properties of galaxies on scales larger than the size of dark-matter halos. Furthermore, P(N;M) provides strong constraints on models of galaxy formation, since it depends upon the merger history of dark-matter halos and the galaxy-galaxy merger rate within halos. We show that measurements from a combination of the 2MASS and SDSS or 2dFGRS datasets will allow P(N;M) averaged over halos occupied by bright galaxies to be accurately measured for N=0-2.

We compare the statistical properties of galaxies found in two different models of hierarchical galaxy formation: the semi-analytic model of Cole et al. and the smoothed particle hydrodynamics (SPH) simulations of Pearce et al. Using a `stripped-down' version of the semi-analytic model which mimics the resolution of the SPH simulations and excludes physical processes not included in them, we find that the two models produce an ensemble of galaxies with remarkably similar properties, although there are some differences in the gas cooling rates and in the number of galaxies that populate halos of different mass. The full semi-analytic model, which has effectively no resolution limit and includes a treatment of star formation and supernovae feedback, produces somewhat different (but readily understandable) results. Agreement is particularly good for the present-day global fractions of hot gas, cold dense (i.e. galactic) gas and uncollapsed gas, for which the SPH and stripped-down semi-analytic calculations differ by at most 25%. In the most massive halos, the stripped-down semi-analytic model predicts, on the whole, up to 50% less gas in galaxies than is seen in the SPH simulations. The two techniques apportion this cold gas somewhat differently amongst galaxies in a given halo. This difference can be tracked down to the greater cooling rate in massive halos in the SPH simulation compared to the semi-analytic model. (abridged)

We present predictions for the reionization of the intergalactic medium by stars in high-redshift galaxies. We calculate ionizing luminosities of galaxies, including the effects of absorption by interstellar gas and dust on the escape fraction, in order to calculate the filling factor of ionized hydrogen in the IGM. For a $\Lambda{\rm CDM}$ cosmology and a physical calculation of the escape fraction, we find that the hydrogen in the IGM will be reionized at redshift z=6.1 if the IGM has uniform density, but only by z=4.5 if the IGM is clumped. We combine our model with an N-body simulation of the distribution of dark matter in the universe in order to calculate the secondary temperature anisotropies induced in the cosmic microwave background (CMB) by scattering off free electrons. The models predict a spectrum of secondary anisotropies covering a broad range of angular scales, with fractional temperature fluctuations $\sim 10^{-7}-10^{-6}$ on arcminute scales. The amplitude depends somewhat on the geometry of reionization, with models in which the regions of highest gas density are reionized first giving larger CMB fluctuations than the case where galaxies ionize surrounding spherical regions, and models where low density regions reionize first giving the smallest fluctuations. Measurement of these anisotropies can therefore put important constraints on the reionization process, in particular, the redshift evolution of the filling factor, and should be a primary objective of a next generation submillimeter telescope such as the Atacama Large Millimeter Array. (abrdiged)

We use a combination of N-body simulations of the hierarchical clustering of dark matter and semi-analytic modelling of the physics of galaxy formation to probe the relationship between the galaxy distribution and the mass distribution. The pairwise line-of-sight velocity dispersion of galaxies is sensitive to the number of galaxies present in halos of different mass. The galaxy velocity dispersion is ~40% lower than that of the dark matter and is in reasonable agreement with the values measured in the Las Campanas redshift survey by Jing et al. over two decades in pair separation. The origin of this offset is statistical rather than dynamical, and depends upon the relative efficiency of galaxy formation in dark matter halos of different mass. Although the model galaxies and the dark matter have markedly different correlation functions in real space, such biases conspire to cause the redshift space correlation functions to be remarkably similar to each other. The predicted redshift-space galaxy correlation function agrees well with observations. We find no evidence in the model for a dependence of clustering strength on intrinsic galaxy luminosity, unless extremely bright galaxies, two magnitudes brighter than L*, are considered. However, there are significant differences when model galaxies are selected either by morphology or by colour. Early type or red galaxies show a much stronger clustering amplitude than late type or blue galaxies, particularly on small scales, again in good agreement with observations.

Current theories of galaxy formation predict that spiral galaxies are embedded in a reservoir of hot gas. This gas is able to cool onto the galaxy replenishing cold gas that is consumed by star formation. Estimates of the X-ray luminosity emitted in the cooling region suggest a bolometric luminosity of order 10 x 10^41 ergs/s in massive systems. We have used ROSAT PSPC data to search for extended X-ray emission from the haloes of three nearby, massive, late-type galaxies: NGC 2841, NGC 4594 and NGC 5529. We infer 95 per cent upper limits on the bolometric X-ray luminosities of the haloes of NGC 2841, NGC 4594 and NGC 5529 of 0.4, 1.2 and 3.8 x 10^41 ergs/s respectively. Thus the true luminosity lies well below the straightforward theoretical prediction. We discuss this discrepancy and suggest a number of ways in which the theoretical model can be brought into agreement with the observational results. A possible solution is that the gravitational potentials of the dark matter haloes of these galaxies are weaker than assumed in the current model. Alternatively, the present day accretion may be substantially less than is required on average to build the disk over the Hubble time. Our results are, however, based on only three galaxies, none of which is ideal for this kind of study. A larger data set is required to explore this important problem further.

We have used a combination of high resolution cosmological N-body simulations and semi-analytic modelling of galaxy formation to investigate the processes that determine the spatial distribution of galaxies in cold dark matter (CDM) models. The galaxy distribution depends sensitively on the efficiency with which galaxies form in halos of different mass. In small mass halos, galaxy formation is inhibited by the reheating of cooled gas by feedback processes, whereas in large mass halos, it is inhibited by the long cooling time of the gas. As a result, the mass-to-light ratio of halos has a deep minimum at the halo mass associated with L* galaxies. This leads to a scale dependent bias in the distribution of galaxies relative to the distribution of mass. On large scales, the bias in the galaxy distribution is related in a simple way to the bias in the distribution of massive halos. On small scales, the correlation function is determined by the interplay between various effects including the spatial exclusion of dark matter halos, the distribution function of the number of galaxies occupying a single dark matter halo and, to a lesser extent, dynamical friction. Remarkably, these processes conspire to produce a correlation function in a flat, Omega_0=0.3, CDM model that is close to a power-law over nearly four orders of magnitude in amplitude. This model agrees well with the correlation function of galaxies measured in the APM survey. On small scales, the model galaxies are less strongly clustered than the dark matter whereas on large scales, they trace the occupied halos. Our clustering predictions are robust to changes in the parameters of the galaxy formation model, provided only those models that match the bright end of the galaxy luminosity function are considered. (abridged)

The main ingredients of recent semi-analytic models of galaxy formation are summarised. We present predictions for the galaxy clustering properties of a well specified LCDM model whose parameters are constrained by observed local galaxy properties. We present preliminary predictions for evolution of clustering that can be probed with deep pencil beam surveys.

The formation and evolution of galaxies can be followed in the context of cosmological structure formation using the technique of semi-analytic modelling. We give a brief outline of the features incorporated into the semi-analytic model of Cole etal (1999). We present two examples of model predictions that can be tested using photometric redshift techniques. The first prediction, of the star formation history of the universe, has already been shown to be in broad argeement with the observational estimates. The second prediction, of the evolution of galaxy clustering with redshift, will be addressed with some of the forthcoming deep, multi-filter imaging surveys discussed at this meeting.

Measurements of galaxy clustering are now becoming possible over a range of redshifts out to z=3. We use a semi-analytic model of galaxy formation to compute the expected evolution of the galaxy correlation function with redshift. We illustrate how the degree of clustering evolution is sensitive to the details of the sample selection. For a fixed apparent magnitude limit, galaxies selected at higher redshifts are located in progressively rarer dark matter haloes, compared with the general population of galaxies in place at each redshift. As a result these galaxies are highly biased tracers of the underlying dark matter distribution and exhibit stronger clustering than the dark matter. In general, the correlation length measured in comoving units, decreases at first with increasing redshift, before increasing again at higher redshift. We show that the epsilon-model often used to interpret the angular correlation function of faint galaxies gives an inadequate description of the evolution of clustering, and offers no physical insight into the clustering process. We compare our predictions with those of a simple, popular model in which a one-to-one correspondence between galaxies and dark halos is assumed. Qualitatively, this model reproduces the correct evolutionary behaviour at high redshift, but the quantitative results can be significantly in error. Our theoretical expectations are in good agreement with the high redshift clustering data of Carlberg etal and Postman etal but are higher than the measurements of Le Fevre etal.

We model the two-point correlation function of galaxies in a cold dark matter Universe by combining two powerful theoretical tools - dissipationless N-body simulations of dark matter clustering (specifically the GIF simulations carried out by MPIA and the Virgo Consortium) and semi-analytic modelling of galaxy formation. We construct catalogues of galaxies containing a wide range of information for each galaxy, including magnitudes in various bands, star formation rates, disk and bulge sizes (from the semi-analytic model described by Cole et al 1998) and, importantly, spatial positions and peculiar velocities from the N-body simulation. We then use this information to study the clustering properties of galaxies.

Any realistic theory of galaxy formation must be set in the context of a model for the formation of structure in the universe. We describe a powerful approach -- semi-analytic modelling -- that combines a set of simple rules describing the gas processes involved in galaxy formation with a scheme to follow the hierarchical growth of dark matter haloes. Surprisingly few free parameters are required to specify the model, and these are fixed with reference to a subset of local observational data. The model produces the full star formation history of a galaxy, allowing a wide range of predictions to be made. We review some of the successes of the models, namely the star formation history of the universe and the evolution of galaxy clustering, before focusing our attention on early-type galaxies. We discuss the observational evidence against the classical picture in which early-type galaxies form at some arbitrarily high redshift in a single monolithic collapse and burst of star formation. The alternative scenario in which spheroidal systems are formed by the merger of disk galaxies is outlined. We review some of the predictions of this model, namely the colour-magnitude relation, the faint counts and the evolution of cluster membership.