We study the sensitivity of cosmological observables to the reheating phase following inflation driven by many scalar fields. We describe a method which allows semi-analytic treatment of the impact of perturbative reheating on cosmological perturbations using the sudden decay approximation. Focusing on $\mathcal{N}$-quadratic inflation, we show how the scalar spectral index and tensor-to-scalar ratio are affected by the rates at which the scalar fields decay into radiation. We find that for certain choices of decay rates, reheating following multiple-field inflation can have a significant impact on the prediction of cosmological observables.

The flatness of the inflaton potential and lightness of the Higgs could have the common origin of the breaking of a global symmetry. This scenario provides a unified framework of Goldstone Inflation and Composite Higgs, where the inflaton and the Higgs both have a pseudo--Goldstone boson nature. The inflaton reheats the Universe via decays to the Higgs and subsequent secondary production of other SM particles via the top and massive vector bosons. We find that inflationary predictions and perturbative reheating conditions are consistent with CMB data for sub--Planckian values of the fields, as well as opening up the possibility of inflation at the TeV scale. We explore this exciting possibility, leading to an interplay between collider data and cosmological constraints.

Observations of the microwave background fluctuations suggest a scale-dependent amplitude asymmetry of roughly 2.5 sigma significance. Inflationary explanations for this 'anomaly' require non-Gaussian fluctuations which couple observable modes to those on much larger scales. In this Letter we describe an analysis of such scenarios which significantly extends previous treatments. We identify the non-Gaussian 'response function' which characterizes the asymmetry, and show that it is non-trivial to construct a model which yields a sufficient amplitude: many independent fine tunings are required, often making such models appear less likely than the anomaly they seek to explain. We present an explicit model satisfying observational constraints and determine for the first time how large its bispectrum would appear to a Planck-like experiment. Although this model is merely illustrative, we expect it is a good proxy for the bispectrum in a sizeable class of models which generate a scale-dependent response using a large eta parameter.

If the primordial bispectrum is sufficiently large then the CMB hemispherical asymmetry may be explained by a large-scale mode of exceptional amplitude which perturbs the zeta two-point function. We extend previous calculations, which were restricted to one- or two-source scenarios, by providing a method to compute the response of the two-point function in any model yielding a 'local-like' bispectrum. In general, this shows that it is not the reduced bispectrum fNL which sources the amplitude and scale-dependence of the mode coupling but rather a combination of 'response functions'. We discuss why it is difficult to construct successful scenarios and enumerate the fine-tunings which seem to be required. Finally, we exhibit a concrete model which can be contrived to match the observational constraints and show that to a Planck-like experiment it would appear to have |fNL-local| ~ |fNL-equi| ~ |fNL-ortho| ~ 1. Therefore, contrary to previous analyses, we conclude that it is possible to generate the asymmetry while respecting observational constraints on the bispectrum and low-ell multipoles even without tuning our location on the long-wavelength mode.

We introduce an alternative parametrisation for the scale dependence of the non-linearity parameter $f_{\rm NL}$ in quasi-local models of non-Gaussianity. Our parametrisation remains valid when $f_{\rm NL}$ changes sign, unlike the commonly adopted power law ansatz $f_{\rm NL}(k) \propto k^{ n_{f_{\rm NL}} }$. We motivate our alternative parametrisation by appealing to the self-interacting curvaton scenario, and as an application, we apply it to the CMB power asymmetry. Explaining the power asymmetry requires a strongly scale dependent non-Gaussianity. We show that regimes of model parameter space where $f_{\rm NL}$ is strongly scale dependent are typically associated with a large $g_{\rm NL}$ and quadrupolar power asymmetry, which can be ruled out by existing observational constraints.

We study the impact of perturbative reheating on primordial observables in models of multiple-field inflation. By performing a sudden decay calculation, we derive analytic expressions for the local-type non-linearity parameter $f_{\rm NL}^{\rm local}$, the scalar spectral index $n_\zeta$, and the tensor-to-scalar ratio $r_T$ as functions of the decay rates of the inflationary fields. We compare our analytic results to a fully numerical classical field theory simulation, finding excellent agreement. We find that the sensitivity of $f_{\rm NL}$, $n_\zeta$, and $r_T$ to the reheating phase depends heavily on the underlying inflationary model. We quantify this sensitivity, and discuss conditions that must be satisfied if observable predictions are to be insensitive to the dynamics of reheating. We demonstrate that upon completion of reheating, all observable quantities take values within finite ranges, the limits of which are determined completely by the conditions during inflation. Furthermore, fluctuations in both fields play an important role in determining the full dependence of the observables on the dynamics of reheating. By applying our formalism to two concrete examples, we demonstrate that variations in $f_{\rm NL}$, $n_\zeta$, and $r_T$ caused by changes in reheating dynamics are well within the sensitivity of Planck, and as such the impact of reheating must be accounted for when making predictions for models of multiple-field inflation. Our final expressions are very general, encompassing a wide range of two-field inflationary models, including the standard curvaton scenario. Our results allow a much more unified approach to studying two-field inflation including the effects of perturbative reheating. As such, entire classes of models can be studied together, allowing a more systematic approach to gaining insight into the physics of the early universe through observation.

We propose a completely new parametrisation of the dark energy equation of state, which uses the dark energy density, $\Omega_e$ as a cosmic clock. We expand the equation of state in a series of orthogonal polynomials, with $\Omega_e$ as the expansion parameter and determine the expansion coefficients by fitting to SNIa and $H(z)$ data. Assuming that $\Omega_e$ is a monotonic function of time, we show that our parametrisation performs better than the popular Chevallier--Polarski--Linder (CPL) and Gerke and Efstathiou (GE) parametrisations, and we demonstrate that it is robust to the choice of prior. Expanding in orthogonal polynomials allows us to relate models of dark energy directly to our parametrisation, which we illustrate by placing constraints on the expansion coefficients extracted from two popular quintessence models. Finally, we comment on how this parametrisation could be modified to accommodate high redshift data, where any non--monotonicity of $\Omega_e$ would need to be accounted for.

Many models have been studied that contain more than one species of dark matter and some of these couple the Cold Dark Matter (CDM) to a light scalar field. In doing this we introduce additional long range forces, which in turn can significantly affect our estimates of cosmological parameters if not properly accounted for. It is, therefore, important to study these models and their resulting cosmological implications. We present a model in which a fraction of the total cold dark matter density is coupled to a scalar field. We study the background and perturbation evolution and calculate the resulting Cosmic Microwave Background anisotropy spectra. The greater the fraction of dark matter coupled to the scalar field and the stronger the coupling strength, the greater the deviation of the background evolution from LCDM. Previous work, with a single coupled dark matter species, has found an upper limit on the coupling strength of order O(0.1). We find that with a coupling of this magnitude more than half the dark matter can be coupled to a scalar field without producing any significant deviations from LCDM.

We study the evolution of the non-linear curvature perturbation during perturbative reheating, and hence how observables evolve to their final values which we may compare against observations. Our study includes the evolution of the two trispectrum parameters, $\gnl$ and $\taunl$, as well as the scale dependence of both $\fnl$ and $\taunl$. In general the evolution is significant and must be taken into account, which means that models of multifield inflation cannot be compared to observations without specifying how the subsequent reheating takes place. If the trispectrum is large at the end of inflation, it normally remains large at the end of reheating. In the classes of models we study, it is very hard to generate $\taunl\gg\fnl^2$, regardless of the decay rates of the fields. Similarly, for the classes of models in which $\gnl\simeq\taunl$ during slow--roll inflation, we find the relation typically remains valid during reheating. Therefore it is possible to observationally test such classes of models without specifying the parameters of reheating, even though the individual observables are sensitive to the details of reheating. It is hard to generate an observably large $\gnl$ however. The runnings, $\nfnl$ and $\ntaunl$, tend to satisfy a consistency relation $\ntaunl=(3/2)\nfnl$, but are in general too small to be observed for the class of models considered regardless of reheating timescale.

This paper is concerned with string cosmology and the dynamics of multiple scalar fields in potentials that can become negative, and their features as (Early) Dark Energy models. Our point of departure is the "String Axiverse", a scenario that motivates the existence of cosmologically light axion fields as a generic consequence of string theory. We couple such an axion to its corresponding modulus. We give a detailed presentation of the rich cosmology of such a model, ranging from the setting of initial conditions on the fields during inflation, to the asymptotic future. We present some simplifying assumptions based on the fixing of the axion decay constant $f_a$, and on the effective field theory when the modulus trajectory is adiabatic, and find the conditions under which these assumptions break down. As a by-product of our analysis, we find that relaxing the assumption of fixed $f_a$ leads to the appearance of a new meta-stable de-Sitter region for the modulus without the need for uplifting by an additional constant. A dynamical systems analysis reveals the existence of many fixed point attractors, repellers and saddle points, which we analyse in detail. We also provide geometric interpretations of the phase space. The fixed points can be used to bound the couplings in the model. A systematic scan of certain regions of parameter space reveals that the future evolution of the universe in this model can be rich, containing multiple epochs of accelerated expansion.

We study the effects of perturbative reheating on the evolution of the curvature perturbation \zeta, in two-field inflation models. We use numerical methods to explore the sensitivity of f_NL, n_s and r to the reheating process, and present simple qualitative arguments to explain our results. In general, if a large non-Gaussian signal exists at the start of reheating, it will remain non zero at the end of reheating. Unless all isocurvature modes have completely decayed before the start of reheating, we find that the non-linearity parameter, f_NL, can be sensitive to the reheating timescale, and that this dependence is most appreciable for `runaway' inflationary potentials that only have a minimum in one direction. For potentials with a minimum in both directions, f_NL can also be sensitive to reheating if a mild hierarchy exists between the decay rates of each field. Within the class of models studied, we find that the spectral index n_s, is fairly insensitive to large changes in the field decay rates, indicating that n_s is a more robust inflationary observable, unlike the non-linearity parameter f_NL. Our results imply that the statistics of \zeta, especially f_NL, can only be reliably used to discriminate between models of two-field inflation if the physics of reheating are properly accounted for.

A sufficiently light scalar field slowly evolving in a potential can account for the dark energy that presently dominates the universe. This quintessence field is expected to couple directly to matter components, unless some symmetry of a more fundamental theory protects or suppresses it. Such a coupling would leave distinctive signatures in the background expansion history of the universe and on cosmic structure formation, particularly at galaxy cluster scales. Using semi--analytic expressions for the CDM halo mass function, we make predictions for halo abundance in models where the quintessence scalar field is coupled to cold dark matter, for a variety of quintessence potentials. We evaluate the linearly extrapolated density contrast at the redshift of collapse using the spherical collapse model and we compare this result to the corresponding prediction obtained from the non--linear perturbation equations in the Newtonian limit. For all the models considered in this work, if there is a continuous flow of energy from the quintessence scalar field to the CDM component, then the predicted number of CDM haloes can only lie below that of $\Lambda$CDM, when each model shares the same cosmological parameters today. In the last stage of our analysis we perform a global MCMC fit to data to find the best fit values for the cosmological model parameters. We find that for some forms of the quintessence potential, coupled dark energy models can offer a viable alternative to $\Lambda$CDM in light of the recent detections of massive high--$z$ galaxy clusters, while other models of coupled quintessence predict a smaller number of massive clusters at high redshift compared to $\Lambda$CDM.