We study the dynamics of a spectator Higgs field which stochastically evolves during inflation onto near-critical trajectories on the edge of a runaway instability. We show that its fluctuations do not produce primordial black holes (PBHs) in sufficient abundance to be the dark matter, nor do they produce significant second-order gravitational waves. First we show that the Higgs produces larger fluctuations on CMB scales than on PBH scales, itself a no-go for a viable PBH scenario. Then we track the superhorizon perturbations nonlinearly through reheating using the delta N formalism to show that they are not converted to large curvature fluctuations. Our conclusions hold regardless of any fine-tuning of the Higgs field for both the Standard Model Higgs and for Higgs potentials modified to prevent unbounded runaway.

We extend the effective field theory of inflation to a general Lagrangian constructed from Arnowitt-Deser-Misner variables that encompasses the most general interactions with up to second derivatives of the scalar field whose background breaks temporal diffeomorphism invariance. Degeneracy conditions, corresponding to 8 distinct types -- only one of which corresponds to known degenerate higher-order scalar-tensor models -- provide necessary conditions for eliminating the Ostrogradsky ghost in a covariant theory at the level of the quadratic action in unitary gauge. Novel implications of the degenerate higher-order system for the Cauchy problem are illustrated with the phase space portrait of an explicit inflationary example: not all field configurations lead to physical solutions for the metric even for positive potentials; solutions are unique for a given configuration only up to a branch choice; solutions on one branch can apparently end at nonsingular points of the metric and their continuation on alternate branches lead to nonsingular bouncing solutions; unitary gauge perturbations can go unstable even when degenerate terms in the Lagrangian are infinitesimal. The attractor solution leads to an inflationary scenario where slow-roll parameters vary and running of the tilt can be large even with no explicit features in the potential far from the end of inflation, requiring the optimized slow-roll approach for predicting observables.

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

We generalize the notion of constant-roll inflation earlier introduced in General Relativity (GR) and $f(R)$ gravity to inflationary models in more general scalar-tensor gravity. A number of novel exact analytic solutions for a FLRW spatially flat cosmological background is found for this case. All forms of the scalar field potential and its coupling to gravity producing the exact de Sitter solution, while the scalar field is varying, are presented. In the particular cases of induced gravity and GR with a non-minimally coupled scalar field, all constant-roll inflationary solutions are found. In the former case they represent power-law inflation, while in the latter case the solution is novel and more complicated. Comparison of scalar perturbations generated during such inflation in induced gravity with observational data shows that the constant-roll parameter should be small, similar to constant-roll inflation in GR and $f(R)$ gravity. Then the solution reduces to the standard slow-roll one with small corrections.

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

We find various exact black hole solutions in the shift-symmetric subclass of the quadratic degenerate higher-order scalar-tensor (DHOST) theories with linearly time dependent scalar field whose kinetic term is constant. The exact solutions are the Schwarzschild and Schwarzschild-(anti-)de Sitter solutions, and the Schwarzschild-type solution with a deficit solid angle, which are accompanied by nontrivial scalar field regular at the black hole event horizon. We derive the conditions for the coupling functions in the DHOST Lagrangian that allow the exact solutions, clarify their compatibility with the degeneracy conditions, and provide general form of coupling functions as well as simple models that satisfy the conditions.

Primordial black holes (PBHs) cannot be produced abundantly enough to be the dark matter in canonical single-field inflation under slow roll. This conclusion is robust to local non-Gaussian correlations between long- and short-wavelength curvature modes, which we show have no effect in slow roll on local primordial black hole abundances. For the prototypical model which evades this no go, ultra-slow roll (USR), these squeezed non-Gaussian correlations have at most an order unity effect on the variance of PBH-producing curvature fluctuations for models that would otherwise fail to form sufficient PBHs. Moreover, the transition out of USR, which is necessary for a successful model, suppresses even this small enhancement unless it causes a large increase in the inflaton kinetic energy in a fraction of an e-fold, which we call a large and fast transition. Along the way we apply the in-in formalism, the delta N formalism, and gauge transformations to compute non-Gaussianities and illuminate different aspects of the physical origin of these results. Local non-Gaussianity in the squeezed limit does not weaken the Gaussian conclusion that PBHs as dark matter in canonical single-field inflation require a complicated and fine-tuned potential shape with an epoch where slow roll is transiently violated.

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

We provide a compact and unified treatment of power spectrum observables for the effective field theory (EFT) of inflation with the complete set of operators that lead to second-order equations of motion in metric perturbations in both space and time derivatives, including Horndeski and GLPV theories. We relate the EFT operators in ADM form to the four additional free functions of time in the scalar and tensor equations. Using the generalized slow roll formalism, we show that each power spectrum can be described by an integral over a single source that is a function of its respective sound horizon. With this correspondence, existing model independent constraints on the source function can be simply reinterpreted in the more general inflationary context. By expanding these sources around an optimized freeze-out epoch, we also provide characterizations of these spectra in terms of five slow-roll hierarchies whose leading order forms are compact and accurate as long as EFT coefficients vary only on timescales greater than an efold. We also clarify the relationship between the unitary gauge observables employed in the EFT and the comoving gauge observables of the post-inflationary universe.

We perform a fully relativistic analysis of even-parity linear perturbations around a static and spherically symmetric solution in the most general scalar-tensor theory with second-order field equations. This paper is a sequel to Kobayashi {\em et al.} (2012), in which the linear perturbation analysis for the odd-parity modes is presented. Expanding the Horndeski action to second order in perturbations and eliminating auxiliary variables, we derive the quadratic action for even-parity perturbations written solely in terms of two dynamical variables. The two perturbations can be interpreted as the gravitational and scalar waves. Correspondingly, we obtain two conditions to evade ghosts and two conditions for the absence of gradient instabilities. Only one in each pair of conditions yields a new stability criterion, as the conditions derived from the stability of the gravitational-wave degree of freedom coincide with those in the odd-parity sector. Similarly, the propagation speed of one of the two modes is the same as that for the odd-parity mode, while the other differs in general from them. Our result is applicable to all the theories of gravitation with an extra single scalar degree of freedom such as the Brans-Dicke theory, $f(R)$ models, and Galileon gravity.

Recent gravitational wave observations of binary black hole mergers and a binary neutron star merger by LIGO and Virgo Collaborations associated with its optical counterpart constrain deviation from General Relativity (GR) both on strong-field regime and cosmological scales with high accuracy, and further strong constraints are expected by near-future observations. Thus, it is important to identify theories of modified gravity that intrinsically possess the same solutions as in GR among a huge number of theories. We clarify the three conditions for theories of modified gravity to allow GR solutions, i.e., solutions with the metric satisfying the Einstein equations in GR and the constant profile of the scalar fields. Our analysis is quite general, as it applies a wide class of single-/multi-field scalar-tensor theories of modified gravity in the presence of matter component, and any spacetime geometry including cosmological background as well as spacetime around black hole and neutron star, for the latter of which these conditions provide a necessary condition for no-hair theorem. The three conditions will be useful for further constraints on modified gravity theories as they classify general theories of modified gravity into three classes, each of which possesses i) unique GR solutions (i.e., no-hair cases), ii) only hairy solutions (except the cases that GR solutions are realized by cancellation between singular coupling functions in the Euler-Lagrange equations), and iii) both GR and hairy solutions, for the last of which one of the two solutions may be selected dynamically.

As the first step to extend our understanding of higher-derivative theories, within the framework of analytic mechanics of point particles, we construct a ghost-free theory involving third-order time derivatives in Lagrangian. While eliminating linear momentum terms in the Hamiltonian is necessary and sufficient to kill the ghosts associated with higher derivatives for Lagrangian with at most second-order derivatives, we find that this is necessary but not sufficient for the Lagrangian with higher than second-order derivatives. We clarify a set of ghost-free conditions under which we show that the Hamiltonian is bounded, and that equations of motion are reducible into a second-order system.

The simple $m^2\phi^2$ potential as an inflationary model is coming under increasing tension with limits on the tensor-to-scalar ratio $r$ and measurements of the scalar spectral index $n_s$. Cubic Galileon interactions in the context of the Horndeski action can potentially reconcile the observables. However, we show that this cannot be achieved with only a constant Galileon mass scale because the interactions turn off too slowly, leading also to gradient instabilities after inflation ends. Allowing for a more rapid transition can reconcile the observables but moderately breaks the slow-roll approximation leading to a relatively large and negative running of the tilt $\alpha_s$ that can be of order $n_s-1$. We show that the observables on CMB and large scale structure scales can be predicted accurately using the optimized slow-roll approach instead of the traditional slow-roll expansion. Upper limits on $|\alpha_s|$ place a lower bound of $r\gtrsim 0.005$ and conversely a given $r$ places a lower bound on $|\alpha_s|$, both of which are potentially observable with next generation CMB and large scale structure surveys.

We consider an inflationary scenario where the rate of inflaton roll defined by $\ddot\phi/H\dot \phi$ remains constant. The rate of roll is small for slow-roll inflation, while a generic rate of roll leads to the interesting case of `constant-roll' inflation. We find a general exact solution for the inflaton potential required for such inflaton behaviour. In this model, due to non-slow evolution of background, the would-be decaying mode of linear scalar (curvature) perturbations may not be neglected. It can even grow for some values of the model parameter, while the other mode always remains constant. However, this always occurs for unstable solutions which are not attractors for the given potential. The most interesting particular cases of constant-roll inflation remaining viable with the most recent observational data are quadratic hilltop inflation (with cutoff) and natural inflation (with an additional negative cosmological constant). In these cases even-order slow-roll parameters approach non-negligible constants while the odd ones are asymptotically vanishing in the quasi-de Sitter regime.

We perform a fully relativistic analysis of odd-type linear perturbations around a static and spherically symmetric solution in the most general scalar-tensor theory with second-order field equations in four-dimensional spacetime. It is shown that, as in the case of general relativity, the quadratic action for the perturbations reduces to the one having only a single dynamical variable, from which concise formulas for no-ghost and no-gradient instability conditions are derived. Our result is applicable to all the theories of gravity with an extra scalar degree of freedom. We demonstrate how the generic formulas can be applied to some particular examples such as the Brans-Dicke theory, $f(R)$ models, and Galileon gravity.

For primordial black holes (PBH) to be the dark matter in single-field inflation, the slow-roll approximation must be violated by at least ${\cal O}(1)$ in order to enhance the curvature power spectrum within the required number of efolds between CMB scales and PBH mass scales. Power spectrum predictions which rely on the inflaton remaining on the slow-roll attractor can fail dramatically leading to qualitatively incorrect conclusions in models like an inflection potential and misestimate the mass scale in a running mass model. We show that an optimized temporal evaluation of the Hubble slow-roll parameters to second order remains a good description for a wide range of PBH formation models where up to a $10^7$ amplification of power occurs in $10$ efolds or more.

The previously introduced class of two-parametric phenomenological inflationary models in General Relativity in which the slow-roll assumption is replaced by the more general, constant-roll condition is generalized to the case of $f(R)$ gravity. A simple constant-roll condition is defined in the original Jordan frame, and exact expressions for a scalaron potential in the Einstein frame, for a function $f(R)$ (in the parametric form) and for inflationary dynamics are obtained. The region of the model parameters permitted by the latest observational constraints on the scalar spectral index and the tensor-to-scalar ratio of primordial metric perturbations generated during inflation is determined.

An invertible field transformation is such that the old field variables correspond one-to-one to the new variables. As such, one may think that two systems that are related by an invertible transformation are physically equivalent. However, if the transformation depends on field derivatives, the equivalence between the two systems is nontrivial due to the appearance of higher derivative terms in the equations of motion. To address this problem, we prove the following theorem on the relation between an invertible transformation and Euler-Lagrange equations: If the field transformation is invertible, then any solution of the original set of Euler-Lagrange equations is mapped to a solution of the new set of Euler-Lagrange equations, and vice versa. We also present applications of the theorem to scalar-tensor theories.

The previously proposed class of phenomenological inflationary models in which the assumption of inflaton slow-roll is replaced by the more general, constant-roll condition is compared with the most recent cosmological observational data, mainly the Planck ones. Models in this two-parametric class which remain viable appear to be close to the slow-roll ones, and their inflaton potentials are close to (but still different from) that of the natural inflation model. Permitted regions for the two model parameters are presented.

We study the development of caustics in shift-symmetric scalar field theories by focusing on simple waves with an $SO(p)$-symmetry in an arbitrary number of space dimensions. We show that the pure Galileon, the DBI-Galileon, and the extreme-relativistic Galileon naturally emerge as the unique set of caustic-free theories, highlighting a link between the caustic-free condition for simple $SO(p)$-waves and the existence of either a global Galilean symmetry or a global (extreme-)relativistic Galilean symmetry.

Regardless of the long history of gauge theories, it is not well recognized under which condition gauge fixing at the action level is legitimate. We address this issue from the Lagrangian point of view, and prove the following theorem on the relation between gauge fixing and Euler-Lagrange equations: In any gauge theory, if a gauge fixing is complete, i.e., the gauge functions are determined uniquely by the gauge conditions, the Euler-Lagrange equations derived from the gauge-fixed action are equivalent to those derived from the original action supplemented with the gauge conditions. Otherwise, it is not appropriate to impose the gauge conditions before deriving Euler-Lagrange equations as it may in general lead to inconsistent results. The criterion to check whether a gauge fixing is complete or not is further investigated. We also provide applications of the theorem to scalar-tensor theories and make comments on recent relevant papers on theories of modified gravity, in which there are confusions on gauge fixing and counting physical degrees of freedom.

In the context of classical mechanics, we study the conditions under which higher-order derivative theories can evade the so-called Ostrogradsky instability. More precisely, we consider general Lagrangians with second order time derivatives, of the form $L(\ddot\phi^a,\dot\phi^a,\phi^a;\dot q^i,q^i)$ with $a = 1,\cdots, n$ and $i = 1,\cdots, m$. For $n=1$, assuming that the $q^i$'s form a nondegenerate subsystem, we confirm that the degeneracy of the kinetic matrix eliminates the Ostrogradsky instability. The degeneracy implies, in the Hamiltonian formulation of the theory, the existence of a primary constraint, which generates a secondary constraint, thus eliminating the Ostrogradsky ghost. For $n>1$, we show that, in addition to the degeneracy of the kinetic matrix, one needs to impose extra conditions to ensure the presence of a sufficient number of secondary constraints that can eliminate all the Ostrogradsky ghosts. When these conditions that ensure the disappearance of the Ostrogradsky instability are satisfied, we show that the Euler-Lagrange equations, which involve a priori higher order derivatives, can be reduced to a second order system.

Self-accelerating backgrounds in massive gravity provide an arena to explore the Cauchy problem for derivatively coupled fields that obey complex constraints which reduce the phase space degrees of freedom. We present here an algorithm based on the Kronecker form of a matrix pencil that finds all hidden constraints, for example those associated with derivatives of the equations of motion, and characteristic curves for any 1+1 dimensional system of linear partial differential equations. With the Regge-Wheeler-Zerilli decomposition of metric perturbations into angular momentum and parity states, this technique applies to fully 3+1 dimensional perturbations of massive gravity around any isotropic self-accelerating background. Five spin modes of the massive graviton propagate once the constraints are imposed: two spin-2 modes with luminal characteristics present in the massless theory as well as two spin-1 modes and one spin-0 mode. Although the new modes all possess the same - typically spacelike - characteristic curves, the spin-1 modes are parabolic while the spin-0 modes are hyperbolic. The joint system, which remains coupled by non-derivative terms, cannot be solved as a simple Cauchy problem from a single non-characteristic surface. We also illustrate the generality of the algorithm with other cases where derivative constraints reduce the number of propagating degrees of freedom or order of the equations.

We show that under a general disformal transformation the linear comoving curvature perturbation is not identically invariant, but is invariant on superhorizon scales for any theory that is disformally related to Horndeski's theory. The difference between disformally related curvature perturbations is found to be given in terms of the comoving density perturbation associated with a single canonical scalar field. In General Relativity it is well-known that this quantity vanishes on superhorizon scales through the Poisson equation that is obtained on combining the Hamiltonian and momentum constraints, and we confirm that a similar result holds for any theory that is disformally related to Horndeski's scalar-tensor theory so long as the invertibility condition for the disformal transformation is satisfied. We also consider the curvature perturbation at full nonlinear order in the unitary gauge, and find that it is invariant under a general disformal transformation if we assume that an attractor regime has been reached. Finally, we also discuss the counting of degrees of freedom in theories disformally related to Horndeski's.

We develop methods to calculate the curvature power spectrum in models where features in the inflaton potential nonlinearly excite modes and generate high frequency features in the spectrum. The first nontrivial effect of excitations generating further excitations arises at third order in deviations from slow roll. If these further excitations are contemporaneous, the series can be resummed, showing the exponential sensitivity of the curvature spectrum to potential features. More generally, this exponential approximation provides a power spectrum template which nonlinearly obeys relations between excitation coefficients and whose parameters may be appropriately adjusted. For a large sharp step in the potential, it greatly improves the analytic power spectrum template and its dependence on potential parameters. For axionic oscillations in the potential, it corrects the mapping between the potential and the amplitude, phase and zero point of the curvature oscillations, which might otherwise cause erroneous inferences in for example the tensor-scalar ratio, formally even when that amplitude is $10^3$ times larger than the slow roll power spectrum. It also estimates when terms that produce double frequency oscillations that are usually omitted when analyzing data should be included. These techniques should allow future studies of high frequency features in the CMB and large scale structure to extend to higher amplitude and/or higher precision.

In models like axion monodromy, temporal features during inflation which are not associated with its ending can produce scalar, and to a lesser extent, tensor power spectra where deviations from scale-free power law spectra can be as large as the deviations from scale invariance itself. Here the standard slow-roll approach breaks down since its parameters evolve on an efolding scale $\Delta N$ much smaller than the efolds to the end of inflation. Using the generalized slow-roll approach, we show that the expansion of observables in a hierarchy of potential or Hubble evolution parameters comes from a Taylor expansion of the features around an evaluation point that can be optimized. Optimization of the leading-order expression provides a sufficiently accurate approximation for current data as long as the power spectrum can be described over the well-observed few efolds by the local tilt and running. Standard second-order approaches, often used in the literature, ironically are worse than leading-order approaches due to inconsistent evaluation of observables. We develop a new optimized next-order approach which predicts observables to $10^{-3}$ even for $\Delta N\sim 1$ where all parameters in the infinite hierarchy are of comparable magnitude. For models with $\Delta N \ll 1$, the generalized slow-roll approach provides integral expressions that are accurate to second order in the deviation from scale invariance. Their evaluation in the monodromy model provides highly accurate explicit relations between the running oscillation amplitude, frequency and phase in the curvature spectrum and parameters of the potential.

We consider the stability of self-accelerating solutions to extended quasidilaton massive gravity in the presence of matter. By making a second or extended fiducial metric dynamical in this model, matter can cause it to evolve from a Lorentzian to Euclidean signature, triggering a ghost instability. We study this possibility with scalar field matter as it can model a wide range of cosmological expansion histories. For the $\Lambda$CDM expansion history, stability considerations substantially limit the available parameter space while for a kinetic energy dominated expansion, no choice of quasidilaton parameters is stable. More generally these results show that there is no mechanism intrinsic to the theory to forbid such pathologies from developing from stable initial conditions and that stability can only be guaranteed for particular choices for the matter configuration.

Self-accelerating solutions in massive gravity provide explicit, calculable examples that exhibit the general interplay between superluminality, the well-posedness of the Cauchy problem, and strong coupling. For three particular classes of vacuum solutions, one of which is new to this work, we construct the conformal diagram for the characteristic surfaces on which isotropic stress-energy perturbations propagate. With one exception, all solutions necessarily possess spacelike characteristics, indicating perturbative superluminality. Foliating the spacetime with these surfaces gives a pathological frame where kinetic terms of the perturbations vanish, confusing the Hamiltonian counting of degrees of freedom. This frame dependence distinguishes the vanishing of kinetic terms from strong coupling of perturbations or an ill-posed Cauchy problem. We give examples where spacelike characteristics do and do not originate from a point where perturbation theory breaks down and where spacelike surfaces do or do not intersect all characteristics in the past light cone of a given observer. The global structure of spacetime also reveals issues that are unique to theories with two metrics: in all three classes of solutions, the Minkowski fiducial space fails to cover the entire de Sitter spacetime allowing worldlines of observers to end in finite proper time at determinant singularities. Characteristics run tangent to these surfaces requiring {\it ad hoc} rules to establish continuity across singularities.

We consider $R^p$ inflation with $p \approx 2$, allowing small deviation from $R^2$ inflation. Using the inflaton potential in the Einstein frame, we construct a consistency relation between the scalar spectral index, the tensor-to-scalar ratio, as well as the running of the scalar spectral index, which will be useful to constrain a deviation from $R^2$ inflation in future observations.

Within the framework of the standard quantum electrodynamics, we compute contribution of vacuum polarization at one-loop order to the power spectrum of the magnetic field on inflationary (de Sitter) background. It is found that the one-loop term exhibits the infrared secular growth that is proportional to the number of $e$-folds. The use of the dynamical renormalization group method, which amounts to partial resummation of higher loop diagrams, shows that the resummed power spectrum is free from the secular growth and the loop effect only changes the power from $(k/a)^4$ at the tree level to ${(k/a)}^{4-\nu}$, where $\nu$ represents the contribution from vacuum polarization. The parameter $\nu$, being proportional to the square of the gauge coupling constant as well as the number of fermion species, is a simple function of a ratio of fermion mass to the Hubble parameter and is positive irrespective of the fermion mass. Thus, the loop effect always enhances the infrared magnetic field strength. We find that $\nu \simeq 5 \times 10^{-3}$ is the possible maximum contribution to $\nu$ from a single fermion. This estimate suggests that either large number of fermion species or large coupling constant is a necessary condition for the loop effect to be responsible for the seed of the cosmic magnetic fields.

Recently, a combined model of the primordial inflation and the present cosmic acceleration has been proposed in the context of f(R) gravity. This model is composed of the late-time acceleration term and an R2 term, which enables the model to avoid high curvature singularity and describe a quasi-de Sitter inflationary phase in the early Universe. An interesting feature of this model is that the reheating dynamics after the inflation is significantly modified, in contrast to the original R2 model, and affects the shape of a gravitational wave background (GWB) spectrum. Here we investigate the production of a GWB during the inflation and reheating eras in the R2-corrected f(R) model and compute a GWB spectrum. We found that interesting region of the model parameters has already been excluded by the cosmological limit on abundance of GWs.

We study the dynamics of the oscillating gauged scalar field in a thermal bath. A Langevin type equation of motion of the scalar field, which contains both dissipation and fluctuation terms, is derived by using the real-time finite temperature effective action approach. The existence of the quantum fluctuation-dissipation relation between the non-local dissipation term and the Gaussian stochastic noise terms is verified. We find the noise variables are anti-correlated at equal-time. The dissipation rate for the each mode is also studied, which turns out to depend on the wavenumber.

Effects invoked by the passage of the cosmic string through the Earth are investigated. The cosmic string induces global oscillations of the Earth whose amplitude and acceleration both linearly depend on the string line density. For the line density maximally allowed by the cosmic observations, the oscillations are perceivable even to human beings and may cause serious damages to the environment. Use of the sophisticated accelerograph enables us to detect the string of a line density down to ten orders of magnitude smaller than the cosmologically relevant value.

It is shown that the tension between recent neutrino oscillation experiments, favoring sterile neutrinos with masses of the order of 1 eV, and cosmological data which impose stringent constraints on neutrino masses from the free streaming suppression of density fluctuations, can be resolved in models of the present accelerated expansion of the Universe based on f(R) gravity.

Ultra slow-roll inflation has recently been used to challenge the non-Gaussianity consistency relation. We show that this inflationary scenario belongs to a one parameter class of models and we study its properties and observational predictions. We demonstrate that the power spectrum remains scale-invariant and that the bi-spectrum is of the local type with fnl=5(3-ns)/4 which, indeed, represents a modification of the consistency relation. However, we also show that the system is unstable and suffers from many physical problems among which is the difficulty to correctly WMAP normalize the model. We conclude that ultra slow-roll inflation remains a very peculiar case, the physical relevance of which is probably not sufficient to call into question the validity of the consistency relation.

We provide a general formalism to calculate the infrared correlators of multiple interacting scalar fields in the de Sitter space by means of the stochastic approach. These scalar fields are treated as test fields and hence our result is applicable to the models such as the curvaton scenario where the fields that yield initially isocurvature modes do not contribute to the cosmic energy density during inflationary expansion. The stochastic formalism combined with the argument of conformal invariance of the correlators reflecting the de Sitter isometries allows us to fix the form and amplitude of the three-point functions completely and partially for the four-point functions in terms of calculable quantities. It turns out that naive scaling argument employed in the previous literature does not necessarily hold and we derive the necessary and sufficient condition for the correlator to obey the naive scaling. We also find that correlation functions can in principle exhibit more complicated structure than argued in the literature.

Within the framework of the recently proposed ghost-free massive gravity, a cosmological constant-type self-accelerating solution has been obtained for Minkowski and de Sitter reference metrics. We ease the assumption on the reference metric and find the self-accelerating solution for the reference metric respecting only isotropy, thus considerably extending the range of known solutions.

Chern-Simons gravitational theories are extensions of general relativity in which the parity is violated due to the Chern-Simons term. We study linear perturbations on the static and spherically symmetric background spacetime both for nondynamical and dynamical Chern-Simons theories. We do not make an assumption that the background Chern-Simons scalar field vanishes, which has been adopted in the literature. By eliminating nondynamical variables using their constraint equations, we derive the reduced second order action from which a set of closed evolution equations containing only dynamical variables are immediately obtained and therefore the number of propagating degrees of freedom as well. It is found that ghost is present both for the nondynamical case and for the dynamical case unless the background Chern-Simons scalar field vanishes. It is also found that if the background scalar field vanishes, ghost degrees of freedom are killed and all the modes propagate at the speed of light.

The reheating dynamics after the inflation induced by $R^2$-corrected $f(R)$ model is considered. To avoid the complexity of solving the fourth order equations, we analyze the inflationary and reheating dynamics in the Einstein frame and its analytical solutions are derived. We also perform numerical calculation including the backreaction from the particle creation and compare the results with the analytical solutions. Based on the results, observational constraints on the model are discussed.

We study linear perturbations around the static and spherically symmetric spacetime for the gravitational theories whose Lagrangian depends on Ricci scalar and the parity violating Chern-Simons term. By an explicit construction, we show that Hamiltonian for the perturbation variables is not bounded from below in general, suggesting that such a background spacetime is unstable against perturbations. This gives a strong limit on a phenomenological gravitational model which violates parity. We also provide a necessary and sufficient condition for the theory to belong to a special class in which no such instability occurs. For such theories, the number of propagating modes for $\ell \ge 2$ is three, one from the odd and the other two from the even. Unlike in the case of $f(R)$ theories, those modes are coupled each other, which can be used as a distinctive feature to test the parity violating theories from observations. All the modes propagate at the speed of light. No-ghost condition and no-tachyon condition are the same as those in $f(R)$ theories. For the dipole perturbations, the odd and the even modes completely decouple. The odd mode gives a slowly-rotating BH solution whose metric is linearized in its angular momentum. We provide an integral expression of such a solution. On the other hand, the even mode propagates at the speed of light. For the monopole perturbation, in addition to a mode which just shifts the mass of the background BH, there is also one even mode that propagates at the speed of light.

It is known that scalar-tensor theory of gravity admits regular crossing of the phantom divide line $w_{\DE}=-1$ for dark energy, and existing viable models of present dark energy for its particular case -- $f(R)$ gravity -- possess one such crossing in the recent past, after the end of the matter dominated stage. It was recently noted that during the future evolution of these models the dark energy equation of state $w_{\DE}$ may oscillate with an arbitrary number of phantom divide crossings. In this paper we prove that the number of crossings can be infinite, present an analytical condition for the existence of this effect and investigate it numerically. With the increase of the present mass of the scalaron (a scalar particle appearing in $f(R)$ gravity) beyond the boundary of the appearance of such oscillations, their amplitude is shown to decrease very fast. As a result, the effect quickly becomes small and its beginning is shifted to the remote future.

We have investigated the evolution of a homogeneous isotropic background of the Universe and inhomogeneous subhorizon matter density perturbations in viable $f(R)$ models of present dark energy and cosmic acceleration analytically and numerically. It is found that viable $f(R)$ models generically exhibit recent crossing of the phantom boundary $w_{\rm DE}=-1$. Furthermore, it is shown that the growth index of perturbations depends both on time and wavenumber. This anomalous growth may explain properties of the observational matter power spectrum from the SDSS data and can also partially counteract the spectrum suppression by massive neutrinos making larger values of the total sum of neutrino rest masses possible.

The effect of massive neutrinos on matter power spectrum is discussed in the context of $f(R)$ gravity. It is shown that the anomalous growth of density fluctuations on small scales due to the scalaron force can be compensated by free streaming of neutrinos. As a result, models which predict observable deviation of the equation-of-state parameter $w_\DE$ from $w_\DE=-1$ can be reconciled with observations of matter clustering if the total neutrino mass is $O(0.5 \eV)$.

Evolution of a background space-time metric and sub-horizon matter density perturbations in the Universe is numerically analyzed in viable $f(R)$ models of present dark energy and cosmic acceleration. It is found that viable models generically exhibit recent crossing of the phantom boundary $w_{\rm DE}=-1$. Furthermore, it is shown that, as a consequence of the anomalous growth of density perturbations during the end of the matter-dominated stage, their growth index evolves non-monotonically with time and may even become negative temporarily.

We present numerical calculation of the evolution of a background space-time metric and sub-horizon matter density perturbations in viable $f(R)$ gravity models of present dark energy and cosmic acceleration. We found that viable models generically exhibit recent crossing of the phantom boundary $w_{DE}=-1$. Moreover, as a consequence of the anomalous growth of density perturbations during the end of the matter-dominated stage, their growth index evolves non-monotonically with time and may even become negative.

For a class of viable cosmological models in $f(R)$ gravity which deviation from the Einstein gravity decreases as a inverse power law of the Ricci scalar $R$ for large $R$, an analytic solution for density perturbations in the matter component during the matter dominated stage is obtained in terms of hypergeometric functions. An analytical expression for the matter transfer function at scales much less than the present Hubble scale is also obtained.