The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged)

We argue that the $\Lambda$CDM tensions of the Hubble-Lemaitre expansion rate $H_0$ and the clustering normalization $\sigma_8$ can be eased, at least in principle, by considering an interaction between dark energy and dark matter in such a way to induce a small and positive early effective equation of state and a weaker gravity. For a dark energy scalar field $\phi$ interacting with dark matter through an exchange of both energy and momentum, we derive a general form of the Lagrangian allowing for the presence of scaling solutions. In a subclass of such interacting theories, we show the existence of a scaling $\phi$-matter-dominated-era ($\phi$MDE) which can potentially alleviate the $H_0$ tension by generating an effective high-redshift equation of state. We also study the evolution of perturbations for a model with $\phi$MDE followed by cosmic acceleration and find that the effective gravitational coupling relevant to the linear growth of large-scale structures can be smaller than the Newton gravitational constant $G$ at low redshifts. The momentum exchange between dark energy and dark matter plays a crucial role for realizing weak gravity, while the energy transfer is also required for the existence of $\phi$MDE.

For a scalar field $\phi$ coupled to cold dark matter (CDM), we provide a general framework for studying the background and perturbation dynamics on the isotropic cosmological background. The dark energy sector is described by a Horndeski Lagrangian with the speed of gravitational waves equivalent to that of light, whereas CDM is dealt as a perfect fluid characterized by the number density $n_c$ and four-velocity $u_c^\mu$. For a very general interacting Lagrangian $f(n_c, \phi, X, Z)$, where $f$ depends on $n_c$, $\phi$, $X=-\partial^{\mu} \phi \partial_{\mu} \phi/2$, and $Z=u_c^{\mu} \partial_{\mu} \phi$, we derive the full linear perturbation equations of motion without fixing any gauge conditions. To realize a vanishing CDM sound speed for the successful structure formation, the interacting function needs to be of the form $f=-f_1(\phi, X, Z)n_c+f_2(\phi, X, Z)$. Employing a quasi-static approximation for the modes deep inside the sound horizon, we obtain analytic formulas for the effective gravitational couplings of CDM and baryon density perturbations as well as gravitational and weak lensing potentials. We apply our general formulas to several interacting theories and show that, in many cases, the CDM gravitational coupling around the quasi de-Sitter background can be smaller than the Newton constant $G$ due to a momentum transfer induced by the $Z$-dependence in $f_2$.

In dark energy models where a scalar field $\phi$ is coupled to the Ricci scalar $R$ of the form $e^{-2Q (\phi-\phi_0)/M_{\rm pl}}R$, where $Q$ is a coupling constant, $\phi_0$ is today's value of $\phi$, and $M_{\rm pl}$ is the reduced Planck mass, we study how the recent Lunar Laser Ranging (LLR) experiment places constraints on the nonminimal coupling from the time variation of gravitational coupling. Besides a potential of the light scalar responsible for cosmic acceleration, we take a cubic Galileon term into account to suppress fifth forces in over-density regions of the Universe. Even if the scalar-matter interaction is screened by the Vainshtein mechanism, the time variation of gravitational coupling induced by the cosmological background field $\phi$ survives in the solar system. For a small Galileon coupling constant $\beta_3$, there exists a kinetically driven $\phi$-matter-dominated-epoch ($\phi$MDE) prior to cosmic acceleration. In this case, we obtain the stringent upper limit $Q \le 3.4 \times 10^{-3}$ from the LLR constraint. For a large $\beta_3$ without the $\phi$MDE, the coupling $Q$ is not particularly bounded from above, but the cosmological Vainshtein screening strongly suppresses the time variation of $\phi$ such that the dark energy equation of state $w_{\rm DE}$ reaches the value close to $-1$ at high redshifts. We study the modified gravitational wave propagation induced by the nonminimal coupling to gravity and show that, under the LLR bound, the difference between the gravitational wave and luminosity distances does not exceed the order $10^{-5}$ over the redshift range $0<z<100$. In dark energy models where the Vainshtein mechanism is at work through scalar derivative self-interactions, it is difficult to probe the signature of nonminimal couplings from the observations of standard sirens.

We study a coupled dark energy scenario in which a massive vector field $A_{\mu}$ with broken $U(1)$ gauge symmetry interacts with the four-velocity $u_c^{\mu}$ of cold dark matter (CDM) through the scalar product $Z=-u_c^{\mu} A_{\mu}$. This new coupling corresponds to the momentum transfer, so that the background vector and CDM continuity equations do not have explicit interacting terms analogous to the energy exchange. Hence the observational preference of uncoupled generalized Proca theories over the $\Lambda$CDM model can be still maintained at the background level. Meanwhile, the same coupling strongly affects the evolution of cosmological perturbations. While the effective sound speed of CDM vanishes, the propagation speed and no-ghost condition of a longitudinal scalar of $A_{\mu}$ and the CDM no-ghost condition are subject to nontrivial modifications by the $Z$ dependence in the Lagrangian. We propose a concrete dark energy model and show that the gravitational interaction on scales relevant to the linear growth of large-scale structures can be smaller than the Newton constant at low redshifts. This leads to the suppression of growth rates of both CDM and total matter density perturbations, so our model allows an interesting possibility for reducing the tension of matter density contrast $\sigma_8$ between high- and low-redshift measurements.

We investigate a dark energy scenario in which a canonical scalar field $\phi$ is coupled to the four velocity $u_{c}^{\mu}$ of cold dark matter (CDM) through a derivative interaction $u_{c}^{\mu} \partial_{\mu} \phi$. The coupling is described by an interacting Lagrangian $f(X, Z)$, where $f$ depends on $X=-\partial^{\mu} \phi \partial_{\mu} \phi/2$ and $Z=u_{c}^{\mu} \partial_{\mu} \phi$. We derive stability conditions of linear scalar perturbations for the wavelength deep inside the Hubble radius and show that the effective CDM sound speed is close to 0 as in the standard uncoupled case, while the scalar-field propagation speed is affected by the interacting term $f$. Under a quasi-static approximation, we also obtain a general expression of the effective gravitational coupling felt by the CDM perturbation. We study the late-time cosmological dynamics for the coupling $f \propto X^{(2-m)/2}Z^m$ and show that the gravitational coupling weaker than the Newton constant can be naturally realized for $m>0$ on scales relevant to the growth of large-scale structures. This allows the possibility for alleviating the tension of $\sigma_8$ between low- and high-redshift measurements.

We provide a general framework for studying the dark energy cosmology in which a scalar field $\phi$ is nonminimally and kinetically coupled to Cold Dark Matter (CDM). The scalar-graviton sector is described by the action of Horndeski theories with the speed of gravitational waves equivalent to that of light, whereas CDM is treated as a perfect fluid given by a Schutz-Sorkin action. We consider two interacting Lagrangians of the forms $f_1(\phi,X)\rho_c (n_c)$ and $f_2 (n_c, \phi,X) J_c^{\mu} \partial_{\mu}\phi$, where $X=-\partial^{\mu} \phi \partial_{\mu} \phi/2$, $\rho_c$ and $n_c$ are the energy density and number density of CDM respectively, and $J_c^{\mu}$ is a vector field related to the CDM four velocity. We derive the scalar perturbation equations of motion without choosing any special gauges and identify conditions for the absence of ghosts and Laplacian instabilities on scales deep inside the sound horizon. Applying a quasi-static approximation in a gauge-invariant manner, we also obtain the effective gravitational couplings felt by CDM and baryons for the modes relevant to the linear growth of large-scale structures. In particular, the $n_c$ dependence in the coupling $f_2$ gives rise to an interesting possibility for realizing the gravitational coupling with CDM weaker than the Newton gravitational constant $G$.

In a class of generalized Proca theories, we study the existence of neutron star solutions with a nonvanishing temporal component of the vector field $A_\mu$ approaching 0 toward spatial infinity, as they may be the endpoints of tachyonic instabilities of neutron star solutions in general relativity with $A_{\mu}=0$. Such a phenomenon is called spontaneous vectorization, which is analogous to spontaneous scalarization in scalar-tensor theories with nonminimal couplings to the curvature or matter. For the nonminimal coupling $\beta X R$, where $\beta$ is a coupling constant and $X=-A_{\mu}A^{\mu}/2$, we show that there exist both 0-node and 1-node vector-field solutions, irrespective of the choice of the equations of state of nuclear matter. The 0-node solution, which is present only for $\beta=-{\cal O}(0.1)$, may be induced by some nonlinear effects such as the selected choice of initial conditions. The 1-node solution exists for $\beta=-{\cal O}(1)$, which suddenly emerges above a critical central density of star and approaches the general relativistic branch with the increasing central density. We compute the mass $M$ and radius $r_s$ of neutron stars for some realistic equations of state and show that the $M$-$r_s$ relations of 0-node and 1-node solutions exhibit notable difference from those of scalarized solutions in scalar-tensor theories. Finally, we discuss the possible endpoints of tachyonic instabilities.

We provide a general framework for studying the evolution of background and cosmological perturbations in the presence of a vector field $A_{\mu}$ coupled to cold dark matter (CDM). We consider an interacting Lagrangian of the form $Q f(X) T_c$, where $Q$ is a coupling constant, $f$ is an arbitrary function of $X=-A_{\mu}A^{\mu}/2$, and $T_c$ is a trace of the CDM energy-momentum tensor. The matter coupling affects the no-ghost condition and sound speed of linear scalar perturbations deep inside the sound horizon, while those of tensor and vector perturbations are not subject to modifications. The existence of interactions also modifies the no-ghost condition of CDM density perturbations. We propose a concrete model of coupled vector dark energy with the tensor propagation speed equivalent to that of light. In comparison to the $Q=0$ case, we show that the decay of CDM to the vector field leads to the phantom dark energy equation of state $w_{\rm DE}$ closer to $-1$. This alleviates the problem of observational incompatibility of uncoupled models in which $w_{\rm DE}$ significantly deviates from $-1$. The maximum values of $w_{\rm DE}$ reached during the matter era are bounded from the CDM no-ghost condition of future de Sitter solutions.

In $f(R)$ gravity and Brans-Dicke theory with scalar potentials, we study the structure of neutron stars on a spherically symmetric and static background for two equations of state: SLy and FPS. In massless BD theory, the presence of a scalar coupling $Q$ with matter works to change the star radius in comparison to General Relativity, while the maximum allowed mass of neutron stars is hardly modified for both SLy and FPS equations of state. In Brans-Dicke theory with the massive potential $V(\phi)=m^2 \phi^2/2$, where $m^2$ is a positive constant, we show the difficulty of realizing neutron star solutions with a stable field profile due to the existence of an exponentially growing mode outside the star. As in $f(R)$ gravity with the $R^2$ term, this property is related to the requirement of extra boundary conditions of the field at the surface of star. For the self-coupling potential $V(\phi)=\lambda \phi^4/4$, this problem can be circumvented by the fact that the second derivative $V_{,\phi \phi}=3\lambda\phi^2$ approaches 0 at spatial infinity. In this case, we numerically show the existence of neutron star solutions for both SLy and FPS equations of state and discuss how the mass-radius relation is modified as compared to General Relativity.

We place observational constraints on the Galileon ghost condensate model, a dark energy proposal in cubic-order Horndeski theories consistent with the gravitational-wave event GW170817. The model extends the covariant Galileon by taking an additional higher-order field derivative $X^2$ into account. This allows for the dark energy equation of state $w_{\rm DE}$ to access the region $-2<w_{\rm DE}<-1$ without ghosts. Indeed, this peculiar evolution of $w_{\rm DE}$ is favored over that of the cosmological constant $\Lambda$ from the joint data analysis of cosmic microwave background (CMB) radiation, baryonic acoustic oscillations (BAOs), supernovae type Ia (SNIa) and redshift-space distortions (RSDs). Furthermore, our model exhibits a better compatibility with the CMB data over the $\Lambda$-cold-dark-matter ($\Lambda$CDM) model by suppressing large-scale temperature anisotropies. The CMB temperature and polarization data lead to an estimation for today's Hubble parameter $H_0$ consistent with its direct measurements at 2$\sigma$. We perform a model selection analysis by using several methods and find a statistically significant preference of the Galileon ghost condensate model over $\Lambda$CDM.

We study observational constraints on a specific dark energy model in the framework of Gleyzes-Langlois-Piazza-Vernizzi theories, which extends the Galileon ghost condensate (GGC) to the domain of beyond Horndeski theories. In this model, we show that the Planck cosmic microwave background (CMB) data, combined with datasets of baryon acoustic oscillations, supernovae type Ia, and redshift-space distortions, give the tight upper bound $|\alpha_{\rm H}^{(0)}| \le {\cal O}(10^{-6})$ on today's beyond-Horndeski (BH) parameter $\alpha_{\rm H}$. This is mostly attributed to the shift of CMB acoustic peaks induced by the early-time changes of cosmological background and perturbations arising from the dominance of $\alpha_{\rm H}$ in the dark energy density. In comparison to the $\Lambda$-cold-dark-matter ($\Lambda$CDM) model, our BH model suppresses the large-scale integrated-Sachs-Wolfe (ISW) tail of CMB temperature anisotropies due to the existence of cubic Galileons, and it modifies the small-scale CMB power spectrum because of the different background evolution. We find that the BH model considered fits the data better than $\Lambda$CDM according to the $\chi^2$ statistics, yet the deviance information criterion (DIC) slightly favors the latter. Given the fact that our BH model with $\alpha_{\rm H}=0$ (i.e., the GGC model) is favored over $\Lambda$CDM even by the DIC, there are no particular signatures for the departure from Horndeski theories in current observations.

The gravitational-wave event GW170817 from a binary neutron star merger together with the electromagnetic counterpart showed that the speed of gravitational waves $c_t$ is very close to that of light for the redshift $z<0.009$. This places tight constraints on dark energy models constructed in the framework of modified gravitational theories. We review models of the late-time cosmic acceleration in scalar-tensor theories with second-order equations of motion (dubbed Horndeski theories) by paying particular attention to the evolution of dark energy equation of state and observables relevant to the cosmic growth history. We provide a gauge-ready formulation of scalar perturbations in full Horndeski theories and estimate observables associated with the evolution of large-scale structures, cosmic microwave background, and weak lensing by employing a so-called quasi-static approximation for the modes deep inside the sound horizon. In light of the recent observational bound of $c_t$, we also classify surviving dark energy models into four classes depending on different structure-formation patterns and discuss how they can be observationally distinguished from each other. In particular, the nonminimally coupled theories in which the scalar field $\phi$ has a coupling with the Ricci scalar $R$ of the form $G_4(\phi) R$, including $f(R)$ gravity, can be tightly constrained not only from the cosmic expansion and growth histories but also from the variation of screened gravitational couplings. The cross correlation of integrated Sachs-Wolfe signal with galaxy distributions can be a key observable for placing bounds on the relative ratio of cubic Galileon density to total dark energy density. The dawn of gravitational-wave astronomy will open up a new window to constrain nonminimally coupled theories further by the modified luminosity distance of tensor perturbations.

The gravitational-wave event GW170817, together with the electromagnetic counterpart, shows that the speed of tensor perturbations $c_T$ on the cosmological background is very close to that of light $c$ for the redshift $z<0.009$. In generalized Proca theories, the Lagrangians compatible with the condition $c_T=c$ are constrained to be derivative interactions up to cubic order, besides those corresponding to intrinsic vector modes. We place observational constraints on a dark energy model in cubic-order generalized Proca theories with intrinsic vector modes by running the Markov chain Monte Carlo (MCMC) code. We use the cross-correlation data of the integrated Sachs-Wolfe (ISW) signal and galaxy distributions in addition to the data sets of cosmic microwave background, baryon acoustic oscillations, type Ia supernovae, local measurements of the Hubble expansion rate, and redshift-space distortions. We show that, unlike cubic-order scalar-tensor theories, the existence of intrinsic vector modes allows the possibility for evading the ISW-galaxy anticorrelation incompatible with the current observational data. As a result, we find that the dark energy model in cubic-order generalized Proca theories exhibits a better fit to the data than the cosmological constant, even by including the ISW-galaxy correlation data in the MCMC analysis.

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

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

In quadratic-order degenerate higher-order scalar-tensor (DHOST) theories compatible with gravitational-wave constraints, we derive the most general Lagrangian allowing for tracker solutions characterized by $\dot{\phi}/H^p={\rm constant}$, where $\dot{\phi}$ is the time derivative of a scalar field $\phi$, $H$ is the Hubble expansion rate, and $p$ is a constant. While the tracker is present up to the cubic-order Horndeski Lagrangian $L=c_2X-c_3X^{(p-1)/(2p)} \square \phi$, where $c_2, c_3$ are constants and $X$ is the kinetic energy of $\phi$, the DHOST interaction breaks this structure for $p \neq 1$. Even in the latter case, however, there exists an approximate tracker solution in the early cosmological epoch with the nearly constant field equation of state $w_{\phi}=-1-2p\dot{H}/(3H^2)$. The scaling solution, which corresponds to $p=1$, is the unique case in which all the terms in the field density $\rho_{\phi}$ and the pressure $P_{\phi}$ obey the scaling relation $\rho_{\phi} \propto P_{\phi} \propto H^2$. Extending the analysis to the coupled DHOST theories with the field-dependent coupling $Q(\phi)$ between the scalar field and matter, we show that the scaling solution exists for $Q(\phi)=1/(\mu_1 \phi+\mu_2)$, where $\mu_1$ and $\mu_2$ are constants. For the constant $Q$, i.e., $\mu_1=0$, we derive fixed points of the dynamical system by using the general Lagrangian with scaling solutions. This result can be applied to the model construction of late-time cosmic acceleration preceded by the scaling $\phi$-matter-dominated epoch.

In cubic-order Horndeski theories where a scalar field $\phi$ is coupled to nonrelativistic matter with a field-dependent coupling $Q(\phi)$, we derive the most general Lagrangian having scaling solutions on the isotropic and homogenous cosmological background. For constant $Q$ including the case of vanishing coupling, the corresponding Lagrangian reduces to the form $L=Xg_2(Y)-g_3(Y)\square \phi$, where $X=-\partial_{\mu}\phi\partial^{\mu}\phi/2$ and $g_2, g_3$ are arbitrary functions of $Y=Xe^{\lambda \phi}$ with constant $\lambda$. We obtain the fixed points of the scaling Lagrangian for constant $Q$ and show that the $\phi$-matter-dominated-epoch ($\phi$MDE) is present for the cubic coupling $g_3(Y)$ containing inverse power-law functions of $Y$. The stability analysis around the fixed points indicates that the $\phi$MDE can be followed by a stable critical point responsible for the cosmic acceleration. We propose a concrete dark energy model allowing for such a cosmological sequence and show that the ghost and Laplacian instabilities can be avoided even in the presence of the cubic coupling.

In scalar-vector-tensor (SVT) theories with parity invariance, we perform a gauge-ready formulation of cosmological perturbations on the flat Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW) background by taking into account a matter perfect fluid. We derive the second-order action of scalar perturbations and resulting linear perturbation equations of motion without fixing any gauge conditions. Depending on physical problems at hand, most convenient gauges can be chosen to study the development of inhomogeneities in the presence of scalar and vector fields coupled to gravity. This versatile framework, which encompasses Horndeski and generalized Proca theories as special cases, is applicable to a wide variety of cosmological phenomena including nonsingular cosmology, inflation, and dark energy. By deriving conditions for the absence of ghost and Laplacian instabilities in several different gauges, we show that, unlike Horndeski theories, it is possible to evade no-go arguments for the absence of stable nonsingular bouncing/genesis solutions in both generalized Proca and SVT theories. We also apply our framework to the case in which scalar and vector fields are responsible for dark energy and find that the separation of observables relevant to the evolution of matter perturbations into tensor, vector, and scalar sectors is transparent in the unitary gauge. Unlike the flat gauge chosen in the literature, this result is convenient to confront SVT theories with observations associated with the cosmic growth history.

In the framework of effective field theories with prominent helicity-0 and helicity-1 fields coupled to each other via a dimension-3 operator, we study the dynamics of inflation driven by the helicity-0 mode, with a given potential energy, as well as the evolution of cosmological perturbations, influenced by the presence of a mixing term between both helicities. In this scenario, the temporal component of the helicity-1 mode is an auxiliary field and can be integrated out in terms of the time derivative of the helicity-0 mode, so that the background dynamics effectively reduces to that in single-field inflation modulated by a parameter $\beta$ associated to the coupling between helicity-0 and helicity-1 modes. We discuss the evolution of a longitudinal scalar perturbation $\psi$ and an inflaton fluctuation $\delta \phi$, and explicitly show that a particular combination of these two, which corresponds to an isocurvature mode, is subject to exponential suppression by the vector mass comparable to the Hubble expansion rate during inflation. Furthermore, we find that the effective single-field description corrected by $\beta$ also holds for the power spectrum of curvature perturbations generated during inflation. We compute the standard inflationary observables such as the scalar spectral index $n_s$ and the tensor-to-scalar ratio $r$ and confront several inflaton potentials with the recent observational data provided by Planck 2018. Our results show that the coupling between helicity-0 and helicity-1 modes can lead to a smaller value of the tensor-to-scalar ratio especially for small-field inflationary models, so our scenario exhibits even better compatibility with the current observational data.

The scalar-vector-tensor theories with second-order equations of motion can accommodate both Horndeski and generalized Proca theories as specific cases. In the presence of a perfect fluid, we study the cosmology in such a most general scheme of scalar-vector-tensor theories with parity invariance by paying particular attention to the application to dark energy. We obtain a closed-form expression of the background equations of motion by using coefficients appearing in the second-order action of scalar perturbations. We also derive conditions for the absence of ghost and Laplacian instabilities of tensor and vector perturbations and show that the existence of matter does not substantially modify the stabilities of dynamical degrees of freedom in the small-scale limit. On the other hand, the sound speed of scalar perturbations is affected by the presence of matter. Employing the quasi-static approximation for scalar perturbations deep inside the sound horizon, we derive analytic expressions of Newtonian and weak lensing gravitational potentials as well as two scalar perturbations arising from the scalar and vector fields. We apply our general framework to dark energy theories with the tensor propagation speed equivalent to the speed of light and show that the observables associated with the growth of matter perturbations and weak lensing potentials are generally affected by intrinsic vector modes and by interactions between scalar and vector fields.

We propose a viable dark energy scenario in the presence of cubic Horndeski interactions and a standard scalar-field kinetic term with two exponential potentials. We show the existence of new scaling solutions along which the cubic coupling $G_3$ provides an important contribution to the field density that scales in the same way as the background fluid density. The solutions finally exit to the epoch of cosmic acceleration driven by a scalar-field dominated fixed point arising from the second exponential potential. We clarify the viable parameter space in which all the theoretically consistent conditions including those for the absence of ghost and Laplacian instabilities are satisfied on scaling and scalar-field dominated critical points. In comparison to Quintessence with the same scalar potential, we find that the cubic coupling gives rise to some novel features: (i) the allowed model parameter space is wider in that a steeper potential can drive the cosmic acceleration; (ii) the dark energy equation of state $w_{\phi}$ today can be closer to $-1$ relative to Quintessence; (iii) even if the density associated with the cubic coupling dominates over the standard field density in the scaling era, the former contribution tends to be suppressed at low redshifts. We also compute quantities associated with the growth of matter perturbations and weak lensing potentials under the quasi-static approximation in the sub-horizon limit and show that the cubic coupling leads to the modified evolution of perturbations which can be distinguished from Quintessence.

We study the cosmology on the Friedmann-Lemaitre-Robertson-Walker background in scalar-vector-tensor theories with a broken $U(1)$ gauge symmetry. For parity-invariant interactions arising in scalar-vector-tensor theories with second-order equations of motion, we derive conditions for the absence of ghosts and Laplacian instabilities associated with tensor, vector, and scalar perturbations at linear order. This general result is applied to the computation of the primordial tensor power spectrum generated during inflation as well as to the speed of gravity relevant to dark energy. We also construct a concrete inflationary model in which a temporal vector component $A_0$ contributes to the dynamics of cosmic acceleration besides a scalar field $\phi$ through their kinetic mixings. In this model, we show that all the stability conditions of perturbations can be consistently satisfied during inflation and subsequent reheating.

In scalar-vector-tensor theories with $U(1)$ gauge invariance, it was recently shown that there exists a new type of hairy black hole (BH) solutions induced by a cubic-order scalar-vector interaction. In this paper, we derive conditions for the absence of ghosts and Laplacian instabilities against odd-parity perturbations on a static and spherically symmetric background for most general $U(1)$ gauge-invariant scalar-vector-tensor theories with second-order equations of motion. We apply those conditions to hairy BH solutions arising from the cubic-order coupling and show that the odd-parity stability in the gravity sector is always ensured outside the event horizon with the speed of gravity equivalent to that of light. We also study the case in which quartic-order interactions are present in addition to the cubic coupling and obtain conditions under which black holes are stable against odd-parity perturbations.

The generalized Proca theories with second-order equations of motion can be healthily extended to a more general framework in which the number of propagating degrees of freedom remains unchanged. In the presence of a quartic-order nonminimal coupling to gravity arising in beyond-generalized Proca theories, the speed of gravitational waves $c_t$ on the Friedmann-Lemaitre-Robertson-Walker (FLRW) cosmological background can be equal to that of light $c$ under a certain condition. By using this condition alone, we show that the speed of gravitational waves in the vicinity of static and spherically symmetric black holes is also equivalent to $c$ for the propagation of odd-parity perturbations along both radial and angular directions. As a by-product, the black holes arising in our beyond-generalized Proca theories are plagued by neither ghost nor Laplacian instabilities against odd-parity perturbations. We show the existence of both exact and numerical black hole solutions endowed with vector hairs induced by the quartic-order coupling.

The Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories up to quartic order are the general scheme of scalar-tensor theories allowing the possibility for realizing the tensor propagation speed $c_t$ equivalent to 1 on the isotropic cosmological background. We propose a dark energy model in which the late-time cosmic acceleration occurs by a simple k-essence Lagrangian analogous to the ghost condensate with cubic and quartic Galileons in the framework of GLPV theories. We show that a wide variety of the variation of the dark energy equation of state $w_{\rm DE}$ including the entry to the region $w_{\rm DE}<-1$ can be realized without violating conditions for the absence of ghosts and Laplacian instabilities. The approach to the tracker equation of state $w_{\rm DE}=-2$ during the matter era, which is disfavored by observational data, can be avoided by the existence of a quadratic k-essence Lagrangian $X^2$. We study the evolution of nonrelativistic matter perturbations for the model $c_t^2=1$ and show that the two quantities $\mu$ and $\Sigma$, which are related to the Newtonian and weak lensing gravitational potentials respectively, are practically equivalent to each other, such that $\mu \simeq \Sigma>1$. For the case in which the deviation of $w_{\rm DE}$ from $-1$ is significant at a later cosmological epoch, the values of $\mu$ and $\Sigma$ tend to be larger at low redshifts. We also find that our dark energy model can be consistent with the bounds on the deviation parameter $\alpha_{\rm H}$ from Horndeski theories arising from the modification of gravitational law inside massive objects.

We study relativistic star solutions in second-order generalized Proca theories characterized by a $U(1)$-breaking vector field with derivative couplings. In the models with cubic and quartic derivative coupling, the mass and radius of stars become larger than those in general relativity for negative derivative coupling constants. This phenomenon is mostly attributed to the increase of star radius induced by a slower decrease of the matter pressure compared to general relativity. There is a tendency that the relativistic star with a smaller mass is not gravitationally bound for a low central density and hence dynamically unstable, but that with a larger mass is gravitationally bound. On the other hand, we show that the intrinsic vector-mode couplings give rise to general relativistic solutions with a trivial field profile, so the mass and radius are not modified from those in general relativity.

In generalized Proca theories with vector-field derivative couplings, a bunch of hairy black hole solutions have been derived on a static and spherically symmetric background. In this paper, we formulate the odd-parity black hole perturbations in generalized Proca theories by expanding the corresponding action up to second order and investigate whether or not black holes with vector hair suffer ghost or Laplacian instabilities. We show that the models with cubic couplings $G_3(X)$, where $X=-A_{\mu}A^{\mu}/2$ with a vector field $A_{\mu}$, do not provide any additional stability condition as in General Relativity. On the other hand, the exact charged stealth Schwarzschild solution with a nonvanishing longitudinal vector component $A_1$, which originates from the coupling to the Einstein tensor $G^{\mu\nu}A_\mu A_\nu$ equivalent to the quartic coupling $G_4(X)$ containing a linear function of $X$, is unstable in the vicinity of the event horizon. The same instability problem also persists for hairy black holes arising from general quartic power-law couplings $G_4(X) \supset \beta_4 X^n$ with the nonvanishing $A_1$, while the other branch with $A_1=0$ can be consistent with conditions for the absence of ghost and Laplacian instabilities. We also discuss the case of other exact and numerical black hole solutions associated with intrinsic vector-field derivative couplings and show that there exists a wide range of parameter spaces in which the solutions suffer neither ghost nor Laplacian instabilities against odd-parity perturbations.

We study the evolution of Gravitational Waves (GWs) during and after inflation as well as the resulting observational consequences in a Lorentz-violating massive gravity theory with one scalar (inflaton) and two tensor degrees of freedom. We consider two explicit examples of the tensor mass $m_g$ that depends either on the inflaton field $\phi$ or on its time derivative $\dot{\phi}$, both of which lead to parametric excitations of GWs during reheating after inflation. The first example is Starobinsky's $R^2$ inflation model with a $\phi$-dependent $m_g$ and the second is a low-energy-scale inflation model with a $\dot{\phi}$-dependent $m_g$. We compute the energy density spectrum $\Omega_{\rm GW}(k)$ today of the GW background. In the Starobinsky's model, we show that the GWs can be amplified up to the detectable ranges of both CMB and DECIGO, but the bound from the big bang nucleosynthesis is quite tight to limit the growth. In low-scale inflation with a fast transition to the reheating stage driven by the potential $V(\phi)=M^2 \phi^2/2$ around $\phi \approx M_{\rm pl}$ (where $M_{\rm pl}$ is the reduced Planck mass), we find that the peak position of $\Omega_{\rm GW}(k)$ induced by the parametric resonance can reach the sensitivity region of advanced LIGO for the Hubble parameter of order 1 GeV at the end of inflation. Thus, our massive gravity scenario offers exciting possibilities for probing the physics of primordial GWs at various different frequencies.

The recent detection of the gravitational wave signal GW170817 together with an electromagnetic counterpart GRB 170817A from the merger of two neutron stars puts a stringent bound on the tensor propagation speed. This constraint can be automatically satisfied in the framework of massive gravity. In this work we consider a general $SO(3)$-invariant massive gravity with five propagating degrees of freedom and derive the conditions for the absence of ghosts and Laplacian instabilities in the presence of a matter perfect fluid on the flat Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW) cosmological background. The graviton potential containing the dependence of three-dimensional metrics and a fiducial metric coupled to a temporal scalar field gives rise to a scenario of the late-time cosmic acceleration in which the dark energy equation of state $w_{\rm DE}$ is equivalent to $-1$ or varies in time. We find that the deviation from the value $w_{\rm DE}=-1$ provides important contributions to the quantities associated with the stability conditions of tensor, vector, and scalar perturbations. In concrete models, we study the dynamics of dark energy arising from the graviton potential and show that there exist viable parameter spaces in which neither ghosts nor Laplacian instabilities are present for both $w_{\rm DE}>-1$ and $w_{\rm DE}<-1$. We also generally obtain the effective gravitational coupling $G_{\rm eff}$ with non-relativistic matter as well as the gravitational slip parameter $\eta_s$ associated with the observations of large-scale structures and weak lensing. We show that, apart from a specific case, the two quantities $G_{\rm eff}$ and $\eta_s$ are similar to those in general relativity for scalar perturbations deep inside the sound horizon.

We present a family of exact black-hole solutions on a static spherically symmetric background in second-order generalized Proca theories with derivative vector-field interactions coupled to gravity. We also derive non-exact solutions in power-law coupling models including vector Galileons and numerically show the existence of regular black holes with a primary hair associated with the longitudinal propagation. The intrinsic vector-field derivative interactions generally give rise to a secondary hair induced by non-trivial field profiles. The deviation from General Relativity is most significant around the horizon and hence there is a golden opportunity for probing the Proca hair by the measurements of gravitational waves (GWs) in the regime of strong gravity.

We derive the profile of a vector field coupled to matter on a static and spherically symmetric background in the context of generalized Proca theories. The cubic Galileon self-interaction leads to the suppression of a longitudinal vector component due to the operation of the Vainshtein mechanism. For quartic and sixth-order derivative interactions, the solutions consistent with those in the continuous limit of small derivative couplings correspond to the branch with the vanishing longitudinal mode. We compute the corrections to gravitational potentials outside a compact body induced by the vector field in the presence of cubic, quartic, and sixth-order derivative couplings, and show that the models can be consistent with local gravity constraints under mild bounds on the temporal vector component. The quintic Galileon interaction does not allow regular solutions of the longitudinal mode for a rapidly decreasing matter density outside the body.

In a model of the late-time cosmic acceleration within the framework of generalized Proca theories, there exists a de Sitter attractor preceded by the dark energy equation of state $w_{\rm DE}=-1-s$, where $s$ is a positive constant. We run the Markov-Chain-Monte-Carlo code to confront the model with the observational data of Cosmic Microwave Background (CMB), baryon acoustic oscillations, supernovae type Ia, and local measurements of the Hubble expansion rate for the background cosmological solutions and obtain the bound $s=0.254^{{}+ 0.118}_{{}-0.097}$ at 95% confidence level (CL). Existence of the additional parameter $s$ to those in the $\Lambda$-Cold-Dark-Matter ($\Lambda$CDM) model allows to reduce tensions of the Hubble constant $H_0$ between the CMB and the low-redshift measurements. Including the cosmic growth data of redshift-space distortions in the galaxy power spectrum and taking into account no-ghost and stability conditions of cosmological perturbations, we find that the bound on $s$ is shifted to $s=0.16^{+0.08}_{-0.08}$ (95 % CL) and hence the model with $s>0$ is still favored over the $\Lambda$CDM model. Apart from the quantities $s, H_0$ and the today's matter density parameter $\Omega_{m0}$, the constraints on other model parameters associated with perturbations are less stringent, reflecting the fact that there are different sets of parameters that give rise to a similar cosmic expansion and growth history.

We study static and spherically symmetric black hole (BH) solutions in second-order generalized Proca theories with nonminimal vector field derivative couplings to the Ricci scalar, the Einstein tensor, and the double dual Riemann tensor. We find concrete Lagrangians which give rise to exact BH solutions by imposing two conditions of the two identical metric components and the constant norm of the vector field. These exact solutions are described by either Reissner-Nordstr\"{o}m (RN), stealth Schwarzschild, or extremal RN solutions with a non-trivial longitudinal mode of the vector field. We then numerically construct BH solutions without imposing these conditions. For cubic and quartic Lagrangians with power-law couplings which encompass vector Galileons as the specific cases, we show the existence of BH solutions with the difference between two non-trivial metric components. The quintic-order power-law couplings do not give rise to non-trivial BH solutions regular throughout the horizon exterior. The sixth-order and intrinsic vector-mode couplings can lead to BH solutions with a secondary hair. For all the solutions, the vector field is regular at least at the future or past horizon. The deviation from General Relativity induced by the Proca hair can be potentially tested by future measurements of gravitational waves in the nonlinear regime of gravity.

The beyond-generalized Proca theories are the extension of second-order massive vector-tensor theories (dubbed generalized Proca theories) with two transverse vector modes and one longitudinal scalar besides two tensor polarizations. Even with this extension, the propagating degrees of freedom remain unchanged on the isotropic cosmological background without an Ostrogradski instability. We study the cosmology in beyond-generalized Proca theories by paying particular attention to the dynamics of late-time cosmic acceleration and resulting observational consequences. We derive conditions for avoiding ghosts and instabilities of tensor, vector, and scalar perturbations and discuss viable parameter spaces in concrete models allowing the dark energy equation of state smaller than $-1$. The propagation speeds of those perturbations are subject to modifications beyond the domain of generalized Proca theories. There is a mixing between scalar and matter sound speeds, but such a mixing is suppressed during most of the cosmic expansion history without causing a new instability. On the other hand, we find that derivative interactions arising in beyond-generalized Proca theories give rise to important modifications to the cosmic growth history. The growth rate of matter perturbations can be compatible with the redshift-space distortion data due to the realization of gravitational interaction weaker than that in generalized Proca theories. Thus, it is possible to distinguish the dark energy model in beyond-generalized Proca theories from the counterpart in generalized Proca theories as well as from the $\Lambda$CDM model.

A non-abelian $SU(2)$ gauge field with a non-minimal Horndeski coupling to gravity gives rise to a de Sitter solution followed by a graceful exit to a radiation-dominated epoch. In this Horndeski Yang-Mills (HYM) theory we derive the second-order action for tensor perturbations on the homogeneous and isotropic quasi de Sitter background. We find that the presence of the Horndeski non-minimal coupling to the gauge field inevitably introduces ghost instabilities in the tensor sector during inflation. Moreover, we also find Laplacian instabilities for the tensor perturbations deep inside the Hubble radius during inflation. Thus, we conclude that the HYM theory does not provide a consistent inflationary framework due to the presence of ghosts and Laplacian instabilities.

In Gleyzes-Langlois-Piazza-Vernizzi (GLPV) scalar-tensor theories, which are outside the domain of second-order Horndeski theories, it is known that there exists a solid angle deficit singularity in the case where the parameter $\alpha_{\rm H}$ characterizing the deviation from Horndeski theories approaches a non-vanishing constant at the center of a spherically symmetric body. Meanwhile, it was recently shown that second-order generalized Proca theories with a massive vector field $A^{\mu}$ can be consistently extended to beyond-generalized Proca theories, which recover shift-symmetric GLPV theories in the scalar limit $A^{\mu} \to \nabla^{\mu} \chi$. In beyond-generalized Proca theories up to quartic-order Lagrangians, we show that solid angle deficit singularities are generally absent due to the existence of a temporal vector component. We also derive the vector-field profiles around a compact object and show that the success of the Vainshtein mechanism operated by vector Galileons is not prevented by new interactions in beyond generalized Proca theories.

We use FIRAS and Planck 2015 data to place observational bounds on inflationary scenarios in multi-fractional spacetimes with $q$-derivatives. While a power-law expansion in the geometric time coordinate is subject to the usual constraints from the tensor-to-scalar ratio, model-independent best fits of the black-body and scalar spectra yield upper limits on the free parameters of the multi-fractal measure of the theory. When the measure describing the fractal spacetime geometry is non-oscillating, information on the CMB black-body spectrum places constraints on the theory independent from but weaker than those obtained from the Standard Model, astrophysical gravitational waves and gamma-ray bursts (GRBs). When log oscillations are included and the measure describes a discrete fractal spacetime at microscopic scales, we obtain the first observational constraints on the amplitudes of such oscillations and find, in general, strong constraints on the multi-scale geometry and on the dimension of space. These results complete the scan and reduction of the parameter space of the theory. Black-body bounds are obtained also for the theory with weighted derivatives.

We consider higher-order derivative interactions beyond second-order generalized Proca theories that propagate only the three desired polarizations of a massive vector field besides the two tensor polarizations from gravity. These new interactions follow the similar construction criteria to those arising in the extension of scalar-tensor Horndeski theories to Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories. On the isotropic cosmological background, we show the existence of a constraint with a vanishing Hamiltonian that removes the would-be Ostrogradski ghost. We study the behavior of linear perturbations on top of the isotropic cosmological background in the presence of a matter perfect fluid and find the same number of propagating degrees of freedom as in generalized Proca theories (two tensor polarizations, two transverse vector modes, and two scalar modes). Moreover, we obtain the conditions for the avoidance of ghosts and Laplacian instabilities of tensor, vector, and scalar perturbations. We observe key differences in the scalar sound speed, which is mixed with the matter sound speed outside the domain of generalized Proca theories.

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.

In beyond-generalized Proca theories including the extension to theories higher than second order, we study the role of a spatial component $v$ of a massive vector field on the anisotropic cosmological background. We show that, as in the case of the isotropic cosmological background, there is no additional ghostly degrees of freedom associated with the Ostrogradski instability. In second-order generalized Proca theories we find the existence of anisotropic solutions on which the ratio between the anisotropic expansion rate $\Sigma$ and the isotropic expansion rate $H$ remains nearly constant in the radiation-dominated epoch. In the regime where $\Sigma/H$ is constant, the spatial vector component $v$ works as a dark radiation with the equation of state close to $1/3$. During the matter era, the ratio $\Sigma/H$ decreases with the decrease of $v$. As long as the conditions $|\Sigma| \ll H$ and $v^2 \ll \phi^2$ are satisfied around the onset of late-time cosmic acceleration, where $\phi$ is the temporal vector component, we find that the solutions approach the isotropic de Sitter fixed point ($\Sigma=0=v$) in accordance with the cosmic no-hair conjecture. In the presence of $v$ and $\Sigma$ the early evolution of the dark energy equation of state $w_{\rm DE}$ in the radiation era is different from that in the isotropic case, but the approach to the isotropic value $w_{\rm DE}^{{\rm (iso)}}$ typically occurs at redshifts $z$ much larger than 1. Thus, apart from the existence of dark radiation, the anisotropic cosmological dynamics at low redshifts is similar to that in isotropic generalized Proca theories. In beyond-generalized Proca theories the only consistent solution to avoid the divergence of a determinant of the dynamical system corresponds to $v=0$, so $\Sigma$ always decreases in time.

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.

For a massive vector field with derivative self-interactions, the breaking of the gauge invariance allows the propagation of a longitudinal mode in addition to the two transverse modes. We consider generalized Proca theories with second-order equations of motion in a curved space-time and study how the longitudinal scalar mode of the vector field gravitates on a spherically symmetric background. We show explicitly that cubic-order self-interactions lead to the suppression of the longitudinal mode through the Vainshtein mechanism. Provided that the dimensionless coupling of the interaction is not negligible, this screening mechanism is sufficiently efficient to give rise to tiny corrections to gravitational potentials consistent with solar-system tests of gravity. We also study the quartic interactions with the presence of non-minimal derivative coupling with the Ricci scalar and find the existence of solutions where the longitudinal mode completely vanishes. Finally, we discuss the case in which the effect of the quartic interactions dominates over the cubic one and show that local gravity constraints can be satisfied under a mild bound on the parameters of the theory.

We measure the redshift-space correlation function from a spectroscopic sample of 2783 emission line galaxies from the FastSound survey. The survey, which uses the Subaru Telescope and covers the redshift ranges of $1.19<z<1.55$, is the first cosmological study at such high redshifts. We detect clear anisotropy due to redshift-space distortions (RSD) both in the correlation function as a function of separations parallel and perpendicular to the line of sight and its quadrupole moment. RSD has been extensively used to test general relativity on cosmological scales at $z<1$. Adopting a LCDM cosmology with the fixed expansion history and no velocity dispersion $\sigma_{\rm v}=0$, and using the RSD measurements on scales above 8Mpc/h, we obtain the first constraint on the growth rate at the redshift, $f(z)\sigma_8(z)=0.482\pm 0.116$ at $z\sim 1.4$ after marginalizing over the galaxy bias parameter $b(z)\sigma_8(z)$. This corresponds to $4.2\sigma$ detection of RSD. Our constraint is consistent with the prediction of general relativity $f\sigma_8\sim 0.392$ within the $1-\sigma$ confidence level. When we allow $\sigma_{\rm v}$ to vary and marginalize it over, the growth rate constraint becomes $f\sigma_8=0.494^{+0.126}_{-0.120}$. We also demonstrate that by combining with the low-z constraints on $f\sigma_8$, high-z galaxy surveys like the FastSound can be useful to distinguish modified gravity models without relying on CMB anisotropy experiments.

In Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories, it is known that the conical singularity arises at the center of a spherically symmetric body ($r=0$) in the case where the parameter $\alpha_{{\rm H}4}$ characterizing the deviation from the Horndeski Lagrangian $L_4$ approaches a non-zero constant as $r \to 0$. We derive spherically symmetric solutions around the center in full GLPV theories and show that the GLPV Lagrangian $L_5$ does not modify the divergent property of the Ricci scalar $R$ induced by the non-zero $\alpha_{{\rm H}4}$. Provided that $\alpha_{{\rm H}4}=0$, curvature scalar quantities can remain finite at $r=0$ even in the presence of $L_5$ beyond the Horndeski domain. For the theories in which the scalar field $\phi$ is directly coupled to $R$, we also obtain spherically symmetric solutions inside/outside the body to study whether the fifth force mediated by $\phi$ can be screened by non-linear field self-interactions. We find that there is one specific model of GLPV theories in which the effect of $L_5$ vanishes in the equations of motion. We also show that, depending on the sign of a $L_5$-dependent term in the field equation, the model can be compatible with solar-system constraints under the Vainshtein mechanism or it is plagued by the problem of a divergence of the field derivative in high-density regions.

We propose a dark energy model where a scalar field $\phi$ has nonlinear self-interactions in the presence of a dilatonic coupling with the Ricci scalar. This belongs to a sub-class of theories beyond Horndeski, which accommodates covariant Galileons and Brans-Dicke theories as specific cases. We derive conditions under which the scalar sound speed squared $c_{\rm s}^2$ is positive from the radiation era to today. Since $c_{\rm s}^2$ remains to be smaller than the order of 1, the deviation from Horndeski theories does not cause heavy oscillations of gauge-invariant gravitational potentials. In this case, the evolution of matter perturbations at low redshifts is similar to that in the coupled dark energy scenario with an enhanced gravitational interaction. On the spherically symmetric background with a matter source, the existence of field self-interactions suppresses the propagation of fifth force inside a Vainshtein radius. We estimate an allowed parameter space in which the model can be compatible with solar system constraints while satisfying conditions for the cosmological viability of background and perturbations.

In a class of Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories, we derive both vacuum and interior Schwarzschild solutions under the condition that the derivatives of a scalar field $\phi$ with respect to the radius $r$ vanish. If the parameter $\alpha_{\rm H}$ characterizing the deviation from Horndeski theories approaches a non-zero constant at the center of a spherically symmetric body, we find that the conical singularity arises at $r=0$ with the Ricci scalar given by $R=-2\alpha_{\rm H}/r^2$. This originates from violation of the geometrical structure of four-dimensional curvature quantities. The conical singularity can disappear for the models in which the parameter $\alpha_{\rm H}$ vanishes in the limit that $r \to 0$. We propose explicit models without the conical singularity by properly designing the classical Lagrangian in such a way that the main contribution to $\alpha_{\rm H}$ comes from the field derivative $\phi'(r)$ around $r=0$. We show that the extension of covariant Galileons with a diatonic coupling allows for the recovery of general relativistic behavior inside a so-called Vainshtein radius. In this case, both the propagation of a fifth force and the deviation from Horndeski theories are suppressed outside a compact body in such a way that the model is compatible with local gravity experiments inside the solar system.

The disformal transformation of metric $g_{\mu \nu} \to \Omega^2 (\phi)g_{\mu \nu}+\Gamma(\phi,X) \partial_{\mu}\phi \partial_{\nu}\phi$, where $\phi$ is a scalar field with the kinetic energy $X= \partial_{\mu}\phi \partial^{\mu}\phi/2$, preserves the Lagrangian structure of Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories (which is the minimum extension of Horndeski theories). In the presence of matter, this transformation gives rise to a kinetic-type coupling between the scalar field $\phi$ and matter. We consider the Einstein frame in which the second-order action of tensor perturbations on the isotropic cosmological background is of the same form as that in General Relativity and study the role of couplings at the levels of both background and linear perturbations. We show that the effective gravitational potential felt by matter perturbations in the Einstein frame can be conveniently expressed in terms of the sum of a General Relativistic contribution and couplings induced by the modification of gravity. For the theories in which the transformed action belongs to a class of Horndeski theories, there is no anisotropic stress between two gravitational potentials in the Einstein frame due to a gravitational de-mixing. We propose a concrete dark energy model encompassing Brans-Dicke theories as well as theories with the tensor propagation speed $c_{\rm t}$ different from 1. We clarify the correspondence between physical quantities in the Jordan/Einstein frames and study the evolution of gravitational potentials and matter perturbations from the matter-dominated epoch to today in both analytic and numerical approaches.

We study the possibility of realizing a growth rate of matter density perturbations lower than that in General Relativity. Using the approach of the effective field theory of modified gravity encompassing theories beyond Horndeski, we derive the effective gravitational coupling $G_{\rm eff}$ and the gravitational slip parameter $\eta$ for perturbations deep inside the Hubble radius. In Horndeski theories we derive a necessary condition for achieving weak gravity associated with tensor perturbations, but this is not a sufficient condition due to the presence of a scalar-matter interaction that always enhances $G_{\rm eff}$. Beyond the Horndeski domain it is possible to realize $G_{\rm eff}$ smaller than Newton's gravitational constant $G$, while the scalar and tensor perturbations satisfy no-ghost and stability conditions. We present a concrete dark energy scenario with varying $c_{\rm t}$ and numerically study the evolution of perturbations to confront the model with the observations of redshift-space distortions and weak lensing.

We show that the universal $\alpha$-attractor models of inflation can be realized by including an auxiliary vector field $A_{\mu}$ for the Starobinsky model with the Lagrangian $f(R)=R+R^2/(6M^2)$. If the same procedure is applied to general modified $f(R)$ theories in which the Ricci scalar $R$ is replaced by $R+A_{\mu} A^{\mu}+\beta \nabla_{\mu}A^{\mu}$ with constant $\beta$, we obtain the Brans-Dicke theory with a scalar potential and the Brans-Dicke parameter $\omega_{\rm BD}=\beta^2/4$. We also place observational constraints on inflationary models based on auxiliary vector modified $f(R)$ theories from the latest Planck measurements of the Cosmic Microwave Background (CMB) anisotropies in both temperature and polarization. In the modified Starobinsky model, we find that the parameter $\beta$ is constrained to be $\beta<25$ (68\,\%\,confidence level) from the bounds of the scalar spectral index and the tensor-to-scalar ratio.

In the approach of the effective field theory of modified gravity, we derive the equations of motion for linear perturbations in the presence of a barotropic perfect fluid on the flat isotropic cosmological background. In a simple version of Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories, which is the minimum extension of Horndeski theories, we show that a slight deviation of the tensor propagation speed squared $c_{\rm t}^2$ from 1 generally leads to the large modification to the propagation speed squared $c_{\rm s}^2$ of a scalar degree of freedom $\phi$. This problem persists whenever the kinetic energy $\rho_X$ of the field $\phi$ is much smaller than the background energy density $\rho_m$, which is the case for most of dark energy models in the asymptotic past. Since the scaling solution characterized by the constant ratio $\rho_X/\rho_m$ is one way out for avoiding such a problem, we study the evolution of perturbations for a scaling dark energy model in the framework of GLPV theories in the Jordan frame. Provided the oscillating mode of scalar perturbations is fine-tuned so that it is initially suppressed, the anisotropic parameter $\eta=-\Phi/\Psi$ between the two gravitational potentials $\Psi$ and $\Phi$ significantly deviates from 1 for $c_{\rm t}^2$ away from 1. For other general initial conditions, the deviation of $c_{\rm t}^2$ from 1 gives rise to the large oscillation of $\Psi$ with the frequency related to $c_{\rm s}^2$. In both cases, the model can leave distinct imprints for the observations of CMB and weak lensing.

In the approach of the effective field theory of modified gravity, we derive the second-order action and the equation of motion for tensor perturbations on the flat isotropic cosmological background. This analysis accommodates a wide range of gravitational theories including Horndeski theories, its generalization, and the theories with spatial derivatives higher than second order (e.g., Horava-Lifshitz gravity). We obtain the inflationary power spectrum of tensor modes by taking into account corrections induced by higher-order spatial derivatives and slow-roll corrections to the de Sitter background. We also show that the leading-order spectrum in concrete modified gravitational theories can be mapped on to that in General Relativity under a disformal transformation. Our general formula will be useful to constrain inflationary models from the future precise measurement of the B-mode polarization in the cosmic microwave background.

It is known that Horndeski theories can be transformed to a sub-class of Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories under the disformal transformation of the metric $g_{\mu \nu} \to \Omega^2(\phi)g_{\mu \nu}+\Gamma (\phi,X) \nabla_{\mu} \phi \nabla_{\nu} \phi$, where $\Omega$ is a function of a scalar field $\phi$ and $\Gamma$ is another function depending on both $\phi$ and $X=g^{\mu \nu}\nabla_{\mu} \phi \nabla_{\nu} \phi$. We show that, with the choice of unitary gauge, both curvature and tensor perturbations on the flat isotropic cosmological background are generally invariant under the disformal transformation. By means of the effective field theories encompassing Horndeski and GLPV theories, we obtain the second-order actions of scalar/tensor perturbations and present the relations for physical quantities between the two frames. The invariance of the inflationary power spectra under the disformal transformation is explicitly proved up to next-to-leading order in slow-roll. In particular, we identify the existence of the Einstein frame in which the tensor power spectrum is of the same form as that in General Relativity and derive the condition under which the spectrum of gravitational waves in GLPV theories is red-tilted.

We review the effective field theory of modified gravity in which the Lagrangian involves three dimensional geometric quantities appearing in the 3+1 decomposition of space-time. On the flat isotropic cosmological background we expand a general action up to second order in the perturbations of geometric scalars, by taking into account spatial derivatives higher than two. Our analysis covers a wide range of gravitational theories-- including Horndeski theory/its recent generalizations and the projectable/non-projectable versions of Ho\v{r}ava-Lifshitz gravity. We derive the equations of motion for linear cosmological perturbations and apply them to the calculations of inflationary power spectra as well as the dark energy dynamics in Galileon theories. We also show that our general results conveniently recover stability conditions of Ho\v{r}ava-Lifshitz gravity already derived in the literature.

We consider perturbations of a static and spherically symmetric background endowed with a metric tensor and a scalar field in the framework of the effective field theory of modified gravity. We employ the previously developed 2+1+1 canonical formalism of a double Arnowitt-Deser-Misner (ADM) decomposition of space-time, which singles out both time and radial directions. Our building block is a general gravitational action that depends on scalar quantities constructed from the 2+1+1 canonical variables and the lapse. Variation of the action up to first-order in perturbations gives rise to three independent background equations of motion, as expected from spherical symmetry. The dynamical equations of linear perturbations follow from the second-order Lagrangian after a suitable gauge fixing. We derive conditions for the avoidance of ghosts and Laplacian instabilities for the odd-type perturbations. We show that our results not only incorporates those derived in the most general scalar-tensor theories with second-order equations of motion (the Horndeski theories) but they can be applied to more generic theories beyond Horndeski.

We mini-review the role of fundamental spin-0 bosons as bosonic coherent motion (BCM) in the Universe. The fundamental spin-0 bosons have the potential to account for the baryon number generation, cold dark matter (CDM) via BCM, dark energy, and inflation. Among these, here we focus on the CDM possibility because it can be experimentally tested with the current experimental techniques. We also comment briefly on the panoply of the other roles of spin-0 bosons.

The prospects for direct measurements of inflationary gravitational waves by next generation interferometric detectors inferred from the possible detection of B-mode polarization of the cosmic microwave background are studied. We compute the spectra of the gravitational wave background and the signal-to-noise ratios by two interferometric detectors (DECIGO and BBO) for large-field inflationary models in which the tensor-to-scalar ratio is greater than the order of 0.01. If the reheating temperature $T_{\rm RH}$ of chaotic inflation with the quadratic potential is high ($T_{\rm RH}>7.9\times10^6$ GeV for upgraded DECIGO and $T_{\rm RH}> 1.8\times 10^{6}$ GeV for BBO), it will be possible to reach the sensitivity of the gravitational background in future experiments at $3\sigma$ confidence level. The direct detection is also possible for natural inflation with the potential $V(\phi)=\Lambda^4 [1-\cos(\phi/f)]$, provided that $f>4.2 M_{\rm pl}$ (upgraded DECIGO) and $f>3.6 M_{\rm pl}$ (BBO) with $T_{\rm RH}$ higher than $10^8$ GeV. The quartic potential $V(\phi)=\lambda \phi^4/4$ with a non-minimal coupling $\xi$ between the inflaton field $\phi$ and the Ricci scalar $R$ gives rise to a detectable level of gravitational waves for $|\xi|$ smaller than the order of 0.01, irrespective of the reheating temperature.

The effective field theory (EFT) of cosmological perturbations is a useful framework to deal with the low-energy degrees of freedom present for inflation and dark energy. We review the EFT for modified gravitational theories by starting from the most general action in unitary gauge that involves the lapse function and the three-dimensional geometric scalar quantities appearing in the Arnowitt-Deser-Misner (ADM) formalism. Expanding the action up to quadratic order in the perturbations and imposing conditions for the elimination of spatial derivatives higher than second order, we obtain the Lagrangian of curvature perturbations and gravitational waves with a single scalar degree of freedom. The resulting second-order Lagrangian is exploited for computing the scalar and tensor power spectra generated during inflation. We also show that the most general scalar-tensor theory with second-order equations of motion-Horndeski theory-belongs to the action of our general EFT framework and that the background equations of motion in Horndeski theory can be conveniently expressed in terms of three EFT parameters. Finally we study the equations of matter density perturbations and the effective gravitational coupling for dark energy models based on Horndeski theory, to confront the models with the observations of large-scale structures and weak lensing.

We study the cosmology of an extended version of Horndeski theories with second-order equations of motion on the flat Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW) background. In addition to a dark energy field $\chi$ associated with the gravitational sector, we take into account multiple scalar fields $\phi_I$ ($I=1,2\cdots,N-1$) characterized by the Lagrangians $P^{(I)}(X_I)$ with $X_I=\partial_{\mu}\phi_I\partial^{\mu}\phi_I$. These additional scalar fields can model the perfect fluids of radiation and non-relativistic matter. We derive propagation speeds of scalar and tensor perturbations as well as conditions for the absence of ghosts. The theories beyond Horndeski induce non-trivial modifications to all the propagation speeds of $N$ scalar fields, but the modifications to those for the matter fields $\phi_I$ are generally suppressed relative to that for the dark energy field $\chi$. We apply our results to the covariantized Galileon with an Einstein-Hilbert term in which partial derivatives of the Minkowski Galileon are replaced by covariant derivatives. Unlike the covariant Galileon with second-order equations of motion in general space-time, the scalar propagation speed square $c_{s1}^2$ associated with the field $\chi$ becomes negative during the matter era for late-time tracking solutions, so the two Galileon theories can be clearly distinguished at the level of linear cosmological perturbations.

Nonlocal massive gravity can provide an interesting explanation for the late-time cosmic acceleration, with a dark energy equation of state $w_{\rm DE}$ smaller than $-1$ in the past. We derive the equations of linear cosmological perturbations to confront such models with the observations of large-scale structures. The effective gravitational coupling to nonrelativistic matter associated with galaxy clusterings is close to Newton's gravitational constant $G$ for a mass scale $m$ slightly smaller than today's Hubble parameter $H_0$. Taking into account the background expansion history as well as the evolution of matter perturbations $\delta_m$, we test for these models with Type Ia Supernovae (SnIa) from Union 2.1, the cosmic microwave background (CMB) measurements from Planck, a collection of baryon acoustic oscillations (BAO), and the growth rate data of $\delta_m$. Using a higher value of $H_0$ derived from its direct measurement ($H_0 \gtrsim 70$ km s$^{-1}$ Mpc$^{-1}$) the data strongly support the nonlocal massive gravity model ($-1.1 \lesssim w_{\rm DE} \lesssim -1.04$ in the past) over the $\Lambda$CDM model ($w_{\rm DE}=-1$), whereas for a lower prior (67 km s$^{-1}$ Mpc$^{-1}$ $\lesssim$ $H_0 \lesssim 70$ km s$^{-1}$ Mpc$^{-1}$) the two models are statistically comparable.

The general k-essence Lagrangian for the existence of cosmological scaling solutions is derived in the presence of multiple scalar fields coupled to a barotropic perfect fluid. In addition to the scaling fixed point associated with the dynamics during the radiation and matter eras, we also obtain a scalar-field dominated solution relevant to dark energy and discuss the stability of them in the two-field set-up. We apply our general results to a model of two canonical fields with coupled exponential potentials arising in string theory. Depending on model parameters and initial conditions, we show that the scaling matter-dominated epochs followed by an attractor with cosmic acceleration can be realized with/without the couplings to scalar fields. The different types of scaling solutions can be distinguished from each other by the evolution of the dark energy equation of state from high-redshifts to today.

In this paper, inflationary cosmology is reviewed, paying particular attention to its observational signatures associated with large-scale density perturbations generated from quantum fluctuations. In the most general scalar-tensor theories with second-order equations of motion, we derive the scalar spectral index $n_s$, the tensor-to-scalar ratio $r$, and the nonlinear estimator $f_{\rm NL}$ of primordial non-Gaussianities to confront models with observations of Cosmic Microwave Background (CMB) temperature anisotropies. Our analysis includes models such as potential-driven slow-roll inflation, k-inflation, Starobinsky inflation, and Higgs inflation with non-minimal/derivative/Galileon couplings. We constrain a host of inflationary models by using the Planck data combined with other measurements to find models most favored observationally in the current literature. We also study anisotropic inflation based on a scalar coupling with a vector (or, two-form) field and discuss its observational signatures appearing in the two-point and three-point correlation functions of scalar and tensor perturbations.

We present a framework for discussing the cosmology of dark energy and dark matter based on two scalar degrees of freedom. An effective field theory of cosmological perturbations is employed. A unitary gauge choice renders the dark energy field into the gravitational sector, for which we adopt a generic Lagrangian depending on three-dimensional geometrical scalar quantities arising in the ADM decomposition. We add to this dark-energy associated gravitational sector a scalar field $\phi$ and its kinetic energy $X$ as dark matter variables. Compared to the single-field case, we find that there are additional conditions to obey in order to keep the equations of motion for linear cosmological perturbations at second order. For such a second-order multi-field theory we derive conditions under which ghosts and Laplacian instabilities of the scalar and tensor perturbations are absent. We apply our general results to models with dark energy emerging in the framework of the Horndeski theory and dark matter described by a k-essence Lagrangian $P(\phi,X)$. We derive the effective coupling between such an imperfect-fluid dark matter and the gravitational sector under the quasi-static approximation on sub-horizon scales. By considering the purely kinetic Lagrangian $P(X)$ as a particular case, the formalism is verified to reproduce the gravitational coupling of a perfect-fluid dark matter.

We place observational likelihood constraints on braneworld and non-commutative inflation for a number of inflaton potentials, using Planck, WMAP polarization and BAO data. Both braneworld and non-commutative scenarios of the kind considered here are limited by the most recent data even more severely than standard general-relativity models. At more than 95 % confidence level, the monomial potential $V(\phi) \propto \phi^p$ is ruled out for $p \geq 2$ in the Randall--Sundrum (RS) braneworld cosmology and, for $p>0$, also in the high-curvature limit of the Gauss--Bonnet (GB) braneworld and in the infrared limit of non-commutative inflation, due to a large scalar spectral index. Some parameter values for natural inflation, small-varying inflaton models and Starobinsky inflation are allowed in all scenarios, although some tuning is required for natural inflation in a non-commutative spacetime.

The Javalambre-Physics of the Accelerated Universe Astrophysical Survey (J-PAS) is a narrow band, very wide field Cosmological Survey to be carried out from the Javalambre Observatory in Spain with a purpose-built, dedicated 2.5m telescope and a 4.7 sq.deg. camera with 1.2Gpix. Starting in late 2015, J-PAS will observe 8500sq.deg. of Northern Sky and measure $0.003(1+z)$ photo-z for $9\times10^7$ LRG and ELG galaxies plus several million QSOs, sampling an effective volume of $\sim 14$ Gpc$^3$ up to $z=1.3$ and becoming the first radial BAO experiment to reach Stage IV. J-PAS will detect $7\times 10^5$ galaxy clusters and groups, setting constrains on Dark Energy which rival those obtained from its BAO measurements. Thanks to the superb characteristics of the site (seeing ~0.7 arcsec), J-PAS is expected to obtain a deep, sub-arcsec image of the Northern sky, which combined with its unique photo-z precision will produce one of the most powerful cosmological lensing surveys before the arrival of Euclid. J-PAS unprecedented spectral time domain information will enable a self-contained SN survey that, without the need for external spectroscopic follow-up, will detect, classify and measure $\sigma_z\sim 0.5\%$ redshifts for $\sim 4000$ SNeIa and $\sim 900$ core-collapse SNe. The key to the J-PAS potential is its innovative approach: a contiguous system of 54 filters with $145\AA$ width, placed $100\AA$ apart over a multi-degree FoV is a powerful "redshift machine", with the survey speed of a 4000 multiplexing low resolution spectrograph, but many times cheaper and much faster to build. The J-PAS camera is equivalent to a 4.7 sq.deg. "IFU" and it will produce a time-resolved, 3D image of the Northern Sky with a very wide range of Astrophysical applications in Galaxy Evolution, the nearby Universe and the study of resolved stellar populations.

We study observational signatures of two classes of anisotropic inflationary models in which an inflaton field couples to (i) a vector kinetic term F_{mu nu}F^{mu nu} and (ii) a two-form kinetic term H_{mu nu lambda}H^{mu nu lambda}. We compute the corrections from the anisotropic sources to the power spectrum of gravitational waves as well as the two-point cross correlation between scalar and tensor perturbations. The signs of the anisotropic parameter g_* are different depending on the vector and the two-form models, but the statistical anisotropies generally lead to a suppressed tensor-to-scalar ratio r and a smaller scalar spectral index n_s in both models. In the light of the recent Planck bounds of n_s and r, we place observational constraints on several different inflaton potentials such as those in chaotic and natural inflation in the presence of anisotropic interactions. In the two-form model we also find that there is no cross correlation between scalar and tensor perturbations, while in the vector model the cross correlation does not vanish. The non-linear estimator f_{NL} of scalar non-Gaussianities in the two-form model is generally smaller than that in the vector model for the same orders of |g_*|, so that the former is easier to be compatible with observational bounds of non-Gaussianities than the latter.

It is known that power-law k-inflation can be realized for the Lagrangian $P=Xg(Y)$, where $X=-(\partial \phi)^2/2$ is the kinetic energy of a scalar field $\phi$ and $g$ is an arbitrary function in terms of $Y=Xe^{\lambda \phi/M_{pl}}$ ($\lambda$ is a constant and $M_{pl}$ is the reduced Planck mass). In the presence of a vector field coupled to the inflaton with an exponential coupling $f(\phi) \propto e^{\mu \phi/M_{pl}}$, we show that the models with the Lagrangian $P=Xg(Y)$ generally give rise to anisotropic inflationary solutions with $\Sigma/H=constant$, where $\Sigma$ is an anisotropic shear and $H$ is an isotropic expansion rate. Provided these anisotropic solutions exist in the regime where the ratio $\Sigma/H$ is much smaller than 1, they are stable attractors irrespective of the forms of $g(Y)$. We apply our results to concrete models of k-inflation such as the generalized dilatonic ghost condensate/the DBI model and we numerically show that the solutions with different initial conditions converge to the anisotropic power-law inflationary attractors. Even in the de Sitter limit ($\lambda \to 0$) such solutions can exist, but in this case the null energy condition is generally violated. The latter property is consistent with the Wald's cosmic conjecture stating that the anisotropic hair does not survive on the de Sitter background in the presence of matter respecting the dominant/strong energy conditions.

We present a general covariant action for massive gravity merging together a class of "non-polynomial" and super-renormalizable or finite theories of gravity with the non-local theory of gravity recently proposed by Jaccard, Maggiore and Mitsou (Phys. Rev. D 88 (2013) 044033). Our diffeomorphism invariant action gives rise to the equations of motion appearing in non-local massive massive gravity plus quadratic curvature terms. Not only the massive graviton propagator reduces smoothly to the massless one without a vDVZ discontinuity, but also our finite theory of gravity is unitary at tree level around the Minkowski background. We also show that, as long as the graviton mass $m$ is much smaller the today's Hubble parameter $H_0$, a late-time cosmic acceleration can be realized without a dark energy component due to the growth of a scalar degree of freedom. In the presence of the cosmological constant $\Lambda$, the dominance of the non-local mass term leads to a kind of "degravitation" for $\Lambda$ at the late cosmological epoch.

Quintessence is a canonical scalar field introduced to explain the late-time cosmic acceleration. The cosmological dynamics of quintessence is reviewed, paying particular attention to the evolution of the dark energy equation of state w. For the field potentials having tracking and thawing properties, the evolution of w can be known analytically in terms of a few model parameters. Using the analytic expression of w, we constrain quintessence models from the observations of supernovae type Ia, cosmic microwave background, and baryon acoustic oscillations. The tracking freezing models are hardly distinguishable from the LCDM model, whereas in thawing models the today's field equation of state is constrained to be w_0<-0.7 (95 % CL). We also derive an analytic formula for the growth rate of matter density perturbations in dynamical dark energy models, which allows a possibility to put further bounds on w from the measurement of redshift-space distortions in the galaxy power spectrum. Finally we review particle physics models of quintessence- such as those motivated by supersymmetric theories. The field potentials of thawing models based on a pseudo-Nambu-Goldstone boson or on extended supergravity theories have a nice property that a tiny mass of quintessence can be protected against radiative corrections.

We study how the Vainshtein mechanism operates in the most general scalar-tensor theories with second-order equations of motion. The field equations of motion, which can be also applicable to most of other screening scenarios proposed in literature, are generally derived in a spherically symmetric space-time with a matter source. In the presence of a field coupling to the Ricci scalar, we clarify conditions under which the Vainshtein mechanism is at work in a weak gravitational background. We also obtain the solutions of the field equation inside a spherically symmetric body and show how they can be connected to exterior solutions that accommodate the Vainshtein mechanism. We apply our general results to a number of concrete models such as the covariant/extended Galileons and the DBI Galileons with Gauss-Bonnet and other terms. In these models the fifth force can be suppressed to be compatible with solar-system constraints, provided that non-linear field kinetic terms coupled to the Einstein tensor do not dominate over other non-linear field self-interactions.

We place observational constraints on slow-variation single-field inflationary models by carrying out the cosmological Monte Carlo simulation with the recent data of Planck combined with the WMAP large-angle polarization, baryon acoustic oscillations, and ACT/SPT temperature data. Our analysis covers a wide variety of models with second-order equations of motion-- including potential-driven slow-roll inflation, non-minimally coupled models, running kinetic couplings, Brans-Dicke theories, potential-driven Galileon inflation, field-derivative couplings to the Einstein tensor, and k-inflation. In the presence of running kinetic exponential couplings, covariant Galileon terms, and field-derivative couplings, the tensor-to-scalar ratio of the self-coupling potential gets smaller relative to that in standard slow-roll inflation, but the models lie outside the 68 % CL observational contour. We also show that k-inflation models can be tightly constrained by adding the bounds from the scalar non-Gaussianities. The small-field inflationary models with asymptotic flat Einstein-frame potentials in the regime phi >> M_{pl} generally fit the data very well. These include the models such as Kahler-moduli inflation, non-minimally coupled Higgs inflation, and inflation in Brans-Dicke theories in the presence of the potential V(phi)=3M^2 (phi-M_{pl})^2/4 with the Brans-Dicke parameter omega_{BD}<O(1) (which covers the Starobinsky's model f(R)=R+R^2/(6M^2) as a special case).

The growth rate of matter density perturbations has been measured from redshift-space distortion (RSD) in the galaxy power spectrum. We constrain the model parameter space for representative modified gravity models to explain the dark energy problem, by using the recent data of f_m(z)sigma_8(z) at the redshifts z = 0.06--0.8 measured by WiggleZ, SDSS LRG, BOSS, and 6dFGRS. We first test the Hu-Sawicki's f(R) dark energy model, and find that only the parameter region close to the standard Lambda Cold Dark Matter (Lambda-CDM) model is allowed (lambda > 12 and 5 for n = 1.5 and 2, respectively, at 95% CL). We then investigate the covariant Galileon model and show that the parameter space consistent with the background expansion history is excluded by the RSD data at more than 10 sigma because of the too large growth rate predicted by the theory. Finally, we consider the extended Galileon scenario, and we find that, in contrast to the covariant Galileon, there is a model parameter space for a tracker solution that is consistent with the RSD data within a 2 sigma level.

We study an inflationary scenario with a two-form field to which an inflaton couples non-trivially. First, we show that anisotropic inflation can be realized as an attractor solution and that the two-form hair remains during inflation. A statistical anisotropy can be developed because of a cumulative anisotropic interaction induced by the background two-form field. The power spectrum of curvature perturbations has a prolate-type anisotropy, in contrast to the vector models having an oblate-type anisotropy. We also evaluate the bispectrum and trispectrum of curvature perturbations by employing the in-in formalism based on the interacting Hamiltonians. We find that the non-linear estimators $f_{NL}$ and $\tau_{NL}$ are correlated with the amplitude $g_*$ of the statistical anisotropy in the power spectrum. Unlike the vector models, both $f_{NL}$ and $\tau_{NL}$ vanish in the squeezed limit. However, the estimator $f_{NL}$ can reach the order of 10 in the equilateral and enfolded limits. These results are consistent with the latest bounds on $f_{NL}$ constrained by Planck.

For two types of quintessence models having thawing and tracking properties, there exist analytic solutions for the dark energy equation of state w expressed in terms of several free parameters. We put observational bounds on the parameters in such scenarios by using the recent data of Supernovae type Ia, Cosmic Microwave Background, and Baryon Acoustic Oscillations. The observational constraints are quite different depending on whether or not the recent BAO data from BOSS are taken into account. With the BOSS data the upper bounds of today's values of w (=w_0) in thawing models is very close to -1, whereas without this data the values of w_0 away from -1 can be still allowed. The tracker equation of state w_{(0)} during the matter era is constrained to be w_{(0)}<-0.949 at 95 % confidence level even without the BOSS data, so that the tracker models with w away from -1 are severely disfavored. We also study observational constraints on scaling models in which w starts to evolve from 0 in the deep matter era and show that the transition to the equation of state close to w=-1 needs to occur at an early cosmological epoch. In the three classes of quintessence models studied in this paper, the past evolution of the Hubble parameters in the best-fit models shows only less than the 2.5 % difference compared to the LCDM.

In the Horndeski's most general scalar-tensor theories, we derive the three-point correlation function of scalar non-Gaussianities generated during single-field inflation in the presence of slow-variation corrections to the leading-order term. Unlike previous works, the resulting bispectrum is valid for any shape of non-Gaussianities. In the squeezed limit, for example, this gives rise to the same consistency relation as that derived by Maldacena in standard single-field slow-roll inflation. We estimate the shape close to the squeezed one at which the effect of the term inversely proportional to the scalar propagation speed squared begins to contribute to the bispectrum. We also show that the leading-order bispectrum can be expressed by the linear combination of two convenient bases whose shapes are highly correlated with equilateral and orthogonal types respectively. We present concrete models in which the orthogonal and enfolded shapes can dominate over the equilateral one.

The red-shift space distortions in the galaxy power spectrum can be used to measure the growth rate of matter density perturbations delta_m. For dynamical dark energy models in General Relativity we provide a convenient analytic formula of f(z) sigma_8(z) written as a function of the redshift z, where f=d ln delta_m/d ln a (a is the cosmological scale factor) and sigma_8 is the rms amplitude of over-density at the scale 8 h^{-1} Mpc. Our formula can be applied to the models of imperfect fluids, quintessence, and k-essence, provided that the dark energy equation of state w does not vary significantly and that the sound speed is not much smaller than 1. We also place observational constraints on dark energy models of constant w and tracking quintessence from the recent data of red-shift space distortions.

For the models of inflation driven by the potential energy of an inflaton field $\phi$, the covariant Galileon Lagrangian $(\partial\phi)^2\Box \phi$ generally works to slow down the evolution of the field. On the other hand, if the Galileon self-interaction is dominant relative to the standard kinetic term, we show that there is no oscillatory regime of inflaton after the end of inflation. This is typically accompanied by the appearance of the negative propagation speed squared $c_s^2$ of a scalar mode, which leads to the instability of small-scale perturbations. For chaotic inflation and natural inflation we clarify the parameter space in which inflaton oscillates coherently during reheating. Using the WMAP constraints of the scalar spectral index and the tensor-to-scalar ratio as well, we find that the self coupling $\lambda$ of the potential $V(\phi)=\lambda \phi^4/4$ is constrained to be very much smaller than 1 and that the symmetry breaking scale $f$ of natural inflation cannot be less than the reduced Planck mass $M_{\rm pl}$. We also show that, in the presence of other covariant Galileon Lagrangians, there are some cases in which inflaton oscillates coherently even for the self coupling $\lambda$ of the order of 0.1, but still the instability associated with negative $c_s^2$ is generally present.

We study inflation and preheating in F(R) supergravity characterized by two mass scales of a scalar degree of freedom (scalaron): M (associated with the inflationary era) and m (associated with the preheating era). The allowed values of the masses M and m are derived from the amplitude of the CMB temperature anisotropies. We show that our model is consistent with the joint observational constraints of WMAP and other measurements in the regime where a sufficient amount of inflation (with the number of e-foldings larger than 50) is realized. In the low-energy regime relevant to preheating, we derive the effective scalar potential in the presence of a pseudo-scalar field chi coupled to the inflaton (scalaron) field phi. If m is much larger than M, we find that there exists the preheating stage in which the field perturbations delta chi and delta phi rapidly grow by a broad parametric resonance.

We place observational constraints on models with the late-time cosmic acceleration based on a number of parametrizations allowing fast transitions for the equation of state of dark energy. In addition to the model of Linder and Huterer where the dark energy equation of state $w$ monotonically grows or decreases in time, we propose two new parametrizations in which $w$ has an extremum. We carry out the likelihood analysis with the three parametrizations by using the observational data of supernovae type Ia, cosmic microwave background, and baryon acoustic oscillations. Although the transient cosmic acceleration models with fast transitions can give rise to the total chi square smaller than that in the $\Lambda$-Cold-Dark-Matter ($\Lambda$CDM) model, these models are not favored over $\Lambda$CDM when one uses the Akaike information criterion which penalizes the extra degrees of freedom present in the parametrizations.

A field kinetic coupling with the Einstein tensor leads to a gravitationally enhanced friction during inflation, by which even steep potentials with theoretically natural model parameters can drive cosmic acceleration. In the presence of this non-minimal derivative coupling we place observational constraints on a number of representative inflationary models such as chaotic inflation, inflation with exponential potentials, natural inflation, and hybrid inflation. We show that most of the models can be made compatible with the current observational data mainly due to the suppressed tensor-to-scalar ratio.

In second-order scalar-tensor theories we study how the Vainshtein mechanism works in a spherically symmetric background with a matter source. In the presence of the field coupling $F(\phi)=e^{-2Q\phi}$ with the Ricci scalar $R$ we generally derive the Vainshtein radius within which the General Relativistic behavior is recovered even for the coupling $Q$ of the order of unity. Our analysis covers the models such as the extended Galileon and Brans-Dicke theories with a dilatonic field self-interaction. We show that, if these models are responsible for the cosmic acceleration today, the corrections to gravitational potentials are generally small enough to be compatible with local gravity constraints.

The extended Galileon models possess tracker solutions with de Sitter attractors along which the dark energy equation of state is constant during the matter-dominated epoch, i.e. w_DE = -1-s, where s is a positive constant. Even with this phantom equation of state there are viable parameter spaces in which the ghosts and Laplacian instabilities are absent. Using the observational data of the supernovae type Ia, the cosmic microwave background (CMB), and baryon acoustic oscillations, we place constraints on the tracker solutions at the background level and find that the parameter s is constrained to be s=0.034 (-0.034,+0.327) (95% CL) in the flat Universe. In order to break the degeneracy between the models we also study the evolution of cosmological density perturbations relevant to the large-scale structure (LSS) and the Integrated-Sachs-Wolfe (ISW) effect in CMB. We show that, depending on the model parameters, the LSS and the ISW effect is either positively or negatively correlated. It is then possible to constrain viable parameter spaces further from the observational data of the ISW-LSS cross-correlation as well as from the matter power spectrum.

We study in detail the power spectra of scalar and tensor perturbations generated during inflation in loop quantum cosmology (LQC). After clarifying in a novel quantitative way how inverse-volume corrections arise in inhomogeneous settings, we show that they can generate large running spectral indices, which generally lead to an enhancement of power at large scales. We provide explicit formulae for the scalar/tensor power spectra under the slow-roll approximation, by taking into account corrections of order higher than the runnings. Via a standard analysis, we place observational bounds on the inverse-volume quantum correction \delta ~ a^{- \sigma} (\sigma >0, $a$ is the scale factor) and the slow-roll parameter \epsilon_V for power-law potentials as well as exponential potentials by using the data of WMAP 7yr combined with other observations. We derive the constraints on \delta for two pivot wavenumbers k_0 for several values of \delta. The quadratic potential can be compatible with the data even in the presence of the LQC corrections, but the quartic potential is in tension with observations. We also find that the upper bounds on \delta (k_0) for given \sigma and k_0 are insensitive to the choice of the inflaton potentials.

In the Horndeski's most general scalar-tensor theories the equations of scalar density perturbations are derived in the presence of non-relativistic matter minimally coupled to gravity. Under a quasi-static approximation on sub-horizon scales we obtain the effective gravitational coupling $G_{eff}$ associated with the growth rate of matter perturbations as well as the effective gravitational potential $\Phi_{eff}$ relevant to the deviation of light rays. We then apply our formulas to a number of modified gravitational models of dark energy--such as those based on f(R) theories, Brans-Dicke theories, kinetic gravity braidings, covariant Galileons, and field derivative couplings with the Einstein tensor. Our results are useful to test the large-distance modification of gravity from the future high-precision observations of large-scale structure, weak lensing, and cosmic microwave background.

In the inflationary scenario of loop quantum cosmology in the presence of inverse-volume corrections, we give analytic formulas for the power spectra of scalar and tensor perturbations convenient to compare with observations. Since inverse-volume corrections can provide strong contributions to the running spectral indices, inclusion of terms higher than the second-order runnings in the power spectra is crucially important. Using the recent data of cosmic microwave background and other cosmological experiments, we place bounds on the quantum corrections.

In the Horndeski's most general scalar-tensor theories with second-order field equations, we derive the conditions for the avoidance of ghosts and Laplacian instabilities associated with scalar, tensor, and vector perturbations in the presence of two perfect fluids on the flat Friedmann-Lemaitre-Robertson-Walker (FLRW) background. Our general results are useful for the construction of theoretically consistent models of dark energy. We apply our formulas to extended Galileon models in which a tracker solution with an equation of state smaller than -1 is present. We clarify the allowed parameter space in which the ghosts and Laplacian instabilities are absent and we numerically confirm that such models are indeed cosmologically viable.

We study the cosmology of a generalized Galileon field $\phi$ with five covariant Lagrangians in which $\phi$ is replaced by general scalar functions $f_{i}(\phi)$ (i=1,...,5). For these theories, the equations of motion remain at second-order in time derivatives. We restrict the functional forms of $f_{i}(\phi)$ from the demand to obtain de Sitter solutions responsible for dark energy. There are two possible choices for power-law functions $f_{i}(\phi)$, depending on whether the coupling $F(\phi)$ with the Ricci scalar $R$ is independent of $\phi$ or depends on $\phi$. The former corresponds to the covariant Galileon theory that respects the Galilean symmetry in the Minkowski space-time. For generalized Galileon theories we derive the conditions for the avoidance of ghosts and Laplacian instabilities associated with scalar and tensor perturbations as well as the condition for the stability of de Sitter solutions. We also carry out detailed analytic and numerical study for the cosmological dynamics in those theories.

For very general scalar-field theories in which the equations of motion are at second-order, we evaluate the three-point correlation function of primordial scalar perturbations generated during inflation. We show that the shape of non-Gaussianities is well approximated by the equilateral type. The equilateral non-linear parameter f_NL^equil is derived on the quasi de Sitter background where the slow-variation parameters are much smaller than unity. We apply our formula for f_NL^equil to a number of single-field models of inflation--such as k-inflation, k-inflation with Galileon terms, potential-driven Galileon inflation, nonminimal coupling models (including field-derivative coupling models), and Gauss-Bonnet gravity.

We study chaotic inflation in the context of modified gravitational theories. Our analysis covers models based on (i) a field coupling $\omega(\phi)$ with the kinetic energy $X$ and a nonmimimal coupling $\zeta \phi^{2} R/2$ with a Ricci scalar $R$, (ii) Brans-Dicke (BD) theories, (iii) Gauss-Bonnet (GB) gravity, and (iv) gravity with a Galileon correction. Dilatonic coupling with the kinetic energy and/or negative nonminimal coupling are shown to lead to compatibility with observations of the Cosmic Microwave Background (CMB) temperature anisotropies for the self-coupling inflaton potential $V(\phi)=\lambda \phi^{4}/4$. BD theory with a quadratic inflaton potential, which covers Starobinsky's $f(R)$ model $f(R)=R+R^{2}/(6M^{2})$ with the BD parameter $\omega_{BD}=0$, gives rise to a smaller tensor-to-scalar ratio for decreasing $\omega_{BD}$. In the presence of a GB term coupled to the field $\phi$, we express the scalar/tensor spectral indices $n_{s}$ and $n_{t}$ as well as the tensor-to-scalar ratio $r$ in terms of two slow-roll parameters and place bounds on the strength of the GB coupling from the joint data analysis of WMAP 7yr combined with other observations. We also study the Galileon-like self-interaction $\Phi(\phi) X \square\phi$ with exponential coupling $\Phi(\phi) \propto e^{\mu\phi}$. Using a CMB likelihood analysis we put bounds on the strength of the Galileon coupling and show that the self coupling potential can in fact be made compatible with observations in the presence of the exponential coupling with $\mu>0$.

We compute the three-point correlation function of primordial scalar density perturbations in a general single-field inflationary scenario, where a scalar field phi has a direct coupling with the Ricci scalar R and the Gauss-Bonnet term GB. Our analysis also covers the models in which the Lagrangian includes a function non-linear in the field kinetic energy X=-(nabla phi)^2/2, and a Galileon-type field self-interaction G(phi, X)*(Box phi), where G is a function of phi and X. We provide a general analytic formula for the equilateral non-Gaussianity parameter f_{NL}^{equil} associated with the bispectrum of curvature perturbations. A quasi de Sitter approximation in terms of slow-variation parameters allows us to derive a simplified form of f_{NL}^{equil} convenient to constrain various inflation models observationally. If the propagation speed of the scalar perturbations is much smaller than the speed of light, the Gauss-Bonnet term as well as the Galileon-type field self-interaction can give rise to large non-Gaussianities testable in future observations. We also show that, in Brans-Dicke theory with a field potential (including f(R) gravity), f_{NL}^{equil} is of the order of slow-roll parameters as in standard inflation driven by a minimally coupled scalar field.

We study observational constraints on the assisted k-inflation models in which multiple scalar fields join an attractor characterized by an effective single field $\phi$. This effective single-field system is described by the Lagrangian $P=Xg(Y)$, where $X$ is the kinetic energy of $\phi$, $\lambda$ is a constant, and $g$ is an arbitrary function in terms of $Y=X e^{\lambda \phi}$. Our analysis covers a wide variety of k-inflation models such as dilatonic ghost condensate, Dirac-Born-Infeld (DBI) field, tachyon, as well as the canonical field with an exponential potential. We place observational bounds on the parameters of each model from the WMAP 7yr data combined with Baryon Acoustic Oscillations (BAO) and the Hubble constant measurement. Using the observational constraints of the equilateral non-Gaussianity parameter $f_{NL}^{equil}$, we further restrict the allowed parameter space of dilatonic ghost condensate and DBI models. We extend the analysis to more general models with several different choices of $g(Y)$ and show that the models such as $g(Y)=c_0+c_p Y^p$ ($p \ge 3$) are excluded by the joint data analysis of the scalar/tensor spectra and primordial non-Gaussianities.

We study the evolution of matter density perturbations in Galileon cosmology where the late-time cosmic acceleration can be realized by a field kinetic energy. We obtain full perturbation equations at linear order in the presence of five covariant Lagrangians ${cal L}_i$ ($i=1,...,5$) satisfying a Galilean symmetry in the flat space-time. The equations for a matter perturbation as well as an effective gravitational potential are derived under a quasi-static approximation on sub-horizon scales. This approximation can reproduce full numerical solutions with high accuracy for the wavelengths relevant to large-scale structures. For the model parameters constrained by the background expansion history of the Universe the growth rate of matter perturbations is larger than that in the LCDM model, with the growth index $gamma$ today typically smaller than 0.4. We also find that, even on very large scales associated with the Integrated-Sachs-Wolfe (ISW) effect in Cosmic Microwave Background (CMB) temperature anisotropies, the effective gravitational potential exhibits a temporal growth during the transition from the matter era to the epoch of cosmic acceleration. These properties are useful to distinguish the Galileon model from the LCDM in future high-precision observations.

We demonstrate that there does exist an equilibrium description of thermodynamics on the apparent horizon in the expanding cosmological background for a wide class of modified gravity theories with the Lagrangian density $f(R, \phi, X)$, where $R$ is the Ricci scalar and $X$ is the kinetic energy of a scalar field $\phi$. This comes from a suitable definition of an energy momentum tensor of the "dark" component obeying the local energy conservation law in the Jordan frame. It is shown that the equilibrium description in terms of the horizon entropy $S$ is convenient because it takes into account the contribution of the horizon entropy $\hat{S}$ in non-equilibrium thermodynamics as well as an entropy production term.

We review recent progress of modified gravity models of dark energy--based on f(R) gravity, scalar-tensor theories, braneworld gravity, Galileon gravity, and other theories. In f(R) gravity and Brans-Dicke theory it is possible to design viable models consistent with local gravity constraints under a chameleon mechanism, while satisfying conditions for the cosmological viability. The Dvali-Gabadazde-Porrati braneworld model can be compatible with local gravity constraints through a nonlinear field self-interaction arising from a brane-bending mode, but the self-accelerating solution contains a ghost mode in addition to the tension with observational data about the cosmic expansion history. The extension of the field self-interaction to more general forms satisfying a Galilean symmetry in the flat space-time allows a possibility to avoid the appearance of ghosts and Laplacian instabilities, while the late-time cosmic acceleration can be realized by the field kinetic energy. We study the evolution of cosmological perturbations in those models to place constraints on model parameters from the observations of large-scale structure, cosmic microwave background, and weak lensing. We also briefly review other modified gravitational models of dark energy-- such as those based on Gauss-Bonnet gravity and Lorentz-violating theories.

We study the cosmology of a covariant Galileon field with five covariant Lagrangians and confront this theory with the most recent cosmological probes: the type Ia supernovae data (Constitution and Union2 sets), cosmic microwave background (WMAP7) and the baryon acoustic oscillations (SDSS7). In the Galileon cosmology with a late-time de Sitter attractor, there is a tracker that attracts solutions with different initial conditions to a common trajectory. Including the cosmic curvature K, we place observational constraints on two distinct cases: (i) the tracker, and (ii) the generic solutions to the equations of motion. We find that the tracker solution can be consistent with the individual observational data, but it is disfavored by the combined data analysis. The generic solutions fare quite well when a non-zero curvature parameter is taken into account, but the Akaike and Bayesian information criteria show that they are not particularly favored over the LCDM model.

In chameleon dark energy models, local gravity constraints tend to rule out parameters in which observable cosmological signatures can be found. We study viable chameleon potentials consistent with a number of recent observational and experimental bounds. A novel chameleon field potential, motivated by f(R) gravity, is constructed where observable cosmological signatures are present both at the background evolution and in the growth-rate of the perturbations. We study the evolution of matter density perturbations on low redshifts for this potential and show that the growth index today gamma_0 can have significant dispersion on scales relevant for large scale structures. The values of gamma_0 can be even smaller than 0.2 with large variations of gamma on very low redshifts for the model parameters constrained by local gravity tests. This gives a possibility to clearly distinguish these chameleon models from the Lambda-Cold-Dark-Matter model in future high-precision observations.

We study the cosmology of a covariant scalar field respecting a Galilean symmetry in flat space-time. We show the existence of a tracker solution that finally approaches a de Sitter fixed point responsible for cosmic acceleration today. The viable region of model parameters is clarified by deriving conditions under which ghosts and Laplacian instabilities of scalar and tensor perturbations are absent. The field equation of state exhibits a peculiar phantom-like behavior along the tracker, which allows a possibility to observationally distinguish the Galileon gravity from the Lambda-CDM model.

In Brans-Dicke theory a non-linear self interaction of a scalar field allows a possibility of realizing the late-time cosmic acceleration, while recovering the General Relativistic behavior at early cosmological epochs. We extend this to more general modified gravitational theories in which a de Sitter solution for dark energy exists without using a field potential. We derive a condition for the stability of the de Sitter point and study the background cosmological dynamics of such theories. We also restrict the allowed region of model parameters from the demand for the avoidance of ghosts and instabilities. A peculiar evolution of the field propagation speed allows us to distinguish those theories from the LCDM model.

Over the past decade, f(R) theories have been extensively studied as one of the simplest modifications to General Relativity. In this article we review various applications of f(R) theories to cosmology and gravity - such as inflation, dark energy, local gravity constraints, cosmological perturbations, and spherically symmetric solutions in weak and strong gravitational backgrounds. We present a number of ways to distinguish those theories from General Relativity observationally and experimentally. We also discuss the extension to other modified gravity theories such as Brans-Dicke theory and Gauss-Bonnet gravity, and address models that can satisfy both cosmological and local gravity constraints.

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.

We try to address quantitatively the question whether a new mass is needed to fit current supernovae data. For this purpose, we consider an infra-red modification of gravity that does not contain any new mass scale but systematic subleading corrections proportional to the curvature. The modifications are of the same type as the one recently derived by enforcing the "Ultra Strong Equivalence Principle" (USEP) upon a Friedmann-Lemaitre-Robertson-Walker universe in the presence of a scalar field. The distance between two comoving observers is altered by these corrections and the observations at high redshift affected at any time during the cosmic evolution. While the specific values of the parameters predicted by USEP are ruled out, there are regions of parameter space that fit SnIa data very well. This allows an interesting possibility to explain the apparent cosmic acceleration today without introducing either a dark energy component or a new mass scale.

We show that it is possible to obtain a picture of equilibrium thermodynamics on the apparent horizon in the expanding cosmological background for a wide class of modified gravity theories with the Lagrangian density $f(R, \phi, X)$, where $R$ is the Ricci scalar and $X$ is the kinetic energy of a scalar field $\phi$. This comes from a suitable definition of an energy momentum tensor of the "dark" component that respects to a local energy conservation in the Jordan frame. In this framework the horizon entropy $S$ corresponding to equilibrium thermodynamics is equal to a quarter of the horizon area $A$ in units of gravitational constant $G$, as in Einstein gravity. For a flat cosmological background with a decreasing Hubble parameter, $S$ globally increases with time, as it happens for viable $f(R)$ inflation and dark energy models. We also show that the equilibrium description in terms of the horizon entropy $S$ is convenient because it takes into account the contribution of both the horizon entropy $\hat{S}$ in non-equilibrium thermodynamics and an entropy production term.

Constantly accumulating observational data continue to confirm that about 70% of the energy density today consists of dark energy responsible for the accelerated expansion of the Universe. We present recent observational bounds on dark energy constrained by the type Ia supernovae, cosmic microwave background, and baryon acoustic oscillations. We review a number of theoretical approaches that have been adopted so far to explain the origin of dark energy. This includes the cosmological constant, modified matter models (such as quintessence, k-essence, coupled dark energy, unified models of dark energy and dark matter), modified gravity models (such as f(R) gravity, scalar-tensor theories, braneworlds), and inhomogeneous models. We also discuss observational and experimental constraints on those models and clarify which models are favored or ruled out in current observations.

We study the evolution of cosmological perturbations in f(G) gravity, where the Lagrangian is the sum of a Ricci scalar R and an arbitrary function f in terms of a Gauss-Bonnet term G. We derive the equations for perturbations assuming matter to be described by a perfect fluid with a constant equation of state w. We show that density perturbations in perfect fluids exhibit negative instabilities during both the radiation and the matter domination, irrespective of the form of f(G). This growth of perturbations gets stronger on smaller scales, which is difficult to be compatible with the observed galaxy spectrum unless the deviation from General Relativity is very small. Thus f(G) cosmological models are effectively ruled out from this Ultra-Violet instability, even though they can be compatible with the late-time cosmic acceleration and local gravity constraints.

Cosmological scaling solutions, which give rise to a scalar-field density proportional to a background fluid density during radiation and matter eras, are attractive to alleviate the energy scale problem of dark energy. In the presence of multiple scalar fields the scaling solution can exit to the epoch of cosmic acceleration through the so-called assisted inflation mechanism. We study cosmological dynamics of a multi-field system in details with a general Lagrangian density p=sum_{i=1}^n X_i g(X_i e^{lambda_i phi_i}), where X_i=-(nabla phi_i)^2/2 is the kinetic energy of the i-th field phi_i, lambda_i is a constant, and g is an arbitrary function in terms of Y_i=X_i e^{lambda_i phi_i}. This covers most of the scalar-field models of dark energy proposed in literature that possess scaling solutions. Using the bound coming from Big-Bang-Nucleosynthesis and the condition under which the each field cannot drive inflation as a single component of the universe, we find the following features: (i) a transient or eternal cosmic acceleration can be realized after the scaling matter era, (ii) a "thawing" property of assisting scalar fields is crucial to determine the evolution of the field equation of state w_{phi}, and (iii) the field equation of state today can be consistent with the observational bound w_{phi}<-0.8 in the presence of multiple scalar fields.

We study the growth of matter density perturbations delta_m for a number of viable f(R) gravity models that satisfy both cosmological and local gravity constraints, where the Lagrangian density f is a function of the Ricci scalar R. If the parameter m=Rf_{,RR}/f_{,R} today is larger than the order of 10^{-6}, linear perturbations relevant to the matter power spectrum evolve with a growth rate s=d (ln delta_m)/d (ln a) (a is the scale factor) that is larger than in the LCDM model. We find the window in the free parameter space of our models for which spatial dispersion of the growth index gamma_0= gamma(z=0) (z is the redshift) appears in the range of values 0.40< gamma_0<0.55, as well as the region in parameter space for which there is essentially no dispersion and gamma_0 converges to values around 0.40<gamma_0<0.43. These latter values are much lower than in the LCDM model. We show that these unusual dispersed or converged spectra are present in most of the viable f(R) models with m(z=0) larger than the order of 10^{-6}. These properties will be essential in the quest for f(R) modified gravity models using future high-precision observations and they confirm the possibility to distinguish clearly most of these models from the LCDM model.

We discuss solar system constraints on f(G) gravity models, where f is a function of the Gauss-Bonnet term G. We focus on cosmologically viable f(G) models that can be responsible for late-time cosmic acceleration. These models generally give rise to corrections of the form epsilon*(r/rs)^p to the vacuum Schwarzschild solution, where epsilon = H^2 rs^2 << 1, rs is the Schwarzschild radius of Sun, and H is the Hubble parameter today. We generally estimate the strength of modifications to General Relativity in order to confront models with a number of experiments such as the deflection of light and the perihelion shift. We show that cosmologically viable f(G) models can be consistent with solar system constraints for a wide range of model parameters.

We study the field profile of a scalar field $\phi$ that couples to a matter fluid (dubbed a chameleon field) in the relativistic gravitational background of a spherically symmetric spacetime. Employing a linear expansion in terms of the gravitational potential $\Phi_c$ at the surface of a compact object with a constant density, we derive the thin-shell field profile both inside and outside the object, as well as the resulting effective coupling with matter, analytically. We also carry out numerical simulations for the class of inverse power-law potentials $V(\phi)=M^{4+n} \phi^{-n}$ by employing the information provided by our analytical solutions to set the boundary conditions around the centre of the object and show that thin-shell solutions in fact exist if the gravitational potential $\Phi_c$ is smaller than 0.3, which marginally covers the case of neutron stars. Thus the chameleon mechanism is present in the relativistic gravitational backgrounds, capable of reducing the effective coupling. Since thin-shell solutions are sensitive to the choice of boundary conditions, our analytic field profile is very helpful to provide appropriate boundary conditions for $\Phi_c \lesssim O(0.1)$.

We derive conditions under which f(G) gravity models, whose Lagrangian densities f are written in terms of a Gauss-Bonnet term G, are cosmologically viable. The most crucial condition to be satisfied is that f_GG, the second derivative of f with respect to G, must be positive, which is required to ensure the stability of a late-time de-Sitter solution as well as the existence of standard radiation/matter dominated epochs. We present a number of explicit f(G) models in which a cosmic acceleration is followed by the matter era. We find that the equation of state of dark energy can cross the phantom divide before reaching the present Universe. The viable models have asymptotic behavior f_GG goes to +0 when |G| goes to infinity, in which case a rapid oscillation of perturbations occurs unless such an oscillating degree of freedom is suppressed relative to a homogeneous mode in the early universe. We also introduce an iterative method to avoid numerical instabilities associated with a large mass of the oscillating mode.

We study second-order cosmological perturbations in scalar-tensor models of dark energy that satisfy local gravity constraints, including f(R) gravity. We derive equations for matter fluctuations under a sub-horizon approximation and clarify conditions under which first-order perturbations in the scalar field can be neglected relative to second-order matter and velocity perturbations. We also compute the skewness of the matter density distribution and find that the difference from the LCDM model is only less than a few percent even if the growth rate of first-order perturbations is significantly different from that in the LCDM model. This shows that the skewness provides a model-independent test for the picture of gravitational instability from Gaussian initial perturbations including scalar-tensor modified gravity models.

We derive analytic solutions of a chameleon scalar field $\phi$ that couples to a non-relativistic matter in the weak gravitational background of a spherically symmetric body, paying particular attention to a field mass $m_A$ inside of the body. The standard thin-shell field profile is recovered by taking the limit $m_A*r_c \to \infty$, where $r_c$ is a radius of the body. We show the existence of "no-shell" solutions where the field is nearly frozen in the whole interior of the body, which does not necessarily correspond to the "zero-shell" limit of thin-shell solutions. In the no-shell case, under the condition $m_A*r_c \gg 1$, the effective coupling of $\phi$ with matter takes the same asymptotic form as that in the thin-shell case. We study experimental bounds coming from the violation of equivalence principle as well as solar-system tests for a number of models including $f(R)$ gravity and find that the field is in either the thin-shell or the no-shell regime under such constraints, depending on the shape of scalar-field potentials. We also show that, for the consistency with local gravity constraints, the field at the center of the body needs to be extremely close to the value $\phi_A$ at the extremum of an effective potential induced by the matter coupling.

We study the effect of modified gravity on weak lensing in a class of scalar-tensor theory that includes $f(R)$ gravity as a special case. These models are designed to satisfy local gravity constraints by having a large scalar-field mass in a region of high curvature. Matter density perturbations in these models are enhanced at small redshifts because of the presence of a coupling $Q$ that characterizes the strength between dark energy and non-relativistic matter. We compute a convergence power spectrum of weak lensing numerically and show that the spectral index and the amplitude of the spectrum in the linear regime can be significantly modified compared to the $\Lambda$CDM model for large values of $|Q|$ of the order of unity. Thus weak lensing provides a powerful tool to constrain such large coupling scalar-tensor models including $f(R)$ gravity.

We estimate the strength of large-scale magnetic fields produced during inflation in the framework of Dirac-Born-Infeld (DBI) theories. This analysis is sufficiently general in the sense that it covers most of conformal symmetry breaking theories in which the electromagnetic field is coupled to a scalar field. In DBI theories there is an additional factor associated with the speed of sound, which allows a possibility to lead to an extra amplification of the magnetic field in a ultra-relativistic region. We clarify the conditions under which seed magnetic fields to feed the galactic dynamo mechanism at a decoupling epoch as well as present magnetic fields on galactic scales are sufficiently generated to satisfy observational bounds.

We study constraints on f(R) dark energy models from solar system experiments combined with experiments on the violation of equivalence principle. When the mass of an equivalent scalar field degree of freedom is heavy in a region with high density, a spherically symmetric body has a thin-shell so that an effective coupling of the fifth force is suppressed through a chameleon mechanism. We place experimental bounds on the cosmologically viable models recently proposed in literature which have an asymptotic form f(R)=R-lambda R_c [1-(R_c/R)^{2n}] in the regime R >> R_c. From the solar-system constraints on the post-Newtonian parameter gamma, we derive the bound n>0.5, whereas the constraints from the violations of weak and strong equivalence principles give the bound n>0.9. This allows a possibility to find the deviation from the LambdaCDM cosmological model. For the model f(R)=R-lambda R_c(R/R_c)^p with 0<p<1 the severest constraint is found to be p<10^{-10}, which shows that this model is hardly distinguishable from the LambdaCDM cosmology.

We construct a family of viable scalar-tensor models of dark energy (DE) which possess a phase of late-time acceleration preceded by a standard matter era, while at the same time satisfying the local gravity constraints (LGC). The coupling Q between the scalar field and the non-relativistic matter in the Einstein frame is assumed to be constant in our scenario, which is a generalization of f(R) gravity theories corresponding to the coupling Q=-1/sqrt{6}. We find that these models can be made compatible with local gravity constraints even when |Q| is of the order of unity through a chameleon mechanism, if the scalar-field potential is chosen to have a sufficiently large mass in the high-curvature regions. We show that these models generally lead to the divergence of the equation of state of DE, which occur at smaller redshifts as the deviation from the LambdaCDM model become more significant. We also study the evolution of matter density perturbations and employ them to place bounds on the coupling |Q| as well as model parameters of the field potential from observations of the matter power spectrum and the CMB anisotropies. We find that, as long as |Q| is smaller than the order of unity, there exist allowed parameter regions that are consistent with both observational and local gravity constraints.

We make a detailed study of matter density perturbations in both metric and Palatini formalisms in theories whose Lagrangian density is a general function, f(R), of the Ricci scalar. We derive these equations in a number of gauges. We show that for viable models that satisfy cosmological and local gravity constraints (LGC), matter perturbation equations derived under a sub-horizon approximation are valid even for super-Hubble scales provided the oscillating mode (scalaron) does not dominate over the matter-induced mode. Such approximate equations are especially reliable in the Palatini formalism because of the absence of scalarons. Using these equations we make a comparative study of the behaviour of density perturbations as well as gravitational potentials for a number of classes of theories. In the metric formalism the parameter m=Rf_{,RR}/f_{,R} characterising the deviation from the Lambda CDM model is constrained to be very small during the matter era in order to ensure compatibility with LGC, but the models in which m grows to the order of 10^{-1} around the present epoch can be allowed. These models also suffer from an additional fine tuning due to the presence of scalaron modes which are absent in the Palatini case. In Palatini formalism LGC and background cosmological constraints provide only weak bounds on |m| by constraining it to be smaller than ~ 0.1. This is in contrast to matter density perturbations which, on galactic scales, place far more stringent constraints on the present deviation parameter m of the order of |m| < 10^{-5} - 10^{-4}. This is due to the peculiar evolution of matter perturbations in the Palatini case which exhibits a rapid growth or a damped oscillation depending on the sign of m.

In a recent paper [1] (PRL 98,131302,2007) we have shown that f(R)=R + mu R^{n} modified gravity dark energy models are not cosmologically viable because during the matter era that precedes the accelerated stage the cosmic expansion is given by a sim t^{1/2} rather than a sim t^{2/3}, where a is a scale factor and t is the cosmic time. A recent work [2] (PLB 639,135-143,2006) by Capozziello et al. criticised our results presenting some apparent counter-examples to our claim in f(R)= mu R^{n} models. We show here that those particular R^{n} models can produce an expansion as a sim t^{2/3} but this does not connect to a late-time acceleration. Hence, though acceptable f(R) dark energy models may exist, the R^{n} models presented in Capozziello et al. are not cosmologically viable, confirming our previous results in Ref. [1].

We identify the class of f(R) dark energy models which have a viable cosmology, i.e. a matter dominated epoch followed by a late-time acceleration. The deviation from a LambdaCDM model (f=R-Lambda) is quantified by the function m=Rf_{,RR}/f_{,R}. The matter epoch corresponds to m(r=-1) simeq +0 (where r=-Rf_{,R}/f) while the accelerated attractor exists in the region 0<m<1. We find that the equation of state w_DE of all such ``viable'' f(R) models exhibits two features: w_DE diverges at some redshift z_{c} and crosses the cosmological constant boundary (``phantom crossing'') at a redshift z_{b} smaller than z_{c}. Using the observational data of Supernova Ia and Cosmic Microwave Background, we obtain the constraint m<O(0.1) and we find that the phantom crossing could occur at z_{b}>1, i.e. within reach of observations. If we add local gravity constraints, the bound on m becomes very stringent, with m several orders of magnitude smaller than unity in the region whose density is much larger than the present cosmological density. The representative models that satisfy both cosmological and local gravity constraints take the asymptotic form m(r)=C(-r-1)^p with p>1 as r approaches -1.

We discuss observational consequences of f(R) dark energy scenarios that satisfy local gravity constraints (LGC) as well as conditions of the cosmological viability. The model we study is given by m(r)=C(-r-1)^p (C>0, p>1) with m=Rf_{,RR}/f_{,R} and r=-Rf_{,R}/f, which cover viable f(R) models proposed so far in a high-curvature region designed to be compatible with LGC. The equation of state of dark energy exhibits a divergence at a redshift z_c that can be as close as a few while satisfying sound horizon constraints of Cosmic Microwave Background (CMB). We study the evolution of matter density perturbations in details and place constraints on model parameters from the difference of spectral indices of power spectra between CMB and galaxy clustering. The models with p>5 can be consistent with those observational constraints as well as LGC. We also discuss the evolution of perturbations in the Ricci scalar R and show that an oscillating mode (scalaron) can easily dominate over a matter-induced mode as we go back to the past. This violates the stability of cosmological solutions, thus posing a problem about how the over-production of scalarons should be avoided in the early universe.

We derive the equation of matter density perturbations on sub-horizon scales for a general Lagrangian density f(R, phi, X) that is a function of a Ricci scalar R, a scalar field phi and a kinetic term X=-(nabla phi)^2/2. This is useful to constrain modified gravity dark energy models from observations of large-scale structure and weak lensing. We obtain the solutions for the matter perturbation delta_m as well as the gravitational potential Phi for some analytically solvable models. In a f(R) dark energy model with the Lagrangian density f(R)=alpha R^{1+m}-Lambda, the growth rates of perturbations exhibit notable differences from those in the standard Einstein gravity unless m is very close to 0. In scalar-tensor models with the Lagrangian density f=F(phi)R+2p(phi,X) we relate the models with coupled dark energy scenarios in the Einstein frame and reproduce the equations of perturbations known in the current literature by making a conformal transformation. We also estimate the evolution of perturbations in both Jordan and Einstein frames when the energy fraction of dark energy is constant during the matter-dominated epoch.

We study D-brane inflation in a warped conifold background that includes brane-position dependent corrections for the nonperturbative superpotential. Instead of stabilizing the volume modulus chi at instantaneous minima of the potential and studying the inflation dynamics with an effective single field (radial distance between a brane and an anti-brane) phi, we investigate the multi-field inflation scenario involving these two fields. The two-field dynamics with the potential V(phi,chi) in this model is significantly different from the effective single-field description in terms of the field phi when the field chi is integrated out. The latter picture underestimates the total number of e-foldings even by one order of magnitude. We show that a correct single-field description is provided by a field psi obtained from a rotation in the two-field space along the background trajectory. This model can give a large number of e-foldings required to solve flatness and horizon problems at the expense of fine-tunings of model parameters. We also estimate the spectra of density perturbations and show that the slow-roll parameter eta_{psi psi}=M_{pl}^2 V_{,psi psi}/V in terms of the rotated field psi determines the spectral index of scalar metric perturbations. We find that it is generally difficult to satisfy, simultaneously, both constraints of the spectral index and the COBE normalization, while the tensor to scalar ratio is sufficiently small to match with observations.

Phantom fields with negative kinetic energy are often plagued by the vacuum quantum instability in the ultraviolet region. We present a Lorentz-violating dark energy model free from this problem and show that the crossing of the cosmological constant boundary w=-1 to the phantom equation of state is realized before reaching a de Sitter attractor. Another interesting feature is a peculiar time-dependence of the effective Newton's constant; the magnitude of this effect is naturally small but may be close to experimental limits. We also derive momentum scales of instabilities at which tachyons or ghosts appear in the infrared region around the present Hubble scale and clarify the conditions under which tachyonic instabilities do not spoil homogeneity of the present/future Universe.

We place observational constraints on a coupling between dark energy and dark matter by using 71 Type Ia supernovae (SNe Ia) from the first year of the five-year Supernova Legacy Survey (SNLS), the cosmic microwave background (CMB) shift parameter from the three-year Wilkinson Microwave Anisotropy Probe (WMAP), and the baryon acoustic oscillation (BAO) peak found in the Sloan Digital Sky Survey (SDSS). The interactions we study are (i) constant coupling delta and (ii) varying coupling delta(z) that depends on a redshift z, both of which have simple parametrizations of the Hubble parameter to confront with observational data. We find that the combination of the three databases marginalized over a present dark energy density gives stringent constraints on the coupling, -0.08 < delta < 0.03 (95% CL) in the constant coupling model and -0.4 < delta_0 < 0.1 (95% CL) in the varying coupling model, where delta_0 is a present value. The uncoupled LambdaCDM model (w_X = -1 and delta = 0) still remains a good fit to the data, but the negative coupling (delta < 0) with the equation of state of dark energy w_X < -1 is slightly favoured over the LambdaCDM model.

We derive the conditions under which dark energy models whose Lagrangian densities f are written in terms of the Ricci scalar R are cosmologically viable. We show that the cosmological behavior of f(R) models can be understood by a geometrical approach consisting in studying the m(r) curve on the (r, m) plane, where m=Rf_{,RR}/f_{,R} and r=-Rf_{,R}/f with f_{,R}=df/dR. This allows us to classify the f(R) models into four general classes, depending on the existence of a standard matter epoch and on the final accelerated stage. The existence of a viable matter dominated epoch prior to a late-time acceleration requires that the variable m satisfies the conditions m(r) approx+0 and dm/dr>-1 at r approx-1. For the existence of a viable late-time acceleration we require instead either (i) m=-r-1, (sqrt{3}-1)/2<m<1 and dm/dr<-1 or (ii) 0<m<1 at r=-2. These conditions identify two regions in the (r, m) space, one for the matter era and the other for the acceleration. Only models with a m(r) curve that connects these regions and satisfy the requirements above lead to an acceptable cosmology. The models of the type f(R)=alpha R^{-n} and f=R+alpha R^{-n} do not satisfy these conditions for any n>0 and n<-1 and are thus cosmologically unacceptable. Similar conclusions can be reached for many other examples discussed in the text. In most cases the standard matter era is replaced by a cosmic expansion with scale factor a propto t^{1/2}. We also find that f(R) models can have a strongly phantom attractor but in this case there is no acceptable matter era.

All $f(R)$ modified gravity theories are conformally identical to models of quintessence in which matter is coupled to dark energy with a strong coupling. This coupling induces a cosmological evolution radically different from standard cosmology. We find that in all $f(R)$ theories that behave as a power of $R$ at large or small $R$ (which include most of those proposed so far in the literature) the scale factor during the matter phase grows as $t^{1/2}$ instead of the standard law $t^{2/3}$. This behaviour is grossly inconsistent with cosmological observations (e.g. WMAP), thereby ruling out these models even if they pass the supernovae test and can escape the local gravity constraints.

We make a systematic study of the cosmological dynamics for a number of f(R) gravity theories in Palatini formalism. We find a number of interesting results: (i) models based on theories of the type (a) f(R)=R-beta R^n and (b) f(R)=R+alpha ln R -beta, unlike the metric formalism, are capable of producing the sequence of radiation-dominated, matter-dominated and de-Sitter periods, and (ii) models based on theories of the type (c) f(R)=R+alpha R^m -beta R^n can produce early as well as late accelerating phases. However for the classes of models considered here, we have been unable to find the sequence of all four dynamical epochs required to account for the complete cosmological dynamics, even though three out of four phases are possible. We also place observational constraints on these models using the recently released supernovae data (SNLS) as well as the baryon acoustic oscillation peak and the CMB shift parameter. The best-fit values are found to be n=0.027, alpha=4.63 for the models based on (a) and alpha=0.11, beta=4.62 for the models based on (b), neither of which are significantly preferred over the LCDM model. The models based on (c) are also consistent with the data with suitable choices of their parameters.

We discuss the realization of inflation and resulting cosmological perturbations in the low-energy effective string theory. In order to obtain nearly scale-invariant spectra of density perturbations and a suppressed tensor-to-scalar ratio, it is generally necessary that the dilaton field $\phi$ is effectively decoupled from gravity together with the existence of a slowly varying dilaton potential. We also study the effect of second-order corrections to the tree-level action which are the sum of a Gauss-Bonnet term coupled to $\phi$ and a kinetic term $(\nabla \phi)^4$. We find that it is possible to realize observationally supported spectra of scalar and tensor perturbations provided that the correction is dominated by the $(\nabla \phi)^4$ term even in the absence of the dilaton potential. When the Gauss-Bonnet term is dominant, tensor perturbations exhibit violent negative instabilities on small-scales about a de Sitter background in spite of the fact that scale-invariant scalar perturbations can be achieved.

The Gauss-Bonnet (GB) curvature invariant coupled to a scalar field $\phi$ can lead to an exit from a scaling matter-dominated epoch to a late-time accelerated expansion, which is attractive to alleviate the coincident problem of dark energy. We derive the condition for the existence of cosmological scaling solutions in the presence of the GB coupling for a general scalar-field Lagrangian density $p(\phi, X)$, where $X=-(1/2)(\nabla \phi)^2$ is a kinematic term of the scalar field. The GB coupling and the Lagrangian density are restricted to be in the form $f(\phi) \propto e^{\lambda \phi}$ and $p=Xg (Xe^{\lambda \phi})$, respectively, where $\lambda$ is a constant and $g$ is an arbitrary function. We also derive fixed points for such a scaling Lagrangian with a GB coupling $f(\phi) \propto e^{\mu \phi}$ and clarify the conditions under which the scaling matter era is followed by a de-Sitter solution which can appear in the presence of the GB coupling. Among scaling models proposed in the current literature, we find that the models which allow such a cosmological evolution are an ordinary scalar field with an exponential potential and a tachyon field with an inverse square potential, although the latter requires a coupling between dark energy and dark matter.

We study cosmologies based on low-energy effective string theory with higher-order string corrections to a tree-level action and with a modulus scalar field (dilaton or compactification modulus). In the presence of such corrections it is possible to construct nonsingular cosmological solutions in the context of Pre-Big-Bang and Ekpyrotic universes. We review the construction of nonsingular bouncing solutions and resulting density perturbations in Pre-Big-Bang and Ekpyrotic models. We also discuss the effect of higher-order string corrections on dark energy universe and show several interesting possibilities of the avoidance of future singularities.

In this paper we review in detail a number of approaches that have been adopted to try and explain the remarkable observation of our accelerating Universe. In particular we discuss the arguments for and recent progress made towards understanding the nature of dark energy. We review the observational evidence for the current accelerated expansion of the universe and present a number of dark energy models in addition to the conventional cosmological constant, paying particular attention to scalar field models such as quintessence, K-essence, tachyon, phantom and dilatonic models. The importance of cosmological scaling solutions is emphasized when studying the dynamical system of scalar fields including coupled dark energy. We study the evolution of cosmological perturbations allowing us to confront them with the observation of the Cosmic Microwave Background and Large Scale Structure and demonstrate how it is possible in principle to reconstruct the equation of state of dark energy by also using Supernovae Ia observational data. We also discuss in detail the nature of tracking solutions in cosmology, particle physics and braneworld models of dark energy, the nature of possible future singularities, the effect of higher order curvature terms to avoid a Big Rip singularity, and approaches to modifying gravity which leads to a late-time accelerated expansion without recourse to a new form of dark energy.

We consider the fate of future singularities in the effective dynamics of loop quantum cosmology. Non-perturbative quantum geometric effects which lead to $\rho^2$ modification of the Friedmann equation at high energies result in generic resolution of singularities whenever energy density $\rho$ diverges at future singularities of Friedmann dynamics. Such quantum effects lead to the avoidance of a Big Rip, which is followed by a recollapsing universe stable against perturbations. Resolution of sudden singularity, the case when pressure diverges but energy density approaches a finite value depends on the ratio of the latter to a critical energy density of the order of Planck. If the value of this ratio is greater than unity, the universe escapes the sudden future singularity and becomes oscillatory.

A cosmological model that aims at solving the coincidence problem should show that dark energy and dark matter follow the same scaling solution from some time onward. At the same time, the model should contain a sufficiently long matter-dominated epoch that takes place before acceleration in order to guarantee a decelerated epoch and structure formation. So a successful cosmological model requires the occurrence of a sequence of epochs, namely a radiation era, a matter-dominated era and a final accelerated scaling attractor with $\Omega_{\phi} \simeq 0.7$. In this paper we derive the generic form of a scalar-field Lagrangian that possesses scaling solutions in the case where the coupling $Q$ between dark energy and dark matter is a free function of the field $\phi$. We then show, rather surprisingly, that the aforementioned sequence of epochs cannot occur for a vast class of generalized coupled scalar field Lagrangians that includes, to our knowledge, all scaling models in the current literature.

We study general properties of attractors for scalar-field dark energy scenarios which possess cosmological scaling solutions. In all such models there exists a scalar-field dominant solution with an energy fraction \Omega_{\phi}=1 together with a scaling solution. A general analytic formula is given to derive fixed points relevant to dark energy coupled to dark matter. We investigate the stability of fixed points without specifying the models of dark energy in the presence of non-relativistic dark matter and provide a general proof that a non-phantom scalar-field dominant solution is unstable when a stable scaling solution exists in the region \Omega_{\phi}<1. A phantom scalar-field dominant fixed point is found to be classically stable. We also generalize the analysis to the case of multiple scalar fields and show that for a non-phantom scalar field assisted acceleration always occurs for all scalar-field models which have scaling solutions. For a phantom field the equation of state approaches that of cosmological constant as we add more scalar fields.

We study inflation arising from the motion of a BPS D3-brane in the background of a stack of k parallel D5-branes. There are two scalar fields in this set up-- (i) the radion field R, a real scalar field, and (ii) a complex tachyonic scalar field chi living on the world volume of the open string stretched between the D3 and D5 branes. We find that inflation is realized by the potential of the radion field, which satisfies observational constraints coming from the Cosmic Microwave Background. After the radion becomes of order the string length scale l_s, the dynamics is governed by the potential of the complex scalar field. Since this field has a standard kinematic term, reheating can be successfully realized by the mechanism of tachyonic preheating with spontaneous symmetry breaking.

We review the theory of inflation with single and multiple fields paying particular attention to the dynamics of adiabatic and entropy/isocurvature perturbations which provide the primary means of testing inflationary models. We review the theory and phenomenology of reheating and preheating after inflation providing a unified discussion of both the gravitational and nongravitational features of multi-field inflation. In addition we cover inflation in theories with extra dimensions and models such as the curvaton scenario and modulated reheating which provide alternative ways of generating large-scale density perturbations. Finally we discuss the interesting observational implications that can result from adiabatic-isocurvature correlations and non-Gaussianity.

We study chaotic dynamics in preheating after inflation in which an inflaton $\phi$ is coupled to another scalar field $\chi$ through an interaction $(1/2)g^2\phi^2\chi^2$. We first estimate the size of the quasi-homogeneous field $\chi$ at the beginning of reheating for large-field inflaton potentials $V(\phi)=V_0\phi^n$ by evaluating the amplitude of the $\chi$ fluctuations on scales larger than the Hubble radius at the end of inflation. Parametric excitations of the field $\chi$ during preheating can give rise to chaos between two dynamical scalar fields. For the quartic potential ($n=4$, $V_0=\lambda/4$) chaos actually occurs for $g^2/\lambda <{\cal O}(10)$ in a linear regime before which the backreaction of created particles becomes important. This analysis is supported by several different criteria for the existence of chaos. For the quadratic potential ($n=2$) the signature of chaos is not found by the time at which the backreaction begins to work, similar to the case of the quartic potential with $g^2/\lambda \gg 1$.

The reconstruction of scalar-field dark energy models is studied for a general Lagrangian density $p(\phi, X)$, where $X$ is a kinematic term of a scalar field $\phi$. We implement the coupling $Q$ between dark energy and dark matter and express reconstruction equations using two observables: the Hubble parameter $H$ and the matter density perturbation $\delta_m$. This allows us to determine the structure of corresponding theoretical Lagrangian together with the coupling $Q$ from observations. We apply our formula to several forms of Lagrangian and present concrete examples of reconstruction by using the recent Gold dataset of supernovae measurements. This analysis includes a generalized ghost condensate model as a way to cross a cosmological-constant boundary even for a single-field case.

We study the motion of a BPS D3-brane in the NS5-brane ring background. The radion field becomes tachyonic in this geometrical set up. We investigate the potential of this geometrical tachyon in the cosmological scenario for inflation as well as dark energy. We evaluate the spectra of scalar and tensor perturbations generated during tachyon inflation and show that this model is compatible with recent observations of Cosmic Microwave Background (CMB) due to an extra freedom of the number of NS5-branes. It is not possible to explain the origin of both inflation and dark energy by using a single tachyon field, since the energy density at the potential minimum is not negligibly small because of the amplitude of scalar perturbations set by CMB anisotropies. However geometrical tachyon can account for dark energy when the number of NS5-branes is large, provided that inflation is realized by another scalar field.

We study the cosmological evolution based upon a $D$-dimensional action in low-energy effective string theory in the presence of second-order curvature corrections and a modulus scalar field (dilaton or compactification modulus). A barotropic perfect fluid coupled to the scalar field is also allowed. Phase space analysis and the stability of asymptotic solutions are performed for a number of models which include ($i$) fixed scalar field, ($ii$) linear dilaton in string frame, and ($iii$) logarithmic modulus in Einstein frame. We confront analytical solutions with observational constraints for deceleration parameter and show that Gauss-Bonnet gravity (with no matter fields) may not explain the current acceleration of the universe. We also study the future evolution of the universe using the GB parametrization and find that big rip singularities can be avoided even in the presence of a phantom fluid because of the balance between the fluid and curvature corrections. A non-minimal coupling between the fluid and the modulus field also opens up the interesting possibility to avoid big rip regardless of the details of the fluid equation of state.

We explore the dynamics of assisted quintessence, where more than one scalar field is present with the same potential. For potentials with tracking solutions, the fields naturally approach the same values; in the context of inflation this leads to the assisted inflation phenomenon where several fields can cooperate to drive a period of inflation though none is able to individually. For exponential potentials, we study the fixed points and their stability confirming results already in the literature, and then carry out a numerical analysis to show how assisted quintessence is realized. For inverse power-law potentials, we find by contrast that there is no assisted behaviour; indeed those are the unique (monotonic) potentials where several fields together behave just as a single field in the same potential. More generally, we provide an algorithm for generating a single-field potential giving equivalent dynamics to multi-field assisted quintessence.

We study the evolution of (phantom) dark energy universe by taking into account the higher-order string corrections to Einstein-Hilbert action with a fixed dilaton. While the presence of a cosmological constant gives stable de-Sitter fixed points in the cases of heterotic and bosonic strings, no stable de-Sitter solutions exist when a phantom fluid is present. We find that the universe can exhibit a Big Crunch singularity with a finite time for type II string, whereas it reaches a Big Rip singularity for heterotic and bosonic strings. Thus the fate of dark energy universe crucially depends upon the type of string theory under consideration.

Cosmological scaling solutions are particularly important in solving the coincidence problem of dark energy. We derive the equations of sub-Hubble linear matter perturbations for a general scalar-field Lagrangian--including quintessence, tachyon, dilatonic ghost condensate and k-essence--and solve them analytically for scaling solutions. We find that matter perturbations are always damped if a phantom field is coupled to dark matter and identify the cases in which the gravitational potential is constant. This provides an interesting possibility to place stringent observational constraints on scaling dark energy models.

In dark energy models of scalar-field coupled to a barotropic perfect fluid, the existence of cosmological scaling solutions restricts the Lagrangian of the field $\vp$ to $p=X g(Xe^{\lambda \vp})$, where $X=-g^{\mu\nu} \partial_\mu \vp \partial_\nu \vp /2$, $\lambda$ is a constant and $g$ is an arbitrary function. We derive general evolution equations in an autonomous form for this Lagrangian and investigate the stability of fixed points for several different dark energy models--(i) ordinary (phantom) field, (ii) dilatonic ghost condensate, and (iii) (phantom) tachyon. We find the existence of scalar-field dominant fixed points ($\Omega_\vp=1$) with an accelerated expansion in all models irrespective of the presence of the coupling $Q$ between dark energy and dark matter. These fixed points are always classically stable for a phantom field, implying that the universe is eventually dominated by the energy density of a scalar field if phantom is responsible for dark energy. When the equation of state $w_\vp$ for the field $\vp$ is larger than -1, we find that scaling solutions are stable if the scalar-field dominant solution is unstable, and vice versa. Therefore in this case the final attractor is either a scaling solution with constant $\Omega_\vp$ satisfying $0<\Omega_\vp<1$ or a scalar-field dominant solution with $\Omega_\vp=1$.

The properties of future singularities are investigated in the universe dominated by dark energy including the phantom-type fluid. We classify the finite-time singularities into four classes and explicitly present the models which give rise to these singularities by assuming the form of the equation of state of dark energy. We show the existence of a stable fixed point with an equation of state $w<-1$ and numerically confirm that this is actually a late-time attractor in the phantom-dominated universe. We also construct a phantom dark energy scenario coupled to dark matter that reproduces singular behaviors of the Big Rip type for the energy density and the curvature of the universe. The effect of quantum corrections coming from conformal anomaly can be important when the curvature grows large, which typically moderates the finite-time singularities.

We study a dark energy scenario in the presence of a tachyon field $\phi$ with potential $V(\phi)$ and a barotropic perfect fluid. The cosmological dynamics crucially depends on the asymptotic behavior of the quantity $\lambda=-M_pV_\phi/V^{3/2}$. If $\lambda$ is a constant, which corresponds to an inverse square potential $V(\phi) \propto \phi^{-2}$, there exists one stable critical point that gives an acceleration of the universe at late times. When $\lambda \to 0$ asymptotically, we can have a viable dark energy scenario in which the system approaches an ``instantaneous'' critical point that dynamically changes with $\lambda$. If $|\lambda|$ approaches infinity asymptotically, the universe does not exhibit an acceleration at late times. In this case, however, we find an interesting possibility that a transient acceleration occurs in a regime where $|\lambda|$ is smaller than of order unity.

We study the variation of the effective fine structure constant alpha for Dirac-Born-Infeld (DBI) type dark energy models. The DBI action based on string theory naturally gives rise to a coupling between gauge fields and a scalar field responsible for accelerated expansion of the universe. This leads to the change of alpha due to a dynamical evolution of the scalar field, which can be compatible with the recently observed cosmological data around the redshift $\tilde{z} \lesssim 3$. We place constraints on several different DBI models including exponential, inverse power-law and rolling massive scalar potentials. We find that these models can satisfy the varying alpha constraint provided that mass scales of the potentials are fine-tuned. When we adopt the mass scales which are motivated by string theory, both exponential and inverse power-law potentials give unacceptably large change of alpha, thus ruled out from observations. On the other hand the rolling massive scalar potential is compatible with the observationally allowed variation of alpha. Therefore the information of varying alpha provides a powerful way to distinguish between a number of string-inspired DBI dark energy models.

We study constraints on a number of patch inflationary models in noncommutative spacetime using a compilation of recent high-precision observational data. In particular, the four-dimensional General Relativistic (GR) case, the Randall-Sundrum (RS) and Gauss-Bonnet (GB) braneworld scenarios are investigated by extending previous commutative analyses to the infrared limit of a maximally symmetric realization of the stringy uncertainty principle. The effect of spacetime noncommutativity modifies the standard consistency relation between the tensor spectral index and the tensor-to-scalar ratio. We perform likelihood analyses in terms of inflationary observables using new consistency relations and confront them with large-field inflationary models with potential V \propto \vp^p in two classes of noncommutative scenarios. We find a number of interesting results: (i) the quartic potential (p=4) is rescued from marginal rejection in the class 2 GR case, and (ii) steep inflation driven by an exponential potential (p \to \infty) is allowed in the class 1 RS case. Spacetime noncommutativity can lead to blue-tilted scalar and tensor spectra even for monomial potentials, thus opening up a possibility to explain the loss of power observed in the cosmic microwave background anisotropies. We also explore patch inflation with a Dirac-Born-Infeld tachyon field and explicitly show that the associated likelihood analysis is equivalent to the one in the ordinary scalar field case by using horizon-flow parameters. It turns out that tachyon inflation is compatible with observations in all patch cosmologies even for large p.

We study the generation of electromagnetic fields in a string-inspired scenario associated with a rolling massive scalar field $\phi$ on the anti-D3 branes of KKLT de Sitter vacua. The 4-dimensional DBI type effective action naturally gives rise to the coupling between the gauge fields and the inflaton $\phi$, which leads to the production of cosmological magnetic fields during inflation due to the breaking of conformal invariance. We find that the amplitude of magnetic fields at decoupling epoch can be larger than the limiting seed value required for the galactic dynamo. We also discuss the mechanism of reheating in our scenario and show that gauge fields are sufficiently enhanced for the modes deep inside the Hubble radius with an energy density greater than that of the inflaton.

In loop quantum cosmology, the universe avoids a big bang singularity and undergoes an early and short super-inflation phase. During super-inflation, non-perturbative quantum corrections to the dynamics drive an inflaton field up its potential hill, thus setting the initial conditions for standard inflation. We show that this effect can raise the inflaton high enough to achieve sufficient e-foldings in the standard inflation era. We analyze the cosmological perturbations generated when slow-roll is violated after super-inflation, and show that loop quantum effects can in principle leave an indirect signature on the largest scales in the CMB, with some loss of power and running of the spectral index.

Our ignorance about the source of cosmic acceleration has stimulated study of a wide range of models and modifications to gravity. Cosmological scaling solutions in any of these theories are privileged because they represent natural backgrounds relevant to dark energy. We study scaling solutions in a generalized background $H^2 \propto \rho_T^n$ in the presence of a scalar field $\vp$ and a barotropic perfect fluid, where $H$ is a Hubble rate and $\rho_T$ is a total energy density. The condition for the existence of scaling solutions restricts the form of Lagrangian to be $p=X^{1/n}g(Xe^{n\lambda \vp})$, where $X=-g^{\mu\nu} \partial_\mu \vp \partial_\nu \vp /2$ and $g$ is an arbitrary function. This is very useful to find out scaling solutions and corresponding scalar-field potentials in a broad class of dark energy models including (coupled)-quintessence, ghost-type scalar field, tachyon and k-essence. We analytically derive the scalar-field equation of state $w_\vp$ and the fractional density $\Omega_\vp$ and apply it to a number of dark energy models.

With the WMAP data we can now begin to test realistic models of inflation involving multiple scalar fields. These naturally lead to correlated adiabatic and isocurvature (entropy) perturbations with a running spectral index. We present the first full (9 parameter) likelihood analysis of double inflation with WMAP data and find that despite the extra freedom, supersymmetric hybrid potentials are strongly constrained with less than 7% correlated isocurvature component allowed when standard priors are imposed on the cosomological parameters. As a result we also find that Akaike & Bayesian model selection criteria rather strongly prefer single-field inflation, just as equivalent analysis prefers a cosmological constant over dynamical dark energy in the late universe. It appears that simplicity is the best guide to our universe.

High-energy modifications to general relativity introduce changes to the perturbations generated during inflation, and the latest high-precision cosmological data can be used to place constraints on such modified inflation models. Recently it was shown that Randall-Sundrum type braneworld inflation leads to tighter constraints on quadratic and quartic potentials than in general relativity. We investigate how this changes with a Gauss-Bonnet correction term, which can be motivated by string theory. Randall-Sundrum models preserve the standard consistency relation between the tensor spectral index and the tensor-to-scalar ratio. The Gauss-Bonnet term breaks this relation, and also modifies the dynamics and perturbation amplitudes at high energies. We find that the Gauss-Bonnet term tends to soften the Randall-Sundrum constraints. The observational compatibility of the quadratic potential is strongly improved. For a broad range of energy scales, the quartic potential is rescued from marginal rejection. Steep inflation driven by an exponential potential is excluded in the Randall-Sundrum case, but the Gauss-Bonnet term leads to marginal compatibility for sufficient e-folds.

We explore a dark energy model with a ghost scalar field in the context of the runaway dilaton scenario in low-energy effective string theory. We address the problem of vacuum stability by implementing higher-order derivative terms and show that a cosmologically viable model of ``phantomized'' dark energy can be constructed without violating the stability of quantum fluctuations. We also analytically derive the condition under which cosmological scaling solutions exist starting from a general Lagrangian including the phantom type scalar field. We apply this method to the case where the dilaton is coupled to non-relativistic dark matter and find that the system tends to become quantum mechanically unstable when a constant coupling is always present. Nevertheless, it is possible to obtain a viable cosmological solution in which the energy density of the dilaton eventually approaches the present value of dark energy provided that the coupling rapidly grows during the transition to the scalar field dominated era.

We investigate a string-inspired scenario associated with a rolling massive scalar field on D-branes and discuss its cosmological implications. In particular, we discuss cosmological evolution of the massive scalar field on the ant-D3 brane of KKLT vacua. Unlike the case of tachyon field, because of the warp factor of the anti-D3 brane, it is possible to obtain the required level of amplitude of density perturbations. We study the spectra of scalar and tensor perturbations generated during the rolling scalar inflation and show that our scenario satisfies the observational constraint coming from the Cosmic Microwave Background anisotropies and other observational data. We also implement the negative cosmological constant arising from the stabilization of the modulus fields in the KKLT vacua and find that this leads to a successful reheating in which the energy density of the scalar field effectively scales as a pressureless dust. The present dark energy can be also explained in our scenario provided that the potential energy of the massive rolling scalar does not exactly cancel with the amplitude of the negative cosmological constant at the potential minimum.

We study cosmological perturbations in generalized Einstein scenarios and show the equivalence of inflationary observables both in the Jordan frame and the Einstein frame. In particular the consistency relation relating the tensor-to-scalar ratio with the spectral index of tensor perturbations coincides with the one in Einstein gravity, which leads to the same likelihood results in terms of inflationary observables. We apply this formalism to nonminimally coupled chaotic inflationary scenarios with potential $V=c\phi^p$ and place constraints on the strength of nonminimal couplings using a compilation of latest observational data. In the case of the quadratic potential ($p=2$), the nonminimal coupling is constrained to be $\xi>-7.0 \times 10^{-3}$ for negative $\xi$ from the $1\sigma$ observational contour bound. Although the quartic potential ($p=4$) is under a strong observational pressure for $\xi=0$, this property is relaxed by taking into account negative nonminimal couplings. We find that inflationary observables are within the $1\sigma$ contour bound as long as $\xi<-1.7 \times 10^{-3}$. We also show that the $p \ge 6$ cases are disfavoured even in the presence of nonminimal couplings.

We obtain observational constraints on Randall--Sundrum type II braneworld inflation using a compilation of data including WMAP, the 2dF and latest SDSS galaxy redshift surveys. We place constraints on three classes of inflation models (large-field, small-field and hybrid models) in the high-energy regime, which exhibit different behaviour compared to the low-energy case. The quartic potential is outside the $2\sigma$ observational contour bound for a number of $e$-folds less than 60, and steep inflation driven by an exponential potential is excluded because of its high tensor-to-scalar ratio. It is more difficult to strongly constrain small-field and hybrid models due to additional freedoms associated with the potentials, but we obtain upper bounds for the energy scale of inflation and the model parameters in certain cases. We also discuss possible ways to break the degeneracy of consistency relations and inflationary observables.

We construct a string-inspired nonsingular cosmological scenario in which an inflaton field is driven up the potential before the graceful exit by employing a low-energy string effective action with an orbifold compactification. This sets up an initial condition for the inflaton to lead to a sufficient amount of $e$-foldings and to generate a nearly scale-invariant primordial density perturbation during slow-roll inflation. Our scenario provides an interesting possibility to explain a suppressed power spectrum at low multipoles due to the presence of the modulus-driven phase prior to slow-roll inflation and thus can leave a strong imprint of extra dimensions in observed CMB anisotropies.

Non-commutative inflation is a modification of standard general relativity inflation which takes into account some effects of the space-time uncertainty principle motivated by ideas from string theory. The corrections to the primordial power spectrum which arise in a model of power-law inflation lead to a suppression of power on large scales, and produce a spectral index that is blue on large scales and red on small scales. This suppression and running of the spectral index are not imposed ad hoc, but arise from an early-Universe stringy phenomenology. We show that it can account for some loss of power on the largest scales that may be indicated by recent WMAP data. Cosmic microwave background anisotropies carry a signature of these very early Universe corrections, and can be used to place constraints on the parameters appearing in the non-commutative model. Applying a likelihood analysis to the WMAP data, we find the best-fit value for the critical wavenumber k_* (which involves the string scale) and for the exponent p (which determines the power-law inflationary expansion). The best-fit value corresponds to a string length of L_s \sim 10^{-28} cm.

The density perturbations generated when the inflaton decay rate is perturbed by a light scalar field $\chi$ are studied. By explicitly solving the perturbation equations for the system of two scalar fields and radiation, we show that even in low energy-scale inflation nearly scale-invariant spectra of scalar perturbations with an amplitude set by observations are obtained through the conversion of $\chi$ fluctuations into adiabatic density perturbations. We demonstrate that the spectra depend on the average decay rate of the inflaton & on the inflaton fluctuations. We then apply this new mechanism to string cosmologies & generalized Einstein theories and discuss the conditions under which scale-invariant spectra are possible.

These lecture notes provide an introduction to cosmic inflation. In particular I will review the basic concepts of inflation, generation of density perturbations, and reheating after inflation.

We study nonsingular cosmological scenarios in a general $D$-dimensional string effective action of the dilaton-modulus-axion system in the presence of the matter source. In the standard dilatonic Brans-Dicke parameter ($\omega=-1$) with radiation, we analytically obtain singularity-free bouncing solutions where the universe starts out in a state with a finite curvature and evolves toward the weakly coupled regime. We apply this analytic method to the string-gas cosmology including the massive state in addition to the leading massless state (radiation), with and without the axion. We numerically find bouncing solutions which asymptotically approach an almost radiation-dominant universe with a decreasing curvature irrespective of the presence of the axion, implying that inclusion of the matter source is crucial for the existence of such solutions for $\omega=-1$. In the theories with $\omega \ne -1$, it is possible to obtain complete regular bouncing solutions with a finite dilaton and curvature in both past and future asymptotics for the general dimension, $D$. We also discuss the case where dilatonic higher-order corrections are involved to the tree-level effective action and demonstrate that the presence of axion/modulus fields and the matter source does not significantly affect the dynamics of the dilaton-driven inflation and the subsequent graceful exit.

We study string-gas cosmology in dilaton gravity, inspired by the fact that it naturally arises in a string theory context. Our main interest is the thermodynamical treatment of the string-gas and the resulting implications for the cosmology. Within an adiabatic approximation, thermodynamical equilibrium and a small, toroidal universe as initial conditions, we numerically solve the corresponding equations of motions in two different regimes describing the string-gas thermodynamics: (i) the Hagedorn regime, with a single scale factor, and (ii) an almost-radiation dominated regime, which includes the leading corrections due to the lightest Kaluza Klein and winding modes, with two scale factors. The scale factor in the Hagedorn regime exhibits very slow time evolution with nearly constant energy and negligible pressure. By contrast, in case (ii) we find interesting cosmological solutions where the large dimensions continue to expand and the small ones are kept undetectably small.

We show explicitly some exciting features of double-inflation: (i) it can often lead to strongly correlated adiabatic and entropy (isocurvature) power spectra. (ii) The two-field slow-roll consistency relations can be violated when the correlation is large at Hubble crossing. (iii) The spectra of adiabatic and entropy perturbations can be strongly scale-dependent and tilted toward either the red or blue. These effects are typically due to a light or time-dependent entropy mass and a non-negligible angular velocity in field space during inflation. They are illustrated via a multi-parameter numerical search for correlations in two concrete models. The correlation is found to be particularly strong in a supersymmetric scenario due to rapid growth of entropy perturbations in the tachyonic region separating the two inflationary stages. Our analysis suggests that realistic double-inflation models will provide a rich and fruitful arena for the application of future cosmic data sets and new approximation schemes which go beyond slow-roll.

We consider the construction of nonsingular Pre-Big-Bang and Ekpyrotic type cosmological models realized by the addition to the action of specific higher-order terms stemming from quantum corrections. We study models involving general relativity coupled to a single scalar field with a potential motivated by the Ekpyrotic scenario. We find that the inclusion of the string loop and quantum correction terms in the string frame makes it possible to obtain solutions of the variational equations which are nonsingular and bouncing in the Einstein frame, even when a negative exponential potential is present, as is the case in the Ekpyrotic scenario. We analyze the spectra of perturbations produced during the bouncing phase and find that the spectrum of curvature fluctuations in the model proposed originally to implement the Ekpyrotic scenario has a large blue tilt ($n_R= 3$). Except for instabilities introduced on small scales, the result agrees with what is obtained by imposing continuity of the induced metric and of the extrinsic curvature across a constant scalar field (up to $k^2$ corrections equal to the constant energy density) matching surface between the contracting and the expanding Einstein Universes. We also discuss nonsingular cosmological solutions obtained when a Gauss-Bonnet term with a coefficient suitably dependent on the scalar matter field is added to the action in the Einstein frame with a potential for the scalar field present.

We discuss the principles governing the selection of inflationary models for which preheating can affect the CMB. This is a (fairly small) subset of those models which have non-negligible entropy/isocurvature perturbations on large scales during inflation. We study new models which belong to this class - two-field inflation with negative nonminimal coupling and hybrid/double/supernatural inflation models where the tachyonic growth of entropy perturbations can lead to the variation of the curvature perturbation, ${\cal R}$, on super-Hubble scales. Finally we present evidence against recent claims for the variation of ${\cal R}$ in the absence of substantial super-Hubble entropy perturbations.

We study nature of singularities in anisotropic string-inspired cosmological models in the presence of a Gauss-Bonnet term. We analyze two string gravity models-- dilaton-driven and modulus-driven cases-- in the Bianchi type-I background without an axion field. In both scenarios singularities can be classified in two ways- the determinant singularity where the main determinant of the system vanishes and the ordinary singularity where at least one of the anisotropic expansion rates of the Universe diverges. In the dilaton case, either of these singularities inevitably appears during the evolution of the system. In the modulus case, nonsingular cosmological solutions exist both in asymptotic past and future with determinant $D=+\infty$ and D=2, respectively. In both scenarios nonsingular trajectories in either future or past typically meet the determinant singularity in past/future when the solutions are singular, apart from the exceptional case where the sign of the time-derivative of dilaton is negative. This implies that the determinant singularity may play a crucial role to lead to singular solutions in an anisotropic background.

We study density perturbations in several cosmological models motivated by string theory. The evolution and the spectra of curvature perturbations ${\cal R}$ are analyzed in the Ekpyrotic scenario and nonsingular string cosmologies. We find that these string-inspired models generally exhibit blue spectra in contrast to standard slow-roll inflationary scenarios. We also clarify the parameter range where ${\cal R}$ is enhanced on superhorizon scales.

We study gravitational creation of light fermions in the presence of classical scalar metric perturbations about a flat Friedmann-Lemaitre- Robertson-Walker (FLRW) background. These perturbations can be large during preheating, breaking the conformal flatness of the background spacetime. We compute numerically the total number of particles generated by the modes of the metric perturbations which have grown sufficiently to become classical. In the absence of inhomogeneities massless fermions are not gravitationally produced, and then this effect may be relevant for abundance estimates of light gravitational relics.

We study cosmological perturbations generated from quantum fluctuations in multi-field inflationary scenarios in generalized Einstein theories, taking both adiabatic and isocurvature modes into account. In the slow-roll approximation, explicit closed-form long-wave solutions for field and metric perturbations are obtained by the analysis in the Einstein frame. Since the evolution of fluctuations depends on specific gravity theories, we make detailed investigations based on analytic and numerical approaches in four generalized Einstein theories: the Jordan-Brans-Dicke (JBD) theory, the Einstein gravity with a non-minimally coupled scalar field, the higher-dimensional Kaluza-Klein theory, and the $R+R^2$ theory with a non-minimally coupled scalar field. We find that solutions obtained in the slow-roll approximation show good agreement with full numerical results except around the end of inflation. Due to the presence of isocurvature perturbations, the gravitational potential $\Phi$ and the curvature perturbations ${\cal R}$ and $\zeta$ do not remain constant on super-horizon scales. In particular, we find that negative non-minimal coupling can lead to strong enhancement of ${\cal R}$ in both the Einstein and higher derivative gravity, in which case it is difficult to unambiguously decompose scalar perturbations into adiabatic and isocurvature modes during the whole stage of inflation.

We study massive spin-1/2 fermion production in nonsingular superstring cosmology, taking into account one-loop quantum corrections to a superstring effective action with dilaton and modulus fields. While no creation occurs in the massless limit, massive fermions can be produced by the existence of a time-dependent frequency. Due to the increase of the Hubble expansion rate during the modulus-driven phase, the occupation of number of fermions continues to grow until the point of the graceful exit, after which fermion creation ceases with the decrease of the Hubble rate.

We study the amplification of large-scale magnetic fields during preheating and inflation in several different models. Preheating can resonantly amplify seed fields on cosmological scales. In the presence of conductivity, however, the effect of resonance is typically weakened while the amplitude of produced magnetic fields depends sensitively on the evolution of the conductivity during the preheating and thermalisation phases. In addition we discuss geometric magnetisation, where amplification of magnetic fields occurs through coupling to curvature invariants. This can be efficient during inflation due to a negative coupling instability. Finally we discuss the breaking of the conformal flatness of the background metric whereby magnetic fields can be stimulated through the growth of scalar metric perturbations during metric preheating.

Preheating after inflation may over-produce primordial black holes (PBH's) in many regions of parameter space. As an example we study two-field models with a massless self-interacting inflaton, taking into account second order field and metric backreaction effects as spatial averages. We find that a complex quilt of parameter regions above the Gaussian PBH over-production threshold emerges due to the enhancement of curvature perturbations on all scales. It should be possible to constrain realistic models of inflation through PBH over-production although many issues, such as rescattering and non-Gaussianity, remain unsolved or unexplored.

We study brane-world preheating in massive chaotic inflationary scenario where scalar fields are confined on the 3-brane. Assuming that quadratic contribution in energy densities dominates the Hubble expansion rate during preheating, the amplitude of inflaton decreases slowly relative to the standard dust-dominated case. This leads to an efficient production of $\chi$ particles via nonperturbative decay of inflaton even if its coupling is of order $g=10^{-5}$. We also discuss massive particle creation heavier than inflaton, which may play important roles for the baryo- and lepto-genesis scenarios.

We study the dynamics and perturbations during inflation and reheating in a multi-field model where a second scalar field $\chi$ is nonminimally coupled to the scalar curvature $(\frac12 \xi R\chi^2$). When $\xi$ is positive, the usual inflationary prediction for large-scale anisotropies is hardly altered while the $\chi$ fluctuation in sub-Hubble modes can be amplified during preheating for large $\xi$. For negative values of $\xi$, however, long-wave modes of the $\chi$ fluctuation exhibit exponential increase during inflation, leading to the strong enhancement of super-Hubble metric perturbations even when $|\xi|$ is less than unity. This is because the effective $\chi$ mass becomes negative during inflation. We constrain the strength of $\xi$ and the initial $\chi$ by the amplitude of produced density perturbations. One way to avoid nonadiabatic growth of super-Hubble curvature perturbations is to stabilize the $\chi$ mass through a coupling to the inflaton. Preheating may thus be necessary in these models to protect the stability of the inflationary phase.

We investigate preheating in a higher-dimensional generalized Kaluza-Klein theory with a quadratic inflaton potential $V(\phi)=\frac12 m^2\phi^2$ including metric perturbations explicitly. The system we consider is the multi-field model where there exists a dilaton field $\sigma$ which corresponds to the scale of compactifications and another scalar field $\chi$ coupled to inflaton with the interaction $\frac12 g^2\phi^2\chi^2+\tilde{g}^2\phi^3\chi$. In the case of $\tilde{g}=0$, we find that the perturbation of dilaton does not undergo parametric amplification while the $\chi$ field fluctuation can be enhanced in the usual manner by parametric resonance. In the presence of the $\tilde{g}^2\phi^3\chi$ coupling, the dilaton fluctuation in sub-Hubble scales is modestly amplified by the growth of metric perturbations for the large coupling $\tilde{g}$. In super-Hubble scales, the enhancement of the dilaton fluctuation as well as metric perturbations is weak, taking into account the backreaction effect of created $\chi$ particles. We argue that not only is it possible to predict the ordinary inflationary spectrum in large scales but extra dimensions can be held static during preheating in our scenario.

Fermion creation during preheating in the presence of multiple scalar fields exhibits a range of interesting behaviour relevant to estimating post-inflation gravitino abundances. We present non-perturbative analysis of this phenomenon over a 6-dimensional parameter space in an expanding background paying particular attention to the interplay between instant and direct fermion preheating. In the broad resonance regime we find that instant fermion production is sensitive to suppression of the long wavelength scalar field modes during inflation. Further, the standard scenario of resonant fermionic preheating through inflaton decay can be significantly modified by instant preheating, and may even lead to a decrease in the number of fermions produced. We explicitly include the effects of metric perturbations and demonstrate that they are important at small coupling but not at strong coupling, due to the rapid saturation of the Pauli bound.

Metric perturbations typically strengthen field resonances during preheating. In contrast we present a model in which the super-Hubble field resonances are completely {\em suppressed} when metric perturbations are included. The model is the nonminimal Fakir-Unruh scenario which is exactly solvable in the long-wavelength limit when metric perturbations are included, but exhibits exponential growth of super-Hubble modes in their absence. This gravitationally enhanced integrability is exceptional, both for its rarity and for the power with which it illustrates the importance of including metric perturbations in consistent studies of preheating. We conjecture a no-go result - there exists no {\em single-field} model with growth of cosmologically-relevant metric perturbations during preheating.

We consider the dynamics of power-law inflation with a nonminimally coupled scalar field $\phi$. It is well known that multiple scalar fields with exponential potentials $V(\phi)=V_0 {\rm exp}(-\sqrt{16\pi/p m_{\rm pl}^2} \phi)$ lead to an inflationary solution even if the each scalar field is not capable to sustain inflation. In this paper, we show that inflation can be assisted even in the one-field case by the effect of nonminimal coupling. When $\xi$ is positive, since an effective potential which arises by a conformal transformation becomes flatter compared with the case of $\xi=0$ for $\phi>0$, we have an inflationary solution even when the universe evolves as non-inflationary in the minimally coupled case. For the negative $\xi$, the assisted inflation can take place when $\phi$ evolves in the region of $\phi<0$ \.

We investigate the particle production of a scalar field $\chi$ coupled to an inflaton field $\phi$ ($g^2\phi^2\chi^2/2$) in the {\it oscillating inflation} model, which was recently proposed by Damour and Mukhanov. Although the fluctuation of the $\phi$ field can be effectively enhanced during a stage of the oscillating inflation, the maximum fluctuation is suppressed as the critical value $\phi_c$ which indicates the scale of the core part of the inflaton potential decreases, in taking into account the back reaction effect of created particles. As for the $\chi$ particle production, we find that larger values of the coupling constant $g$ are required to lead to an efficient parametric resonance with the decrease of $\phi_c$, because an effective mass of inflaton around the minimum of its potential becomes larger. However, it is possible to generate the superheavy $\chi$ particle whose mass is greater than $10^{14}$ GeV, which would result in an important consequence for the GUT baryogenesis.

Recently, Cormier and Holman pointed out that fluctuations of an inflaton field $\phi$ are significantly enhanced in the model of {\it spinodal inflation} with a potential $V(\phi)$ for which the second derivative $V^{(2)}(\phi)$ changes sign. As an application of this model, we investigate particle production in the natural inflation model with a potential $V(\phi)=m^4 [1+\cos (\phi/f)]$ by making use of the Hartree approximation. For typical mass scales $f \sim m_{\rm pl} \sim 10^{19}$GeV, and $m \sim m_{\rm GUT} \sim 10^{16}$GeV, we find that growth of fluctuations relevantly occurs for the initial value of inflaton $\phi(0)<0.1 m_{\rm pl}$. Especially for $\phi(0)<10^{-6}m_{\rm pl}$, maximum fluctuations are so large that secondary inflation takes place by produced fluctuations. In this case, the achieved number of $e$-folding becomes much larger than in the case where an effect of spinodal instability is neglected.

We study parametric amplification of Kaluza-Klein (KK) modes in a higher $D$-dimensional generalized Kaluza-Klein theory, which was originally considered by Mukohyama in the narrow resonance case. It was suggested that KK modes can be enhanced by an oscillation of a scale of compactification by the $d$-dimensional sphere $S^d~(d=D-4)$ and by the direct product $S^{d_1}\times S^{d_2}~(d_1+d_2=D-4)$. We extend this past work to the more general case where initial values of the scale of compactification and the quantum number of the angular momentum $l$ of KK modes are not small. We perform analytic approaches based on the Mathieu equation as well as numerical calculations, and find that the expansion of the universe rapidly makes the KK field deviate from instability bands. As a result, KK modes are not enhanced sufficiently in an expanding universe in these two classes of models.

We investigate preheating of an inflaton field $\phi$ coupled nonminimally to a spacetime curvature. In the case of a self-coupling inflaton potential $V(\phi)=\lambda \phi^4/4$, the dynamics of preheating changes by the effect of the negative $\xi$. We find that the nonminimal coupling works in two ways. First, since the initial value of inflaton field for reheating becomes smaller with the increase of $|\xi|$, the evolution of the inflaton quanta is delayed for fixed $\lambda$. Second, the oscillation of the inflaton field is modified and the nonadiabatic change around $\phi=0$ occurs significantly. That makes the resonant band of the fluctuation field wider. Especially for strong coupling regimes $|\xi| \gg 1$, the growth of the inflaton flutuation is dominated by the resonance due to the nonminimal coupling, which leads to the significant enhancement of low momentum modes. Although the final variance of the inflaton fluctuation does notchange significantly compared with the minimally coupled case, we have found that the energy transfer from the homogeneous inflaton to created particles efficiently occurs for $\xi<-60$.

In higher-curvature inflation models ($R+\alpha_n R^n$), we study a parametric preheating of a scalar field $\chi$ coupled non-minimally to a spacetime curvature $R$ ($\xi R \chi^2$). In the case of $R^2$-inflation model, efficient preheating becomes possible for rather small values of $\xi$, i.e. $|\xi|< several. Although the maximal fluctuation $\sqrt{< \chi^2 >}_{max} \approx 2 \times10^{17}$ GeV for $\xi \approx -4$ is almost the same as the chaotic inflation model with a non-minimally coupled $\chi$ field, the growth rate of the fluctuation becomes much larger and efficient preheating is realized. We also investigate preheating for $R^4$ model and find that the maximal fluctuation is $\sqrt{< \chi^2 >}_{max} \approx 8 \times 10^{16}$ GeV for $\xi \approx -35$.

We investigate a resonant particle production of a scalar field $\chi$ coupled non-minimally to a spacetime curvature $R$ ($\xi R \chi^2$) as well as to an inflaton field $\phi$ ($g^2\phi^2\chi^2$). In the case of $g < 3 \times 10^{-4}$, $\xi$ effect assists $g$-resonance in certain parameter regimes. However, for $g > 3 \times 10^{-4}$, $g$-resonance is not enhanced by $\xi$ effect because of $\xi$ suppression effect as well as a back reaction effect. If $\xi \approx -4$, the maximal fluctuation of produced $\chi$-particle is $\sqrt{< \chi^2 >}_{max} \approx 2 \times 10^{17}$ GeV for $g < 1 \times 10^{-5}$, which is larger than the minimally coupled case with $g \approx 1 \times 10^{-3}$.