The distributions of the pairwise line-of-sight velocity between galaxies and their host clusters are segregated according to the galaxy's colour and morphology. We investigate the velocity distribution of red-spiral galaxies, which represents a rare population within galaxy clusters. We find that the probability distribution function of the pairwise line-of-sight velocity $v_{\rm{los}}$ between red-spiral galaxies and galaxy clusters has a dip at $v_{\rm{los}} = 0$, which is a very odd feature, at 93\% confidence level. To understand its origin, we construct a model of the phase space distribution of galaxies surrounding galaxy clusters in three-dimensional space by using cosmological $N$-body simulations. We adopt a two component model that consists of the infall component, which corresponds to galaxies that are now falling into galaxy clusters, and the splashback component, which corresponds to galaxies that are on their first (or more) orbit after falling into galaxy clusters. We find that we can reproduce the distribution of the line-of-sight velocity of red-spiral galaxies with the dip with a very simple assumption that red-spiral galaxies reside predominantly in the infall component, regardless of the choice of the functional form of their spatial distribution. Our results constrain the quenching timescale of red-spiral galaxies to a few Gyrs, and the radius where the morphological transformation is effective as $r \sim 0.2 h^{-1} \rm{Mpc}$.

We present multi-wavelength studies of optically-defined merging clusters, based on the Hyper Suprime-Cam Subaru Strategic Program. Luminous red galaxies, tracing cluster mass distributions, enable to identify cluster subhalos at various merging stages, and thus make a homogeneous sample of cluster mergers, which is unbiased with respect to the merger boost of the intracluster medium (ICM). We define, using a peak-finding method, merging clusters with multiple-peaks and single clusters with single-peaks from the CAMIRA cluster catalog. Stacked weak-lensing analysis indicates that our sample of the merging clusters is categorized into major mergers. The average halo concentration for the merging clusters is $\sim70\%$ smaller than that of the single-peak clusters, which agrees well with predictions of numerical simulations. The spatial distribution of subhalos is less centrally concentrated than the mass distribution of the main halo. The fractions of red galaxies in the merging clusters are not higher than those of the single-peak clusters. We find a signature of the merger boost of the ICM from stacked Planck Sunyaev-Zeldovich effect and ROSAT X-ray luminosity, but not in optical richness. The stacked X-ray surface brightness distribution, aligned with the main-subhalo pairs of low redshift and massive clusters, shows that the central gas core is elongated along the merger axis and overall gas distribution is misaligned by $\sim60$ deg. The homogeneous, unbiased sample of cluster mergers and multi-wavelength follow-up studies provide a unique opportunity to make a complete picture of merger physics over the whole process.

Velocity dispersions have been employed as a method to measure masses of clusters. To complement this conventional method, we explore the possibility of constraining cluster masses from the stacked phase space distribution of galaxies at larger radii, where infall velocities are expected to have a sensitivity to cluster masses. First, we construct a two component model of the three-dimensional phase space distribution of haloes surrounding clusters up to 50 $h^{-1}$Mpc from cluster centres based on $N$-body simulations. We find that the three-dimensional phase space distribution shows a clear cluster mass dependence up to the largest scale examined. We then calculate the probability distribution function of pairwise line-of-sight velocities between clusters and haloes by projecting the three-dimensional phase space distribution along the line-of-sight with the effect of the Hubble flow. We find that this projected phase space distribution, which can directly be compared with observations, shows a complex mass dependence due to the interplay between infall velocities and the Hubble flow. Using this model, we estimate the accuracy of dynamical mass measurements from the projected phase space distribution at the transverse distance from cluster centres larger than $2h^{-1}$Mpc. We estimate that, by using $1.5\times 10^5$ spectroscopic galaxies, we can constrain the mean cluster masses with an accuracy of 14.5\% if we fully take account of the systematic error coming from the inaccuracy of our model. This can be improved down to 5.7\% by improving the accuracy of the model.