Gravitational redshift as a relativistic effect in cosmological objects is investigated. Possible signatures of the gravitational redshift in measurements of satellite galaxies in clusters of galaxies, intracluster gas, as well as galaxies associated with voids are investigated by developing simple theoretical models. In the analysis of the gravitational redshift of satellite galaxies, we develop a very simple analytic model for satellite galaxies virialised in halos, which enables us to evaluate the signals depending on the properties of the halo occupation distribution of galaxies. We obtain results consistent with recent previous results, though our results are restricted to the satellite galaxies inside the virial radius. In the analysis of intracluster gas, we develop a simple analytic model including the effect of random motions of gases, which are assumed to generate nonthermal pressure. We demonstrate a possible contribution of the random motions of gases to the gravitational redshift. We also investigate a possible signature of the gravitational redshift in measurements of galaxies associated with voids, for the first time as far as we know, by utilising a simple analytic model. We show that the second-order Hubble term, which appears in the expansion of the scale factor around the centre of a void, may make a significant contribution depending on the way the galaxy samples are analysed.

We obtain a constraint on the parameters of a generalized cubic Galileon gravity model exhibiting the Vainshtein mechanism by using multi-wavelength observations of the Coma Cluster. The generalized cubic Galileon model is characterized by three parameters of the turning scale associated with the Vainshtein mechanism, and the amplitude of modifying a gravitational potential and a lensing potential. X-ray and Sunyaev-Zeldovich (SZ) observations of the intra-cluster medium are sensitive to the gravitational potential, while the weak-lensing (WL) measurement is specified by the lensing potential. A joint fit of a complementary multi-wavelength dataset of X-ray, SZ and WL measurements enables us to simultaneously constrain these three parameters of the generalized cubic Galileon model for the first time. We also find a degeneracy between the cluster mass parameters and the gravitational modification parameters, which is influential in the limit of the weak screening of the fifth force.

The chameleon gravity model postulates the existence of a scalar field that couples with matter to mediate a fifth force. If it exists, this fifth force would influence the hot X-ray emitting gas filling the potential wells of galaxy clusters. However, it would not influence the clusters' weak lensing signal. Therefore, by comparing X-ray and weak lensing profiles, one can place upper limits on the strength of a fifth force. This technique has been attempted before using a single, nearby cluster (Coma, $z=0.02$). Here we apply the technique to the stacked profiles of 58 clusters at higher redshifts ($0.1<z<1.2$), including 12 new to the literature, using X-ray data from the XMM Cluster Survey (XCS) and weak lensing data from the Canada France Hawaii Telescope Lensing Survey (CFHTLenS). Using a multi-parameter MCMC analysis, we constrain the two chameleon gravity parameters ($\beta$ and $\phi_{\infty}$). Our fits are consistent with general relativity, not requiring a fifth force. In the special case of $f(R)$ gravity (where $\beta = \sqrt{1/6}$), we set an upper limit on the background field amplitude today of $|f_{\rm{R0}}| < 6 \times 10^{-5}$ (95% CL). This is one of the strongest constraints to date on $|f_{\rm{R0}}|$ on cosmological scales. We hope to improve this constraint in future by extending the study to hundreds of clusters using data from the Dark Energy Survey.

We study the third order solutions of the cosmological density perturbations in the Horndeski's most general scalar-tensor theory under the condition that the Vainshtein mechanism is at work. In this work, we thoroughly investigate the independence property of the functions describing the nonlinear mode-couplings, which is also useful for models within the general relativity. Then, we find that the solutions of the density contrast and the velocity divergence up to the third order ones are characterized by 6 parameters. Furthermore, the 1-loop order power spectra obtained with the third order solutions are described by 4 parameters. We exemplify the behavior of the 1-loop order power spectra assuming the kinetic gravity braiding model, which demonstrates that the effect of the modified gravity appears more significantly in the power spectrum of the velocity divergence than the density contrast.

We investigate a potential of the higher multipole power spectra of the galaxy distribution in redshift space as a cosmological probe on halo scales. Based on the fact that a halo model explains well the multipole power spectra of the luminous red galaxy (LRG) sample in the Sloan Digital Sky Survey (SDSS), we focus our investigation on the random motions of the satellite LRGs that determine the higher multipole spectra at large wavenumbers. We show that our theoretical model fits the higher multipole spectra at large wave numbers from N-body numerical simulations and we apply these results for testing the gravity theory and the velocity structure of galaxies on the halo scales. In this analysis, we use the multipole spectra P_4(k) and P_6(k) on the small scales of the range of wavenumber 0.3<k/[h{Mpc}^{-1}]<0.6, which is in contrast to the usual method of testing gravity by targeting the linear growth rate on very large scales. We demonstrate that our method could be useful for testing gravity on the halo scales.

We propose a novel method to test the gravitational interactions in the outskirts of galaxy clusters. When gravity is modified, this is typically accompanied by the introduction of an additional scalar degree of freedom, which mediates an attractive fifth force. The presence of an extra gravitational coupling, however, is tightly constrained by local measurements. In chameleon modifications of gravity, local tests can be evaded by employing a screening mechanism that suppresses the fifth force in dense environments. While the chameleon field may be screened in the interior of the cluster, its outer region can still be affected by the extra force, introducing a deviation between the hydrostatic and lensing mass of the cluster. Thus, the chameleon modification can be tested by combining the gas and lensing measurements of the cluster. We demonstrate the operability of our method with the Coma cluster, for which both a lensing measurement and gas observations from the X-ray surface brightness, the X-ray temperature, and the Sunyaev-Zel'dovich effect are available. Using the joint observational data set, we perform a Markov chain Monte Carlo analysis of the parameter space describing the different profiles in both the Newtonian and chameleon scenarios. We report competitive constraints on the chameleon field amplitude and its coupling strength to matter. In the case of f(R) gravity, corresponding to a specific choice of the coupling, we find an upper bound on the background field amplitude of |f_{R0}|<6*10^{-5}, which is currently the tightest constraint on cosmological scales.

We study the bispectrum of the matter density perturbations induced by the large scale structure formation in the most general second-order scalar-tensor theory that may possess the Vainshtein mechanism as a screening mechanism. On the basis of the standard perturbation theory, we derive the bispectrum being expressed by a kernel of the second order of the density perturbations. We find that the kernel at the leading order is characterized by one parameter, which is determined by the solutions of the linear density perturbations, the Hubble parameter and the other function specifying nonlinear interactions. This does not allow for varied behavior in the bispectrum of the matter density perturbations in the most general second-order scalar-tensor theory equipped with the Vainshtein mechanism. We exemplify the typical behavior of the bispectrum in a kinetic gravity braiding model.

We investigate the gas density, temperature, and pressure profiles in a dark matter halo under the influence of the chameleon force. We solve the hydrostatic equilibrium equation for the gas coupled with the chameleon field in an analytic manner, using an approximate solution for the chameleon field equation with the source term, with a generalized Navarro-Frenk-White universal density profile. We find that the gas distribution becomes compact because a larger pressure gradient is necessary due to the additional chameleon force. By confronting the theoretical prediction with the data of the temperature profile of the Hydra A cluster according to Suzaku x-ray observations out to the virial radius, we demonstrate that a useful constraint on a model parameter can be obtained depending on the value of the coupling constant. For example, the upper bound of the background value of chameleon field, $\phi_\infty<10^{-4}\mpl$, is obtained in the case $\beta=1$, where $\beta$ is the coupling constant between the chameleon field and the matter, and $\mpl$ is the Planck mass. However, the error of the present data is not so small that we obtain a useful constraint in the case $\beta=1/\sqrt{6}$, which corresponds to an $f(R)$ model.