We present a detailed study of the static spherically symmetric solutions in de Rham-Gabadadze-Tolley (dRGT) theory. Since the diffeomorphism invariance can be restored by introducing the St\"{u}ckelberg fields $\phi^a$, there is new invariant $I^{ab}=g^{\mu\nu}\partial_{\mu}\phi^a\partial_\nu\phi^b$ in the massive gravity, which adds to the ones usually encountered in general relativity (GR). In the unitary gauge $\phi^a=x^\mu\delta_\mu^a$, any inverse metric $g^{\mu\nu}$ that has divergence including the coordinate singularity in GR would exhibit a singularity in the invariant $I^{ab}$. Therefore, there is no conventional Schwarzschild metric if we choose unitary gauge. In this paper, we obtain a self-consistent static spherically symmetric ansatz in the nonunitary gauge. Under this ansatz, we find that there are seven solutions including the Schwarzschild solution, Reissner-Nordstr\"{o}m solution and five other solutions. These solutions may possess an event horizon depending upon the physical parameters (Schwarzschild radius $r_s$, scalar charge $S$ and/or electric charge $Q$). If these solutions possess an event horizon, we show that the singularity of $I^{ab}$ is absent at the horizon. Therefore, these solutions may become candidates for black holes in dRGT.

In this paper, we shall study the dynamical behavior of the universe accelerated by the so called Veneziano ghost dark energy component locally and globally by using the linearization and nullcline method developed in this paper. The energy density is generalized to be proportional to the Hawking temperature defined on the trapping horizon instead of Hubble horizon of the Friedmann-Robertson-Walker (FRW) universe. We also give a prediction of the fate of the universe and present the bifurcation phenomenon of the dynamical system of the universe. It seems that the universe could be dominated by dark energy at present in some region of the parameter space.

Depending on the choice of its potential, the scalar phantom field $\phi$ (the equation of state parameter $w<-1$) leads to various catastrophic fates of the universe including big rip, little rip and other future singularity. For example, big rip results from the evolution of the phantom field with an exponential potential and little rip stems from a quadratic potential in general relativity (GR). By choosing the same potential as in GR, we suggest a new mechanism to avoid these unexpected fates (big and little rip) in the inverse-\textit{R} gravity. As a pedagogical illustration, we give an exact solution where phantom field leads to a power-law evolution of the scale factor in an exponential type potential. We also find the sufficient condition for a universe in which the equation of state parameter crosses $w=-1$ divide. The phantom field with different potentials, including quadratic, cubic, quantic, exponential and logarithmic potentials are studied via numerical calculation in the inverse-\textit{R} gravity with $R^{2}$ correction. The singularity is avoidable under all these potentials. Hence, we conclude that the avoidance of big or little rip is hardly dependent on special potential.

A new cosmological model based on the de Sitter gravity is investigated by dynamical analysis and numerical discussions. Via some transformations, the evolution equations of this model can form an autonomous system with 8 physical critical points. Among these critical points there exist one positive attractor and one negative attractor. The positive attractor describes the asymptotic behavior of late-time universe, which indicates that the universe will enter the exponential expansion phase, finally. Some numerical calculations are also carried out, which convince us of this conclusion derived from the dynamical analysis.

We present new static spherically-symmetric solutions of Einstein equations with the quintessence-like matter surrounding a global monopole. These new solutions of the coupling scalar-Einstein equations are more complicated, which depend on the parameter of equation of state $-1 < w_{q} <-{1/3}$. A gravitating global monopole produces a gravitational field of de Sitter kind outside the core in addition to a solid angular deficit. In the $w_{q} = -{1/3}$ case, we have proved that the solution cannot exist since the density of quintessence-like tends to zero if $w_{q} \to -{1/3}$. As a new feature, these monopoles have the outer horizon depending on both Goldstone field and quintessence-like. Since current observations constrain $-1.14 < w_{q} < -0.93$, new global monopoles have interesting astrophysical applications.