Studies of the red clump giant population in the inner Milky Way suggest the Galactic bulge/bar has a boxy/peanut/X-shaped structure as predicted by its formation via a disc buckling instability. We used a non-parametric method of estimating the Galactic bulge morphology that is based on maximum entropy regularisation. This enabled us to extract the three-dimensional distribution of the red giant stars in the bulge from deep photometric catalogues of the VISTA Variables in the Via Lactea (VVV) survey. Our high-resolution reconstruction confirms the well-known boxy/peanut/X-shaped structure of the bulge. We also find spiral arm structures that extend to around three kpc in front of and behind the bulge and are on different sides of the bulge major axis. However, we show that the details of the features behind and in front of the bulge are sensitive to the amplitudes of the red giant branch bump and asymptotic giant branch bump we used in our luminosity function.

The abundance and narrow magnitude dispersion of Red Clump (RC) stars make them a popular candidate for mapping the morphology of the bulge region of the Milky Way. Using an estimate of the RC's intrinsic luminosity function, we extracted the three-dimensional density distribution of the RC from deep photometric catalogues of the VISTA Variables in the Via Lactea (VVV) survey. We used maximum entropy based deconvolution to extract the spatial distribution of the bulge from Ks-band star counts. We obtained our extrapolated non-parametric model of the bulge over the inner 40 by 40 degrees squared region of the Galactic centre. Our reconstruction also naturally matches onto a parametric fit to the bulge outside the VVV region and inpaints overcrowded and high extinction regions. We found a range of bulge properties consistent with other recent investigations based on the VVV data. In particular, we estimated the bulge mass to be in the range 13 to 17 billion solar masses, the X-component to be between 18% and 25% of the bulge mass, and the bulge angle with respect to the Sun-Galactic centre line to be between 18 and 32 degrees. Studies of the Fermi Large Area Telescope (LAT) gamma-ray Galactic centre excess suggests that the excess may be traced by Galactic bulge distributed sources. We applied our deconvolved density in a template fitting analysis of this Fermi-LAT GeV excess and found an improvement in the fit compared to previous parametric based templates.

There is growing evidence that the Galactic Center Excess identified in the $\textit{Fermi}$-LAT gamma-ray data arises from a population of faint astrophysical sources. We provide compelling supporting evidence by showing that the morphology of the excess traces the stellar over-density of the Galactic bulge. By adopting a template of the bulge stars obtained from a triaxial 3D fit to the diffuse near-infrared emission, we show that it is detected at high significance. The significance deteriorates when either the position or the orientation of the template is artificially shifted, supporting the correlation of the gamma-ray data with the Galactic bulge. In deriving these results, we have used more sophisticated templates at low-latitudes for the $\textit{Fermi}$ bubbles compared to previous work and the three-dimensional Inverse Compton (IC) maps recently released by the ${\tt GALPROP}$ team. Our results provide strong constraints on Millisecond Pulsar (MSP) formation scenarios proposed to explain the excess. We find that an $\textit{admixture formation}$ scenario, in which some of the relevant binaries are $\textit{primordial}$ and the rest are formed $\textit{dynamically}$, is preferred over a primordial-only formation scenario at $7.6\sigma$ confidence level. Our detailed morphological analysis also disfavors models of the disrupted globular clusters scenario that predict a spherically symmetric distribution of MSPs in the Galactic bulge. For the first time, we report evidence of a high energy tail in the nuclear bulge spectrum that could be the result of IC emission from electrons and positrons injected by a population of MSPs and star formation activity from the same site.

An anomalous gamma-ray excess emission has been found in Fermi Large Area Telescope data covering the centre of the Galaxy. Several theories have been proposed for this `Galactic Centre Excess'. They include self-annihilation of dark matter particles, an unresolved population of millisecond pulsars, an unresolved population of young pulsars, or a series of burst events. Here we report on a new analysis that exploits hydrodynamical modelling to register the position of interstellar gas associated with diffuse Galactic gamma-ray emission. We find evidence that the Galactic Centre Excess gamma rays are statistically better described by the stellar over-density in the Galactic bulge and the nuclear stellar bulge, rather than a spherical excess. Given its non-spherical nature, we argue that the Galactic Centre Excess is not a dark matter phenomenon but rather associated with the stellar population of the Galactic bulge and nuclear bulge.

Fermi Large Area Telescope data reveal an excess of GeV gamma rays from the direction of the Galactic Center and bulge. Several explanations have been proposed for this excess including an unresolved population of millisecond pulsars (MSPs) and self-annihilating dark matter. It has been claimed that a key discriminant for or against the MSP explanation can be extracted from the properties of the luminosity function describing this source population. Specifically, is the luminosity function of the putative MSPs in the Galactic Center consistent with that characterizing the resolved MSPs in the Galactic disk? To investigate this we have used a Bayesian Markov Chain Monte Carlo to evaluate the posterior distribution of the parameters of the MSP luminosity function describing both resolved MSPs and the Galactic Center excess. At variance with some other claims, our analysis reveals that, within current uncertainties, both data sets can be well fit with the same luminosity function.

Excess GeV gamma rays from the Galactic Center (GC) have been measured with the Fermi Large Area Telescope (LAT). The presence of the GC excess (GCE) appears to be robust with respect to changes in the diffuse galactic background modeling. The three main proposals for the GCE are an unresolved population of millisecond pulsars (MSPs), outbursts of cosmic rays from the GC region, and self-annihilating dark matter (DM). The injection of secondary electrons and positrons into the interstellar medium (ISM) by an unresolved population of MSPs or DM annihilations can lead to observable gamma-ray emission via inverse Compton scattering or bremsstrahlung. Here we investigate how to determine whether secondaries are important in a model for the GCE. We develop a method of testing model fit which accounts for the different spatial morphologies of the secondary emission. We examine several models which give secondary emission and illustrate a case where a broadband analysis is not sufficient to determine the need for secondary emission.

The High-Energy Stereoscopic System (HESS) has detected diffuse TeV emission correlated with the distribution of molecular gas along the Ridge at the Galactic Center. Diffuse, non-thermal emission is also seen by the Fermi large area telescope (Fermi-LAT) in the GeV range and by radio telescopes in the GHz range. Additionally, there is a distinct, spherically symmetric excess of gamma rays seen by Fermi-LAT in the GeV range. A cosmic ray flare, occurring in the Galactic Center, $10^4$ years ago has been proposed to explain the TeV Ridge. An alternative, steady-state model explaining all three data sets (TeV, GeV, and radio) invokes purely leptonic processes. We show that the flare model from the Galactic Center also provides an acceptable fit to the GeV and radio data, provided the diffusion coefficient is energy independent. However, if Kolmogorov-type turbulence is assumed for the diffusion coefficient, we find that two flares are needed, one for the TeV data (occurring approximately $10^4 $ years ago) and an older one for the GeV data (approximately $10^5$ years old). We find that the flare models we investigate do not fit the spherically symmetric GeV excess as well as the usual generalized Navarro-Frenk-White spatial profile, but are better suited to explaining the Ridge. We also show that a range of single-zone, steady-state models are able to explain all three spectral data sets. Large gas densities equal to the volumetric average in the region can be accommodated by an energy independent diffusion or streaming based steady-state model. Additionally, we investigate how the flare and steady-state models may be distinguished with future gamma-ray data looking for a spatial dependence of the gamma-ray spectral index.

The Fermi Large Area Telescope data appear to have an excess of gamma rays from the inner 150 pc of the Galactic Center (GC). The main explanations proposed for this are: an unresolved population of millisecond pulsars (MSPs), dark matter (DM) annihilation, and cosmic rays interacting with molecular clouds at the GC. In this conference proceeding article we highlight some of the cosmic-ray results of [1]. The MSPs and DM explanations were modeled with spatial templates well fitted by the square of a generalized Navarro-Frenk-White (NFW) profile with inner slope $\gamma=1.2$. The cosmic-ray option was modeled with a 20-cm continuum emission Galactic Ridge template. A template based on the HESS residuals were shown to give similar results. The gamma-ray excess was found to be best fit by a combination of the generalized NFW squared template and a Galactic Ridge template. We also found the spectra of each template was not significantly affected in the combined fit and is consistent with previous single template fits. Due to the ridge like morphology of the molecular gas at the GC, the cosmic-ray option was not able to explain all the excess gamma rays, indicating that DM annihilation and/or an unresolved population of MSPs is also needed.

The Fermi-LAT data appear to have an excess of gamma rays from the inner 150 pc of the Galactic Center. The main explanations proposed for this are: an unresolved population of millisecond pulsars (MSPs), dark matter (DM) annihilation, and nonthermal bremsstrahlung produced by a population of electrons interacting with neutral gas in molecular clouds. The first two options have spatial templates well fitted by the square of a generalized Navarro-Frenk-White (NFW) profile with inner slope $\gamma=1.2$. We model the third option with a 20-cm continuum emission Galactic Ridge template. A template based on the HESS residuals is shown to give similar results. The gamma-ray excess is found to be best fit by a combination of the generalized NFW squared template and a Galactic Ridge template. We also find the spectra of each template is not significantly affected in the combined fit and is consistent with previous single template fits. That is the generalized NFW squared spectrum can be fit by either of order 1000 unresolved MSPs or DM with mass around 30 GeV, a thermal cross section, and mainly annihilating to $b \bar{b}$ quarks. While the Galactic Ridge continues to have a spectrum consistent with a population of nonthermal electrons whose spectrum also provides a good fit to synchrotron emission measurements. We also show that the current DM fit may be hard to test, even with 10 years of Fermi-LAT data, especially if there is a mixture of DM and MSPs contributing to the signal, in which case the implied DM cross section will be suppressed.

Employing Fermi-LAT gamma ray observations, several independent groups have found excess extended gamma ray emission at the Galactic center (GC). Both, annihilating dark matter (DM) or a population of $\sim 10^3$ unresolved millisecond pulsars (MSPs) are regarded as well motivated possible explanations. However, there is significant uncertainties in the diffuse galactic background at the GC. We have performed a revaluation of these two models for the extended gamma ray source at the GC by accounting for the systematic uncertainties of the Galactic diffuse emission model. We also marginalize over point source and diffuse background parameters in the region of interest. We show that the excess emission is significantly more extended than a point source. We find that the DM (or pulsars population) signal is larger than the systematic errors and therefore proceed to determine the sectors of parameter space that provide an acceptable fit to the data. We found that a population of order a 1000 MSPs with parameters consistent with the average spectral shape of Fermi-LAT measured MSPs was able to fit the GC excess emission. For DM, we found that a pure $\tau^+\tau^-$ annihilation channel is not a good fit to the data. But a mixture of $\tau^+\tau^-$ and $b\bar{b}$ with a $\left<\sigma v\right>$ of order the thermal relic value and a DM mass of around 20 to 60 GeV provides an adequate fit.

Current data are in good agreement with the predictions of single field inflation. However, the hemispherical asymmetry seen in the cosmic microwave background data, may hint at a potential problem. Generalizing to multi-field models may provide one possible explanation. A useful way of modeling perturbations in multi-field inflation is to investigate the projection of the perturbation along and perpendicular to the background fields' trajectory. These correspond to the adiabatic and isocurvature perturbations. However, it is important to note that in general there are no corresponding adiabatic and isocurvature fields. The purpose of this article is to highlight the distinction between a field redefinition and a perturbation projection. We provide a detailed derivation of the evolution of the isocurvature perturbation to show that no assumption of an adiabatic or isocurvature field is needed. We also show how this evolution equation is consistent with the field covariant evolution equations for the isocurvature perturbation in the flat field space limit.

Based on three years of Fermi Large Area Telescope (LAT) gamma-ray data of the Virgo cluster, evidence for an extended emission associated with dark matter pair annihilation in the b b-bar channel has been reported by Han et al. [1]. After an in depth spatial and temporal analysis, we argue that the tentative evidence for a gamma-ray excess from the Virgo cluster is mainly due to the appearance of a population of previously unresolved gamma-ray point sources in the region of interest. These point sources are not part of the LAT second source catalogue (2FGL), but are found to be above the standard detection significance threshold when three or more years of LAT data is included.

Number counts of massive high-redshift clusters provide a window to study primordial non-Gaussianity. The current quality of data, however, forces the statistical analysis to probe a region of parameter space -- the extreme tail of the mass function -- which is neither accessible in any of the currently available theoretical prescriptions for calculating the mass function, nor calibrated in N-body simulations. In this work we present a new analytical prescription for calculating a "resummed" non-Gaussian halo mass function, which is constructed to remain stable in the extreme tail. We show that the prescription works well in the parameter regime that has been currently explored in simulations. We then use Fisher matrix techniques to compare our prescription with an extrapolated fit to N-body simulations, which has recently been used to obtain constraints from data collected by the South Pole Telecope. We show that for the current data, both prescriptions would lead to statistically consistent constraints. As the data improve, however, there is a possibility of introducing a statistically significant bias in the constraints due to the choice of prescription, especially if non-Gaussianity is scale dependent and becomes relatively large on cluster scales. It would then be necessary to test the accuracy of the prescriptions in N-body simulations that can probe clusters with high masses and redshifts in the presence of large non-Gaussianity.

Bursts of particle production during inflation provide a well-motivated mechanism for creating bump like features in the primordial power spectrum. Current data constrains these features to be less than about 5% the size of the featureless primordial power spectrum at wavenumbers of about 0.1 h Mpc^{-1}. We forecast that the Planck cosmic microwave background experiment will be able to strengthen this constraint to the 0.5% level. We also predict that adding data from a square kilometer array (SKA) galaxy redshift survey would improve the constraint to about the 0.1% level. For features at larger wave-numbers, Planck will be limited by Silk damping and foregrounds. While, SKA will be limited by non-linear effects. We forecast for a Cosmic Inflation Probe (CIP) galaxy redshift survey, similar constraints can be achieved up to about a wavenumber of 1 h Mpc^{-1}.

A large bulk flow, which is in tension with the Lambda Cold Dark Matter ($\Lambda$CDM) cosmological model, has been observed. In this paper, we provide a physically plausible explanation of this bulk flow, based on the assumption that some fraction of the observed dipole in the cosmic microwave background is due to an intrinsic fluctuation, so that the subtraction of the observed dipole leads to a mismatch between the cosmic microwave background (CMB) defined rest frame and the matter rest frame. We investigate a model that takes into account the relative velocity (hereafter the tilted velocity) between the two frames, and develop a Bayesian statistic to explore the likelihood of this tilted velocity. By studying various independent peculiar velocity catalogs, we find that: (1) the magnitude of the tilted velocity $u$ is around 400 km/s, and its direction is close to what is found from previous bulk flow analyses; for most catalogs analysed, u=0 is excluded at about the $2.5 \sigma$ level;(2) constraints on the magnitude of the tilted velocity can result in constraints on the duration of inflation, due to the fact that inflation can neither be too long (no dipole effect) nor too short (very large dipole effect); (3) Under the assumption of a super-horizon isocurvature fluctuation, the constraints on the tilted velocity require that inflation lasts at least 6 e-folds longer (at the 95% confidence interval) than that required to solve the horizon problem. This opens a new window for testing inflation and models of the early Universe from observations of large scale structure.

Very massive high redshift clusters can be used to constrain and test the Lambda CDM model. Taking into account the observational constraints of Jee et al. (2009) we have calculated the probability for the most massive cluster to be found in the range (5.2-7.6)e14 M_sun, between redshifts 1.4-2.2, with a sky area of 11 sqdeg and under non-Gaussian initial conditions. Clusters constrain the non-Gaussianity on smaller scales than current cosmic microwave background or halo bias data and so can be used to test for running of the non-Gaussianity parameter f_NL.

We study how direct detection of the inflationary gravitational wave background constrains inflationary parameters and complements CMB polarization measurements. The error ellipsoids calculated using the Fisher information matrix approach with Planck and the direct detection experiment, BBO (Big Bang Observer), show different directions of parameter degeneracy, and the degeneracy is broken when they are combined. For a slow-roll parameterization, we show that BBO could significantly improve the constraints on the tensor-to-scalar ratio compared with Planck alone. We also look at a quadratic and a natural inflation model. In both cases, if the temperature of reheating is also treated as a free parameter, then the addition of BBO can significantly improve the error bars. In the case of natural inflation, we find that the addition of BBO could even partially improve the error bars of a cosmic variance-limited CMB experiment.

A "compensated" isocurvature perturbation consists of an overdensity (or underdensity) in the cold dark matter which is completely cancelled out by a corresponding underdensity (or overdensity) in the baryons. Such a configuration may be generated by a curvaton model of inflation if the cold dark matter is created before curvaton decay and the baryon number is created by the curvaton decay (or vice-versa). Compensated isocurvature perturbations, at the level producible by the curvaton model, have no observable effect on cosmic microwave background anisotropies or on galaxy surveys. They can be detected through their effect on the distribution of neutral hydrogen between redshifts 30 to 300 using 21 cm absorption observations. However, to obtain a good signal to noise ratio, very large observing arrays are needed. We estimate that a fast Fourier transform telescope would need a total collecting area of about 20 square kilometers to detect a curvaton generated compensated isocurvature perturbation at more than 5 sigma significance.

We investigate how well galaxy cluster number counts can constrain the primordial power spectrum. Measurements of the primary anisotropies in the cosmic microwave background (CMB) may be limited, by the presence of foregrounds from secondary sources, to probing the primordial power spectrum at wave numbers less than about 0.30 h Mpc^{-1}. We break up the primordial power spectrum into a number of nodes and interpolate linearly between each node. This allows us to show that cluster number counts could then extend the constraints on the form of the primordial power spectrum up to wave numbers of about 0.45 h Mpc^{-1}. We estimate combinations of constraints from PLANCK and SPT primary CMB and their respective SZ surveys. We find that their constraining ability is limited by uncertainties in the mass scaling relations. We also estimate the constraint from clusters detected from a SNAP like gravitational lensing survey. As there is an unambiguous and simple relationship between the filtered shear of the lensing survey and the cluster mass, it may be possible to obtain much tighter constraints on the primordial power spectrum in this case.

We estimate the solar system motion relative to the cosmic microwave background using type Ia supernovae (SNe) measurements. We take into account the correlations in the error bars of the SNe measurements arising from correlated peculiar velocities. Without accounting for correlations in the peculiar velocities, the SNe data we use appear to detect the peculiar velocity of the solar system at about the 3.5 sigma level. However, when the correlations are correctly accounted for, the SNe data only detects the solar system peculiar velocity at about the 2.5 sigma level. We forecast that the solar system peculiar velocity will be detected at the 9 sigma level by GAIA and the 11 sigma level by the LSST. For these surveys we find the correlations are much less important as most of the signal comes from higher redshifts where the number density of SNe is insufficient for the correlations to be important.

Bayesian model selection provides a formal method of determining the level of support for new parameters in a model. However, if there is not a specific enough underlying physical motivation for the new parameters it can be hard to assign them meaningful priors, an essential ingredient of Bayesian model selection. Here we look at methods maximizing the prior so as to work out what is the maximum support the data could give for the new parameters. If the maximum support is not high enough then one can confidently conclude that the new parameters are unnecessary without needing to worry that some other prior may make them significant. We discuss a computationally efficient means of doing this which involves mapping p-values onto upper bounds of the Bayes factor (or odds) for the new parameters. A p-value of 0.05 ($1.96\sigma$) corresponds to odds less than or equal to 5:2 which is below the `weak' support at best threshold. A p-value of 0.0003 ($3.6\sigma$) corresponds to odds of less than or equal to 150:1 which is the `strong' support at best threshold. Applying this method we find that the odds on the scalar spectral index being different from one are 49:1 at best. We also find that the odds that there is primordial hemispherical asymmetry in the cosmic microwave background are 9:1 at best.

We detect the correlated peculiar velocities of nearby type Ia supernovae (SNe), while highlighting an error in some of the literature. We find sigma_8=0.79 +/- 0.22 from SNe, and examine the potential of this method to constrain cosmological parameters in the future. We demonstrate that a survey of 300 low-z SNe (such as the nearby SNfactory) will underestimate the errors on w by about 35% if the coherent peculiar velocities are not included.

A linear modulation of the primordial perturbations is proposed as an explanation for the observed asymmetry between the northern and southern hemispheres of the Wilkinson Microwave Anisotropy Probe (WMAP) data. A cut sky, reduced resolution third year "Internal Linear Combination" (ILC) map was used to estimate the modulation parameters. A foreground template and a modulated plus unmodulated monopole and dipole were projected out of the likelihood. The effective chi squared is reduced by nine for three extra parameters. The mean Galactic colatitude and longitude, of the modulation, with 68%, 95% and 99.7% confidence intervals were 56^{+17 +36 +65}_{-17 -35 -51} and 63^{+28 +59 +105}_{-26 -58 -213}. The mean percentage change of the variance, across the pole's of the modulation, was 62^{+18 +35 +57}_{-18 -35 -47}. Implications of these results and possible generating mechanisms are discussed.

We introduce a class of models in which statistical isotropy is broken spontaneously in the CMB by a non-linear response to long-wavelength fluctuations in a mediating field. These fluctuations appear as a gradient locally and pick out a single preferred direction. The non-linear response imprints this direction in a range of multipole moments. We consider two manifestations of isotropy breaking: additive contributions and multiplicative modulation of the intrinsic anisotropy. Since WMAP exhibits an alignment of power deficits, an additive contribution is less likely to produce the observed alignments than the usual isotropic fluctuations, a fact which we illustrate with an explicit cosmological model of long-wavelength quintessence fluctuations. This problem applies to other models involving foregrounds or background anisotropy that seek to restore power to the CMB. Additive models that account directly for the observed power exacerbate the low power of the intrinsic fluctuations. Multiplicative models can overcome these difficulties. We construct a proof of principle model that significantly improves the likelihood and generates stronger alignments than WMAP in 30-45% of realizations.

There is some observational evidence that the dark energy may not be smooth on large scales. This makes it worth while to try and get as simple and as intuitive a picture of how dark energy perturbations behave so as to be able to better constrain possible models of dark energy and the generation of large scale perturbations. The separate Universe method provides an easy way to evaluate cosmological perturbations, as all that is required is an understanding of the background behavior. Here, this method is used to show how the size of the dark energy perturbations, preferred by observations, is larger than would be expected, and so some mechanism may be required to amplify them.

We propose a model which produces dark energy perturbations large enough to explain the lack of power seen at the quadrupole scale in the cosmic microwave background. If the dark energy is frozen from horizon exit during inflation until dark energy domination, then it is not possible to have perturbations in the dark energy which are large enough. We propose using a tachyonic amplification mechanism to overcome this. The dark energy is taken to be a complex scalar field, where the radial field has a Mexican hat potential. During inflation, the radial component is trapped near the maximum of its potential. At the end of inflation, it rolls down to the minimum. The dark energy today is taken to be a pseudo-Nambu-Goldstone boson. The perturbations generated during inflation are amplified by the rolling of the radial field. We also examine the use of the variable decay mechanism in order to generate an anti-correlation between the dark energy perturbations and the curvature perturbation. We show that using this mechanism then constrains the properties of the dark energy and its evolution from redshift one until today.

We explicate the origin of the temperature quadrupole in the adiabatic dark energy model and explore the mechanism by which scale invariant isocurvature dark energy perturbations can lead to its sharp suppression. The model requires anticorrelated curvature and isocurvature fluctuations and is favored by the WMAP data at about the 95% confidence level in a flat scale invariant model. In an inflationary context, the anticorrelation may be established if the curvature fluctuations originate from a variable decay rate of the inflaton; such models however tend to overpredict gravitational waves. This isocurvature model can in the future be distinguished from alternatives involving a reduction in large scale power or modifications to the sound speed of the dark energy through the polarization and its cross correlation with the temperature. The isocurvature model retains the same polarization fluctuations as its adiabatic counterpart but reduces the correlated temperature fluctuations. We present a pedagogical discussion of dark energy fluctuations in a quintessence and k-essence context in the Appendix.

In the curvaton model of inflation, where a second scalar field, the "curvaton", is responsible for the observed inhomogeneity, a non-zero neutrino degeneracy may lead to a characteristic pattern of isocurvature perturbations in the neutrino, cold dark matter and baryon components. We find the current data can only place upper limits on the level of isocurvature perturbations. These can be translated into upper limits on the neutrino degeneracy parameter. In the case that lepton number is created before curvaton decay, we find that the limit on the neutrino degeneracy parameter is comparable with that obtained from Big-bang nucleosynthesis. For the case that lepton number is created by curvaton decay we find that the absolute value of the non-Gaussianity parameter, |f_nl|, must be less than 10 (95% confidence interval).

Simple curvaton models can generate a mixture of of correlated primordial adiabatic and isocurvature perturbations. The baryon and cold dark matter isocurvature modes differ only by an observationally null mode in which the two perturbations almost exactly compensate, and therefore have proportional effects at linear order. We discuss the CMB anisotropy in general mixed models, and give a simple approximate analytic result for the large scale CMB anisotropy. Working numerically we use the latest WMAP observations and a variety of other data to constrain the curvaton model. We find that models with an isocurvature contribution are not favored relative to simple purely adiabatic models. However a significant primordial totally correlated baryon isocurvature perturbation is not ruled out. Certain classes of curvaton model are thereby ruled out, other classes predict enough non-Gaussianity to be detectable by the Planck satellite. In the appendices we review the relevant equations in the covariant formulation and give series solutions for the radiation dominated era.

We investigate the effects on cosmological density perturbations of dark radiation in a Randall-Sundrum 2 type brane-world. Dark radiation in the background is limited by observational constraints to be a small fraction of the radiation energy density, but it has an interesting qualitative effect in the radiation era. On large scales, it serves to slightly suppress the radiation density perturbations at late times, while boosting the perturbations in dark radiation. In a kinetic (stiff) era, the suppression is much stronger, and drives the density perturbations to zero.

The ekpyrotic and cyclic universe scenarios have revived the idea that the density perturbations apparent in today's universe could have been generated in a `pre-singularity' epoch before the big bang. These scenarios provide explicit mechanisms whereby a scale invariant spectrum of adiabatic perturbations may be generated without the need for cosmic inflation, albeit in a phase preceding the hot big bang singularity. A key question they face is whether there exists a unique prescription for following perturbations through the bounce, an issue which is not yet definitively settled. This goal of this paper is more modest, namely to study a bouncing Universe model in which neither General Relativity nor the Weak Energy Condition is violated. We show that a perturbation which is pure growing mode before the bounce does not match to a pure decaying mode perturbation after the bounce. Analytical estimates of when the comoving curvature perturbation varies around the bounce are given. It is found that in general it is necessary to evaluate the evolution of the perturbation through the bounce in detail rather than using matching conditions.

We compare the latest cosmic microwave background data with theoretical predictions including correlated adiabatic and CDM isocurvature perturbations with a simple power-law dependence. We find that there is a degeneracy between the amplitude of correlated isocurvature perturbations and the spectral tilt. A negative (red) tilt is found to be compatible with a larger isocurvature contribution. Estimates of the baryon and CDM densities are found to be almost independent of the isocurvature amplitude. The main result is that current microwave background data do not exclude a dominant contribution from CDM isocurvature fluctuations on large scales.

This PhD thesis presents a study of the effect and generation of non-adiabatic perturbations in Cosmology. We study adiabatic (curvature) and entropy (isocurvature) perturbations produced during a period of cosmological inflation that is driven by multiple scalar fields with an arbitrary interaction potential. A local rotation in field space is performed to separate out the adiabatic and entropy modes. The resulting field equations show explicitly how on large scales entropy perturbations can source adiabatic perturbations if the background solution follows a curved trajectory in field space, and how adiabatic perturbations cannot source entropy perturbations in the long-wavelength limit. We show why one in general expects the adiabatic and entropy perturbations to be correlated at the end of inflation. Then, we consider two-field preheating after inflation, examining conditions under which entropy perturbations can alter the large-scale curvature perturbation and showing how our new formalism has advantages in numerical stability when the background solution follows a non-trivial trajectory in field space. Then we compare the latest cosmic microwave background data with theoretical predictions including correlated adiabatic and CDM isocurvature perturbations. We find that there is a degeneracy between the amplitude of correlated isocurvature perturbations and the spectral tilt. We then study perturbations in Randall-Sundrum-type brane-world cosmologies. The density perturbations generate Weyl curvature in the bulk, which in turn backreacts on the brane via stress-energy perturbations. Bulk effects produce a non-adiabatic mode, even when the matter perturbations are adiabatic, and alter the background dynamics.

The initial condition problem of inflation is examined from the perspective of both spacetime embedding and scalar field dynamics. The spacetime embedding problem is solved for arbitrary initial spatial curvature Omega, which generalizes previous works that primarily treat the flat case Omega=1. Scalar field dynamics that is consistent with the embedding constraints are examined, with the additional treatment of damping effects. The effects of inhomogeneities on the embedding problem also are considered. A category of initial conditions are identified that are not acausal and can develop into an inflationary regime.

In Randall-Sundrum-type brane-world cosmologies, density perturbations generate Weyl curvature in the bulk, which in turn backreacts on the brane via stress-energy perturbations. On large scales, the perturbation equations contain a closed system on the brane, which may be solved without solving for the bulk perturbations. Bulk effects produce a non-adiabatic mode, even when the matter perturbations are adiabatic, and alter the background dynamics. As a consequence, the standard evolution of large-scale fluctuations in general relativity is modified. The metric perturbation on large-scales is not constant during high-energy inflation. It is constant during the radiation era, except at most during the very beginning, if the energy is high enough.

We study adiabatic (curvature) and entropy (isocurvature) perturbations produced during a period of cosmological inflation that is driven by multiple scalar fields with an arbitrary interaction potential. A local rotation in field space is performed to separate out the adiabatic and entropy modes. The resulting field equations show explicitly how on large scales entropy perturbations can source adiabatic perturbations if the background solution follows a curved trajectory in field space, and how adiabatic perturbations cannot source entropy perturbations in the long-wavelength limit. It is the effective mass of the entropy field that determines the amplitude of entropy perturbations during inflation. We present two applications of the equations. First, we show why one in general expects the adiabatic and entropy perturbations to be correlated at the end of inflation, and calculate the cross-correlation in the context of a double inflation model with two non-interacting fields. Second, we consider two-field preheating after inflation, examining conditions under which entropy perturbations can alter the large-scale curvature perturbation and showing how our new formalism has advantages in numerical stability when the background solution follows a non-trivial trajectory in field space.

The relative growth of field and metric perturbations during preheating is sensitive to initial conditions set in the preceding inflationary phase. Recent work suggests this may protect super-Hubble metric perturbations from resonant amplification during preheating. We show that this possibility is fragile and sensitive to the specific form of the interactions between the inflaton and other fields. The suppression is naturally absent in two classes of preheating in which either (1) the vacua of the non-inflaton fields during inflation are deformed away from the origin, or (2) the effective masses of non-inflaton fields during inflation are small but during preheating are large. Unlike the simple toy model of a $g^2 \phi^2 \chi^2$ coupling, most realistic particle physics models contain these other features. Moreover, they generically lead to both adiabatic and isocurvature modes and non-Gaussian scars on super-Hubble scales. Large-scale coherent magnetic fields may also appear naturally.