The analytic "equilibrium model" for galaxy evolution using a mass balance equation is able to reproduce mean observed galaxy scaling relations between stellar mass, halo mass, star formation rate (SFR) and metallicity across the majority of cosmic time with a small number of parameters related to feedback. Here we aim to test this data-constrained model to quantify deviations from the mean relation between stellar mass and SFR, i.e. the star-forming galaxy main sequence (MS). We implement fluctuation in halo accretion rates parameterised from merger-based simulations, and quantify the intrinsic scatter introduced into the MS under the assumption that fluctuations in star formation follow baryonic inflow fluctuations. We predict the 1-sigma MS scatter to be ~ 0.2 - 0.25 dex over the stellar mass range 10^8 Mo to 10^11 Mo and a redshift range 0.5 < z < 3 for SFRs averaged over 100 Myr. The scatter increases modestly at z > 3, as well as by averaging over shorter timescales. The contribution from merger-induced star formation is generally small, around 5% today and 10 - 15% during the peak epoch of cosmic star formation. These results are generally consistent with available observations, suggesting that deviations from the MS primarily reflect stochasticity in the inflow rate owing to halo mergers.

We present a model for the dynamical evolution of subhaloes based on an approach combining numerical and analytical methods. Our method is based on tracking subhaloes in an N-body simulation up to the last point that it can be resolved, and applying an analytic prescription for its merger timescale that takes dynamical friction and tidal disruption into account. When applied to cosmological N-body simulations with mass resolutions that differ by two orders of magnitude, the technique produces halo occupation distributions that agree to within 3%.

We present predictions for the two-point correlation function of galaxy clustering as a function of stellar mass, computed using two new versions of the GALFORM semi-analytic galaxy formation model. These models make use of a high resolution, large volume N-body simulation, set in the WMAP7 cosmology. One model uses a universal stellar initial mass function (IMF), while the other assumes different IMFs for quiescent star formation and bursts. Particular consideration is given to how the assumptions required to estimate the stellar masses of observed galaxies (such as the choice of IMF, stellar population synthesis model and dust extinction) influence the perceived dependence of galaxy clustering on stellar mass. Broad-band spectral energy distribution fitting is carried out to estimate stellar masses for the model galaxies in the same manner as in observational studies. We show clear differences between the clustering signals computed using the true and estimated model stellar masses. As such, we highlight the importance of applying our methodology to compare theoretical models to observations. We introduce an alternative scheme for the calculation of the merger timescales for satellite galaxies in GALFORM, which takes into account the dark matter subhalo information from the simulation. This reduces the amplitude of small-scale clustering. The new merger scheme offers improved or similar agreement with observational clustering measurements, over the redshift range 0 < z < 0.7. We find reasonable agreement with clustering measurements from GAMA, but find larger discrepancies for some stellar mass ranges and separation scales with respect to measurements from SDSS and VIPERS, depending on the GALFORM model used.

We examine the star formation histories (SFHs) of galaxies in smoothed particle hydrodynamics (SPH) simulations, compare them to parametric models that are commonly used in fitting observed galaxy spectral energy distributions, and examine the efficacy of these parametric models as practical tools for recovering the physical parameters of galaxies. The commonly used tau-model, with SFR ~ exp(-t/tau), provides a poor match to the SFH of our SPH galaxies, with a mismatch between early and late star formation that leads to systematic errors in predicting colours and stellar mass-to-light ratios. A one-parameter lin-exp model, with SFR ~ t*exp(-t/tau), is much more successful on average, but it fails to match the late-time behavior of the bluest, most actively star-forming galaxies and the passive, "red and dead" galaxies. We introduce a 4-parameter model, which transitions from lin-exp to a linear ramp after a transition time, which describes our simulated galaxies very well. We test the ability of these parametrised models to recover (at z=0, 0.5, and 1) the stellar mass-to-light ratios, specific star formation rates, and stellar population ages from the galaxy colours, computed from the full SPH star formation histories using the FSPS code of Conroy et al. (2009). Fits with tau-models systematically overestimate M/L by ~ 0.2 dex, overestimate population ages by ~ 1-2 Gyr, and underestimate sSFR by ~ 0.05 dex. Fits with lin-exp are less biased on average, but the 4-parameter model yields the best results for the full range of galaxies. Marginalizing over the free parameters of the 4-parameter model leads to slightly larger statistical errors than 1-parameter fits but essentially removes all systematic biases, so this is our recommended procedure for fitting real galaxies.

We place constraints on the matter density of the Universe and the amplitude of clustering using measurements of the galaxy two-point correlation function from the Sloan Digital Sky Survey (SDSS). We generate model predictions for different cosmologies by populating rescaled N-body simulations with galaxies using the subhalo abundance matching (SHAM) technique. We find Omega-M = 0.29 +/- 0.03 and sigma-8 = 0.86 +/- 0.04 at 68% confidence from fitting the observed two-point galaxy correlation function of galaxies brighter than M_r = -18 in a volume limited sample of galaxies obtained by the SDSS. We discuss and quantify potential sources of systematic error, and conclude that while there is scope for improving its robustness, the technique presented in this paper provides a powerful low redshift constraint on the cosmological parameters that is complementary to other commonly used methods.

Subhalo abundance matching (SHAM) is a technique for populating simulated dark matter distributions with galaxies, assuming a monotonic relation between a galaxy's stellar mass or luminosity and the mass of its parent dark matter halo or subhalo. We examine the accuracy of SHAM in two cosmological SPH simulations, one of which includes momentum-driven winds. The SPH simulations indeed show a nearly monotonic relation between stellar mass and halo mass provided that, for satellite galaxies, we use the mass of the subhalo at the epoch when it became a satellite. In each simulation, the median relation for central and satellite galaxies is nearly identical, though a somewhat larger fraction of satellites are outliers. SHAM-assigned masses (at z=0-2), luminosities (R-band at z=0), or star formation rates (at z=2) have a 68% scatter of 0.09-0.15 dex relative to the true simulation values. When we apply SHAM to the subhalo population of collisionless N-body simulation with the same initial conditions as the SPH runs, we find generally good agreement for the halo occupation distributions and halo radial profiles of galaxy samples defined by thresholds in stellar mass. However, because a small fraction of SPH galaxies suffer severe stellar mass loss after becoming satellites, SHAM slightly overpopulates high mass halos; this effect is more significant for the wind simulation, which produces galaxies that are less massive and more fragile. SHAM recovers the two-point correlation function of the SPH galaxies in the no-wind simulation to better than 10% at scales 0.1 < r < 10 Mpc/h. For the wind simulation, agreement is better than 15% at r > 2 Mpc/h, but overpopulation of massive halos increases the correlation function by a factor of ~2.5 on small scales.

We examine the accretion and merger histories of central and satellite galaxies in a smoothed particle hydrodynamics (SPH) cosmological simulation that resolves galaxies down to 7e9 M[Sun]. Most friends-of-friends halos in the simulation have a distinct central galaxy, typically two to five times more massive than the most massive satellite. As expected, satellites have systematically higher assembly redshifts than central galaxies of the same baryonic mass, and satellites in more massive halos form earlier. However, contrary to the simplest expectations, satellite galaxies continue to accrete gas and convert it to stars; the gas accretion declines steadily over a period of 0.5-1 Gyr after the satellite halo merges with a larger parent halo. Satellites in a cluster mass halo eventually begin to lose baryonic mass. Since z=1, 27% of central galaxies (above 3e10 M[Sun]) and 22% of present-day satellite galaxies have merged with a smaller system above a 1:4 mass ratio; about half of the satellite mergers occurred after the galaxy became a satellite and half before. In effect, satellite galaxies can remain "central" objects of halo substructures, with continuing accretion and mergers, making the transition in assembly histories and physical properties a gradual one. Implementing such a gradual transformation in semi-analytic models would improve their agreement with the observed colour distributions of satellite galaxies in groups and with the observed colour dependence of galaxy clustering.

We explore constraints on extensions to the standard models of cosmology and particle physics which modify the early-Universe expansion rate S = H'/H (parametrized by the effective number of neutrinos N_nu). The constraints on N_nu and the baryon density parameter (eta_B = n_B/n_gamma = 10^(-10)*eta_10) from BBN at 20 minutes are compared with those from the CMB at 400 kyr and LSS at 14 Gyr. BBN provides the strongest constraint on N_nu (1.6 < N_nu < 3.3 at 95% confidence), but a weaker constraint on eta_B. The CMB/LSS best constrain the baryon density (5.9 < eta_10 < 6.4 at 95% confidence), independent of N_nu, but provide a relatively weak N_nu constraint, consistent with N_nu = 3. Using the best fit values and the allowed ranges of the CMB/LSS-derived parameters to calculate the BBN-predicted primordial abundances yields excellent agreement with the observationally inferred abundance of deuterium and good agreement with 4He, confirming the consistency between the BBN and CMB/LSS results. However, the BBN-predicted abundance of 7Li is high, by a factor of 3 or more. We comment on the value of N_nu and a possible anomaly in the matter power spectrum inferred from observations of the Ly-alpha forest. The good agreement between our BBN and CMB/LSS results permit us to constrain any post-BBN entropy production as well as to limit the production of any non-thermalized relativistic particles and, allow us to combine them finding 95% ranges, 1.8 < N_nu < 3.2 and 5.9 < eta_10 < 6.4.

The relic cosmic background neutrinos accompanying the cosmic microwave background (CMB) photons may hide a universal lepton asymmetry orders of magnitude larger than the universal baryon asymmetry. At present, the only direct way to probe such an asymmetry is through its effect on the abundances of the light elements produced during primordial nucleosynthesis. The relic light element abundances also depend on the baryon asymmetry, parameterized by the baryon density parameter (eta_B = n_B/n_gamma = 10^(-10)*eta_10), and on the early-universe expansion rate, parameterized by the expansion rate factor (S = H'/H) or, equivalently by the effective number of neutrinos (N_nu = 3 + 43(S^2 - 1)/7). We use data from the CMB (and Large Scale Structure: LSS) along with the observationally-inferred relic abundances of deuterium and helium-4 to provide new bounds on the universal lepton asymmetry, finding for eta_L, the analog of eta_B, 0.072 +/- 0.053 if it is assumed that N_nu = 3 and, 0.115 +/- 0.095 along with N_nu = 3.3^{+0.7}_{-0.6}, if N_nu is free to vary.