This work studies wave propagation in the most general scalar-tensor theories, particularly focusing on the causal structure realized in these theories and also the shock formation process induced by nonlinear effects. For these studies we use the Horndeski theory and its generalization to the two scalar field case. We show that propagation speeds of gravitational wave and scalar field wave in these theories may differ from the light speed depending on background field configuration, and find that a Killing horizon becomes a boundary of causal domain if the scalar fields share the symmetry of the background spacetime. About the shock formation, we focus on transport of discontinuity in second derivatives of the metric and scalar field in the shift-symmetric Horndeski theory. We find that amplitude of the discontinuity generically diverges within finite time, which corresponds to shock formation. It turns out that the canonical scalar field and the scalar DBI model, among other theories described by the Horndeski theory, are free from such shock formation even when the background geometry and scalar field configuration are nontrivial. We also observe that gravitational wave is protected against shock formation when the background has some symmetries at least. This fact may indicate that the gravitational wave in this theory is more well-behaved compared to the scalar field, which typically suffers from shock formation.

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

The Horndeski theory is known as the most general scalar-tensor theory with second-order field equations. In this paper, we explore the bi-scalar extension of the Horndeski theory. Following Horndeski's approach, we determine all the possible terms appearing in the second-order field equations of the bi-scalar-tensor theory. We compare the field equations with those of the generalized multi-Galileons, and confirm that our theory contains new terms that are not included in the latter theory. We also discuss the construction of the Lagrangian leading to our most general field equations.

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

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

We derive the cosmological matching conditions for the homogeneous and isotropic background and for linear perturbations in Horndeski's most general second-order scalar-tensor theory. In general relativity, the matching is done in such a way that the extrinsic curvature is continuous across the transition hypersurface. This procedure is generalized so as to incorporate the mixing of scalar and gravity kinetic terms in the field equations of Horndeski's theory. Our matching conditions have a wide range of applications including the galilean genesis and the bounce scenarios, in which stable, null energy condition violating solutions play a central role. We demonstrate how our matching conditions are used in the galilean genesis scenario. In doing so, we extend the previous genesis models and provide a unified description of the theory admitting the solution that starts expanding from the Minkowski spacetime.

We propose a multi-field extension of (generalized) G-inflation, based on covariant multi-galileons and their generalization preserving second-order field equations. We compute the quadratic action for cosmological perturbations. By comparing the formulas for cosmological perturbations, it is highlighted that multi-field DBI galileon inflation is not included in the multi-field version of generalized G-inflation. Our result indicates that the generalized covariant multi-galileon theory is not the most general multi-scalar-tensor theory with second-order field equations.

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

In Randall-Sundrum II (RS-II) braneworld model, it has been conjectured according to the AdS/CFT correspondence that brane-localized black hole (BH) larger than the bulk AdS curvature scale $\ell$ cannot be static, and it is dual to a four dimensional BH emitting the Hawking radiation through some quantum fields. In this scenario, the number of the quantum field species is so large that this radiation changes the orbital evolution of a BH binary. We derived the correction to the gravitational waveform phase due to this effect and estimated the upper bounds on $\ell$ by performing Fisher analyses. We found that DECIGO/BBO can put a stronger constraint than the current table-top result by detecting gravitational waves from small mass BH/BH and BH/neutron star (NS) binaries. Furthermore, DECIGO/BBO is expected to detect 10$^5$ BH/NS binaries per year. Taking this advantage, we found that DECIGO/BBO can actually measure $\ell$ down to $\ell=0.33 \mu$m for 5 year observation if we know that binaries are circular a priori. This is about 40 times smaller than the upper bound obtained from the table-top experiment. On the other hand, when we take eccentricities into binary parameters, the detection limit weakens to $\ell=1.5 \mu$m due to strong degeneracies between $\ell$ and eccentricities. We also derived the upper bound on $\ell$ from the expected detection number of extreme mass ratio inspirals (EMRIs) with LISA and BH/NS binaries with DECIGO/BBO, extending the discussion made recently by McWilliams. We found that these less robust constraints are weaker than the ones from phase differences.

We study the fields of spherically symmetric thin shell sources, a.k.a. Dyson spheres, in a {\it fully nonlinear covariant} theory of gravity with the simplest galileon field. We integrate exactly all the field equations once, reducing them to first order nonlinear equations. For the simplest galileon, static solutions come on {\it six} distinct branches. On one, a Dyson sphere surrounds itself with a galileon hair, which far away looks like a hair of any Brans-Dicke field. The hair changes below the Vainshtein scale, where the extra galileon terms dominate the minimal gradients of the field. Their hair looks more like a fuzz, because the galileon terms are suppressed by the derivative of the volume determinant. It shuts off the `hair bunching' over the `angular' 2-sphere. Hence the fuzz remains dilute even close to the source. This is really why the Vainshtein's suppression of the modifications of gravity works close to the source. On the other five branches, the static solutions are all {\it singular} far from the source, and shuttered off from asymptotic infinity. One of them, however, is really the self-accelerating branch, and the singularity is removed by turning on time dependence. We give examples of regulated solutions, where the Dyson sphere explodes outward, and its self-accelerating side is nonsingular. These constructions may open channels for nonperturbative transitions between branches, which need to be addressed further to determine phenomenological viability of multi-branch gravities.