We present a simple fitting formula for the scale-dependent growth rate of Hu-Sawicki model in $f(R)$ modified gravity. We compare the accuracy of the fitting function against numerical results and report achieving a sub-percent maximum error for all different studied cases. The validity of the fitting function is tested against a wide range of scales, $10^{-4} \le k \le 5 \ [\mathrm{Mpc}^{-1}]$ , for various redshift values in the range of $z \in [0.0,\ 3.0]$ and different $|f_R^0|$ values. This formula is useful in producing predictions for cosmological probes that are sensitive to growth rate such as redshift space distortions, galaxy-galaxy correlation function, peculiar velocity statistics and pairwise velocities of galaxy clusters.

We perform an extensive analysis of optical counterparts of Planck PSZ2 clusters, considering matches with three recent catalogs built from SDSS data: AMF DR9, redMaPPer (v6.3) and Wen et al (WHL). We significantly extend the number of optical counterparts of detected Planck clusters, and characterize the optical properties when multiple identifications in different catalogs exist. For Planck clusters which already possess an external validation, we analyze the redshift assignment for both optical and X--ray determinations. We then analyze the Planck Cosmology sample and comment on redshift determination and potential mass mis-determinations due to alignment issues. Finally, we inspect the reconstructed $y$ map from Planck and reason on the detectability of optical clusters. Overall, AMF DR9 main (extended) finds 485 (511) optical matches, with 45 (55) previously unmatched PSZ2 clusters, to be compared with the 374 optical matches already present in PSZ2. 29 of the 55 previously unmatched clusters do not yet have a followup in the literature. 18 of these are found in more than one SDSS catalog with consistent redshifts. We provide redshift and mass estimates for the newly matched clusters, and discuss the comparison with the follow-ups, when present. We find good agreement between the redMaPPer and AMF DR9 redshift determinations. From the Planck Cosmology sample, we find 14 clusters which merit further investigation, and discuss possible alignment issues for 9 of these clusters. After inspecting the $y$ map, we provide a list of 229 optical clusters not included in the Planck PSZ2 catalog but showing a prominent $y$ signal. We have further investigated the 86 clusters with Planck S/N $>4.5$ using the MMF technique (applied to the Planck HFI maps), and were able to detect 20 new cluster candidates that are not included in the PSZ2 catalog.

We introduce a new estimator for the mean pairwise velocities of galaxy clusters, which is based on the measurement of the clusters' $\textit{transverse}$ velocity components. The Rees-Sciama (RS) effect offers an opportunity to measure transverse peculiar velocities through its distinct dipolar signature around the halo centers in the Cosmic Microwave Background (CMB) temperature map. We exploit this dipolar structure to extract the magnitude and direction of the transverse velocity vectors from CMB maps simulated with the expected characteristics of future surveys like CMB-S4. Although in the presence of lensed CMB and instrumental noise individual velocities are not reliably reconstructed, we demonstrate that the mean pairwise velocity measurement obtained using the estimator yields a signal-to-noise ratio of $5.2$ for $\sim21,000$ halos with $M > 7\times10^{13}\rm M_\odot$ in a $40\times40$ [deg$^2$] patch at $z=0.5$. While the proposed estimator carries promising prospects for measuring pairwise velocities through the RS effect in CMB stage IV experiments, its applications extend to any other potential probe of transverse velocities.

In this work we study the simultaneous effect of primordial non-Gaussianity and the modification of the gravity in $f(R)$ framework on large scale structure observations. We show that non-Gaussianity and modified gravity introduce a scale dependent bias and growth rate functions. The deviation from $\Lambda$CDM in the case of primordial non-Gaussian models is in large scales, while the growth rate deviates from $\Lambda$CDM in small scales for modified gravity theories. We show that the redshift space distortion can be used to distinguish positive and negative $f_{NL}$ in standard background, while in $f(R)$ theories they are not easily distinguishable. The galaxy power spectrum is generally enhanced in presence of non-Gaussianity and modified gravity. We also obtain the scale dependence of this enhancement. Finally we define galaxy growth rate and galaxy growth rate bias as new observational parameters to constrain cosmology.