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.

In this work we revisit the topological mass generation of 2-forms and establish a connection to the unique derivative coupling arising in the quartic Lagrangian of the systematic construction of massive $2-$form interactions, relating in this way BF theories to Galileon-like theories of 2-forms. In terms of a massless $1-$form $A$ and a massless $2-$form $B$, the topological term manifests itself as the interaction $B\wedge F$, where $F = {\rm d} A$ is the field strength of the $1-$form. Such an interaction leads to a mechanism of generation of mass, usually referred to as "topological generation of mass" in which the single degree of freedom propagated by the $2-$form is absorbed by the $1-$form, generating a massive mode for the $1-$form. Using the systematical construction in terms of the Levi-Civita tensor, it was shown that, apart from the quadratic and quartic Lagrangians, Galileon-like derivative self-interactions for the massive 2-form do not exist. A unique quartic Lagrangian $\epsilon^{\mu\nu\rho\sigma}\epsilon^{\alpha\beta\gamma}_{\;\;\;\;\;\;\sigma}\partial_{\mu}B_{\alpha\rho}\partial_{\nu}B_{\beta\gamma}$ arises in this construction in a way that it corresponds to a total derivative on its own but ceases to be so once an overall general function is introduced. We show that it exactly corresponds to the same interaction of topological mass generation. Based on the decoupling limit analysis of the interactions, we bring out supporting arguments for the uniqueness of such a topological mass term and absence of the Galileon-like interactions. Finally, we discuss some preliminary applications in cosmology.

In this paper we consider a model based on interacting $p-$forms and explore some cosmological applications. Restricting to gauge invariant actions, we build a general Lagrangian allowing for arbitrary interactions between the $p-$forms (including interactions with a $0-$form, scalar field) in a given background in $D$ dimensions. For simplicity, we restrict the construction to up to first order derivatives of the fields in the Lagrangian. We discuss with detail the four dimensional case and devote some attention to the mechanism of topological mass generation originated by couplings of the form $B\wedge F$ between a $p-$form and a $(3-p)-$form. As a result, we show the system of the interacting $p-$forms $(p=1,2,3)$ is equivalent to a parity violating, massive, Proca vector field model. Finally, we discuss some cosmological applications. In a first case we study a very minimalistic system composed by a $3-$form coupled to a $0-$form. The $3-$form induces an effective potential which acts as a cosmological constant term suitable to drive the late time accelerated expansion of the universe dominated by dark energy. We study the dynamics of the system and determine its critical points and stability. Additionally, we study a system composed by a scalar field and a $1$-form. This case is interesting because the presence of a coupled $1-$form can generate non vanishing anisotropic signatures during the late time accelerated expansion. We discuss the evolution of cosmological parameters such as the equation of state in this model.

We use conformal symmetry to constrain the shape of inflationary correlators in the presence of long-lived vector field perturbations. Applying conformal Ward identities, we derive general expressions, up to amplitudes and normalization factors, for the two and three point correlators in the presence of vector fields mediated by the interaction $f(\phi)\left(F_{\mu \nu}F^{\mu \nu}+\alpha\tilde{F}_{\mu \nu}F^{\mu \nu}\right)$, where $f(\phi)$ is a suitable coupling function between the scalar and the vector field. The previous interaction allows for isotropy and parity symmetry breaking and is consistent with super horizon conformal symmetry. As an application of the conformal field theory techniques followed here, we evaluate the mixed tensor-scalar $\langle \gamma \zeta \rangle$ and tensor-scalar-scalar $\langle \gamma \zeta \zeta \rangle$ correlators which are interesting to look for parity violating effects related with chiral gravitational waves. Finally, we derive consistency relations for the three point correlators obtained.

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 cosmology in the presence of arbitrary couplings between $p$-forms in 4-dimensional space-time for a general action respecting gauge symmetry and parity invariance. The interaction between 0-form (scalar field $\phi$) and 3-form fields gives rise to an effective potential $V_{\rm eff}(\phi)$ for the former after integrating out the contribution of the latter. We explore the dynamics of inflation on an anisotropic cosmological background for a coupled system of 0-, 1-, and 2-forms. In the absence of interactions between 1- and 2-forms, we derive conditions under which the anisotropic shear endowed with nearly constant energy densities of 1- and 2-forms survives during slow-roll inflation for an arbitrary scalar potential $V_{\rm eff}(\phi)$. If 1- and 2-forms are coupled to each other, we show the existence of a new class of anisotropic inflationary solutions in which the energy density of 2-form is sustained by that of 1-form through their interactions. Our general analytic formulas for the anisotropic shear are also confirmed by the numerical analysis for a concrete inflaton potential.

We study the dynamics of dark energy in the presence of a 2-form field coupled to a canonical scalar field $\phi$. We consider the coupling proportional to $e^{-\mu \phi/M_{\rm pl}} H_{\alpha \beta \gamma}H^{\alpha \beta \gamma}$ and the scalar potential $V(\phi) \propto e^{-\lambda \phi/M_{\rm pl}}$, where $H_{\alpha \beta \gamma}$ is the 2-form field strength, $\mu, \lambda$ are constants, and $M_{\rm pl}$ is the reduced Planck mass. We show the existence of an anisotropic matter-dominated scaling solution followed by a stable accelerated fixed point with a non-vanishing shear. Even if $\lambda \geq {\cal O}(1)$, it is possible to realize the dark energy equation of state $w_{\rm DE}$ close to $-1$ at low redshifts for $\mu \gg \lambda$. The existence of anisotropic hair and the oscillating behavior of $w_{\rm DE}$ are key features for distinguishing our scenario from other dark energy models like quintessence.

In this paper we use conformal field theory techniques to constrain the form of the correlations functions of an inflationary scalar-vector model described by the interaction term $f_1(\phi)F_{\mu \nu}F^{\mu \nu} + f_2(\phi)\tilde{F}_{\mu \nu}F^{\mu \nu}$. We use the fact that the conformal group is the relevant symmetry group acting on super horizon scales in an inflationary de Sitter background. As a result, we find that super horizon conformal symmetry, constraints the form of the coupling functions $f_1, f_2$ to be homogeneous functions of the same degree. We derive the general form of the correlators involving scalar and vector perturbations in this model and determine its squeezed limit scaling behaviour for super horizon scales. The approach followed here is useful to constraint the shape of scalar-vector correlators, and our results agree with recent literature on the subject, but don't allow us to determine amplitude factors of the correlators.

Inflationary models including vector fields have attracted a great deal of attention over the past decade. Such an interest owes to the fact that they might contribute to, or even be fully responsible for, the curvature perturbation imprinted in the CMB. However, the necessary breaking of the vector field's conformal invariance during inflation is not without problems. In recent years it has been realized that a number of instabilities endangering the consistency of the theory arise when the conformal invariance is broken by means of a non-minimal coupling to gravity. In this paper we consider a massive vector field non-minimally coupled to gravity through the Gauss-Bonnet invariant, and investigate whether the vector can obtain a nearly scale-invariant perturbation spectrum while evading the emergence of perturbative instabilities. We find that the strength of the coupling must be extremely small if the vector field is to have a chance to contribute to the total curvature perturbation.

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.

A full-relativistic approach is used to compute the radius of the innermost stable circular orbit (ISCO), the Keplerian, frame-dragging, precession and oscillation frequencies of the radial and vertical motions of neutral test particles orbiting the equatorial plane of a magnetized neutron star. The space-time around the star is modelled by the six parametric solution derived by Pachon et al. It is shown that the inclusion of an intense magnetic field, such as the one of a neutron star, have non-negligible effects on the above physical quantities, and therefore, its inclusion is necessary in order to obtain a more accurate and realistic description of the physical processes occurring in the neighbourhood of this kind of objects such as the dynamics of accretion disk. The results discussed here also suggest that the consideration of strong magnetic fields may introduce non-negligible corrections in, e.g., the relativistic precession model and therefore on the predictions made on the mass of neutron stars.

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.

Whether analytic exact vacuum(electrovacuum) solutions of the Einstein(Einstein-Maxwell) field equations can accurately describe or not the exterior spacetime of compact stars remains still an interesting open question in Relativistic Astrophysics. As an attempt to establish their level of accuracy, the radii of the Innermost Stable Circular Orbits (ISCOs) of test particles given by analytic exterior spacetime geometries have been compared with the ones given by numerical solutions for neutron stars (NSs) obeying a realistic equation of state (EoS). It has been so shown that the six-parametric solution of Pach\'on, Rueda, and Sanabria (2006) (hereafter PRS) is more accurate to describe the NS ISCO radii than other analytic models. We propose here an additional test of accuracy for analytic exterior geometries based on the comparison of orbital frequencies of neutral test particles. We compute the Keplerian, frame-dragging, as well as the precession and oscillation frequencies of the radial and vertical motions of neutral test particles for the Kerr and PRS geometries; then we compare them with the numerical values obtained by Morsink and Stella (1999) for realistic NSs. We identify the role of high-order multipole moments such as the mass quadrupole and current octupole in the determination of the orbital frequencies especially in the rapid rotation regime. The results of this work are relevant to cast a separatrix between black hole (BH) and NS signatures as well as probe the nuclear matter EoS and NS parameters from the Quasi-Periodic Oscillations (QPOs) observed in Low Mass X-Ray Binaries.

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.

The perturbation of a light field might affect preheating and hence generate a contribution to the spectrum and non-gaussianity of the curvature perturbation \zeta. The field might appear directly in the preheating model (curvaton-type preheating) or indirectly through its effect on a mass or coupling (modulated preheating). We give general expressions for \zeta based on the \delta N formula, and apply them to the cases of quadratic and quartic chaotic inflation. For the quadratic case, curvaton-type preheating is ineffective in contributing to \zeta, but modulated preheating can be effective. For quartic inflation, curvaton-type preheating may be effective but the usual \delta N formalism has to be modified. We see under what circumstances the recent numerical simulation of Bond et al. [0903.3407] may be enough to provide a rough estimate for this case. This paper is dedicated to the memory of Lev Kofman who died on 12th November 2009

We study the statistical descriptors for some cosmological inflationary models that allow us to get large levels of non-gaussianity and violations of statistical isotropy. Basically, we study two different class of models: a model that include only scalar field perturbations, specifically a subclass of small-field slow-roll models of inflation with canonical kinetic terms, and models that admit both vector and scalar field perturbations. We study the former to show that it is possible to attain very high, including observable, values for the levels of non-gaussianity f_{NL} and \tao_{NL} in the bispectrum B_\zeta and trispectrum T_\zeta of the primordial curvature perturbation \zeta respectively. Such a result is obtained by taking care of loop corrections in the spectrum P_\zeta, the bispectrum B_\zeta and the trispectrum T_\zeta . Sizeable values for f_{NL} and \tao_{NL} arise even if \zeta is generated during inflation. For the latter we study the spectrum P_\zeta, bispectrum B_\zeta and trispectrum $T_\zeta of the primordial curvature perturbation when \zeta 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. The levels of non-gaussianity f_{NL} and \tao_{NL}, are 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 levels of non-gaussianity may be very high, in some cases exceeding the current observational limit.

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.

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.