We propose a class of theories that can limit scalars constructed from the extrinsic curvature. Applied to cosmology, this framework allows us to control not only the Hubble parameter but also anisotropies without the problem of Ostrogradsky ghost, which is in sharp contrast to the case of limiting spacetime curvature scalars. Our theory can be viewed as a generalization of mimetic and cuscuton theories (thus clarifying their relation), which are known to possess a structure that limits only the Hubble parameter on homogeneous and isotropic backgrounds. As an application of our framework, we construct a model where both anisotropies and the Hubble parameter are kept finite at any stage in the evolution of the universe in the diagonal Bianchi type I setup. The universe starts from a constant-anisotropy phase and recovers Einstein gravity at low energies. We also show that the cosmological solution is stable against essentially all modes of perturbation, provided some reasonable conditions are satisfied.

Finding effective theories of modified gravity that can resolve cosmological singularities and avoid other physical pathologies such as ghost and gradient instabilities has turned out to be a rather difficult task. The concept of limiting curvature, where one bounds a finite number of curvature-invariant functions thanks to constraint equations, is a promising avenue in that direction, but its implementation has only led to mixed results. Cuscuton gravity, which can be defined as a special subclass of $k$-essence theory for instance, is a minimal modification of gravity since it does not introduce any new degrees of freedom on a cosmological background. Importantly, it naturally incorporates the idea of limiting curvature. Accordingly, models of cuscuton gravity are shown to possess non-singular cosmological solutions and those appear stable at first sight. Yet, various subtleties arise in the perturbations such as apparent divergences, e.g., when the Hubble parameter crosses zero. We revisit the cosmological perturbations in various gauges and demonstrate that the stability results are robust even at those crossing points, although certain gauges are better suited to analyze the perturbations. In particular, the spatially-flat gauge is found to be ill defined when $H=0$. Otherwise, the sound speed is confirmed to be generally close to unity in the ultraviolet, and curvature perturbations are shown to remain essentially constant in the infrared throughout a bounce phase. Perturbations for a model of extended cuscuton (as a subclass of Horndeski theory) are also studied and similar conclusions are recovered.

We find that the recently-proposed ghost-free interaction of a 2-form gauge field in four dimensions, which contains derivative couplings in a nonperturbative manner, can be regarded as a resummation of ghostly interaction terms. We investigate the higher derivative structure of this model in a minisuperspace description and demonstrate that the higher derivative terms can be removed by taking appropriate combinations of the Euler-Lagrange equations, while a truncation at a finite order spoils this structure. We also show that this nature is peculiar to four dimensions.

We show that a simple sub-class of Horndeski theory can describe a time crystal Universe. The time crystal Universe can be regarded as a baby Universe nucleated from a flat space, which is mediated by an extension of Giddings-Strominger instanton in a 2-form theory dual to the Horndeski theory. Remarkably, when a cosmological constant is included, de Sitter Universe can be created by tunneling from the time crystal Universe. It gives rise to a past completion of an inflationary Universe.

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

We propose a new class of higher derivative scalar-tensor theories without the Ostrogradsky's ghost instabilities. The construction of our theory is originally motivated by a scalar field with spacelike gradient, which enables us to fix a gauge in which the scalar field appears to be non-dynamical. We dub such a gauge as the spatial gauge. Though the scalar field loses its dynamics, the spatial gauge fixing breaks the time diffeomorphism invariance and thus excites a scalar mode in the gravity sector. We generalize this idea and construct a general class of scalar-tensor theories through a non-dynamical scalar field, which preserves only spatial covariance. We perform a Hamiltonian analysis and confirm that there are at most three (two tensors and one scalar) dynamical degrees of freedom, which ensures the absence of a degree of freedom due to higher derivatives. Our construction opens a new branch of scalar-tensor theories with higher derivatives.

Extendibility of inflationary spacetimes with flat spatial geometry is investigated. We find that the past boundary of an inflationary spacetime becomes a so-called parallely propagated curvature singularity if the ratio $\dot{H}/a^2$ diverges at the boundary, where $\dot{H}$ and $a$ represent the time derivative of the Hubble parameter and the scale factor, respectively. On the other hand, if the ratio $\dot{H}/a^2$ converges, then the past boundary is regular and continuously extendible. We also develop a method to judge the continuous ($C^0$)extendibility of spacetime in the case of slow-roll inflation driven by a canonical scalar field. As applications of this method, we find that Starobinsky inflation has a $C^0$ parallely propagated curvature singularity, but a small field inflation model with a Higgs-like potential does not. We also find that an inflationary solution in a modified gravity theory with limited curvature invariants is free of such a singularity and is smoothly extendible.

We study cosmological applications of extended vector-tensor theories, whose Lagrangians contain up to two derivatives with respect to metric and vector field. We derive background equations under the assumption of homogeneous and isotropic universe and study the nature of cosmological perturbations on top of the background. We found an example of healthy cosmological solutions in non-trivial degenerate theory of gravity for the first time, where those perturbations do not suffer from any instabilities, that is, ghost and gradient instabilities.

We show that the parametric resonance of gravitational waves occurs due to the axion coherent oscillation and the circular polarization of gravitational waves is induced by the Chern-Simons coupling. However, we have never observed these signatures in the data of gravitational waves. Using this fact, we give stringent constraints on the Chern-Simons coupling constant $\ell$ and the abundance of the light string axion. In particular, we improved the current bound $\ell\leq10^{8}\,\,{\rm km}$ by many orders of magnitude.

We investigate inflation models in Jordan frame supergravity, in which an inflaton non-minimally couples to the scalar curvature. By imposing the condition that an inflaton would have the canonical kinetic term in the Jordan frame, we construct inflation models with asymptotically flat potential through pole inflation technique and discuss their relation to the models based on Einstein frame supergravity. We also show that the model proposed by Ferrara et al. has special position and the relation between the K\"ahler potential and the frame function is uniquely determined by requiring that scalars take the canonical kinetic terms in the Jordan frame and that a frame function consists only of a holomorphic term (and its anti-holomorphic counterpart) for symmetry breaking terms. Our case corresponds to relaxing the latter condition.

We investigate static, spherically symmetric solutions in gravitational theories which have limited curvature invariants, aiming to remove the singularity in the Schwarzschild space-time. We find that if we only limit the Gauss-Bonnet term and the Ricci scalar, then the singularity at the origin persists. Moreover we find that the event horizon can develop a curvature singularity. We also investigate a new class of theories in which all components of the Riemann tensor are bounded. We find that the divergence of the quadratic curvature invariants at the event horizon is avoidable in this theory. However, other kinds of singularities due to the dynamics of additional degrees of freedom cannot be removed, and the space-time remains singular.

It is widely believed that axions are ubiquitous in string theory and could be the dark matter. The peculiar features of the axion dark matter are coherent oscillations and a coupling to the electromagnetic field through the Chern-Simons term. In this paper, we study consequences of these two features of the axion with the mass in a range from $10^{-13}\,{\rm eV}$ to $10^{3}\,{\rm eV}$. First, we study the parametric resonance of electromagnetic waves induced by the coherent oscillation of the axion. As a result of the resonance, the amplitude of the electromagnetic waves is enhanced and the circularly polarized monochromatic waves will be generated. Second, we study the velocity of light in the background of the axion dark matter. In the presence of the Chern-Simons term, the dispersion relation is modified and the speed of light will oscillate in time. It turns out that the change of speed of light would be difficult to observe. We argue that the future radio wave observations of the resonance can give rise to a stronger constraint on the coupling constant and/or the density of the axion dark matter.

Memory effects of gravitational waves from astronomical events or primordial universe might have the information of new physics. It is intriguing to observe that the memory effect exists in electrodynamics as a net momentum kick, while the memory effect in gravity appears as a net relatively displacement. In particular, Winicour has shown that the B-mode memory, which characterizes parity odd global distribution of memory, does not exist. We study the memory effect in axion electrodynamics and find that the B-mode memory effect can exist provided the existence of coherently oscillating axion background field. Moreover, we examine the detectability of the axion dark matter using this effect. We also argue the existence of the B-mode gravitational memory effect in the presence of the axion dark matter.

We revisit nonsingular cosmologies in which the limiting curvature hypothesis is realized. We study the cosmological perturbations of the theory and determine the general criteria for stability. For the simplest model, we find generic Ostrogradski instabilities unless the action contains the Weyl tensor squared with the appropriate coefficient. When considering two specific nonsingular cosmological scenarios(one inflationary and one genesis model), we find ghost and gradient instabilities throughout most of the cosmic evolution. Furthermore, we show that the theory is equivalent to a theory of gravity where the action is a general function of the Ricci and Gauss-Bonnet scalars, and this type of theory is known to suffer from instabilities in anisotropic backgrounds. This leads us to construct a new type of curvature invariant scalar function. We show that it does not have Ostrogradski instabilities, and it avoids ghost and gradient instabilities for most of the interesting background inflationary and genesis trajectories. We further show that it does not possess additional new degrees of freedom in an anisotropic spacetime. This opens the door for studying stable alternative nonsingular very early universe cosmologies.

Recently, several extensions of massive vector theory in curved space-time have been proposed in many literatures. In this paper, we consider the most general vector-tensor theories that contain up to two derivatives with respect to metric and vector field. By imposing a degeneracy condition of the Lagrangian in the context of ADM decomposition of space-time to eliminate an unwanted mode, we construct a new class of massive vector theories where five degrees of freedom can propagate, corresponding to three for massive vector modes and two for massless tensor modes. We find that the generalized Proca and the beyond generalized Proca theories up to the quartic Lagrangian, which should be included in this formulation, are degenerate theories even in curved space-time. Finally, introducing new metric and vector field transformations, we investigate the properties of thus obtained theories under such transformations.

We study quasinormal modes of black holes in Lovelock gravity. We formulate the WKB method adapted to Lovelock gravity for the calculation of quasinormal frequencies (QNFs). As a demonstration, we calculate various QNFs of Lovelock black holes in seven and eight dimensions. We find that the QNFs show remarkable features depending on the coefficients of the Lovelock terms, the species of perturbations, and spacetime dimensions. In the case of the scalar field, when we increase the coefficient of the third order Lovelock term, the real part of QNFs increases, but the decay rate becomes small irrespective of the mass of the black hole. For small black holes, the decay rate ceases to depend on the Gauss-Bonnet term. In the case of tensor type perturbations of the metric field, the tendency of the real part of QNFs is opposite to that of the scalar field. The QNFs of vector type perturbations of the metric show no particular behavior. The behavior of QNFs of the scalar type perturbations of the metric field is similar to the vector type. However, available data are rather sparse, which indicates that the WKB method is not applicable to many models for this sector.

The covariant St\"{u}ckelberg formalism is applied to bi-gravity in order to revisit the issue of absence of the Boulware-Deser (BD) ghost. We first confirm that the leading order action in the decoupling limit for helicity-2 modes of metrics and helicity-0 mode of St\"{u}ckelberg perturbations does not lead to higher time derivative in equations of motion, which suggests the absence of the BD ghost. By extending this method, we reconfirm that the BD ghost does not appear for arbitrary order of the perturbations at the decoupling limit in bi-gravity.

We investigate perturbations of a class of spherically symmetric solutions in massive gravity and bi-gravity. The background equations of motion for the particular class of solutions we are interested in reduce to a set of the Einstein equations with a cosmological constant. Thus, the solutions in this class include all the spherically symmetric solutions in general relativity, such as the Friedmann-Lema\^{i}tre-Robertson-Walker solution and the Schwarzschild (-de Sitter) solution, though the one-parameter family of two parameters of the theory admits such a class of solutions. We find that the equations of motion for the perturbations of this class of solutions also reduce to the perturbed Einstein equations at first and second order. Therefore, the stability of the solutions coincides with that of the corresponding solutions in general relativity. In particular, these solutions do not suffer from non-linear instabilities which often appear in the other cosmological solutions in massive gravity and bi-gravity.

The St\"uckelberg analysis of nonlinear massive gravity in the presence of a general fiducial metric is investigated. We develop a "covariant" formalism for the St\"uckelberg expansion by working with a local inertial frame, through which helicity modes can be characterized correctly. Within this covariant approach, an extended $\Lambda_3$ decoupling limit analysis can be consistently performed, which keeps $\bar{R}_{\mu\nu\rho\sigma}/m^2$ fixed with $\bar{R}_{\mu\nu\rho\sigma}$ the Riemann tensor of the fiducial metric. In this extended decoupling limit, the scalar mode $\pi$ acquires self-interactions due to the presence of the curvature of the fiducial metric. However, the equation of motion for $\pi$ remains of second order in derivatives, which extends the understanding of the absence of the Boulware Deser ghost in the case of a flat fiducial metric.

We find new, simple cosmological solutions with flat, open, and closed spatial geometries, contrary to the previous wisdom that only the open model is allowed. The metric and the St\"{u}ckelberg fields are given explicitly, showing nontrivial configurations of the St\"{u}ckelberg in the usual Friedmann-Lema\^{i}tre-Robertson-Walker coordinates. The solutions exhibit self-acceleration, while being free from ghost instabilities. Our solutions can accommodate inhomogeneous dust collapse represented by the Lema\^{i}tre-Tolman-Bondi metric as well. Thus, our results can be used not only to describe homogeneous and isotropic cosmology but also to study gravitational collapse in massive gravity.