In this paper, in the Starobinsky gravity (SG) model, the gravitational wave (GW) luminosity distance dLGW(z) is analytically and numerically investigated by considering the dark matter as the two different kinds of shear viscous fluids in universe. Concretely, on the basis of the propagation equation of the metric perturbation derived by using transverse-traceless gauge for the linear perturbation of the metric in FRW background, we obtain the analytical expressions of the modified friction term δ(z) and the ratio dLGW(z)/dLEM(z) of GW to electromagnetic (EM) wave luminosity distance, and find that they only depend on the SG model parameter α and the shear viscous coefficient η, regardless of the bulk viscous coefficient ζ. Evidently, both of them are affected by α and η. It worth stressing that the results given by us can reduce to ones in documents [13, 22, 30]. It follows that we can extend the previous results as the special cases in our model. Furthermore, by combining the latest observational data, we more tightly constrain the values of related parameters and numerically analyze the influence of both α and η on δ(z) and dLGW(z)/dLEM(z). We observe that, in the range of 0<z≤8.2, the evolutionary curves of δ(z) and dLGW(z)/dLEM(z) exhibit that δ(z) (dLGW(z)/dLEM(z)) monotonically decrease (increase) as redshift rises, moreover there are always δ(z)<0 and dLGW(z)≥dLEM(z), specially dLGW(0)=dLEM(0) in the limit of z→0, which means that detecting GW signals might represent a more efficient tool than detecting EM signals to test modified gravity model on the given scale of cosmic distances. It is worth noting that the numerical analyses to δ(z) and dLGW(z)/dLEM(z) illustrate once again that the GW luminosity distance dLGW(z) is indeed influenced by α and η, in particular more sensitive to α than η in our model. In addition, by comparing the GW and EM luminosity distances of SG model with dL(GR) of Einstein GR, we find that the GW and EM signals in SG1 and SG2 change from stronger to weaker than the ones in GR as redshift increases and more difficultly to be detected in the same distance and detection sensitivity. Meanwhile, the EM signals change from strong to weak earlier than GW signals do. Finally, we confirm that δ(z) and dLGW(z) of our model meet the constraints imposed by the standard siren GW170817 and its electromagnetic companion GRB170817A.
Read full abstract