Motivated by a scenario of magnetogenesis in which a homogeneous magnetic field is generated during inflation, we study the magnetohydrodynamic evolution of the primordial plasma motions for two kinds of initial conditions -- (i) a spatially homogeneous imposed field with an unlimited correlation length, and (ii) a zero flux scale-invariant statistically homogeneous magnetic field. In both cases, we apply, for a short initial time interval, monochromatic forcing at a certain wave number so that the correlation length is finite, but much smaller than the typical length scale of turbulence. In particular, we investigate the decay of non-helical and helical hydromagnetic turbulence. We show that, in the presence of a homogeneous magnetic field, the decay of helical and nonhelical small-scale fields can occur rapidly. This is a special property of a system with a perfectly homogeneous magnetic field, which is sometimes considered as a local approximation to a slowly varying background field. It can never change and acts as an imposed magnetic field. This is in a sharp contrast to the case of a statistically homogeneous magnetic field, where we recover familiar decay properties: a much slower decay of magnetic energy and a faster growth of the correlation length, especially in the case with magnetic helicity. The result suggests that a homogeneous magnetic field, if generated during inflation, should persist under the influence of small-scale fields and could be the origin of the large-scale magnetic field in the universe.

The idea of gauge theories of gravity predicts that there should exist not only the massless graviton but also massive particles carrying the gravitational force. We study the cosmology in a quadratic gravity with dynamical torsion where gravity may be interpreted as a gauge force associated with the Poincar\'{e} group. In addition to the massless spin-2 graviton, the model contains four non-ghost massive particle species: a couple of spin-0, a spin-1 and a spin-2. Supposing the restoration of the local Weyl invariance in the UV limit and the parity invariance, we find the most general minisuperspace action describing a homogeneous and isotropic universe with a flat spatial geometry. We then transform the minisuperspace action to a quasi-Einstein frame in which the field space is a hyperboloid and the field potential is a combination of those of a Starobinsky-like inflation and a natural inflation. Remarkably, thanks to the multi-field dynamics, the Starobinsky-like inflationary trajectory can be realized even if the initial condition is away from the top of the Starobinsky-like potential. We also study linear tensor perturbations and find qualitatively different features than the Starobinsky inflation, spontaneous parity violation and mixing of the massless and massive spin-2 modes, which might reveal the underlying nature of gravity through inflationary observables.

We study the phenomenology of a class of minimally modified gravity theories called $f(\mathcal{H})$ theories, in which the usual general relativistic Hamiltonian constraint is replaced by a free function of it. After reviewing the construction of the theory and a consistent matter coupling, we analyze the dynamics of cosmology at the levels of both background and perturbations, and present a concrete example of the theory with a $3$-parameter family of the function $f$. Finally, we compare this example model to Planck data as well as some later-time probes, showing that such a realization of $f(\mathcal{H})$ theories fits the data significantly better than the standard $\Lambda$CDM model, in particular by modifying gravity at intermediate redshifts, $z\simeq743$.

In a very recent paper [1], we have proposed a novel $4$-dimensional gravitational theory with two dynamical degrees of freedom, which serves as a consistent realization of $D\to4$ Einstein-Gauss-Bonnet gravity with the rescaled Gauss-Bonnet coupling constant $\tilde{\alpha}$. In the present paper, we study cosmological implications of the theory in the presence of a perfect fluid and clarify the similarities and differences between the results obtained from the consistent $4$-dimensional theory and those from the previously considered, naive (and inconsistent) $D\rightarrow 4$ limit. Studying the linear perturbations, we explicitly show that the theory only has tensorial gravitational degrees of freedom (besides the matter degree) and that for $\tilde{\alpha}>0$ and $\dot{H}<0$, perturbations are free of any pathologies so that we can implement the setup to construct early and/or late time cosmological models. Interestingly, a $k^4$ term appears in the dispersion relation of tensor modes which plays significant roles at small scales and makes the theory different than not only general relativity but also many other modified gravity theories as well as the naive (and inconsistent) $D\to 4$ limit. Taking into account the $k^4$ term, the observational constraint on the propagation of gravitational waves yields the bound $\tilde{\alpha} \lesssim {\cal O}(1)\,{\rm eV}^{-2}$. This is the first bound on the only parameter (besides the Newton's constant and the choice of a constraint that stems from a temporal gauge fixing) in the consistent theory of $D\to 4$ Einstein-Gauss-Bonnet gravity.

We propose a modified gravity theory that propagates only two local gravitational degrees of freedom and that does not have an Einstein frame. According to the classification in JCAP 01 (2019) 017 [arXiv:1810.01047 [gr-qc]], this is a type-II minimally modified gravity theory. The theory is characterized by the gravitational constant $G_{\rm N}$ and a function $V(\phi)$ of a non-dynamical auxiliary field $\phi$ that plays the role of dark energy. Once one fixes a homogeneous and isotropic cosmological background, the form of $V(\phi)$ is determined and the theory no longer possesses a free parameter or a free function, besides $G_{\rm N}$. For $V'(\phi) = 0$ the theory reduces to general relativity (GR) with $G_N$ being the Newton's constant and $V=const.$ being the cosmological constant. For $V'(\phi) \ne 0$, it is shown that gravity behaves differently from GR but that GR with $G_{\rm N}$ being the Newton's constant is recovered for weak gravity at distance and time scales sufficiently shorter than the scale associated with $V(\phi)$. Therefore this theory provides the simplest framework of cosmology in which deviations from GR can be tested by observational data.

The constant-roll inflation with small positive value of the constant-roll parameter $\beta\equiv \frac{\ddot\phi}{H\dot\phi}={\rm const.}$ has been known to produce a slightly red-tilted curvature power spectrum compatible with the current observational constraints. In this work, we shed light on the constant-roll inflation with negative $\beta$ and investigate how a stage of constant-roll inflation may realize the growth in the primordial curvature power spectrum necessary to produce a peaked spectrum of primordial black hole abundance. We first review the behavior of constant-roll models in the range of parameters $-\frac32<\beta<0$, which allows for a constant-roll attractor stage generating a blue-tilted curvature power spectrum without superhorizon growth. As a concrete realization, we consider a potential with two slow-roll stages, separated by the constant-roll stage, in a way that satisfies the current constraints on the power spectrum and the primordial black hole abundance. The model can produce primordial black holes as all dark matter, LIGO-Virgo events, or OGLE microlensing events. Due to the range of different scalar tilts allowed by the constant-roll potential, this construction is particularly robust and testable by future observations.

We study canonical transformations of general relativity (GR) to provide a novel matter coupling to gravity. Although the transformed theory is equivalent to GR in vacuum, the equivalence no longer holds if a matter field minimally couples to the canonically transformed gravitational field. We find that a naive matter coupling to the transformed field leads to the appearance of an extra mode in the phase space, rendering the theory inconsistent. We then find a consistent and novel way of matter coupling: after imposing a gauge fixing condition, a matter field can minimally couple to gravity without generating an unwanted extra mode. As a result, the way matter field couples to the gravitational field determines the preferred time direction and the resultant theory has only two gravitational degrees of freedom. We also discuss the cosmological solution and linear perturbations around it, and confirm that their dynamics indeed differ from those in GR. The novel matter coupling can be used for a new framework of modified gravity theories.

The Horndeski gauge-gravity coupling is the leading non-minimal interaction between gravity and gauge bosons, and it preserves all the symmetries and the number of physical degrees of freedom of the standard model of particle physics and general relativity. In this paper we study the effects of the non-minimal interaction in astronomy and cosmology, and obtain upper bounds on the associated dimensionless coupling constant $\lambda$. From the modification of equations of motion of gauge bosons applied to compact astronomical objects, we find upper bounds $|\lambda| \lesssim 10^{88}$, $|\lambda| \lesssim 10^{70}$ and $|\lambda| \lesssim 10^{81}$ from a black hole shadow, neutron stars and white dwarfs, respectively. The bound $|\lambda| \lesssim 10^{70}$ that is deduced from neutron stars is the strongest and provides twenty orders of magnitude improvement of the previously known best bound on this parameter. On the other hand, the effects of this term on modification of the gravitational Poisson equation lead to a weaker bound $|\lambda| \lesssim 10^{98}$. From the propagation of gravitational waves we also find $|\lambda| \lesssim 10^{119}$, which is even weaker.

In scalar-tensor theories we revisit the issue of strong coupling of perturbations around stealth solutions, i.e.\ backgrounds with the same forms of the metric as in General Relativity but with non-trivial configurations of the scalar field. The simplest among them is a stealth Minkowski (or de Sitter) solution with a constant, timelike derivative of the scalar field, i.e.\ ghost condensation. In the decoupling limit the effective field theory (EFT) describing perturbations around the stealth Minkowski (or de Sitter) solution shows the universal dispersion relation of the form $\omega^2 = \alpha k^4/M^2$, where $M$ is a mass scale characterizing the background scalar field and $\alpha$ is a dimensionless constant. Provided that $\alpha$ is positive and of order unity, a simple scaling argument shows that the EFT is weakly coupled all the way up to $M$. On the other hand, if the structure of the underlining theory forces the perturbations to follow second-order equations of motion then $\alpha=0$ and the dispersion relation loses dependence on the spatial momentum. This not only explains the origin of the strong coupling problem that was recently pointed out in a class of degenerate theories but also provides a hint for a possible solution of the problem. We then argue that a controlled detuning of the degeneracy condition, which we call scordatura, renders the perturbations weakly coupled without changing the properties of the stealth solutions of degenerate theories at astrophysical scales.

We provide two derivations of the baryonic equations that can be straightforwardly implemented in existing Einstein--Boltzmann solvers. One of the derivations begins with an action principle, while the other exploits the conservation of the stress-energy tensor. While our result is manifestly covariant and satisfies the Bianchi identities, we point out that this is not the case for the implementation of the seminal work by Ma and Bertschinger and in the existing Boltzmann codes. We also study the tight coupling approximation up to the second order without choosing any gauge using the covariant full baryon equations. We implement the improved baryon equations in a Boltzmann code and investigate the change in the estimate of cosmological parameters by performing an MCMC analysis. With the covariantly correct baryon equations of motion, we find 1% deviation for the best fit values of the cosmological parameters that should be taken into account. While in this paper, we study the Lambda-CDM model only, our baryon equations can be easily implemented in other models and various modified gravity theories.

We study a model of two scalar fields with a hyperbolic field space and show that it reduces to a single-field Dirac-Born-Infeld (DBI) model in the limit where the field space becomes infinitely curved. We apply the de Sitter swampland conjecture to the two-field model and take the same limit. It is shown that in the limit, all quantities appearing in the swampland conjecture remain well-defined within the single-field DBI model. Based on a consistency argument, we then speculate that the condition derived in this way can be considered as the de Sitter swampland conjecture for a DBI scalar field by its own. The condition differs from those proposed in the literature and only the one in the present paper passes the consistency argument. As a byproduct, we also point out that one of the inequalities in the swampland conjecture for a multi-field model with linear kinetic terms should involve the lowest mass squared for scalar perturbations and that this quantity can be significantly different from the lowest eigenvalue of the Hessian of the potential in the local orthonormal frame if the field space is highly curved. Finally, we propose an extension of the de Sitter swampland conjecture to a more general scalar field with the Lagrangian of the form $P(X,\varphi)$, where $X=-(\partial\varphi)^2/2$.

We study the constraint on the inflationary energy scale from the recently proposed Trans-Planckian Censorship Conjecture (TCC). We find a universal upper bound on the inflationary Hubble expansion rate $H_\text{inf}$ which is solely determined by the reheating temperature $T_\text{rh}$: $H_\text{inf}/M_\text{Pl}\lesssim T_0/T_\text{rh}$, where $T_0$ is the photon temperature today. The upper limit can be saturated by a post-inflationary oscillatory stage with the critical equation-of-state parameter $w\approx-1/3$, or by an inflation model with multiple stages. In the lowest reheating temperature required for big bang nucleosynthesis, the upper bound on the tensor-to-scalar ratio $r$ at the CMB scales is $r\lesssim10^{-8}$, which can be realized in many string-inspired inflation models.

We calculate tensor bispectrum in a theory of massive tensor gravitons which predicts blue-tilted and largely amplified primordial gravitational waves. We find that a new 3-point interaction can produce a larger tensor bispectrum than the conventional one from general relativity. The bispectrum peaks at the squeezed configuration and its slope towards the squeezed limit is determined by the graviton mass. This squeezed tensor bispectrum may be observed as the quadrupolar modulation of the tensor power spectrum by interferometers.

String theory suggests a unique and unambiguous modification to General Relativity: the symmetry of $\mathbf{O}(D,D)$ T-duality promotes the entire closed-string massless NS-NS sector to stringy graviton fields. The symmetry fixes the couplings to other matter fields unambiguously and the Einstein field equations are enriched to comprise $D^{2}+1$ components, dubbed recently as the Einstein Double Field Equations. Here we explore the cosmological implications of this `Stringy Gravity'. We derive the most general homogeneous and isotropic ansatzes for both stringy graviton fields and the stringy energy-momentum tensor. Substituting them into the Einstein Double Field Equations, we obtain the $\mathbf{O}(D,D)$ completion of the Friedmann equations along with a generalized continuity equation. We discuss how this gives an enriched and novel framework beyond typical string cosmology, with solutions that may be characterized by two equation-of-state parameters, $w$ (conventional) and $\lambda$ (new). When $\lambda+3w=1$, the dilaton remains constant throughout the cosmological evolution, and one recovers the standard Friedmann equations for generic matter content (i.e. for any $w$), an improvement over conventional string cosmology where this occurs only for a radiation equation of state ($w=1/3$). We further point out that, in contrast to General Relativity, in Stringy Gravity there is no de Sitter solution arising from either an $\mathbf{O}(D,D)$-symmetric cosmological constant or scalar field with positive energy density.

We consider a new form of theories of gravity in which the action is written in terms of the Ricci scalar and its first and second derivatives. Despite the higher derivative nature of the action, the theory is free from ghost under an appropriate choice of the functional form of the Lagrangian. This model possesses $2+2$ physical degrees of freedom, namely $2$ scalar degrees and $2$ tensor degrees. We exhaust all such theories with the Lagrangian of the form $f(R, (\nabla R)^2, \Box R)$, where $R$ is the Ricci scalar, and then show some examples beyond this ansatz. In course of analysis, we prove the equivalence between these examples and generalized bi-Galileon theories.

It is well known that the standard scalar field mimetic cosmology provides a dark matter-like energy density component. Considering $SU(2)$ gauge symmetry, we study the gauge field extension of the mimetic scenario in spatially flat and curved FLRW spacetimes. Because of the mimetic constraint, the standard Yang-Mills term plays the role of the cosmological constant while the mimetic term provides two different contributions: one is the standard radiation scaling like $ a^{-4}$ while the other contribution in energy density scales as $\propto a^{-2}$. Consequently, in the Friedmann equation we have two different energy densities which scale as $\propto a^{-2}$: one is the mimetic spatial curvature-like and the other is the standard spatial curvature which can compete with each other. The degeneracy between these two contributions are disentangled in this scenario since the mimetic spatial curvature-like term shows up only at the dynamical level while the standard spatial curvature term shows up at both dynamical and kinematical levels.

Spontaneous scalarization is an interesting mechanism for modification of gravity by nonminimal coupling of a scalar field to matter or curvature invariants in the context of scalar-tensor theories, and its onset is signaled by linear instability of the scalar field around the corresponding general relativity solution. We thus perform the linear stability analysis of the scalar field about general relativity solutions and highlight a crucial difference between a spherically symmetric profile and a planar symmetric profile. We clarify that the critical value for the instability is sensitive to the morphology and that the spontaneous scalarization occurs much more easily with the planar symmetric shape than with the spherically symmetric shape.

We study cosmology in a class of minimally modified gravity (MMG) with two local gravitational degrees of freedom. We classify modified gravity theories into type-I and type-II: theories of type-I have an Einstein frame and can be recast by change of variables as general relativity (GR) with a non-minimal matter coupling, while theories of type-II have no Einstein frame. Considering a canonical transformation of the lapse, the 3-dimensional induced metric and their conjugate momenta we generate type-I MMG. We then show that phenomenological deviations from GR, such as the speed of gravitational waves $c_T$ and the effective gravitational constant for scalar perturbations $G_{\rm eff}$, are characterized by two functions of an auxiliary variable. We study the phenomenology of several models all having $c_T=1$. We obtain a scenario with $c_T=1$ in which the effective equation-of-state parameter of dark energy is different from $-1$ even though the cosmic acceleration is caused by a bare cosmological constant, and we find that it is possible to reconstruct the theory on choosing a selected time-evolution for the effective dark energy component.

We numerically study the gravitational collapse of a massless scalar field with spherical symmetry in Einstein-aether theory, and show that apparent, spin-0 and dynamical universal horizons (dUHs) can be all formed. The spacetime and the aether field are well-behaved and regular, including regions nearby these horizons (but away from the center of spherical symmetry). The spacetime outside the apparent and spin-0 horizons settles down to a static configuration, and some of such resulting static black holes were already found numerically in the literature. On the other hand, the proper distance of the outermost dUH from the apparent (or spin-0) horizon keeps increasing on $\mbox{\ae}$ther-orthogonal time slices. This indicates that the outermost dUH is evolving into the causal boundary, even for excitations with large speeds of propagation.

We study a theory of massive tensor gravitons which predicts blue-tilted and largely amplified primordial gravitational waves. After inflation, while their mass is significant until it diminishes to a small value, gravitons are diluted as non-relativistic matter and hence their amplitude can be substantially amplified compared to the massless gravitons which decay as radiation. We show that such gravitational waves can be detected by interferometer experiments, even if their signal is not observed on the CMB scales.

We extend the mimetic cosmology to models containing gauge invariant $p$-forms. The $0$-form case reproduces the well-known results of the mimetic dark matter, the $1$-form corresponds to the gauge field mimetic model while the $2$-form model is the Hodge dual of the $0$-form model in $4$ spacetime dimensions. We study the cosmological applications of the new gauge field mimetic model and show that it generates an energy density component which mimics the roles of spatial curvature. In the presence of the Maxwell term, the model also supports the flat, open and closed de Sitter-like cosmological backgrounds while the spatial geometry is flat for all three cases. We perform the cosmological perturbations analysis and show that the model is stable in the case of open de Sitter-like solution while it suffers from ghost instabilities in the case of the closed de Sitter-like solution.

A detailed stability analysis is presented for the de Sitter solution with a homogeneous magnetic field that was recently found in the context of a $U(1)$ gauge theory nonminimally coupled to scalar-tensor gravity. The magnetic field is `stealth' in the sense that the corresponding stress-energy tensor is of the form of an effective cosmological constant and thus is isotropic despite the fact that the magnetic field has a preferred spatial direction. We study the stability of the solution against linear perturbations in the subhorizon and superhorizon limits. We then present some explicit examples that satisfy all stability conditions. The stable de Sitter solution with a homogeneous magnetic field opens up a new possibility for inflationary magnetogenesis, in which magnetic fields in the Universe at all scales may originate from a classical, homogeneous magnetic field sustained during inflation.

Motivated by the GW170817/GRB170817A constraint on the deviation of the speed of gravitational waves from that of photons, we study disformal transformations of the metric in the context of the generalized Proca theory. The constraint restricts the form of the gravity Lagrangian, the way the electromagnetism couples to the gravity sector on cosmological backgrounds, or in general a combination of both. Since different ways of coupling matter to gravity are typically related to each other by disformal transformations, it is important to understand how the structure of the generalized Proca Lagrangian changes under disformal transformations. For disformal transformations with constant coefficients we provide the complete transformation rule of the Lagrangian. We find that additional terms, which were considered as beyond generalized Proca in the literature, are generated by the transformations. Once these additional terms are included, on the other hand, the structure of the gravity Lagrangian is preserved under the transformations. We then derive the transformation rules for the sound speeds of the scalar, vector and tensor perturbations on a homogeneous and isotropic background. We explicitly show that they transform following the natural expectation of metric transformations, that is, according to the transformation of the background lightcone structure. We end by arguing that inhomogeneities due to structures in the universe, e.g. dark matter halos, generically changes the speed of gravitational waves from its cosmological value. We give a rough estimate of the effect of inhomogeneities and find that the fine-tuning should not depend on the background or that the fine-tuned theory has to be further fine-tuned to pass the tight constraint.

The minimal theory of massive gravity (MTMG) has two branches of stable cosmological solutions: a self-accelerating branch, which, except for the mass of tensor modes has exactly the same behavior of linear perturbations as $\Lambda$CDM in general relativity (GR), and a normal branch with nontrivial behavior. We explore the influence of the integrated Sachs-Wolfe-galaxy correlation constraints on the normal branch of MTMG, which, in its simplest implementation, has one free parameter more than $\Lambda$CDM in GR (or the self-accelerating branch of MTMG): $\theta$. This parameter is related to the graviton mass and only affects the behavior of the cosmological linear perturbation dynamics. Using 2d-mass and SDSS data, we check which values of $\theta$ lead to a positive or negative cross-correlation. We find that positive cross-correlation is achieved for a large parameter-space interval. Within this allowed region of parameter space, we perform a $\chi^2$ analysis in terms of the parameter $\theta$, while keeping the other background parameters fixed to the best-fit values of Planck. We then infer that the normal branch of MTMG fits the data well in a nontrivial portion of the parameter space, and future experiments should be able to distinguish such a model from $\Lambda$CDM in GR (or the self-accelerating branch of MTMG).

In this paper, we carry out a systematic analysis of the theoretical and observational constraints on the dimensionless coupling constants $c_i$ ($i=1,2,3,4$) of the Einstein-aether theory, taking into account the events GW170817 and GRB 170817A. The combination of these events restricts the deviation of the speed $c_T$ of the spin-2 graviton to the range, $- 3\times 10^{-15} < c_T -1 < 7\times 10^{-16}$, which for the Einstein-aether theory implies $\left|c_{13}\right| \le 10^{-15}$ with $c_{ij} \equiv c_{i} + c_{j}$. The rest of the constraints are divided into two groups: those on the ($c_1, c_{14}$)-plane and those on the ($c_2, c_{14}$)-plane, except the strong-field constraints. The latter depend on the sensitivities $\sigma_{\ae}$ of neutron stars, which are not known at present in the new ranges of the parameters found in this paper.

We study the emergence of quantum entanglement in multi-field inflation. In this scenario, the perturbations of one field contribute to the observable curvature perturbation, while multi-field dynamics with the other fields affect the curvature perturbation through particle production and entanglement. We develop a general formalism which defines the quantum entanglement between the perturbations of the multiple fields both in the Heisenberg and Schr\"odinger pictures, and show that entanglement between different fields can arise dynamically in the context of multi-field inflationary scenarios. We also present a simple model in which a sudden change in the kinetic matrix of the scalar fields generates entanglement and an oscillatory feature appears in the power spectrum of the inflaton perturbation.

The minimal theory of quasidilaton massive gravity with or without a Horndeski-type kinetic term for the quasidilaton field propagates only three physical modes: the two massive tensor polarizations and one scalar mode. This reduced number of degrees of freedom is realized by a Lorentz symmetry violation at cosmological scales and the presence of appropriate constraints that remove unwanted modes. Vacuum cosmological solutions have been considered in a previous work, and it has been shown that the late-time de Sitter attractor is stable under inhomogeneous perturbations. In this work, we explore the stability of cosmological solutions in the presence of matter fields. Assuming for simplicity that the quasidilaton scalar is on an attractor at the level of the background, we derive stability conditions in the subhorizon limit, and find the scalar sound speeds, as well as the modification with respect to general relativity to the gravitational potential in the quasistatic approximation. We also find that the speed limit of gravitational waves coincides with the speed of light for any homogeneous and isotropic cosmological background, on or away from the attractor.

In Ho\v{r}ava-Lifshitz gravity a scaling isotropic in space but anisotropic in spacetime, often called anisotropic scaling with the dynamical critical exponent z=3, lies at the base of its renormalizability. This scaling also leads to a novel mechanism of generating scale-invariant cosmological perturbations, solving the horizon problem without inflation. In this paper we propose a possible solution to the flatness problem, in which we assume that the initial condition of the Universe is set by a small instanton respecting the same scaling. We argue that the mechanism may be more general than the concrete model presented here, and rely simply on the deformed dispersion relations of the theory, and on equipartition of the various forms of energy at the starting point.

The recently proposed chameleonic extension of bigravity theory, by including a scalar field dependence in the graviton potential, avoids several fine-tunings found to be necessary in usual massive bigravity. In particular it ensures that the Higuchi bound is satisfied at all scales, that no Vainshtein mechanism is needed to satisfy solar-system experiments, and that the strong coupling scale is always above the scale of cosmological interest all the way up to the early universe. This paper extends the previous work by presenting a stable example of cosmology in the chameleon bigravity model. We find a set of initial conditions and parameters such that the derived stability conditions on general flat Friedmann background are satisfied at all times. The evolution goes through radiation dominated, matter dominated, and de Sitter eras. We argue that the parameter space allowing for such a stable evolution may be large enough to encompass an observationally viable evolution. We also argue that our model satisfies all known constraints due to gravitational wave observations so far and thus can be considered as a unique testing ground of gravitational wave phenomenologies in bimetric theories of gravity.

This article extends bimetric formulations of massive gravity to make the mass of the graviton to depend on its environment. This minimal extension offers a novel way to reconcile massive gravity with local tests of general relativity without invoking the Vainshtein mechanism. On cosmological scales, it is argued that the model is stable and that it circumvents the Higuchi bound, hence relaxing the constraints on the parameter space. Moreover, with this extension the strong coupling scale is also environmentally dependent in such a way that it is kept sufficiently higher than the expansion rate all the way up to the very early universe, while the present graviton mass is low enough to be phenomenologically interesting. In this sense the extended bigravity theory serves as a partial UV completion of the standard bigravity theory. This extension is very generic and robust and a simple specific example is described.

We study disformal transformations of cosmological perturbations by vector fields in theories invariant under U(1) gauge transformations. Three types of vector disformal transformations are considered: (i) disformal transformations by a single timelike vector; (ii) disformal transformations by a single spacelike vector; and (iii) disformal transformations by three spacelike vectors. We show that transformations of type (i) do not change either curvature perturbation or gravitational waves; that those of type (ii) do not change curvature perturbation but change gravitational waves; and that those of type (iii) change both curvature perturbation and gravitational waves. Therefore, coupling matter fields to the metric after disformal transformations of type (ii) or (iii) in principle have observable consequences. While the recent multi-messenger observation of binary neutron stars has singled out a proper disformal frame at the present epoch with a high precision, the result of the present paper may thus help distinguishing disformal frames in the early universe.

We propose a scenario that can naturally explain the observed dark matter-baryon ratio in the context of bimetric theory with a chameleon field. We introduce two additional gravitational degrees of freedom, the massive graviton and the chameleon field, corresponding to dark matter and dark energy, respectively. The chameleon field is assumed to be non-minimally coupled to dark matter, i.e., the massive graviton, through the graviton mass terms. We find that the dark matter-baryon ratio is dynamically adjusted to the observed value due to the energy transfer by the chameleon field. As a result, the model can explain the observed dark matter-baryon ratio independently from the initial abundance of them.

We study primordial perturbations from hyperinflation, proposed recently and based on a hyperbolic field-space. In the previous work, it was shown that the field-space angular momentum supported by the negative curvature modifies the background dynamics and enhances fluctuations of the scalar fields qualitatively, assuming that the inflationary background is almost de Sitter. In this work, we confirm and extend the analysis based on the standard approach of cosmological perturbation in multi-field inflation. At the background level, to quantify the deviation from de Sitter, we introduce the slow-varying parameters and show that steep potentials, which usually can not drive inflation, can drive inflation. At the linear perturbation level, we obtain the power spectrum of primordial curvature perturbation and express the spectral tilt and running in terms of the slow-varying parameters. We show that hyperinflation with power-law type potentials has already been excluded by the recent Planck observations, while exponential-type potential with the exponent of order unity can be made consistent with observations as far as the power spectrum is concerned. We also argue that, in the context of a simple $D$-brane inflation, the hyperinflation requires exponentially large hyperbolic extra dimensions but that masses of Kaluza-Klein gravitons can be kept relatively heavy.

The minimal theory of quasidilaton massive gravity allows for a stable self-accelerating de Sitter solution in a wide range of parameters. On the other hand, in order for the theory to be compatible with local gravity tests, the fifth force due to the quasidilaton scalar needs to be screened at local scales. The present paper thus extends the theory by inclusion of a cubic Horndeski term in a way that (i) respects the quasidilaton global symmetry, that (ii) maintains the physical degrees of freedom in the theory being three, that (iii) can accommodate the Vainshtein screening mechanism and that still (iv) allows for a stable self-accelerating de Sitter solution. After adding the Horndeski term (and a k-essence type nonlinear kinetic term as well) to the precursor action, we switch to the Hamiltonian language and find a complete set of independent constraints. We then construct the minimal theory with three physical degrees of freedom by carefully adding a pair of constraints to the total Hamiltonian of the precursor theory. Switching back to the Lagrangian language, we study cosmological solutions and their stability in the minimal theory. In particular, we show that a self-accelerating de Sitter solution is stable for a wide range of parameters. Furthermore, as in the minimal theory of massive gravity, the propagation speed of the massive gravitational waves in the high momentum limit precisely agrees with the speed of light.

We introduce a quasidilaton scalar field to the minimal theory of massive gravity with the Minkowski fiducial metric, in such a way that the quasidilaton global symmetry is maintained and that the theory admits a stable self-accelerating de Sitter solution. We start with a precursor theory that contains three propagating gravitational degrees of freedom without a quasidilaton scalar and introduce St\"uckelberg fields to covariantize its action. This makes it possible for us to formulate the quasidilaton global symmetry that mixes the St\"uckelberg fields and the quasidilaton scalar field. By the Hamiltonian analysis we confirm that the precursor theory with the quasidilaton scalar contains four degrees of freedom, three from the precursor massive gravity and one from the quasidilaton scalar. We further remove one propagating degree of freedom to construct the minimal quasidilaton theory with three propagating degrees of freedom, corresponding to two polarizations of gravitational waves from the minimal theory of massive gravity and one scalar from the quasidilaton field, by carefully introducing two additional constraints to the system in the Hamiltonian language. Switching to the Lagrangian language, we find self-accelerating de Sitter solutions in the minimal quasidilaton theory and analyze their stability. It is found that the self-accelerating de Sitter solution is stable in a wide range of parameters.

We investigate the Hamiltonian structure of a class of gravitational theories whose actions are linear in the lapse function. We derive the necessary and sufficient condition for a theory in this class to have two or less local physical degrees of freedom. As an application we then find several concrete examples of modified gravity theories in which the total number of local physical degrees of freedom in the gravity sector is two.

We investigate the minimal theory of massive gravity (MTMG) recently introduced. After reviewing the original construction based on its Hamiltonian in the vielbein formalism, we reformulate it in terms of its Lagrangian in both the vielbein and the metric formalisms. It then becomes obvious that, unlike previous attempts in the literature of Lorentz-violating massive gravity, not only the potential but also the kinetic structure of the action is modified from the de Rham-Gabadadze-Tolley (dRGT) massive gravity theory. We confirm that the number of physical degrees of freedom in MTMG is two at fully nonlinear level. This proves the absence of various possible pathologies such as superluminality, acausality and strong coupling. Afterwards, we discuss the phenomenology of MTMG in the presence of a dust fluid. We find that on a flat homogeneous and isotropic background we have two branches. One of them (self-accelerating branch) naturally leads to acceleration without the genuine cosmological constant or dark energy. For this branch both the scalar and the vector modes behave exactly as in general relativity (GR). The phenomenology of this branch differs from GR in the tensor modes sector, as the tensor modes acquire a non-zero mass. Hence, MTMG serves as a stable nonlinear completion of the self-accelerating cosmological solution found originally in dRGT theory. The other branch (normal branch) has a dynamics which depends on the time-dependent fiducial metric. For the normal branch, the scalar mode sector, even though as in GR only one scalar mode is present (due to the dust fluid), differs from the one in GR, and, in general, structure formation will follow a different phenomenology. The tensor modes will be massive, whereas the vector modes, for both branches, will have the same phenomenology as in GR.

We present a novel cosmological solution in the framework of extended quasidilaton theory which underwent scrutiny recently. We only consider terms that do not generate the Boulware-Deser degree of freedom, hence the "ghost-free" quasidilaton theory, and show three new branches of cosmological evolution therein. One of the solutions passes the perturbative stability tests. This new solution exhibits a late time self-acceleration and all graviton polarizations acquire masses that converge to a constant in the asymptotic future. Moreover, all modes propagate at the speed of light. We propose that this solution can be used as a benchmark model for future phenomenological studies.

In the extended quasidilaton massive gravity we perform a nonlinear transformation of the shift vector and then calculate the second derivatives of the Hamiltonian density with respect to the lapse function and the (nonlinearly transformed) shift vector. It is then shown that the $4\times 4$ Hessian matrix is invertible, meaning that the equations of motion for the lapse function and the shift vector simply determine themselves. Therefore, there is no primary constraint that removes the Boulware-Deser ghost.

In this paper, we revisit the issue of static hairs of black holes in gravitational theories with broken Lorentz invariance in the case that the speed $c_{\phi}$ of the khronon field becomes infinitely large, $c_{\phi} = \infty$, for which the sound horizon of the khronon field coincides with the universal horizon, and the boundary conditions at the sound horizon reduce to those given normally at the universal horizons. As a result, less boundary conditions are present in this extreme case in comparison with the case $c_{\phi} = $ finite. Then, it would be expected that static hairs might exist. However, we show analytically that even in this case static hairs still cannot exist, based on a decoupling limit analysis. We also consider the cases in which $c_{\phi}$ is finite but with $c_{\phi} \gg 1$, and obtain the same conclusion.

Many models of inflation driven by vector fields alone have been known to be plagued by pathological behaviors, namely ghost and/or gradient instabilities. In this work, we seek a new class of vector-driven inflationary models that evade all of the mentioned instabilities. We build our analysis on the Generalized Proca Theory with an extension to three vector fields to realize isotropic expansion. We obtain the conditions required for quasi de-Sitter solutions to be an attractor analogous to the standard slow-roll one and those for their stability at the level of linearized perturbations. Identifying the remedy to the existing unstable models, we provide a simple example and explicitly show its stability. This significantly broadens our knowledge on vector inflationary scenarios, reviving potential phenomenological interests for this class of models.

Some classes of inflationary models naturally introduce two distinct metrics/frames, and their equivalence in terms of observables has often been put in question. D-brane inflation proposes candidates for an inflaton embedded in the string theory and possesses descriptions on the brane and bulk metrics/frames, which are connected by a conformal/disformal transformation that depends on the inflaton and its derivatives. It has been shown that curvature perturbations generated by the inflaton are identical in both frames, meaning that observables such as the spectrum of cosmic microwave background (CMB) anisotropies are independent of whether matter fields---including those in the standard model of particle physics---minimally couple to the brane or the bulk metric/frame. This is true despite the fact that the observables are eventually measured by the matter fields and that the total action including the matter fields is different in the two cases. In contrast, in curvaton scenarios, the observables depend on the frame to which the curvaton minimally couples. Among all inflationary scenarios, we focus on two models motivated by the KKLMMT fine-tuning problem: a slow-roll inflation with an inflection-point potential and a model of a rapidly rolling inflaton that conformally couples to gravity. In the first model, the difference between the frames in which the curvaton resides is encoded in the spectral index of the curvature perturbations, depicting the nature of the frame transformation. In the second model, the curvaton on the brane induces a spectral index significantly different from that in the bulk and is even falsified by the observations. This work thus demonstrates that two frames connected by a conformal/disformal transformation lead to different physical observables such as CMB anisotropies in curvaton models.

The standard $\Lambda$-CDM predicts a growth of structures which tends to be higher than the values of redshift space distortion (RSD) measurements, if the cosmological parameters are fixed by the CMB data. In this paper we point out that this discrepancy can be resolved/understood if we assume that the graviton has a small but non-zero mass. In the context of the Minimal Theory of Massive Gravity (MTMG), due to infrared Lorentz violations measurable only at present cosmological scales, the graviton acquires a mass without being haunted by unwanted extra degrees of freedom. While the so-called self-accelerating branch of cosmological solutions in MTMG has the same phenomenology for the background as well as the scalar- and vector-type linear perturbations as the $\Lambda$-CDM in General Relativity (GR), it is possible to choose another branch so that the background is the same as that in GR but the evolution of matter perturbations gets modified by the graviton mass. On studying the fit of such modified dynamics to the above-mentioned RSD measurements, we find that the $\Lambda$-CDM model is less probable than MTMG by two orders of magnitude. With the help of the cross correlation between the integrated Sachs-Wolfe (ISW) effect and the large scale structure (LSS), the data also pin-down the graviton mass squared around $\mu^2\approx - (3\times 10^{-33} \rm{eV})^2$, which is consistent with the latest bound $|\mu^2|<(1.2\times 10^{-22} \rm{eV})^2$ set by the recent LIGO observation.

In the context of a U(1) gauge theory non-minimally coupled to scalar-tensor gravity, we find a cosmological attractor solution that represents a de Sitter universe with a homogeneous magnetic field. The solution fully takes into account backreaction of the magnetic field to the geometry and the scalar field. Such a solution is made possible by scaling-type global symmetry and fine-tuning of two parameters of the theory. If the fine-tuning is relaxed then the solution is deformed to an axisymmetric Bianchi type-I universe with constant curvature invariants, a homogeneous magnetic field and a homogeneous electric field. Implications to inflationary magnetogenesis are briefly discussed.

The extended quasidilaton theory is one of the simplest Lorentz-invariant massive gravity theories which can accommodate a stable self-accelerating vacuum solution. In this paper we revisit this theory and study the effect of matter fields. For a matter sector that couples minimally to the physical metric, we find hints of a Jeans type instability in the IR. In the analogue k-essence field set-up, this instability manifests itself as an IR ghost for the scalar field perturbation, but this can be interpreted as a classical instability that becomes relevant below some momentum scale in terms of matter density perturbations. We also consider the effect of the background evolution influenced by matter on the stability of the gravity sector perturbations. In particular, we address the previous claims of ghost instability in the IR around the late time attractor. We show that, although the matter-induced modification of the evolution potentially brings tension to the stability conditions, one goes beyond the regime of validity of the effective theory well before the solutions become unstable. We also draw attention to the fact that the IR stability conditions are also enforced by the existence requirements of consistent background solutions.

In the setup of ghost condensation model the generalized second law of black hole thermodynamics can be respected under a radiatively stable assumption that couplings between the field responsible for ghost condensate and matter fields such as those in the Standard Model are suppressed by the Planck scale. Since not only black holes but also cosmology are expected to play important roles towards our better understanding of gravity, we consider a cosmological setup to test the theory of ghost condensation. In particular we shall show that the de Sitter entropy bound proposed by Arkani-Hamed, et.al. is satisfied if ghost inflation happened in the early epoch of our universe and if there remains a tiny positive cosmological constant in the future infinity. We then propose a notion of cosmological Page time after inflation.

We consider the finite interactions of the generalized Proca theory including the sixth-order Lagrangian and derive the full linear perturbation equations of motion on the flat Friedmann-Lema\^{i}tre-Robertson-Walker background in the presence of a matter perfect fluid. By construction, the propagating degrees of freedom (besides the matter perfect fluid) are two transverse vector perturbations, one longitudinal scalar, and two tensor polarizations. The Lagrangians associated with intrinsic vector modes neither affect the background equations of motion nor the second-order action of tensor perturbations, but they do give rise to non-trivial modifications to the no-ghost condition of vector perturbations and to the propagation speeds of vector and scalar perturbations. We derive the effective gravitational coupling $G_{\rm eff}$ with matter density perturbations under a quasi-static approximation on scales deep inside the sound horizon. We find that the existence of intrinsic vector modes allows a possibility for reducing $G_{\rm eff}$. In fact, within the parameter space, $G_{\rm eff}$ can be even smaller than the Newton gravitational constant $G$ at the late cosmological epoch, with a peculiar phantom dark energy equation of state (without ghosts). The modifications to the slip parameter $\eta$ and the evolution of growth rate $f\sigma_8$ are discussed as well. Thus, dark energy models in the framework of generalized Proca theories can be observationally distinguished from the $\Lambda$CDM model according to both cosmic growth and expansion history. Furthermore, we study the evolution of vector perturbations and show that outside the vector sound horizon the perturbations are nearly frozen and start to decay with oscillations after the horizon entry.

We consider the possibility that the massive graviton is a viable candidate of dark matter in the context of bimetric gravity. We first derive the energy-momentum tensor of the massive graviton and show that it indeed behaves as that of dark matter fluid. We then discuss a production mechanism and the present abundance of massive gravitons as dark matter. Since the metric to which ordinary matter fields couple is a linear combination of the two mass eigenstates of bigravity, production of massive gravitons, i.e. the dark matter particles, is inevitably accompanied by generation of massless gravitons, i.e. the gravitational waves. Therefore, in this scenario some information about dark matter in our universe is encoded in gravitational waves. For instance, if LIGO detects gravitational waves generated by the preheating after inflation then the massive graviton with the mass of $\sim 0.01$ GeV is a candidate of the dark matter.

Caustic singularity formations in shift-symmetric $k$-essence and Horndeski theories on a fixed Minkowski spacetime were recently argued. In $n$ dimensions, this singularity is the $(n-2)$-dimensional plane in spacetime at which second derivatives of a field diverge and the field loses single-valued description for its evolution. This does not necessarily imply a pathological behavior of the system but rather invalidates the effective description. The effective theory would thus have to be replaced by another to describe the evolution thereafter. In this paper, adopting the planar-symmetric $1$+$1$-dimensional approach employed in the original analysis, we seek all $k$-essence theories in which generic simple wave solutions are free from such caustic singularities. Contrary to the previous claim, we find that not only the standard canonical scalar but also the DBI scalar are free from caustics, as far as planar-symmetric simple wave solutions are concerned. Addition of shift-symmetric Horndeski terms does not change the conclusion.

We consider a massive vector field with derivative interactions that propagates only the 3 desired polarizations (besides two tensor polarizations from gravity) with second-order equations of motion in curved space-time. The cosmological implications of such generalized Proca theories are investigated for both the background and the linear perturbation by taking into account the Lagrangian up to quintic order. In the presence of a matter fluid with a temporal component of the vector field, we derive the background equations of motion and show the existence of de Sitter solutions relevant to the late-time cosmic acceleration. We also obtain conditions for the absence of ghosts and Laplacian instabilities of tensor, vector, and scalar perturbations in the small-scale limit. Our results are applied to concrete examples of the general functions in the theory, which encompass vector Galileons as a specific case. In such examples, we show that the de Sitter fixed point is always a stable attractor and study viable parameter spaces in which the no-ghost and stability conditions are satisfied during the cosmic expansion history.

In this work we study the partially constrained vielbein formulation of the new quasidilaton theory of massive gravity which couples to both physical and fiducial metrics simultaneously via a composite effective metric. This formalism improves the new quasidilaton model since the Boulware-Deser ghost is removed fully non-linearly at all scales. This also yields crucial implications in the cosmological applications. We derive the governing cosmological background evolution and study the stability of the attractor solution.

We propose a new theory of massive gravity with only two propagating degrees of freedom. After defining the theory in the unitary gauge in the vielbein language, we shall perform a Hamiltonian analysis to count the number of physical degrees of freedom, and then study some phenomenologies. While the homogeneous and isotropic background cosmology and the tensor linear perturbations around it are described by exactly the same equations as those in the de Rham-Gabadadze-Tolley (dRGT) massive gravity, the scalar and vector gravitational degrees of freedom are absent in the new theory at the fully nonlinear level. Hence the new theory provides a stable nonlinear completion of the self-accelerating cosmological solution that was originally found in the dRGT theory.

We study metric transformations which depend on a scalar field $\phi$ and its first derivatives and confirm that the number of physical degrees of freedom does not change under such transformations, as long as they are not singular. We perform a Hamiltonian analysis of a simple model in the gauge $\phi = t$. In addition, we explicitly show that the transformation and the gauge fixing do commute in transforming the action. We then extend the analysis to more general gravitational theories and transformations in general gauges. We verify that the set of all constraints and the constraint algebra are left unchanged by such transformations and conclude that the number of degrees of freedom is not modified by a regular and invertible generic transformation among two metrics. We also discuss the implications on the recently called "hidden" constraints and on the case of a singular transformation, a.k.a. mimetic gravity.

We consider a consistent linear effective vielbein matter coupling without introducing the Boulware-Deser ghost in ghost-free massive gravity. This is achieved in the partially constrained vielbein formulation. We first introduce the formalism and prove the absence of ghost at all scales. As next we investigate the cosmological application of this coupling in this new formulation. We show that even if the background evolution accords with the metric formulation, the perturbations display important different features in the partially constrained vielbein formulation. We study the cosmological perturbations of the two branches of solutions separately. The tensor perturbations coincide with those in the metric formulation. Concerning the vector and scalar perturbations, the requirement of absence of ghost and gradient instabilities yields slightly different allowed parameter space.

We study the nature of constraints and count the number of degrees of freedom in the nonprojectable version of the $U(1)$ extension of Ho\v{r}ava-Lifshitz gravity, using the standard method of Hamiltonian analysis in the classical field theory. This makes it possible for us to investigate the condition under which the scalar graviton is absent at a fully nonlinear level. We show that the scalar graviton does not exist at the classical level if and only if two specific coupling constants are exactly zero. The operators corresponding to these two coupling constants are marginal for any values of the dynamical critical exponent of the Lifshitz scaling and thus should be generated by quantum corrections even if they are eliminated from the bare action. We thus conclude that the theory in general contains the scalar graviton.

We study the cosmology of bimetric theory with a composite matter coupling. We find two possible branches of background evolution. We investigate the question of stability of cosmological perturbations. For the tensor and vector perturbations, we derive conditions on the absence of ghost and gradient instabilities. For the scalar modes, we obtain conditions for avoiding ghost degrees. In the first branch, we find that one of the scalar modes becomes a ghost at the late stages of the evolution. Conversely, this problem can be avoided in the second branch. However, we also find that the constraint for the second branch prevents the doubly coupled matter fields from being the standard ingredients of cosmology. We thus conclude that a realistic and stable cosmological model requires additional minimally coupled matter fields.

We present a new quasidilaton theory of Poincare invariant massive gravity, based on the recently proposed framework of matter coupling that makes it possible for the kinetic energy of the quasidilaton scalar to couple to both physical and fiducial metrics simultaneously. We find a scaling-type exact solution that expresses a self-accelerating de Sitter universe, and then analyze linear perturbations around it. It is shown that in a range of parameters all physical degrees of freedom have non-vanishing quadratic kinetic terms and are stable in the subhorizon limit, while the effective Newton's constant for the background is kept positive.

We investigate the cosmological perturbations around FLRW solutions to non- linear massive gravity with a new effective coupling to matter proposed recently. Unlike the case with minimal matter coupling, all five degrees of freedom in the gravity sector propagate on generic self-accelerating FLRW backgrounds. We study the stability of the cosmological solutions and put constraints on the parameters of the theory by demanding the correct sign for the kinetic terms for scalar, vector and tensor perturbations.

We investigate the cosmology of SO(3)-invariant massive gravity with 5 degrees of freedom. In contrast with previous studies, we allow for a non-trivial fiducial metric, which can be justified by invoking, for example, a dilaton-like global symmetry. We write the homogeneous and isotropic equations of motion in this more general setup and identify, in particular, de Sitter solutions. We then study the linear perturbations around the homogeneous cosmological solutions, by deriving the quadratic Lagrangian governing the dynamics of scalar, vector and tensor modes. We thus obtain the conditions for the perturbations to be well-behaved. We show that it is possible to find de Sitter solutions whose perturbations are weakly coupled and stable, i.e. without ghost-like or gradient instabilities.

We discuss a mechanism of particle production during inflation that can result in a blue gravitational wave (GW) spectrum, compatible with the BICEP2 result and with the r < 0.11 limit on the tensor-to-scalar ratio at the Planck pivot scale. The mechanism is based on the production of vector quanta from a rolling pseudo-scalar field. Both the vector and the pseudo-scalar are only gravitationally coupled to the inflaton, to keep the production of inflaton quanta at an unobservable level (the overproduction of non-gaussian scalar perturbations is a generic difficulty for mechanisms that aim to generate a visible GW signal from particle production during inflation). This mechanism can produce a detectable amount of GWs for any inflationary energy scale. The produced GWs are chiral and non-gaussian; both these aspects can be tested with large-scale polarization data (starting from Planck). We study how to reconstruct the pseudo-scalar potential from the GW spectrum.

We study cosmological perturbations in bimetric theory with two fluids each of which is coupled to one of the two metrics. Focusing on a healthy branch of background solutions, we clarify the stability of the cosmological perturbations. For this purpose, we extend the condition for the absence of the so-called Higuchi ghost, and show that the condition is guaranteed to be satisfied on the healthy branch. We also calculate the squared propagation speeds of perturbations and derive the conditions for the absence of the gradient instability. To avoid the gradient instability, we find that the model parameters are weakly constrained.

This article investigates the construction of fermions and the formulation of the Standard Model of particle physics in a theory in which the Lorentz signature emerges from an underlying microscopic purely Euclidean $SO(4)$ theory. Couplings to a clock field are responsible for triggering the change of signature of the effective metric in which the standard fields propagate. We demonstrate that Weyl and Majorana fermions can be constructed in this framework. This construction differs from other studies of Euclidean fermions, as the coupling to the clock field allows us to write down an action which flows to the usual action in Minkowski spacetime. We then show how the Standard Model can be obtained in this theory and consider the constraints on non-Standard-Model operators which can appear in the QED sector due to CPT and Lorentz violation.

In this paper, we first propose a universal coupling between the gravity and matter in the framework of the Ho\v{r}ava-Lifshitz theory of gravity with an extra U(1) symmetry for both the projectable and non-projectable cases. Then, using this universal coupling we study the post-Newtonian approximations and obtain the parameterized post-Newtonian (PPN) parameters in terms of the coupling constants of the theory. Contrary to the previous works in which only two PPN parameters were calculated, we obtain {\it all} PPN parameters. We then, for the first time in either projectable or non-projectable case, find that all the solar system tests carried out so far are satisfied in a large region of the parameters space. In particular, the same results obtained in general relativity can be easily realized here. A remarkable feature is that the solar system tests impose no constraint on the parameter $\lambda$ appearing in the kinetic part of the action. As a result, the solar system tests, when combined with the condition for avoidance of strong coupling, do not lead to an upper bound on the energy scale $M_{*}$ that suppresses higher dimensional operators in the theory. This is in sharp contrast to other versions of the HL theory.

We investigate the Euclidean vacuum mode functions of a massive vector field in a spatially open chart of de Sitter spacetime. In the one-bubble open inflationary scenario that naturally predicts a negative spatial curvature after a quantum tunneling, it is known that a light scalar field has the so-called supercurvature mode, i.e. an additional discrete mode which describes fluctuations over scales larger than the spatial curvature scale. If such supercurvature modes exist for a vector field with a sufficiently light mass, then they would decay slower and easily survive the inflationary era. However, the existence of supercurvature mode strongly depends on details of the system. To clarify whether a massive vector field has supercurvature modes, we consider a U(1) gauge field with gauge and conformal invariances spontaneously broken through the Higgs mechanism, and present explicit expressions for the Euclidean vacuum mode functions. We find that, for any values of the vector field mass, there is no supercurvature mode. In the massless limit, the absence of supercurvature modes in the scalar sector stems from the gauge symmetry.

Recently the first example of a unitary theory of Lorentz-invariant massive gravity allowing for stable self-accelerating de Sitter solutions was found, extending the quasidilaton theory. In this paper we further generalize this new action for the quasidilaton field by introducing general Lagrangian terms which are consistent with the quasidilaton symmetry while leading to second order equations of motion. We find that the structure of the theory, compared to the simplest stable example, does not change on introducing these new terms.

We present a review of cosmological solutions in non-linear massive gravity, focusing on the stability of perturbations. Although homogeneous and isotropic solutions have been found, these are now known to suffer from either Higuchi ghost or a new non-linear ghost instability. We discuss two approaches to alleviate this issue. By relaxing the symmetry of the background by e.g. breaking isotropy in the hidden sector, it is possible to accommodate a stable cosmological solution. Alternatively, extending the theory to allow for new dynamical degrees of freedom can also remove the conditions which lead to the instability. As examples for this case, we study the stability of self-accelerating solutions in the quasi-dilatonic extension and generic cosmological solutions in the varying mass extension. While the quasi-dilaton case turns out to be unstable, the varying mass case allows stable regimes of parameters. Viable self-accelerating solutions in the varying mass theory yet remain to be found.

We present the first example of a unitary theory of Lorentz-invariant massive gravity, with all degrees of freedom propagating on a strictly homogeneous and isotropic, self-accelerating de Sitter background. The theory is a simple extension of the quasidilaton theory, respecting the symmetry of the original theory but allowing for a new type of coupling between the massive graviton and the quasidilaton scalar.

We investigate nonlinear stability of two classes of cosmological solutions in massive gravity: isotropic Friedmann-Lemaitre-Robertson-Walker (FLRW) solutions and anisotropic FLRW solutions. For this purpose we construct the linear cosmological perturbation theory around axisymmetric Bianchi type--I backgrounds. We then expand the background around the two classes of solutions, which are fixed points of the background evolution equation, and analyze linear perturbations on top of it. This provides a consistent truncation of nonlinear perturbations around these fixed point solutions and allows us to analyze nonlinear stability in a simple way. In particular, it is shown that isotropic FLRW solutions exhibit nonlinear ghost instability. On the other hand, anisotropic FLRW solutions are shown to be ghost-free for a range of parameters and initial conditions.

We suggest a new scenario of gravitation in which gravity at the fundamental level is described by a Riemannian (i.e. locally Euclidean) theory without the notion of time. The Lorentzian metric structure and the notion of time emerge as effective properties at long distances. On the other hand, at short distances, higher derivative terms compatible with the Riemannian diffeomorphism become important and the system is described by a power-counting renormalizable Riemannian theory.

The Lorentzian metric structure used in any field theory allows one to implement the relativistic notion of causality and to define a notion of time dimension. This article investigates the possibility that at the microscopic level the metric is Riemannian, i.e. locally Euclidean, and that the Lorentzian structure, that we usually consider as fundamental, is in fact an effective property that emerges in some regions of a 4-dimensional space with a positive definite metric. In such a model, there is no dynamics nor signature flip across some hypersurface; instead, all the fields develop a Lorentzian dynamics in these regions because they propagate in an effective metric. It is shown that one can construct a decent classical field theory for scalars, vectors and (Dirac) spinors in flat spacetime. It is then shown that gravity can be included but that the theory for the effective Lorentzian metric is not general relativity but of the covariant Galileon type. The constraints arising from stability, the equivalence principle and the constancy of fundamental constants are detailed and a phenomenological picture of the emergence of the Lorentzian metric is also given. The construction, while restricted to classical fields in this article, offers a new view on the notion of time.

We study cosmological perturbations around self-accelerating solutions to two extensions of nonlinear massive gravity: the quasi-dilaton theory and the mass-varying theory. We examine stability of the cosmological solutions, and the extent to which the vanishing of the kinetic terms for scalar and vector perturbations of self-accelerating solutions in massive gravity is generic when the theory is extended. We find that these kinetic terms are in general non-vanishing in both extensions, though there are constraints on the parameters and background evolution from demanding that they have the correct sign. In particular, the self-accelerating solutions of the quasi-dilaton theory are always unstable to scalar perturbations with wavelength shorter than the Hubble length.

We study spherically symmetric, stationary vacuum configurations in general covariant theory (U(1) extension) of Ho\v{r}ava-Lifshitz gravity with the projectability condition and an arbitrary coupling constant $\lambda$, and obtain all the solutions in closed forms. If the gauge field $A$ and the Newtonian prepotential $\varphi$ do not directly couple to matter fields, the theory is inconsistent with solar system tests for $\lambda\not=1$, no matter how small $|\lambda-1|$ is. This is shown to be true also with the most general ansatz of spherical (but not necessarily stationary) configurations. Therefore, to be consistent with observations, one needs either to find a mechanism to restrict $\lambda$ precisely to $\lambda_{GR}=1$, or to consider $A$ and/or $\varphi$ as parts of the 4-dimensional metric on which matter fields propagate. In the latter, requiring that the line element be invariant not only under the foliation-preserving diffeomorphism but also under the local U(1) transformations, we propose the replacements, $N \rightarrow N - \upsilon(A - {\cal{A}})/c^2$ and $N^i \rightarrow N^i+N\nabla^{i}\varphi$, where $\upsilon$ is a dimensionless coupling constant to be constrained by observations, $N$ and $N^i$ are, respectively, the lapse function and shift vector, and ${\cal{A}} \equiv - \dot{\varphi} + N^i\nabla_{i}\varphi + N(\nabla_{i}\varphi)^2/2$. With this prescription, we show explicitly that the aforementioned solutions are consistent with solar system tests for both $\lambda=1$ and $\lambda\not=1$, provided that $|\upsilon-1|<10^{-5}$. From this result, the physical and geometrical interpretations of the fields $A$ and $\varphi$ become clear. However, it still remains to be understood how to obtain such a prescription from the action principle.

A stochastic gravitational-wave background (SGWB) is expected to arise from the superposition of many independent and unresolved gravitational-wave signals, of either cosmological or astrophysical origin. Some cosmological models (characterized, for instance, by a pseudo-scalar inflaton, or by some modification of gravity) break parity, leading to a polarized SGWB. We present a new technique to measure this parity violation, which we then apply to the recent results from LIGO to produce the first upper limit on parity violation in the SGWB, assuming a generic power-law SGWB spectrum across the LIGO sensitive frequency region. We also estimate sensitivity to parity violation of the future generations of gravitational-wave detectors, both for a power-law spectrum and for a model of axion inflation. This technique offers a new way of differentiating between the cosmological and astrophysical sources of the isotropic SGWB, as astrophysical sources are not expected to produce a polarized SGWB.

We discuss the detectability of gravitational waves with a time dependent mass contribution, by means of the stochastic gravitational wave observations. Such a mass term typically arises in the cosmological solutions of massive gravity theories. We conduct the analysis based on a general quadratic action, and thus the results apply universally to any massive gravity theories in which modification of general relativity appears primarily in the tensor modes. The primary manifestation of the modification in the gravitational wave spectrum is a sharp peak. The position and height of the peak carry information on the present value of the mass term, as well as the duration of the inflationary stage. We also discuss the detectability of such a gravitational wave signal using the future-planned gravitational wave observatories.

We report the strictest observational verification of CPT invariance in the photon sector, as a result of gamma-ray polarization measurement of distant gamma-ray bursts (GRBs), which are brightest stellar-size explosions in the universe. We detected the gamma-ray polarization of three GRBs with high significance, and the source distances may be constrained by a well-known luminosity indicator for GRBs. For the Lorentz- and CPT-violating dispersion relation E_{\pm}^2=p^2 \pm 2\xi p^3/M_{Pl}, where \pm denotes different circular polarization states of the photon, the parameter \xi is constrained as |\xi|<O(10^{-15}). Barring precise cancellation between quantum gravity effects and dark energy effects, the stringent limit on the CPT-violating effect leads to the expectation that quantum gravity presumably respects the CPT invariance.

Recently observational lower bounds on the strength of cosmic magnetic fields were reported, based on gamma-ray flux from distant blazars. If inflation is responsible for the generation of such magnetic fields then the inflation energy scale is bounded from above as rho_{inf}^{1/4} < 2.5 times 10^{-7}M_{Pl} times (B_{obs}/10^{-15}G)^{-2} in a wide class of inflationary magnetogenesis models, where B_{obs} is the observed strength of cosmic magnetic fields. The tensor-to-scalar ratio is correspondingly constrained as r< 10^{-19} times (B_{obs}/10^{-15}G)^{-8}. Therefore, if the reported strength B_{obs} \geq 10^{-15}G is confirmed and if any signatures of gravitational waves from inflation are detected in the near future, then our result indicates some tensions between inflationary magnetogenesis and observations.

We argue that all homogeneous and isotropic solutions in nonlinear massive gravity are unstable. For this purpose, we study the propagating modes on a Bianchi type--I manifold. We analyze their kinetic terms and dispersion relations as the background manifold approaches the homogeneous and isotropic limit. We show that in this limit, at least one ghost always exists and that its frequency tends to vanish for large scales, meaning that it cannot be integrated out from the low energy effective theory. This ghost mode is interpreted as a leading nonlinear perturbation around a homogeneous and isotropic background.

In the scope of the nonlinear massive gravity, we study fixed points of evolution equations for a Bianchi type--I universe. We find a new attractor solution with non-vanishing anisotropy, on which the physical metric is isotropic but the Stuckelberg configuration is anisotropic. As a result, at the background level, the solution describes a homogeneous and isotropic universe, while a statistical anisotropy is expected from perturbations, suppressed by smallness of the graviton mass.

We present a fully nonlinear study of long wavelength cosmological perturbations within the framework of the projectable Horava-Lifshitz gravity, coupled to a single scalar field. Adopting the gradient expansion technique, we explicitly integrate the dynamical equations up to any order of the expansion, then restrict the integration constants by imposing the momentum constraint. While the gradient expansion relies on the long wavelength approximation, amplitudes of perturbations do not have to be small. When the $\lambda\to 1$ limit is taken, the obtained nonlinear solutions exhibit a continuous behavior at any order of the gradient expansion, recovering general relativity in the presence of a scalar field and the "dark matter as an integration constant". This is in sharp contrast to the results in the literature based on the "standard" (and naive) perturbative approach where in the same limit, the perturbative expansion of the action breaks down and the scalar graviton mode appears to be strongly coupled. We carry out a detailed analysis on the source of these apparent pathologies and determine that they originate from an improper application of the perturbative approximation in the momentum constraint. We also show that there is a new branch of solutions, valid in the regime where $|\lambda-1|$ is smaller than the order of perturbations. In the limit $\lambda\to1$, this new branch allows the theory to be continuously connected to general relativity (plus "dark matter").

We study cosmological perturbations of self-accelerating universe solutions in the recently proposed nonlinear theory of massive gravity, with general matter content. While the broken diffeomorphism invariance implies that there generically are 2 tensor, 2 vector and 2 scalar degrees of freedom in the gravity sector, we find that the scalar and vector degrees have vanishing kinetic terms and nonzero mass terms. Depending on their nonlinear behavior, this indicates either nondynamical nature of these degrees or strong couplings. Assuming the former, we integrate out the 2 vector and 2 scalar degrees of freedom. We then find that in the scalar and vector sectors, gauge-invariant variables constructed from metric and matter perturbations have exactly the same quadratic action as in general relativity. The difference from general relativity arises only in the tensor sector, where the graviton mass modifies the dispersion relation of gravitational waves, with a time-dependent effective mass. This may lead to modification of stochastic gravitational wave spectrum.

In the context of a recently proposed nonlinear massive gravity with Lorentz-invariant mass terms, we investigate open Friedmann-Robertson-Walker (FRW) universes driven by arbitrary matter source. While the flat FRW solutions were recently shown to be absent, the proof does not extend to the open universes. We find three independent branches of solutions to the equations of motion for the St\"uckelberg scalars. One of the branches does not allow any nontrivial FRW cosmologies, as in the previous no-go result. On the other hand, both of the other two branches allow general open FRW universes governed by the Friedmann equation with the matter source, the standard curvature term and an effective cosmological constant $\Lambda_{\pm}=c_{\pm}m_g^2$. Here, $m_g$ is the graviton mass, +and - represent the two branches, and $c_{\pm}$ are constants determined by the two dimensionless parameters of the theory. Since an open FRW universe with a sufficiently small curvature constant can approximate a flat FRW universe but there is no exactly flat FRW solution, the theory exhibits a discontinuity at the flat FRW limit.

We perform a fully nonlinear analysis of superhorizon perturbation in Ho\v{r}ava-Lifshitz gravity, based on the gradient expansion method. We present a concrete expression for the solution of gravity equations up to the second order in the gradient expansion, and prove that the solution can be extended to any order. The result provides yet another example for analogue of the Vainshtein effect: the nonlinear solution is regular in the limit $\lambda\to 1$ and recovers general relativity coupled to dark matter at low energy. Finally, we propose a definition of nonlinear curvature perturbation ${\cal R}$ in Ho\v{r}ava-Lifshitz gravity and show that it is conserved up to the first order in the gradient expansion.

In the framework of the power-counting renormalizable theory of gravitation, recently proposed by Ho\v{r}ava, we study the limit $\lambda\to\infty$, which is arguably the most natural candidate for the ultraviolet fixed point of the renormalization group flow. In the projectable version with dynamical critical exponent $z=3$, we analyze the Friedmann-Robertson-Walker background driven by the so-called "dark matter as integration constant" and perturbations around it. We show that amplitudes of quantum fluctuations for both scalar and tensor gravitons remain finite in the limit and that the theory is weakly coupled under a certain condition.

We consider CPT violation in neutrino sector, which is induced by ghost condensation. A model with extra dimension is suggested where ghost condensation occurs at a distant location separated from the SM brane. Right handed neutrinos in the bulk, which are originally introduced to explain small Yukawa couplings, play the role of messenger fields communicating ghost condensation and the standard model sector and lead to a sizable CPT violation in neutrino sector at the leading order. The model provides a resolution to the recent MINOS anomaly without spoiling any experimental constraints and may be able to be tested by observing an interesting phenomenon, twinkling cosmic microwave background radiation, with timescale about O(10-100) minutes at future CMB observations e.g. Planck.

A Lifshitz scalar with the dynamical critical exponent z = 3 obtains scale-invariant, super-horizon field fluctuations without the need of an inflationary era. Since this mechanism is due to the special scaling of the Lifshitz scalar and persists in the presence of unsuppressed self-couplings, the resulting fluctuation spectrum can deviate from a Gaussian distribution. We study the non-Gaussian nature of the Lifshitz scalar's intrinsic field fluctuations, and show that primordial curvature perturbations sourced from such field fluctuations can have large non-Gaussianity of order f_NL = O(100), which will be detected by upcoming CMB observations. We compute the bispectrum and trispectrum of the fluctuations, and discuss their configurations in momentum space. In particular, the bispectrum is found to take various shapes, including the local, equilateral, and orthogonal shapes. Intriguingly, all integrals in the in-in formalism can be performed analytically.

This article reviews basic construction and cosmological implications of a power-counting renormalizable theory of gravitation recently proposed by Horava. We explain that (i) at low energy this theory does not exactly recover general relativity but instead mimic general relativity plus dark matter; that (ii) higher spatial curvature terms allow bouncing and cyclic universes as regular solutions; and that (iii) the anisotropic scaling with the dynamical critical exponent z=3 solves the horizon problem and leads to scale-invariant cosmological perturbations even without inflation. We also comment on issues related to an extra scalar degree of freedom called scalar graviton. In particular, for spherically-symmetric, static, vacuum configurations we prove non-perturbative continuity of the lambda->1+0 limit, where lambda is a parameter in the kinetic action and general relativity has the value lambda=1. We also derive the condition under which linear instability of the scalar graviton does not show up.

We derive the equations of linear cosmological perturbations for the general Lagrangian density $f (R,\phi, X)/2+L_c$, where $R$ is a Ricci scalar, $\phi$ is a scalar field, and $X=-(\nabla \phi)^2/2$ is a field kinetic energy. We take into account a nonlinear self-interaction term $L_c$ recently studied in the context of "Galileon" cosmology, which keeps the field equations at second order. Taking into account a scalar-field mass explicitly, the equations of matter density perturbations and gravitational potentials are obtained under a quasi-static approximation on sub-horizon scales. We also derive conditions for the avoidance of ghosts and Laplacian instabilities associated with propagation speeds. Our analysis includes most of modified gravity models of dark energy proposed in literature and thus it is convenient to test the viability of such models from both theoretical and observational points of view.

In the inflationary universe, there can be light fields other than the inflaton. We explore a possibility that such light fields source the primordial perturbations, while minimally affecting the inflaton dynamics. We show that during inflation, fluctuations of the light fields can be converted to adiabatic curvature perturbations, which accumulate and become significant by the end of the inflationary era. An additional goal of this work is to distinguish between light fields which can/cannot be ignored during inflation. Such criteria become useful for examining cosmological scenarios with multiple fields. As concrete examples, our results are applied to D-brane inflation models. We consider effects from KK modes (oscillation modes) of wrapped branes in monodromy-driven large-field models, and angular directions of throat geometries in warped D-brane inflation.

We revisit ghost dark matter, the possibility that ghost condensation may serve as an alternative to dark matter. In particular, we investigate the Friedmann-Robertson-Walker (FRW) background evolution and the large-scale structure (LSS) in the $\Lambda$GDM universe, i.e. a late-time universe dominated by a cosmological constant and ghost dark matter. The FRW background of the $\Lambda$GDM universe is indistinguishable from that of the standard $\Lambda$CDM universe if $M\gtrsim 1 {\rm eV}$, where $M$ is the scale of spontaneous Lorentz breaking. From the LSS we find a stronger bound: $M\gtrsim 10 {\rm eV}$. For smaller $M$, ghost dark matter would have non-negligible sound speed after the matter-radiation equality, and thus the matter power spectrum would significantly differ from observation. These bounds are compatible with the phenomenological upper bound $M\lesssim 100 {\rm GeV}$ known in the literature.

Ghost inflation predicts almost scale-invariant primordial cosmological perturbations with relatively large non-Gaussianity. The bispectrum is known to have a large contribution at the wavenumbers forming an equilateral triangle and the corresponding nonlinear parameter $f_{NL}^{equil}$ is typically of order $O(10^2)$. In this paper we calculate trispectrum from ghost inflation and show that the corresponding nonlinear parameter $\tau_{NL}$ is typically of order $O(10^4)$. We investigate the shape dependence of the trispectrum and see that it has some features different from DBI inflation. Therefore, our result may be useful as a template to distinguish ghost inflation from other models of inflation by future experiments.

We develop a theory of nonlinear cosmological perturbations on superhorizon scales for a single scalar field with a general kinetic term and a general form of the potential. We employ the ADM formalism and the spatial gradient expansion approach, characterised by $O(\epsilon^m)$, where $\epsilon=1/(HL)$ is a small parameter representing the ratio of the Hubble radius to the characteristic length scale $L$ of perturbations. We obtain the general solution for a full nonlinear version of the curvature perturbation valid up through second-order in $\epsilon$ ($m=2$). We find the solution satisfies a nonlinear second-order differential equation as an extension of the equation for the linear curvature perturbation on the comoving hypersurface. Then we formulate a general method to match a perturbative solution accurate to $n$-th-order in perturbation inside the horizon to our nonlinear solution accurate to second-order ($m=2$) in the gradient expansion on scales slightly greater than the Hubble radius. The formalism developed in this paper allows us to calculate the superhorizon evolution of a primordial non-Gaussianity beyond the so-called $\delta N$ formalism or separate universe approach which is equivalent to leading order ($m=0$) in the gradient expansion. In particular, it can deal with the case when there is a temporary violation of slow-roll conditions. As an application of our formalism, we consider Starobinsky's model, which is a single field model having a temporary non-slow-roll stage due to a sharp change in the potential slope. We find that a large non-Gaussianity can be generated even on superhorizon scales due to this temporary suspension of slow-roll inflation.

In Horava-Lifshitz gravity, regularity of a solution requires smoothness of not only the spacetime geometry but also the foliation. As a result, the regularity condition at the center of a star is more restrictive than in general relativity. Assuming that the energy density is a piecewise-continuous, non-negative function of the pressure and that the pressure at the center is positive, we prove that the momentum conservation law is incompatible with the regularity at the center for any spherically-symmetric, static configurations. The proof is totally insensitive to the structure of higher spatial curvature terms and, thus, holds for any values of the dynamical critical exponent $z$. Therefore, we conclude that a spherically-symmetric star should include a time-dependent region near the center. We also comment on the condition under which linear instability of the scalar graviton does not show up.

In scalar-tensor theories the scalar fields generically couple nontrivially to gravity. We study the observable properties of inflationary models with non-minimally coupled inflaton and Dirac-Born-Infeld (DBI) kinetic term. Within the assumptions of the priors of our Monte-Carlo simulations we find these models can generate new interesting observable signatures. Our discussion focuses on string theory inspired phenomenological models of relativistic D-brane inflation. While successful string theory constructions of ultra-violet DBI brane inflation remain elusive, we show that in suitable regions of the parameter space it is possible to use cosmological observables to probe the non-minimial coupling. Fortunately, the most observationally promising range of parameters include models yielding intermediate levels of non-gaussianity in the range consistent with WMAP 5-year data, and to be constrained further by the Planck satellite.

The origin of large-scale magnetic field in the universe is one of the greatest mysteries in modern cosmology. We present a new mechanism for generation of large-scale magnetic field, based on the power-counting renormalizable theory of gravitation recently proposed by Horava. Contrary to the usual case in general relativity, the U(1) gauge symmetry of a Maxwell action in this theory permits terms breaking conformal invariance in the ultraviolet. Moreover, for high frequency modes, the anisotropic scaling intrinsic to the theory inevitably makes the sound horizon far outside the Hubble horizon. Consequently, non-inflationary cosmic expansion in the early universe naturally generates super-horizon quantum fluctuations of the magnetic field. Specializing our consideration to the case with the dynamical critical exponent $z=3$, we show an explicit set of parameters for which (i) the amplitude of generated magnetic field is large enough as a seed for the dynamo mechanism; (ii) backreaction to the cosmic expansion is small enough; and (iii) the high-energy dispersion relation is consistent with the most recent observational limits from MAGIC and FERMI.

By looking for possible violation of the generalized second law, we might be able to find regions in the space of theories and states that do not allow holographic dual descriptions. We revisit three proposals for violation of the generalized second law in the simplest Higgs phase of gravity called ghost condensate. Two of them, (i) analogue of Penrose process and (ii) semiclassical heat flow, are based on Lorentz breaking effects, by which particles of different species can have different limits of speed. We show that processes in both (i) and (ii) are always slower than accretion of ghost condensate and cannot decrease the total entropy before the accretion increases the entropy. The other proposal is to use (iii) negative energy carried by excitations of ghost condensate. We prove an averaged null energy condition, which we conjecture prevents the proposal (iii) from violating the generalized second law in a coarse-grained sense.

We obtain cosmological constraints on models of inflation which exhibit rapid roll solutions. Rapid roll attractors exist for potentials with large mass terms and are thus of interest for inflationary model building within string theory. We constrain a general ansatz for the power spectrum arising from rapid roll inflation that, in the small field limit, can be associated with tree level hybrid potentials with variable mass terms and nonminimal gravitational coupling $\xi R\phi^2$. We consider perturbations generated through modulated reheating and/or curvaton mechanisms in place of the observationally unacceptable primary spectra generated by inflaton fluctuations in these models. The lack of a hierarchy amongst higher-order $k$-dependencies of the power spectrum results in models with potentially large runnings, allowing us to impose tight constraints on such models using CMB and LSS data. In particular, we find $n_s <1$ and $|\alpha| < 0.01$. We conclude with a concrete realization of rapid roll inflation within warped throat brane inflation that is in good agreement with current data.

Dubovsky and Sibiryakov recently proposed a scenario in which particles of different species propagate with different speeds due to their direct couplings to ghost condensate. It was argued that this extended version of ghost condensate allows a gedanken experiment leading to violation of the generalized second law. However, in the original ghost condensate scenario, difference in propagation speeds is suppressed by M^2/M_{Pl}^2, where M is the order parameter of spontaneous Lorentz breaking and M_{Pl} is the Planck scale. In this case the energy transfer necessary for the gedanken experiment is so slow that the timescale of decrease of entropy, if any, is always longer than the Jeans timescale of ghost condensate. Hence the generalized second law is not violated by the gedanken experiment in the original ghost condensate scenario. This conclusion trivially extends to gauged ghost condensation by taking into account accretion of gauged ghost condensate into a black hole.

There are at least four versions of Horava-Lishitz gravity in the literature. We consider the version without the detailed balance condition with the projectability condition and address one aspect of the theory: avoidance of caustics for constant time hypersurfaces. We show that there is no caustic with plane symmetry in the absence of matter source if \lambda\ne 1. If \lambda=1 is a stable IR fixed point of the renormalization group flow then \lambda is expected to deviate from 1 near would-be caustics, where the extrinsic curvature increases and high-energy corrections become important. Therefore, the absence of caustics with \lambda\ne 1 implies that caustics cannot form with this symmetry in the absence of matter source. We argue that inclusion of matter source will not change the conclusion. We also argue that caustics with codimension higher than one will not form because of repulsive gravity generated by nonlinear higher curvature terms. These arguments support our conjecture that there is no caustic for constant time hypersurfaces. Finally, we discuss implications to the recently proposed scenario of ``dark matter as integration constant''.

In the non-relativistic theory of gravitation recently proposed by Horava, the Hamiltonian constraint is not a local equation satisfied at each spatial point but an equation integrated over a whole space. The global Hamiltonian constraint is less restrictive than its local version, and allows a richer set of solutions than in general relativity. We show that a component which behaves like pressureless dust emerges as an "integration constant" of dynamical equations and momentum constraint equations. Consequently, classical solutions to the infrared limit of Horava-Lifshitz gravity can mimic general relativity plus cold dark matter.

We present a curvaton model from type IIB string theory compactified on a warped throat with approximate isometries. Considering an (anti-)D3-brane sitting at the throat tip as a prototype standard model brane, we show that the brane's position in the isometry directions can play the role of curvatons. The basic picture is that the fluctuations of the (anti-)D3-brane in the angular isometry directions during inflation eventually turns into the primordial curvature perturbations, and subsequently the brane's oscillation excites other open string modes on the brane and reheat the universe. We find in the explicit case of the KS throat that a wide range of parameters allows a consistent curvaton scenario. It is also shown that the oscillations of branes at throat tips are capable of producing large non-Gaussianity, either through curvature or isocurvature perturbations. Since such setups naturally arise in warped (multi-)throat compactifications and are constrained by observational data, the model can provide tests for compactification scenarios. This work gives an explicit example of string theory providing light fields for generating curvature perturbations. Such mechanisms free the inflaton from being responsible for the perturbations, thus open up new possibilities for inflation models.

Based on the renormalizable theory of gravitation recently proposed by Horava, we present a simple scenario to generate almost scale-invariant, super-horizon curvature perturbations. The anisotropic scaling with dynamical critical exponent z=3 implies that the amplitude of quantum fluctuations of a free scalar field generated in the early epoch of the expanding universe is insensitive to the Hubble expansion rate and, thus, scale-invariant. Those fluctuations are later converted to curvature perturbations by the curvaton mechanism or/and the modulated decay of heavy particles/oscillating fields. This scenario works, for example, for power law expansion a\propto t^p with p>1/3 and, thus, does not require inflation. Also, this scenario does not rely on any additional assumptions such as the detailed balance condition.

We show that, assuming the dispersion relation proposed recently by Horava in the context of quantum gravity, radiation energy density exhibits a peculiar dependence on the scale factor; the radiation energy density decreases proportional to a^{-6}. This simple scaling can have an impact on cosmology. As an example, we show that the resultant baryon asymmetry as well as the stochastic gravity waves can be enhanced. We also discuss current observational constraint on the dispersion relation.

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 investigate density perturbations generated through modulated reheating while inflation is driven by a conformally coupled scalar field. A large running of the spectral index is obtained, which reflects the basic nature of conformal inflation that higher-order time derivatives of the Hubble parameter during inflation are not necessarily small. This feature may allow us to distinguish between conformal inflation models and standard minimally coupled ones. We also investigate how the resulting fluctuations are modified when there is a deviation from an exact conformal coupling between the inflaton and gravity. Finally, we apply our results to the warped brane inflation model and see that observational bounds from the WMAP5 data suggest a blue tilted density perturbation spectrum.

We propose a mechanism of reheating after inflation in multi-throat scenarios of warped extra dimensions. Validity of an effective field theory on the standard model (SM) brane requires that the position of the SM brane during inflation be different from the position after inflation. The latter is supposed to be near the tip of the SM throat but the former is not. After inflation, when the Hubble expansion rate becomes sufficiently low, the SM brane starts moving towards the tip and eventually oscillates. The SM fields are excited by the brane motion and the universe is reheated. Since interaction between the brane position modulus and the SM fields is suppressed only by the local string scale, the modulus effectively decays into the SM fields.

We develop a theory of non-linear cosmological perturbations at superhorizon scales for a scalar field with a Lagrangian of the form $P(X,\phi)$, where $X=-\partial^{\mu}\phi\partial_{\mu}\phi$ and $\phi$ is the scalar field. We employ the ADM formalism and the spatial gradient expansion approach to obtain general solutions valid up to the second order in the gradient expansion. This formulation can be applied to, for example, DBI inflation models to investigate superhorizon evolution of non-Gaussianities. With slight modification, we also obtain general solutions valid up to the same order for a perfect fluid with a general equation of state $P=P(\rho)$.

We derive constraints on the tensor to scalar ratio and on the background charge of the warped throat for DBI inflation driven by D5- and D7-branes wrapped over cycles of the throat. It is shown that the background charge well beyond the known maximal value is required in most cases for DBI inflation to generate cosmological observables compatible with the WMAP3 data. Most of the results derived in this paper are insensitive to the details of the inflaton potential, and could be applied to generic warped throats.

Usual inflation is realized with a slow rolling scalar field minimally coupled to gravity. In contrast, we consider dynamics of a scalar with a flat effective potential, conformally coupled to gravity. Surprisingly, it contains an attractor inflationary solution with the rapidly rolling inflaton field. We discuss models with the conformal inflaton with a flat potential (including hybrid inflation). There is no generation of cosmological fluctuations from the conformally coupled inflaton. We consider realizations of modulated (inhomogeneous reheating) or curvaton cosmological fluctuations in these models. We also implement these unusual features for the popular string-theoretic warped inflationary scenario, based on the interacting D3-anti D3 branes. The original warped brane inflation suffers a large inflaton mass due to conformal coupling to 4-dimensional gravity. Instead of considering this as a problem and trying to cure it with extra engineering, we show that warped inflation with the conformally coupled, rapidly rolling inflaton is yet possible with N=37 efoldings, which requires low energy scales 1-100 TeV of inflation. Coincidentally, the same warping numerology can be responsible for the hierarchy. It is shown that the scalars associated with angular isometries of the warped geometry of compact manifold (e.g. S^3 of KS geometry) have solutions identical to conformally coupled modes and also cannot be responsible for cosmological fluctuations. We discuss other possibilities.

In the simplest Higgs phase of gravity called ghost condensation, an accelerating universe with a phantom era (w<-1) can be realized without ghost or any other instabilities. In this paper we show how to reconstruct the potential in the Higgs sector Lagrangian from a given cosmological history (H(t), \rho(t)). This in principle allows us to constrain the potential by geometrical information of the universe such as supernova distance-redshift relation. We also derive the evolution equation for cosmological perturbations in the Higgs phase of gravity by employing a systematic low energy expansion. This formalism is expected to be useful to test the theory by dynamical information of large scale structure in the universe such as cosmic microwave background anisotropy, weak gravitational lensing and galaxy clustering.

We investigate the universal low-energy dynamics of the simplest Higgs phase for gravity, `ghost condensation.' We show that the nonlinear dynamics of the `ghostone' field dominate for all interesting gravitational sources. Away from caustic singularities, the dynamics is equivalent to the irrotational flow of a perfect fluid with equation of state p \propto \rho^2, where the fluid particles can have negative mass. We argue that this theory is free from catastrophic instabilities due to growing modes, even though the null energy condition is violated. Numerical simulations show that solutions generally have singularities in which negative energy regions shrink to zero size. We exhibit partial UV completions of the theory in which these singularities are smoothly resolved, so this does not signal any inconsistency in the effective theory. We also consider the bounds on the symmetry breaking scale M in this theory. We argue that the nonlinear dynamics cuts off the Jeans instability of the linear theory, and allows M \lsim 100MeV.

It is pointed out that in the warped string compactification, motion of anti-D-branes near the bottom of a throat behaves like dark matter. Several scenarios for production of the dark matter are suggested, including one based on the D/anti-D interaction at the late stage of D/anti-D inflation.

We find a simple exact solution of 6-dimensional braneworld which captures some essential features of warped flux compactification, including a warped geometry, compactification, a magnetic flux, and one or two 3-brane(s). In this setup we analyze how the Hubble expansion rate on each brane changes when the brane tension changes. It is shown that effective Newton's constant resulting from this analysis agrees with that inferred by simply integrating extra dimensions out. Based on the result, a general formula for effective Newton's constant is conjectured and its application to cosmology with type IIB warped string compactification is discussed.

A tachyon is considered to be sick in the context of particle mechanics, but in field theory just indicates instability of a background. We consider a similar possibility that a ghost in field theory might be just an indication of instability of a background and that it can condense to form a different background around which there is no ghost. We construct a low energy effective field theory based on the derivative expansion around the stable background. Possible applications are discussed, including dark energy, dark matter, inflation and black hole.

We investigate how the ghost condensate reacts to black holes immersed in it. A ghost condensate defines a hypersurface-orthogonal congruence of timelike curves, each of which has the tangent vector u^\mu=-g^{\mu\nu}\partial_\nu\phi. It is argued that the ghost condensate in this picture approximately corresponds to a congruence of geodesics. In other words, the ghost condensate accretes into a black hole just like a pressure-less dust. Correspondingly, if the energy density of the ghost condensate at large distance is set to an extremely small value by cosmic expansion then the late-time accretion rate of the ghost condensate should be negligible. The accretion rate remains very small even if effects of higher derivative terms are taken into account, provided that the black hole is sufficiently large. It is also discussed how to reconcile the black hole accretion with the possibility that the ghost condensate might behave like dark matter.

One of the most disturbing difficulties in thinking about the cosmological constant is that it is not stable under radiative corrections. The feedback mechanism proposed in [hep-th/0306108] is a dynamical way to protect a zero or small cosmological constant against radiative corrections. Hence, while this by itself does not solve the cosmological constant problem, it can help solving the problem. In the present paper we investigate stability and gravity in this approach and show that the feedback mechanism is both classically and quantum mechanically stable and has self-consistent, stable dynamics.

We consider a dynamical approach to the cosmological constant. There is a scalar field with a potential whose minimum occurs at a generic, but negative, value for the vacuum energy, and it has a non-standard kinetic term whose coefficient diverges at zero curvature as well as the standard kinetic term. Because of the divergent coefficient of the kinetic term, the lowest energy state is never achieved. Instead, the cosmological constant automatically stalls at or near zero. The merit of this model is that it is stable under radiative corrections and leads to stable dynamics, despite the singular kinetic term. The model is not complete, however, in that some reheating is required. Nonetheless, our approach can at the very least reduce fine-tuning by 60 orders of magnitude or provide a new mechanism for sampling possible cosmological constants and implementing the anthropic principle.

We propose a new scenario for early cosmology, where an inflationary de Sitter phase is obtained with a ghost condensate. The transition to radiation dominance is triggered by the ghost itself, without any slow-roll potential. Density perturbations are generated by fluctuations around the ghost condensate and can be reliably computed in the effective field theory. The fluctuations are scale invariant as a consequence of the de Sitter symmetries, however, the size of the perturbations are parametrically different from conventional slow-roll inflation, and the inflation happens at far lower energy scales. The model makes definite predictions that distinguish it from standard inflation, and can be sharply excluded or confirmed by experiments in the near future. The tilt in the scalar spectrum is predicted to vanish (n_s=1), and the gravity wave signal is negligible. The non-Gaussianities in the spectrum are predicted to be observable: the 3-point function is determined up to an overall order 1 constant, and its magnitude is much bigger than in conventional inflation, with an equivalent f_NL ~ 100, not far from the present WMAP bounds.

We propose a theoretically consistent modification of gravity in the infrared, which is compatible with all current experimental observations. This is an analog of Higgs mechanism in general relativity, and can be thought of as arising from ghost condensation--a background where a scalar field \phi has a constant velocity, <\dot\phi> = M^2. The ghost condensate is a new kind of fluid that can fill the universe, which has the same equation of state, \rho = -p, as a cosmological constant, and can hence drive de Sitter expansion of the universe. However, unlike a cosmological constant, it is a physical fluid with a physical scalar excitation, which can be described by a systematic effective field theory at low energies. The excitation has an unusual low-energy dispersion relation \omega^2 \sim k^4 / M^2. If coupled to matter directly, it gives rise to small Lorentz-violating effects and a new long-range 1/r^2 spin dependent force. In the ghost condensate, the energy that gravitates is not the same as the particle physics energy, leading to the possibility of both sources that can gravitate and antigravitate. The Newtonian potential is modified with an oscillatory behavior starting at the distance scale M_{Pl}/M^2 and the time scale M_{Pl}^2/M^3. This theory opens up a number of new avenues for attacking cosmological problems, including inflation, dark matter and dark energy.

We propose a new, self-consistent and dynamical scenario which gives rise to well-defined initial conditions for five-dimensional brane-world cosmologies with radion stabilization. At high energies, the five-dimensional effective theory is assumed to have a scale invariance so that it admits an expanding scaling solution as a future attractor. The system automatically approaches the scaling solution and, hence, the initial condition for the subsequent low-energy brane cosmology is set by the scaling solution. At low energies, the scale invariance is broken and a radion stabilization mechanism drives the dynamics of the brane-world system. We present an exact, analytic scaling solution for a class of scale-invariant effective theories of five-dimensional brane-world models which includes the five-dimensional reduction of the Horava-Witten theory, and provide convincing evidence that the scaling solution is a future attractor.

We investigate light-cone structure on the world-volume of an unstable D-brane with a tachyon decaying inhomogeneously by using a field theoretical description. It is shown that (i) light-cones governing open strings are narrower than those governing closed strings and will eventually collapse inward in all directions except at kinks, where the tachyon remains at the top of its potential; and that (ii) light-cones governing open strings at a kink will be narrowed only in the direction perpendicular to the kink surface. It is also shown that (iii) future-directed light-cones governing open strings near a kink are tilted towards the kink, compared with those governing closed strings. The result (i) implies that open strings except at kinks are redshifted, compared with closed strings, and will eventually cease to be dynamical. On the other hand, the result (ii) shows that open strings on a kink surface can move freely along the kink surface and are dynamical but do not feel the existence of the spatial dimension perpendicular to the kink surface. The result (iii) indicates that open strings near a kink have tendency to move towards the kink. Hence, the light-cone structure vividly illustrates how open strings behave during the dynamical formation of a kink. We also discuss about a possibility that the early universe has a network of various dimensional D-branes, black-branes and tachyon matter. A problem associated with the network and a possible solution to the problem are discussed.

Brane cosmology driven by the tachyon rolling down to its ground state is investigated. We adopt an effective field theoretical description for the tachyon and Randall-Sundrum type brane world scenario. After formulating basic equations, we show that the standard cosmology with a usual scalar field can mimic the low energy behavior of the system near the tachyon ground state. We also investigate qualitative behavior of the system beyond the low energy regime for positive, negative and vanishing 4-dimensional effective cosmological constant $\Lambda_4=\kappa_5^4V(T_0)^2/12-|\Lambda_5|/2$, where $\kappa_5$ and $\Lambda_5$ are 5-dimensional gravitational coupling constant and (negative) cosmological constant, respectively, and $V(T_0)$ is the (positive) tension of the brane in the tachyon ground state. In particular, for $\Lambda_4<0$ the tachyon never settles down to its potential minimum and the universe eventually hits a big-crunch singularity.

Four dimensional gravity in the low energy limit of a higher dimensional theory has been expected to be a (generalized) Brans-Dicke theory. A subtle point in brane world scenarios is that the system of four dimensional effective gravitational equations is not closed due to bulk gravitational waves and bulk scalars. Nonetheless, weak gravity on the brane can be analyzed completely. We revisit the theory of weak brane gravity using gauge-invariant gravitational and scalar perturbations around a background warped geometry with a bulk scalar between two flat branes. We obtain a simple condition for the radion stabilization in terms of the scalar field potentials. We show that for general potentials of the scalar field which provides radion stabilization and a general conformal transformation to a frame in which matter on the branes are minimally coupled to the metric, 4-dimensional Einstein gravity, not BD gravity, is restored at low energies on either brane. In contrast, in RS brane world scenario without a bulk scalar, low energy gravity is BD one. We conjecture that in general brane world scenarios with more than one scalar field, one will again encounter the situation that low energy gravity is not described by the Einstein theory. Equipped with the weak gravity results, we discuss the properties of 4d brane gravitational equations, in particular, the value and sign of 4d Newton's gravitational coupling.

In brane world scenarios with a bulk scalar field between two branes it is known that 4-dimensional Einstein gravity is restored at low energies on either brane. By using a gauge-invariant gravitational and scalar perturbation formalism we extend the theory of weak gravity in the brane world scenarios to higher energies, or shorter distances. We argue that weak gravity on either brane is indistinguishable from 4-dimensional higher derivative gravity, provided that the inter-brane distance (radion) is stabilized, that the background bulk scalar field is changing near the branes and that the background bulk geometry near the branes is warped. This argument holds for a general conformal transformation to a frame in which matter on the branes is minimally coupled to the metric. In particular, Newton's constant and the coefficients of curvature-squared terms in the 4-dimensional effective action are determined up to an ambiguity of adding a Gauss-Bonnet topological term. In other words, we provide the brane-world realization of the so called $R^2$-model without utilizing a quantum theory. We discuss the appearance of composite spin-2 and spin-0 fields in addition to the graviton on the brane and point out a possibility that the spin-0 field may play the role of an effective inflaton to drive brane-world inflation. Finally, we conjecture that the sequence of higher derivative terms is an infinite series and, thus, indicates non-locality in the brane world scenarios.

In this paper we investigate a parametric resonance phenomenon of a Kaluza-Klein mode in a $D$-dimensional generalized Kaluza-Klein theory. As the origin of the parametric resonance we consider a small oscillation of a scale of the compactification around a today's value of it. To make our arguments definite and for simplicity we consider two classes of models of the compactification: those by $S_{d}$ ($d=D-4$) and those by $S_{d_{1}}\times S_{d_{2}}$ ($d_1\ge d_2$, $d_{1}+d_{2}=D-4$). For these models we show that parametric resonance can occur for the Kaluza-Klein mode. After that, we give formulas of a creation rate and a number of created quanta of the Kaluza-Klein mode due to the parametric resonance, taking into account the first and the second resonance band. By using the formulas we calculate those quantities for each model of the compactification. Finally we give conditions for the parametric resonance to be efficient and discuss cosmological implications.

We consider the horizon problem in a homogeneous but anisotropic universe (Bianchi type I). We show that the problem cannot be solved if (1) the matter obeys the strong energy condition with the positive energy density and (2) the Einstein equations hold. The strong energy condition is violated during cosmological inflation.