Unlike the electromagnetic radiation from astrophysical objects, gravitational waves (GWs) from binary star mergers have much longer wavelengths and are coherent. For ground-based GW detectors, when the lens object between the source and the earth has mass $\sim 1-10^5M_\odot$, the diffraction effect should be considered since the chirping wavelengths are comparable to the scale of the barrier (its Schwarzschild radius). The waveform will thus be distorted as the fringes. In this work, we show that signals from the third-generation GW detectors like the Einstein Telescope (ET) would be a smoking gun for probing the nature of compact dark matter (CDM) or primordial black holes. Detection of the lensing effects becomes harder when the lens mass is smaller. ET is more sensitive than LIGO, the constraint is available for CDM mass $>5M_\odot$ while LIGO can only detect the mass $>100M_\odot$. For a null search of the fringes, one-year observation of ET can constrain the CDM density fraction to $\sim10^{-2}-10^{-5}$ in the mass range $M_{\rm{CDM}}=10M_\odot-100M_\odot$.

Testing gravity theory in the strong field region becomes a reality due to the observations of gravitational waves and black hole shadows. In this paper, we discuss how to constrain the possible deviations of the classical general relativity with the image of M87* observed by the Event Horizon Telescope. More precisely, we want to know where is the event horizon for a non-rotating black hole. General relativity predicts the horizon is located at the Schwarzschild radius $r_\textrm{s}$, while other gravity theories may give different predictions. We propose a parameterized Schwarzschild metric (PSM) in which the horizon is located at $r=nr_\textrm{s}$, where $n$ is a real free parameter, and prove general relativity with nonlinear electrodynamics allows $n\neq1$. In the weak field region, the PSM is equivalent to the Schwarzschild metric regardless of the value of $n$. In the strong field region, the difference between the PSM and Schwarzschild metric would leave an imprint on the shadow image. We present detailed calculations and discussions on the black hole shadows with large background light source and accretion disk in the PSM framework. More importantly, we point out that $n\approx2$ can be used to explain why the black hole mass measured by the shadow is a factor of about two larger than the previous gas dynamics measurements. If this explanation is confirmed to be right, then this phenomenon, together with the late-time cosmological acceleration, will be very important to test gravity theories.

The Newtonian approximation with a nonvanishing nonlocal background field is analyzed for the scalar-tensor nonlocal gravity and nonlocal Gauss-Bonnet gravity. For these two theories, our calculations show that the Newtonian gravitational constant $G$ is time-varying and $|\dot{G}/G|=\mathcal{O}(H_0)$ for the general case of cosmological background evolution, which is similar to the results of the Deser-Woodard and Maggiore-Mancarella theories. Therefore, observations about the orbit period of binary star (or star-planet) systems could rule out these theories. One thing worth mentioning is that the nonlocal Gauss-Bonnet gravity gives $\Psi=\Phi$ and a constant $G$ in the de Sitter phase. Our results also highlight the uniqueness of the RT model [M. Maggiore, \href{http://dx.doi.org/10.1103/PhysRevD.89.043008}{Phys. Rev. D {\bf 89}, 043008 (2014)}], which is the only nonlocal gravity theory that can successfully describe the gravitational phenomena from solar system to cosmological scales for now.

Einstein field equations are second-order differential equations. In this paper, we propose a new gravity theory with pure fourth-order field equations, which we call 4G for brevity. We discuss the applications of 4G in cosmology, gravitational waves, and local gravitational systems. 4G predicts the scale factor $a\propto t^{4/3}$ for the matter-dominated universe, and $a\propto t$ for the radiation-dominated universe. The former can explain the late-time acceleration without suffering from the coincidence and fine-tuning problems, while the latter can solve the horizon problem. 4G is a massless gravity, which means the speed of gravitational waves equals to the speed of light. Based on the discussions about exact vacuum solution and weak field approximation, we argue that Schwarzschild metric should be the real physical metric to describe solar system in 4G.

Strong gravitational lensing is an important tool to probe the universe. In the theoretical analysis of gravitational lensing, it is assumed that continuous density profile can correctly describe the lens galaxy. But in fact this assumption has never been rigorously tested. In this paper, we discuss this issue, and point out that if we use discrete density profile to model the lens galaxy, then the position of the images does not change, but the magnification will be increased. Strongly lensed gravitational waves could test this conclusion in the future.

The cosmological constant problem has become one of the most important ones in modern cosmology. In this paper, we try to construct a model that can avoid the cosmological constant problem and have the potential to explain the apparent late-time accelerating expansion of the universe in both luminosity distance and angular diameter distance measurement channels. In our model, the core is to modify the relation between cosmological redshift and scale factor for photons. We point out three ways to test our hypothesis: the supernova time dilation; the gravitational waves and its electromagnetic counterparts emitted by the binary neutron star systems; and the Sandage--Loeb effect. All of this method is feasible now or in the near future.

In this paper, by analyzing the dynamics of the inhomogeneous quintessence dark energy, we find that the gradient energy of dark energy will oscillate and gradually vanish, which indicates the gradient energy of the scalar field present in the early universe does not affect the current dynamics of the universe. Moreover, with the decaying of gradient energy, there exists a possible mutual transformation between kinetic energy and gradient energy. In the framework of interacting dark energy models, we argue that inhomogeneous dark energy may have a significant effect on the evolution of the cosmic background, the investigation of which still requires fully relativistic $N$-body numerical simulations in the future.

In this paper we use a recently compiled data set, which comprises 118 galactic-scale strong gravitational lensing (SGL) systems to constrain the statistic property of SGL system, as well as the curvature of universe without assuming any fiducial cosmological model. Based on the singular isothermal ellipsoid (SIE) model of SGL system, we obtain that the constrained curvature parameter $\Omega_{\rm k}$ is close to zero from the SGL data, which is consistent with the latest result of planck measurement. More interestingly, we find that the parameter $f$ in the SIE model is strongly correlated with the curvature $\Omega_{\rm k}$. Neglecting this correlation in the analysis will significantly overestimate the constraining power of SGL data on the curvature. Furthermore, the obtained constraint on $f$ is different from previous results: $f=1.105\pm0.030$ ($68\%$ C.L.), which means that the standard singular isothermal sphere (SIS) model ($f=1$) is disfavored by the current SGL data at more than $3\sigma$ confidence level. We also divide the whole SGL data into two parts according to the centric stellar velocity dispersion $\sigma_{\rm c}$ and find that the larger value of $\sigma_{\rm c}$ the subsample has, the more favored the standard SIS model is. Finally, we extend the SIE model by assuming the power-law density profiles for the total mass density, $\rho=\rho_0(r/r_0)^{-\alpha}$, and luminosity density, $\nu=\nu_0(r/r_0)^{-\delta}$, and obtain the constraints on the power-law indexes: $\alpha=1.95\pm0.04$ and $\delta=2.40\pm0.13$ at 68\% confidence level. When assuming the power-law index $\alpha=\delta=\gamma$, this scenario is totally disfavored by the current SGL data, $\chi^2_{\rm min,\gamma} - \chi^2_{\rm min,SIE} \simeq 53$.