We investigate the parameter distributions of the viable generalized scalar-tensor theory with conventional dust matter after GW170817 in a model-independent way. We numerically construct the models by computing the time evolution of a scalar field, which leads to a positive definite second-order Hamiltonian and are consistent with the observed Hubble parameter. We show the model parameter distributions in the degenerate higher-order scalar-tensor (DHOST) theory, and its popular subclasses, e.g., Horndeski and GLPV theories, etc.. We find that 1) the Planck mass run rate, $\alpha_M$, is insensitive to distinguish the theories. 2) the kinetic-braiding parameter, $\alpha_B$, clearly discriminates the models from those of the Horndeski theory, 3) the parameters for the higher-order theories, $\alpha_H$ and $\beta_1$, are relatively smaller in magnitude (by several factors) than $\alpha_M$ and $\alpha_B$, but can still be used for discriminating the theories except for the GLPV theory. Based on the above three facts, we propose a minimal set of parameters that sensibly distinguishes the subclasses of DHOST theories, ($\alpha_M$, $\alpha_B-\alpha_M/2$, $\beta_1$).

The properties of gravitational-wave (GW) propagation are modified in alternative theories of gravity and are crucial observables to test gravity at cosmological distance. The propagation speed has already been measured from GW170817 so precisely and pinned down to the speed of light, while other properties of GW propagation have not constrained tightly yet. In this paper, we investigate the measurement precisions of the amplitude damping rate (equivalently, the time variation of the gravitational coupling for GWs) and graviton mass in the generalized framework of GW propagation with the future detectors such as Voyager, Cosmic Explorer, and Einstein Telescope. As a result, we show that the future GW observation can reach 1% error for the amplitude damping. We also study the time variation of the gravitational couplings in Horndeski theory by performing Monte Carlo-based numerical simulations. From the simulation results, we find that the current accelerating Universe prefers the models with less damping of GWs and that the equivalence principle can be tested at the level of 1% by the future GW observation.

Gravitational waves (GW) are generally affected by modification of a gravity theory during propagation in cosmological distance. We numerically perform a quantitative analysis on Horndeski theory at cosmological scale to constrain the Horndeski theory by GW observations in model-independent way. We formulate a parameterization for a numerical simulation based on the Monte Carlo method and obtain the classification of the models that agrees with cosmic accelerating expansion within observational errors of the Hubble parameter. As a result, we find that a large group of the models in the Horndeski theory that mimic cosmic expansion of ${\Lambda}$CDM model can be excluded from the simultaneous detection of a GW and its electromagnetic transient counterpart. Based on our result and the latest detection of GW170817 and GRB170817A, we conclude that the subclass of Horndeski theory including arbitrary functions $G_4$ and $G_5$ can hardly explain cosmic accelerating expansion without fine-tuning.

Instead of Lorentz invariance, gravitational degrees of freedom may obey Lifshitz scaling at high energies, as it happens in Ho\v{r}ava's proposal for quantum gravity. We study consequences of this proposal for the spectra of primordial perturbations generated at inflation. Breaking of 4D diffeomorphism (Diff) invariance down to the foliation-preserving Diff in Ho\v{r}ava-Lifshitz (HL) gravity leads to appearance of a scalar degree of freedom in the gravity sector, khronon, which describes dynamics of the time foliation. One can naively expect that mixing between inflaton and khronon will jeopardize conservation of adiabatic perturbations at super Hubble scales. This indeed happens in the projectable version of the theory. By contrast, we find that in the non-projectable version of HL gravity, khronon acquires an effective mass which is much larger than the Hubble scale well before the Hubble crossing time and decouples from the adiabatic curvature perturbation $\zeta$ sourced by the inflaton fluctuations. As a result, at super Hubble scales the adiabatic perturbation $\zeta$ behaves as in an effectively single field system and its spectrum is conserved in time. Lifshitz scaling is imprinted in the power spectrum of $\zeta$ through the modified dispersion relation of the inflaton. We point out violation of the consistency relation between the tensor-to-scalar ratio and the spectral tilt of primordial gravitational waves and suggest that it can provide a signal of Lorentz violation in inflationary era.

The Einstein Equivalence Principle~(EEP) can be verified by the measurement of the spectral distortions of the Cosmic Microwave Background (CMB). The existence of energy-dependency in the cosmological redshift effect means the EEP violation. Introducing the energy-dependent Friedmann-Robertson-Walker metric motivated by rainbow gravity, we show that the energy-dependent redshift effect causes the CMB spectral distortions. Assuming the simple energy-dependent form of the metric, we evaluate the distortions. From the COBE/FIRAS bound, we find that the deviation degree from the EEP, which is comparable to the difference of the parameterized-post-Newtonian parameter "gamma" in energy, is less than 10^{-9} at the CMB energy scale.