We show that simple scalar field models can give rise to curvature singularities in the effective Friedmann dynamics of Loop Quantum Cosmology (LQC). We find singular solutions for spatially flat Friedmann-Robertson-Walker cosmologies with a canonical scalar field and a negative exponential potential, or with a phantom scalar field and a positive potential. While LQC avoids big bang or big rip type singularities, we find sudden singularities where the Hubble rate is bounded, but the Ricci curvature scalar diverges. We conclude that the effective equations of LQC are not in themselves sufficient to avoid the occurrence of singularities.

We highlight the unexpected impact of nucleosynthesis and other early universe constraints on the detectability of tracking quintessence dynamics at late times, showing that such dynamics may well be invisible until the unveiling of the Stage-IV dark energy experiments (DUNE, JDEM, LSST, SKA). Nucleosynthesis forces |w'(0)| < 0.2 for the models we consider and strongly limits potential deviations from LCDM. Surprisingly, the standard CPL parametrisation, w(z) = w_0 + w_a z/(1+z), cannot match the nucleosynthesis bound for minimally coupled tracking scalar fields. Given that such models are arguably the best-motivated alternatives to a cosmological constant these results may significantly impact future cosmological survey design and imply that dark energy may well be dynamical even if we do not detect any dynamics in the next decade.

We consider linear perturbation equations for long-wavelength scalar metric perturbations in generalised gravity, applicable to non-singular cosmological models including a bounce from collapse to expansion in the very early universe. We present the general form for the perturbation equations which follows from requiring that the inhomogeneous universe on large scales obeys the same local equations as the homogeneous Friedmann-Robertson-Walker background cosmology (the separate universes approach). In a pseudo-longitudinal gauge this becomes a homogeneous second-order differential equation for adiabatic perturbations, which reduces to the usual equation for the longitudinal gauge metric perturbation in general relativity with vanishing anisotropic stress. As an application we show that the scale-invariant spectrum of perturbations in the longitudinal gauge generated during an ekpyrotic collapse are not transferred to the growing mode of adiabatic density perturbations in the expanding phase in a simple bounce model.

We solve for the behaviour of cosmological perturbations in the Dvali-Gabadadze-Porrati (DGP) braneworld model using a new numerical method. Unlike some other approaches in the literature, our method uses no approximations other than linear theory and is valid on large scales. We examine the behaviour of late-universe density perturbations for both the self-accelerating and normal branches of DGP cosmology.

We solve for the behaviour of cosmological perturbations in the Dvali-Gabadadze-Porrati (DGP) braneworld model using a new numerical method. Unlike some other approaches in the literature, our method uses no approximations other than linear theory and is valid on large scales. We examine the behaviour of late-universe density perturbations for both the self-accelerating and normal branches of DGP cosmology. Our numerical results can form the basis of a detailed comparison between the DGP model and cosmological observations.

We study the evolution of scalar perturbations in the radiation-dominated era of Randall-Sundrum braneworld cosmology by numerically solving the coupled bulk and brane master wave equations. We find that density perturbations with wavelengths less than a critical value (set by the bulk curvature length) are amplified during horizon re-entry. Conversely, we explicitly confirm from simulations that the spectrum is identical to GR on large scales. Although this magnification is not relevant for the cosmic microwave background or measurements of large scale structure, it may have some bearing on the formation of primordial black holes in Randall-Sundrum models.

We study the behaviour of scalar perturbations in the radiation-dominated era of Randall-Sundrum braneworld cosmology by numerically solving the coupled bulk and brane master wave equations. We find that density perturbations with wavelengths less than a critical value (set by the bulk curvature length) are amplified during horizon re-entry. This means that the radiation era matter power spectrum will be at least an order of magnitude larger than the predictions of general relativity (GR) on small scales. Conversely, we explicitly confirm from simulations that the spectrum is identical to GR on large scales. Although this magnification is not relevant for the cosmic microwave background or measurements of large scale structure, it will have some bearing on the formation of primordial black holes in Randall-Sundrum models.

We reconsider the original model of quadratic hybrid inflation in light of the WMAP three-year results and study the possibility of obtaining a spectral index of primordial density perturbations, $n_s$, smaller than one from this model. The original hybrid inflation model naturally predicts $n_s\geq1$ in the false vacuum dominated regime but it is also possible to have $n_s<1$ when the quadratic term dominates. We therefore investigate whether there is also an intermediate regime compatible with the latest constraints, where the scalar field value during the last 50 e-folds of inflation is less than the Planck scale.