We investigate gravitational lensing and particle motions around non-asymptotically flat black hole spacetime in non-linear, ghost-free massive gravity theory, called de Rham-Gabadadze-Tolley (dRGT) massive gravity. Deflection angle formulae are derived in terms of perihelion parameter. The deflection angle can be positive, zero or even negative with various perihelion distance. The negative angle reveals repulsive behaviour of gravity from a linear term $\gamma$ in the dRGT black hole solution. We also find an analytically approximated formula of deflection angle in two regimes: large and small $\gamma$ term regimes which are shown to be consistent with direct numerical integration. Null and timelike geodesic motions on equatorial plane are explored. Particle trajectories around the dRGT black hole are plotted and discussed in details.

Recent low-redshift observations give value of the present-time Hubble parameter $H_{0}\simeq 74~\rm{km}/{sec}/\rm{Mpc}$, roughly 10\% higher than the predicted value $H_{0}=67.4~\rm{km}/{sec}/\rm{Mpc}$ from Planck's observations of the CMB and the $\Lambda$CDM model. Phenomenologically, we show that the Friedmann equation requires an extra unknown component $X$ to contribute a negative density to the Universe in order to resolve the Hubble tension without changing the Planck's constraint on the matter and dark energy densities. For the extra negative density to be sufficiently small, its equation-of-state parameter must satisfy $1/3\leq w_{X}\leq1$. We propose a quintom model of two scalar fields that realizes this condition and successfully resolve the Hubble tension. One scalar field acts as a quintessence while another "phantom" scalar conformally couples to matter in such a way that viable cosmological scenario can be achieved. The model depends only on two parameters, $\lambda_{\phi}$ and $\delta$ which represent rolling tendency of the self-interacting potential of the quintessence and the strength of conformal phantom-matter coupling respectively.

We present exact spherical black hole solutions in de Rham, Gabadadze and Tolley (dRGT) massive gravity for a generic choice of the parameters in the theory, and also discuss the thermodynamical and phase structure of the black hole in both the grand canonical and canonical ensembles (for charged case). It turns out that the dGRT black hole solutions includes the known solutions to the Einstein field equations, such as, the monopole-de Sitter-Schwarzschild ones with the coefficients for the third and fourth terms in the potential and the graviton mass in massive gravity naturally generates the cosmological constant and the global monopole term. Furthermore, we compute the mass, temperature, and entropy of dGRT black hole solutions and also perform thermodynamical stability. It turns out that the presence of the graviton mass completely changes the black hole thermodynamics, and it can provide the Hawking-Page phase transition which is also true for the obtained charged black holes. Interestingly, the entropy of a black hole is unaffected and still obeys area law. In particular, our results, in the limit $m_g \rightarrow 0$, reduced exactly to \emph{vis-$\grave{a}$-vis} the general relativity results.

For a large class of mass-varying massive gravity models, the graviton mass cannot provide the late-time cosmic expansion of the universe due to its vanishing at late time. In this work, we propose a new class of mass-varying massive gravity in which the graviton mass varies according to a kinetic term of a k-essence field. By using a more general form of the fiducial metric, we found a solution such that a non-vanishing graviton mass can drive the accelerated expansion of the universe at late time. We also perform dynamical analyses of such model and found that without introducing the k-essence Lagrangian, the graviton mass can be responsible for both dark contents of the universe, namely dark energy that drives the accelerated expansion of the universe and non-relativistic matter that plays the role of dark matter. Moreover, by including the k-essence Lagrangian, we found that it is possible to alleviate the so-called cosmic coincidence problem.