In the context of the dark energy scenario, the Einstein Yang-Mills Higgs model in the SO(3) representation was studied for the first time by M. Rinaldi (see JCAP 1510, 023 (2015)) in a homogeneous and isotropic spacetime. We revisit this model, finding in particular that the interaction between the Higgs field and the gauge fields generates contributions to the momentum density, anisotropic stress and pressures, thus making the model inconsistent with the assumed background. We instead consider a homogeneous but anisotropic Bianchi-I spacetime background in this paper and analyze the corresponding dynamical behaviour of the system. We find that the only attractor point corresponds to an isotropic accelerated expansion dominated by the Higgs potential. However, the model predicts non-negligible anisotropic shear contributions nowadays, i.e. the current universe can have hair although it will loose it in the future. We investigate the evolution of the equation of state for dark energy and highlight some possible consequences of its behaviour related to the process of large-scale structure formation. As a supplement, we propose the "Higgs triad" as a possibility to make the Einstein Yang-Mills Higgs model be consistent with a homogeneous and isotropic spacetime.

We study a new class of vector dark energy models where multi-Proca fields $A_\mu^a$ are coupled to cold dark matter by the term $f(X)\tilde{\mathcal{L}}_{m}$ where $f(X)$ is a general function of $X\equiv -\frac{1}{2}A^\mu_ a A^a_\mu$ and $\tilde{\mathcal{L}}_{m}$ is the cold dark matter Lagrangian. From here, we derive the general covariant form of the novel interaction term sourcing the field equations. This result is quite general in the sense that encompasses Abelian and non-Abelian vector fields. In particular, we investigate the effects of this type of coupling in a simple dark energy model based on three copies of canonical Maxwell fields to realize isotropic expansion. The cosmological background dynamics of the model is examined by means of a dynamical system analysis to determine the stability of the emergent cosmological solutions. As an interesting result, we find that the coupling function leads to the existence of a novel scaling solution during the dark matter dominance. Furthermore, the critical points show an early contribution of the vector field in the form of dark radiation and a stable de Sitter-type attractor at late times mimicking dark energy. The cosmological evolution of the system as well as the aforementioned features are verified by numerical computations. Observational constraints are also discussed to put the model in a more phenomenological context in the light of future observations.

Under the same spirit of the Galileon-Horndeski theories and their more modern extensions, the generalized SU(2) Proca theory was built by demanding that its action may be free of the Ostrogradski's instability. Nevertheless, the theory must also be free of other instability problems in order to ensure its viability. As a first approach to address this issue, we concentrate on a quite general variant of the theory and investigate the general conditions for the absence of ghost and gradient instabilities in the tensor sector without the need for resolving the dynamical background. The phenomenological interest of this approach as well as of the variant investigated lies on the possibility of building cosmological models driven solely by non-Abelian vector fields that may account for a successful description of both the early inflation and the late-time accelerated expansion of the universe.

To date, different alternative theories of gravity, although related, involving Proca fields have been proposed. Unfortunately, the procedure to obtain the relevant terms in some formulations has not been systematic enough or exhaustive, thus resulting in some missing terms or ambiguity in the process carried out. In this paper, we propose a systematic procedure to build the beyond generalized theory for a Proca field in four dimensions containing only the field itself and its first-order derivatives. We examine the validity of our procedure at the fourth level of the generalized Proca theory. In our approach, we employ all the possible Lorentz-invariant Lagrangian pieces made of the Proca field and its first-order derivatives, including those that violate parity, and find the relevant combination that propagates only three degrees of freedom and has healthy dynamics for the longitudinal mode. The key step in our procedure is to retain the flat space-time divergences of the currents in the theory during the covariantization process. In the curved space-time theory, some of the retained terms are no longer current divergences so that they induce the new terms that identify the beyond generalized Proca field theory. The procedure constitutes a systematic method to build general theories for multiple vector fields with or without internal symmetries.

Inspired in the Standard Model of Elementary Particles, the Einstein Yang-Mills Higgs action with the Higgs field in the SU(2) representation was proposed in Class. Quantum Grav. 32 (2015) 045002 as the element responsible for the dark energy phenomenon. We revisit this action emphasizing in a very important aspect not sufficiently explored in the original work and that substantially changes its conclusions. This aspect is the role that the Yang-Mills Higgs interaction plays at fixing the gauge for the Higgs field, in order to sustain a homogeneous and isotropic background, and at driving the late accelerated expansion of the Universe by moving the Higgs field away of the minimum of its potential and holding it towards an asymptotic finite value. We analyse the dynamical behaviour of this system and supplement this analysis with a numerical solution whose initial conditions are in agreement with the current observed values for the density parameters. This scenario represents a step towards a successful merging of cosmology and well-tested particle physics phenomenology.

We study the role that a cosmic triad in the generalized $SU(2)$ Proca theory, specifically in one of the pieces of the Lagrangian that involves the symmetric version $S_{\mu \nu}$ of the gauge field strength tensor $F_{\mu \nu}$, has on dark energy and primordial inflation. Regarding dark energy, the triad behaves asymptotically as a couple of radiation perfect fluids whose energy densities are negative for the $S$ term but positive for the Yang-Mills term. This leads to an interesting dynamical fine-tuning mechanism that gives rise to a combined equation of state parameter $\omega \simeq -1$ and, therefore, to an eternal period of accelerated isotropic expansion for an ample spectrum of initial conditions. Regarding primordial inflation, one of the critical points of the associated dynamical system can describe a prolonged period of isotropic slow-roll inflation sustained by the $S$ term. This period ends up when the Yang-Mills term dominates the energy density leading to the radiation dominated epoch. Unfortunately, in contrast to the dark energy case, the primordial inflation scenario is strongly sensitive to the coupling constants and initial conditions. The whole model, including the other pieces of the Lagrangian that involve $S_{\mu \nu}$, might evade the recent strong constraints coming from the gravitational wave signal GW170817 and its electromagnetic counterpart GRB 170817A.

An alternative for the construction of fundamental theories is the introduction of Galileons. These are fields whose action leads to non higher than second-order equations of motion. As this is a necessary but not sufficient condition to make the Hamiltonian bounded from below, as long as the action is not degenerate, the Galileon construction is a way to avoid pathologies both at the classical and quantum levels. Galileon actions are, therefore, of great interest in many branches of physics, specially in high energy physics and cosmology. This proceedings contribution presents the generalities of the construction of both scalar and vector Galileons following two different but complimentary routes.

Following previous works on generalized Abelian Proca theory, also called vector Galileon, we investigate the massive extension of an SU(2) gauge theory, i.e., the generalized SU(2) Proca model, which could be dubbed non-Abelian vector Galileon. This particular symmetry group permits fruitful applications in cosmology such as inflation driven by gauge fields. Our approach consists in building, in an exhaustive way, all the Lagrangians containing up to six contracted Lorentz indices. For this purpose, and after identifying by group theoretical considerations all the independent Lagrangians which can be written at these orders, we consider the only linear combinations propagating three degrees of freedom and having healthy dynamics for their longitudinal mode, i.e., whose pure St\"uckelberg contribution turns into the SU(2) multi-Galileon dynamics. Finally, and after having considered the curved space-time expansion of these Lagrangians, we discuss the form of the theory at all subsequent orders.

We summarize previous results on the most general Proca theory in 4 dimensions containing only first-order derivatives in the vector field (second-order at most in the associated St\"uckelberg scalar) and having only three propagating degrees of freedom with dynamics controlled by second-order equations of motion. Discussing the Hessian condition used in previous works, we conjecture that, as in the scalar galileon case, the most complete action contains only a finite number of terms with second-order derivatives of the St\"uckelberg field describing the longitudinal mode, which is in agreement with the results of JCAP 1405, 015 (2014) and Phys. Lett. B 757, 405 (2016) and complements those of JCAP 1602, 004 (2016). We also correct and complete the parity violating sector, obtaining an extra term on top of the arbitrary function of the field $A_\mu$, the Faraday tensor $F_{\mu \nu}$ and its Hodge dual $\tilde{F}_{\mu \nu}$.

Gauge-flation model at zeroth order in cosmological perturbation theory offers an interesting scenario for realizing inflation within a particle physics context, allowing us to investigate interesting possible connections between inflation and the subsequent evolution of the Universe. Difficulties, however, arise at the perturbative level, thus motivating a modification of the original model. In order to agree with the latest Planck observations, we modify the model such that the new dynamics can produce a relation between the spectral index $n_{s}$ and the tensor-to-scalar ratio $r$ allowed by the data. By including an identical mass term for each of the fields of the system, we find interesting dynamics leading to slow-roll inflation of the right length. The presence of the mass term has the potential to modify the $n_{s}$ vs. $r$ relation so as to agree with the data. As a first step, we study the model at zeroth order in cosmological perturbation theory, finding the conditions required for slow-roll inflation and the number of e-foldings of inflation. Numerical solutions are used to explore the impact of the mass term. We conclude that the massive version of Gauge-flation offers a viable inflationary model.

We revisit the most general theory for a massive vector field with derivative self-interactions, extending previous works on the subject to account for terms having trivial total derivative interactions for the longitudinal mode. In the flat spacetime (Minkowski) case, we obtain all the possible terms containing products of up to five first-order derivatives of the vector field, and provide a conjecture about higher-order terms. Rendering the metric dynamical, we covariantize the results and add all possible terms implying curvature.

We will expose in this paper our advances towards a proof of the equivalence between FRW background expansion, during some period of time that contains primordial inflation, and the statistical isotropy of the primordial curvature perturbation $\zeta$ at the end of this period of time. Our motivation rests on the growing interest in the existence of a preferred direction in the Universe hinted by the continuous presence of anomalies in the CMB data.

There are many inflationary models that allow the formation of the large-scale structure of the observable universe. The non-gaussianity parameter $f_{NL}$ is a useful tool to discriminate among these cosmological models when comparing the theoretical predictions with the satellite survey results like those from WMAP and Planck. The goal of this proceeding contribution is to review the moment transport equations methodology and the subsequent calculation of the $f_{NL}$ parameter.

We analyze two possible vector-field models using the techniques of dynamical systems. The first model involves a U(1)-vector field and the second a triad of SU(2)-vector fields. Both models include a gauge-fixing term and a power-law potential. A dynamical system is formulated and it is found that one of the critical points, for each model, corresponds to inflation, the origin of these critical points being the respective gauge-fixing terms. The conditions for the existence of an inflationary era which lasts for at least 60 efolds are studied.

We consider cosmological inflationary models in which vector fields play some role in the generation of the primordial curvature perturbation $\zeta$. Such models are interesting because the involved vector fields naturally seed statistical anisotropy in the primordial fluctuations which could eventually leave a measurable imprint on the cosmic microwave background fluctuations. In this article, we estimate the scale and shape dependent effects on the non-Gaussianity (NG) parameters due to the scale dependent statistical anisotropy in the distribution of the fluctuations. For concreteness, we use a power spectrum (PS) of the fluctuations of the quadrupolar form: $P_\zeta(\vec{k})\equiv P_\zeta(k)[1+g_\zeta(k)(\hat{n} \cdot \hat{k})^2 ]$, where $g_{\zeta}(k)$ is the only quantity which parametrizes the level of statistical anisotropy and $\hat{n}$ is a unitary vector which points towards the preferred direction. Then, we evaluate the contribution of the running of $g_{\zeta}(k)$ on the NG parameters by means of the $\delta N$ formalism. We focus specifically on the details for the $f_{\rm NL}$ NG parameter, associated with the bispectrum $B_\zeta$, but the structure of higher order NG parameters is straightforward to generalize. Although the level of statistical anisotropy in the PS is severely constrained by recent observations, the importance of statistical anisotropy signals in higher order correlators remains to be determined, this being the main task that we address here. The precise measurement of the shape and scale dependence (or running) of statistical parameters such as the NG parameters and the statistical anisotropy level could provide relevant elements for model building and for the determination of the presence (or nonpresence) of inflationary vector fields and their role in the inflationary mechanism.

We study the most general contributions due to scalar field perturbations, vector field perturbations, and anisotropic expansion to the generation of statistical anisotropy in the primordial curvature perturbation \zeta. Such a study is done using the \delta N formalism where only linear terms are considered. Here, we consider two specific cases that lead to determine the power spectrum P_\zeta(k) of the primordial curvature perturbation. In the first one, we consider the possibility that the n-point correlators of the field perturbations in real space are invariant under rotations in space (statistical isotropy); as a result, we obtain as many levels of statistical anisotropy as vector fields present and, therefore, several preferred directions. The second possibility arises when we consider anisotropic expansion, which leads us to obtain I+a additional contributions to the generation of statistical anisotropy of \zeta compared with the former case, being I and a the number of scalar and vector fields involved respectively.

We analyze a massive vector field with a non-canonical kinetic term in the action, minimally coupled to gravity, where the mass and kinetic function of the vector field vary as functions of time during inflation. The vector field is introduced following the same idea of a scalar curvaton, which must not affect the inflationary dynamics since its energy density during inflation is negligible compared to the total energy density in the Universe. Using this hypothesis, the vector curvaton will be solely responsible for generating the primordial curvature perturbation \zeta. We have found that the spectra of the vector field perturbations are scale-invariant in superhorizon scales due to the suitable choice of the time dependence of the kinetic function and the effective mass during inflation. The preferred direction, generated by the vector field, makes the spectrum of \zeta depend on the wavevector, i.e. there exists statistical anisotropy in \zeta. This is discussed principally in the case where the mass of the vector field increases with time during inflation, where it is possible to find a heavy field (M >> H) at the end of inflation, making the particle production be practically isotropic; thus, the longitudinal and transverse spectra are nearly the same order which in turn causes that the statistical anisotropy generated by the vector field is within the observational bounds.

We present the different consistency relations that can be seen as variations of the well known Suyama-Yamaguchi (SY) consistency relation \tau_{NL} \geqslant (\frac{6}{5} f_{NL})^2, the latter involving the levels of non-gaussianity f_{NL} and \tau_{NL} in the primordial curvature perturbation \zeta, they being scale-invariant. We explicitly state under which conditions the SY consistency relation has been claimed to hold in its different varieties (implicitly) presented in the literature since its inception back in 2008; as a result, we show for the first time that the variety \tau_{NL} ({\bf k}_1, {\bf k}_1) \geqslant (\frac{6}{5} f_{NL} ({\bf k}_1))^2, which we call "the fifth variety", is always satisfied even when there is strong scale-dependence and high levels of statistical anisotropy as long as statistical homogeneity holds: thus, an observed violation of this specific variety would prevent the comparison between theory and observation, shaking this way the foundations of cosmology as a science.

We present the different consistency relations that can be seen as variations of the well known Suyama-Yamaguchi (SY) consistency relation \tau_{NL} \geqslant ((6/5) f_{NL})^2. It has been claimed that the following variation: \tau_{NL} ({\bf k}_1, {\bf k_3}) \geqslant (6/5)^2 f_{NL} ({\bf k}_1) f_{NL} ({\bf k}_3), which we call "the fourth variety", in the collapsed (for \tau_{NL}) and squeezed (for f_{NL}) limits is always satisfied independently of any physics; however, the proof depends sensitively on the assumption of scale-invariance which only applies for cosmological models involving Lorentz-invariant scalar fields (at least at tree level), leaving room for a strong violation of this variety of the consistency relation when non-trivial degrees of freedom, for instance vector fields, are in charge of the generation of \zeta. With this in mind as a motivation, we explicitly state under which conditions the SY consistency relation has been claimed to hold in its different varieties (implicitly) presented in the literature; as a result, we show for the first time that the variety \tau_{NL} ({\bf k}_1, {\bf k}_1) \geqslant ((6/5) f_{NL} ({\bf k}_1))^2, which we call "the fifth variety", is always satisfied even when there is strong scale-dependence as long as statistical homogeneity holds: thus, an observed violation of this specific variety would prevent the comparison between theory and observation, shaking this way the foundations of cosmology as a science. Later, we concern about the existence of non-trivial degrees of freedom, concretely vector fields for which the levels of non-gaussianity have been calculated for very few models, finding that the fourth variety of the SY consistency relation is indeed strongly violated for some specific wavevector configurations while the fifth variety continues to be well satisfied. (Abridged)

We consider inflationary models in which vector fields are responsible for part or eventually all of the primordial curvature perturbation \zeta. Such models are phenomenologically interesting since they naturally introduce anisotropies in the probability distribution function of the primordial fluctuations that can leave a measurable imprint in the cosmic microwave background. Assuming that non-Gaussianity is generated due to the superhorizon evolution, we use the \delta N formalism to do a complete tree level calculation of the non-Gaussianity parameters f_{NL} and \tau_{NL} in the presence of vector fields. We isolate the isotropic pieces of the non-Gaussianity parameters, which anyway have contributions from the vector fields, and show that they obey the Suyama-Yamaguchi consistency relation \tau^{iso}_{NL}>=(6/5f^{iso}_{NL})^2. Other ways of defining the non-Gaussianity parameters, which could be observationally relevant, are stated and the respective Suyama-Yamaguchi-like consistency relations are obtained.

We argue that the true nature of the renormalizability of Horava-Lifshitz gravity lies in the presence of higher order spatial derivatives and not in the anisotropic Lifshitz scaling of space and time. We discuss the possibility of constructing a higher order spatial derivatives model that has the same renormalization properties of Horava-Lifshitz gravity but that does not make use of the Lifshitz scaling. In addition, the state-of-the-art of the Lorentz symmetry restoration in Horava-Lifshitz-type theories of gravitation is reviewed.

A diagrammatic approach to calculate n-point correlators of the primordial curvature perturbation \zeta was developed a few years ago following the spirit of the Feynman rules in Quantum Field Theory. The methodology is very useful and time-saving, as it is for the case of the Feynman rules in the particle physics context, but, unfortunately, is not very well known by the cosmology community. In the present work, we extend such an approach in order to include not only scalar field perturbations as the generators of \zeta, but also vector field perturbations. The purpose is twofold: first, we would like the diagrammatic approach (which we would call the Feynman-like rules) to become widespread among the cosmology community; second, we intend to give an easy tool to formulate any correlator of \zeta for those cases that involve vector field perturbations and that, therefore, may generate prolonged stages of anisotropic expansion and/or important levels of statistical anisotropy. Indeed, the usual way of formulating such correlators, using the Wick's theorem, may become very clutter and time-consuming.

We use the \delta N formalism to study the trispectrum T_\zeta of the primordial curvature perturbation \zeta when the latter is generated by vector field perturbations, considering the tree-level and one-loop contributions. The order of magnitude of the level of non-gaussianity in the trispectrum, \tau_{NL}, is calculated in this scenario and related to the order of magnitude of the level of non-gaussianity in the bispectrum, f_{NL}, and the level of statistical anisotropy in the power spectrum, g_\zeta. Such consistency relations will put under test this scenario against future observations. Comparison with the expected observational bound on \tau_{NL} from WMAP, for generic inflationary models, is done.

We calculate the trispectrum T_\zeta of the primordial curvature perturbation \zeta, generated during a {\it slow-roll} inflationary epoch by considering a two-field quadratic model of inflation with {\it canonical} kinetic terms. We consider loop contributions as well as tree level terms, and show that it is possible to attain very high, {\it including observable}, values for the level of non-gaussianity \tau_{NL} if T_\zeta is dominated by the one-loop contribution. Special attention is paid to the claim in JCAP {\bf 0902}, 017 (2009) [arXiv:0812.0807 [astro-ph]] that, in the model studied in this paper and for the specific inflationary trajectory we choose, the quantum fluctuations of the fields overwhelm the classical evolution. We argue that such a claim actually does not apply to our model, although more research is needed in order to understand the role of quantum diffusion. We also consider the probability that an observer in an ensemble of realizations of the density field sees a non-gaussian distribution. In that respect, we show that the probability associated to the chosen inflationary trajectory is non-negligible. Finally, the levels of non-gaussianity f_{NL} and \tau_{NL} in the bispectrum B_\zeta and trispectrum T_\zeta of \zeta, respectively, are also studied for the case in which \zeta is not generated during inflation.

We study the spectrum P_\zeta and bispectrum B_\zeta of the primordial curvature perturbation \zeta when the latter is generated by scalar and vector field perturbations. The tree-level and one-loop contributions from vector field perturbations are worked out considering the possibility that the one-loop contributions may be dominant over the tree level terms (both (either) in P_\zeta and (or) in B_\zeta) and viceversa. The level of non-gaussianity in the bispectrum, f_{NL}, is calculated and related to the level of statistical anisotropy in the power spectrum, g_\zeta. For very small amounts of statistical anisotropy in the power spectrum, the level of non-gaussianity may be very high, in some cases exceeding the current observational limit.

The \delta N formula for the primordial curvature perturbation \zeta is extended to include vector as well as scalar fields. Formulas for the tree-level contributions to the spectrum and bispectrum of \zeta are given, exhibiting statistical anisotropy. The one-loop contribution to the spectrum of \zeta is also worked out. We then consider the generation of vector field perturbations from the vacuum, including the longitudinal component that will be present if there is no gauge invariance. Finally, the \delta N formula is applied to the vector curvaton and vector inflation models with the tensor perturbation also evaluated in the latter case.

A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in furthering our understanding of two fundamental aspects of cosmology and astrophysics: (i) the physics of the very early universe that created the primordial seeds for large-scale structures, and (ii) the subsequent growth of structures via gravitational instability and gas physics at later times. To date, observations of fluctuations in the Cosmic Microwave Background (CMB) and the Large-Scale Structure of the Universe (LSS) have focused largely on the Gaussian contribution as measured by the two-point correlations (or the power spectrum) of density fluctuations. However, an even greater amount of information is contained in non-Gaussianity and a large discovery space therefore still remains to be explored. Many observational probes can be used to measure non-Gaussianity, including CMB, LSS, gravitational lensing, Lyman-alpha forest, 21-cm fluctuations, and the abundance of rare objects such as clusters of galaxies and high-redshift galaxies. Not only does the study of non-Gaussianity maximize the science return from a plethora of present and future cosmological experiments and observations, but it also carries great potential for important discoveries in the coming decade.

We show in this paper that it is possible to attain very high, {\it including observable}, values for the level of non-gaussianity f_{NL} associated with the bispectrum B_\zeta of the primordial curvature perturbation \zeta, in a subclass of small-field {\it slow-roll} models of inflation with canonical kinetic terms. Such a result is obtained by taking care of loop corrections both in the spectrum P_\zeta and the bispectrum B_\zeta. Sizeable values for f_{NL} arise even if \zeta is generated during inflation. Five issues are considered when constraining the available parameter space: 1. we must ensure that we are in a perturbative regime so that the \zeta series expansion, and its truncation, are valid. 2. we must apply the correct condition for the (possible) loop dominance in B_\zeta and/or P_\zeta. 3. we must satisfy the spectrum normalisation condition. 4. we must satisfy the spectral tilt constraint. 5. we must have enough inflation to solve the horizon problem.

In this paper we calculate additional contributions to that part of the non-Gaussianity of the curvature perturbation, which come from the three-point correlator of the field perturbations. We estimate this contribution in the following models for its origin: single-component inflation, multi-component chaotic inflation, a two-component "hybrid" inflationary model, and the curvaton scenario. In all of these models, the additional contributions to the primordial non-gaussianity considered here, are too small to be ever detected.

We review some theoretical and statistical aspects of the origin of the large-scale structure in the Universe, in view of the two most widely known and accepted scenarios: the inflaton and the curvaton scenarios. Among the theoretical aspects, we point out the impossibility of having a low inflationary energy scale in the simplest curvaton model. A couple of modifications to the simplest setup are explored, corresponding to the implementation of a second (thermal) inflationary period whose end makes the curvaton field `heavy', triggering either its oscillations or immediate decay. Low scale inflation is then possible to attain with H_\ast being as low as 1 TeV. Among the statistical aspects, we study the bispectrum B_\zeta(k_1,k_2,k_3) of the primordial curvature perturbation \zeta whose normalisation \fnl gives information about the level of non-gaussianity in \zeta. In connection with \fnl, several conserved and/or gauge invariant quantities described as the second-order curvature perturbation have been given in the literature. We review each of these quantities showing how to interpret one in terms of the others, and analyze the respective expected \fnl in both the inflaton and the curvaton scenarios as well as in other less known models for the generation of primordial perturbations and/or non-gaussianities. The \delta N formalism turns out to be a powerful technique to compute \fnl in multi-component slow-roll inflation, as the knowledge of the evolution of some family of unperturbed universes is the only requirement. We present for the first time this formalism and apply it to selected examples.

We extend the \delta N formalism so that it gives all of the stochastic properties of the primordial curvature perturbation \zeta if the initial field perturbations are gaussian. The calculation requires only the knowledge of some family of unperturbed universes. A formula is given for the normalisation \fnl of the bispectrum of \zeta, which is the main signal of non-gaussianity. Examples of the use of the formula are given, and its relation to cosmological perturbation theory is explained.

The end of a thermal inflation era, driven by the rolling of a flaton field coupled to the curvaton, cause a huge increment in the curvaton mass and decay rate while the curvaton is still frozen. It is shown that, if this increment is enough for the curvaton to immediately decay, low scale inflation with Hubble parameter H_\ast \sim 10^3 GeV is achieved for more natural values of the flaton-curvaton coupling constant \lambda (10^{-10} \lsim \lambda \lsim 10^{-4}) and the curvaton bare mass m_\sigma (m_\sigma \lsim 1 GeV).

Several conserved and/or gauge invariant quantities described as the second-order curvature perturbation have been given in the literature. We revisit various scenarios for the generation of second-order non-gaussianity in the primordial curvature perturbation \zeta, employing for the first time a unified notation and focusing on the normalisation f_{NL} of the bispectrum. When the classical curvature perturbation first appears a few Hubble times after horizon exit, |f_{NL}| is much less than 1 and is, therefore, negligible. Thereafter \zeta (and hence f_{NL}) is conserved as long as the pressure is a unique function of energy density (adiabatic pressure). Non-adiabatic pressure comes presumably only from the effect of fields, other than the one pointing along the inflationary trajectory, which are light during inflation (`light non-inflaton fields'). During single-component inflation f_{NL} is constant, but multi-component inflation might generate |f_{NL}| \sim 1 or bigger. Preheating can affect f_{NL} only in atypical scenarios where it involves light non-inflaton fields. The curvaton scenario typically gives f_{NL} \ll -1 or f_{NL} = +5/4. The inhomogeneous reheating scenario can give a wide range of values for f_{NL}. Unless there is a detection, observation can eventually provide a limit |f_{NL}| \lsim 1, at which level it will be crucial to calculate the precise observational limit using second order theory.

The primordial curvature perturbation may be due to a `curvaton' field, which dominates (or almost dominates) the energy density before it decays. In the simplest version of the curvaton model the scale of inflation has to be quite high corresponding to a Hubble parameter H>10^7 GeV. We here explore two modifications of the curvaton model which can instead allow inflation at a low scale. (i) The curvaton is a Pseudo Nambu-Goldstone Boson (PNGB), with a symmetry-breaking phase transition during inflation. (ii) The curvaton mass increases suddenly at some moment after the end of inflation but before the onset of the curvaton oscillations. Both proposals can work but not in a completely natural way. Also, the lower bound on the scale of inflation depends somewhat on the details of the framework used. Nevertheless, we show that inflation with H as low as 1 TeV or lower is possible to be attained.

We show how the mechanism recently proposed by Alexander, Peskin, and Sheikh-Jabbari, in which the observed cosmic matter-antimatter asymmetry comes from the contribution of the cosmological tensor perturbations to the gravitational anomaly in the Standard Model, does not work. This holds when taking into account the commutation relations for the creation-annihilation operators associated to the gravitational waves and the geometrical properties of the polarization tensors.