We find vacuum solutions such that massive gravitons are confined in a local spacetime region by their gravitational energy in asymptotically flat spacetimes in the context of the bigravity theory. We call such self-gravitating objects massive graviton geons. The basic equations can be reduced to the Schr\"odinger-Poisson equations with the tensor "wavefunction" in the Newtonian limit. We obtain a non-spherically symmetric solution with $j=2,\ell=0$ as well as a spherically symmetric solution with $j=0,\ell=2$ in this system where $j$ is the total angular momentum quantum number and $\ell$ is the orbital angular momentum quantum number, respectively. The energy eigenvalue of the Schr\"odinger equation in the non-spherical solution is smaller than that in the spherical solution. We then study the perturbative stability of the spherical solution and find that there is an unstable mode in the quadrupole mode perturbations which may be interpreted as the transition mode to the non-spherical solution. The results suggest that the non-spherically symmetric solution is the ground state of the massive graviton geon. The massive graviton geons may decay in time due to emissions of gravitational waves but this timescale can be quite long when the massive gravitons are non-relativistic and then the geons can be long-lived. We also argue possible prospects of the massive graviton geons: applications to the ultralight dark matter scenario, nonlinear (in)stability of the Minkowski spacetime, and a quantum transition of the spacetime.

We study a vacuum Bianchi IX universe in the context of Ho\v{r}ava-Lifshitz (HL) gravity. In particular, we focus on the classical dynamics of the universe and analyze how anisotropy changes the history of the universe. For small anisotropy, we find an oscillating universe as well as a bounce universe just as the case of the Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetime. However, if the initial anisotropy is large, we find the universe which ends up with a big crunch after oscillations if a cosmological constant $\Lambda$ is zero or negative. For $\Lambda>0$, we find a variety of histories of the universe, that is de Sitter expanding universe after oscillations in addition to the oscillating solution and the previous big crunch solution. This fate of the universe shows sensitive dependence of initial conditions, which is one of the typical properties of a chaotic system. If the initial anisotropy is near the upper bound, we find the universe starting from a big bang and ending up with a big crunch for $\Lambda\leq 0$, while de Sitter expanding universe starting from a big bang for $\Lambda>0$.

We study the dynamics of isotropic and homogeneous universes in the generalized Ho\v{r}ava-Lifshitz gravity, and classify all possible evolutions of vacuum spacetime. In the case without the detailed balance condition, we find a variety of phase structures of vacuum spacetimes depending on the coupling constants as well as the spatial curvature $K$ and a cosmological constant $\Lambda$. A bounce universe solution is obtained for $\Lambda> 0, K=\pm 1$ or $\Lambda= 0, K=- 1$, while an oscillation spacetime is found for $\Lambda\geq 0, K=1$, or $\Lambda< 0, K=\pm 1$. We also propose a quantum tunneling scenario from an oscillating spacetime to an inflationary universe, resulting in a macroscopic cyclic universe.