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 investigate effects of massive graviton on the rotation curves of the Milky Way, spiral galaxies and Low Surface Brightness~(LSB) galaxies. Using a simple de Rham, Gabadadze, and Tolley (dRGT) massive gravity model, we find static spherically symmetric metric and a modified Tolman-Oppenheimer-Volkoff (TOV) equation. The dRGT nonlinear graviton interactions generate density and pressures which behave like a dark energy that can mimic the gravitational effects of a dark matter halo. We found that rotation curves of most galaxies can be fitted well by a single constant-gravity parameter $\gamma \sim m_{g}^{2}C \sim 10^{-28}~{\rm m^{-1}}$ corresponding to the graviton mass in the range $m_g \sim 10^{-21}-10^{-30} {\rm eV}$ depending on the choice of the fiducial metric parameter $C\sim 1-10^{18}~\text{m}$. Fitting rotation curve of the Milky Way puts strong constraint on the Yukawa-type coupling of the massive graviton exchange as a result of the shell effects.

We study the cosmological dynamics of D-BIonic and DBI scalar field, which is coupled to matter fluid. For the exponential potential and the exponential couplings, we find a new analytic scaling solution yielding the accelerated expansion of the Universe. Since it is shown to be an attractor for some range of the coupling parameters, the density parameter of matter fluid can be the observed value, as in the coupled quintessence with a canonical scalar field. Contrary to the usual coupled quintessence, where the value of matter couple giving observed density parameter is too large to satisfy observational constraint from CMB, we show that the D-BIonic theory can give similar solution with much smaller value of matter coupling. As a result, together with the fact that the D-BIonic theory has a screening mechanism, the D-BIonic theory can solve the so-called coincidence problem as well as the dark energy problem.

We investigate the effects of chameleon scalar field to the effective density and pressure of a dark matter halo. The pressure is generated from the chameleonic fifth force on the matter. We demonstrate that the thick-shell non-singular boundary condition which forbids singular point leads to extremely stringent constraint on the matter-chameleon coupling when applied to galaxy. We argue that chameleon profile with central singularity is more likely to develop in general physical situation. The chameleonic fifth force from the chameleon profile with central singularity experienced by the dark matter could significantly modify the rotation curve of galaxies. The chameleonic fifth force could generate steeper cusp to the rotation curves in any dark matter profiles starting from the Navarro-Frenk-White (NFW) to the pseudo-isothermal (ISO) profile. Upper limits on the coupling constant between the chameleon and the dark matter are estimated from observational data of the late-type Low-Surface-Brightness galaxies (LSB). It is in the order of $\beta < 10^{-3}$.