The largest galaxies acquire their mass early on, when the Universe is still youthful. Cold streams violently feed these young galaxies a vast amount of fresh gas, resulting in very efficient star formation. Using a well resolved hydrodynamical simulation of galaxy formation, we demonstrate that these mammoth galaxies are already in place a couple of billion years after the Big Bang. Contrary to local starforming galaxies, where dust re-emits a large part of the stellar ultraviolet (UV) light at infrared and sub-millimetre wavelengths, our self-consistent modelling of dust extinction predicts that a substantial fraction of UV photons should escape from primordial galaxies. Such a model allows us to compute reliably the number of high redshift objects as a function of luminosity, and yields galaxies whose UV luminosities closely match those measured in the deepest observational surveys available. This agreement is remarkably good considering our admittedly still simple modelling of the interstellar medium (ISM) physics. The luminosity functions (LF) of virtual UV luminous galaxies coincide with the existing data over the whole redshift range from 4 to 7, provided cosmological parameters are set to their currently favoured values. Despite their considerable emission at short wavelengths, we anticipate that the counterparts of the brightest UV galaxies will be detected by future sub-millimetre facilities like ALMA

We show how to estimate the covariance of the power spectrum of a statistically homogeneous and isotropic density field from a single periodic simulation, by applying a set of weightings to the density field, and by measuring the scatter in power spectra between different weightings. We recommend a specific set of 52 weightings containing only combinations of fundamental modes, constructed to yield a minimum variance estimate of the covariance of power. Numerical tests reveal that at nonlinear scales the variance of power estimated by the weightings method substantially exceeds that estimated from a simple ensemble method. We argue that the discrepancy is caused by beat-coupling, in which products of closely spaced Fourier modes couple by nonlinear gravitational growth to the beat mode between them. Beat-coupling appears whenever nonlinear power is measured from Fourier modes with a finite spread of wavevector, and is therefore present in the weightings method but not the ensemble method. Beat-coupling inevitably affects real galaxy surveys, whose Fourier modes have finite width. Surprisingly, the beat-coupling contribution dominates the covariance of power at nonlinear scales, so that, counter-intuitively, it is expected that the covariance of nonlinear power in galaxy surveys is dominated not by small scale structure, but rather by beat-coupling to the largest scales of the survey.

We measure the covariance of the non-linear matter power spectrum from N-body simulations using two methods. In the first case, the covariance of power is estimated from the scatter over many random realizations of the density field. In the second, we use a novel technique to measure the covariance matrix from each simulation individually by re-weighting the density field with a carefully chosen set of functions. The two methods agree at linear scales, but unexpectedly they disagree substantially at increasingly non-linear scales. Moreover, the covariance of non-linear power measured using the re-weightings method changes with box size. The numerical results are consistent with an explanation given in a companion paper, which argues that the cause of the discrepancy is beat-coupling, in which products of Fourier modes separated by a small wavevector couple by gravitational growth to the large-scale beat mode between them. We calculate the information content of the non-linear power spectrum (about the amplitude of the initial, linear power spectrum) using both methods and confirm the result of a previous paper, that at translinear scales the power spectrum contains little information over and above that in the linear power spectrum, but that there is a marked increase in information at non-linear scales. We suggest that, in real galaxy surveys, the covariance of power at non-linear scales is likely to be dominated by beat-coupling to the largest scales of the survey and that, as a result, only part of the information potentially available at non-linear scales is actually measurable from real galaxy surveys.

Using the halo model, we investigate the cosmological Fisher information in the non-linear dark-matter power spectrum about the initial amplitude of linear power. We find that there is little information on `translinear' scales (where the one- and two-halo terms are both significant) beyond what is on linear scales, but that additional information is present on small scales, where the one-halo term dominates. This behavior agrees with the surprising results that Rimes & Hamilton (2005, 2006) found using N-body simulations. We argue that the translinear plateau in cumulative information arises largely from fluctuations in the numbers of large haloes in a finite volume. This implies that more information could be extracted on non-linear scales if the masses of the largest haloes in a survey are known.

We use an ensemble of N-body simulations of the currently favoured (concordance) cosmological model to measure the amount of information contained in the non-linear matter power spectrum about the amplitude of the initial power spectrum. Two surprising results emerge from this study: (i) that there is very little independent information in the power spectrum in the translinear regime (k ~ 0.2-0.8 Mpc/h at the present day) over and above the information at linear scales and (ii) that the cumulative information begins to rise sharply again with increasing wavenumber in the non-linear regime. In the fully non-linear regime, the simulations are consistent with no loss of information during translinear and non-linear evolution. If this is indeed the case then the results suggest a picture in which translinear collapse is very rapid, and is followed by a bounce prior to virialization, impelling a wholesale revision of the HKLM-PD formalism.

We present a phenomenological model for understanding the origins and evolution of galaxy morphologies from a phenomenological perspective. The model includes an observationally motivated prescription for star formation in galaxy disks, as well as a merger-driven starburst mode of star formation. We consider the formation of bulges and ellipticals both in mergers and by global instabilities in disks. We use our model to investigate the fundamental properties of disk galaxies and the effects of surface-brightness limits on large galaxy surveys.