Issue |
A&A
Volume 695, March 2025
|
|
---|---|---|
Article Number | A46 | |
Number of page(s) | 9 | |
Section | Cosmology (including clusters of galaxies) | |
DOI | https://doi.org/10.1051/0004-6361/202451722 | |
Published online | 03 March 2025 |
Peculiar velocities of galaxy clusters in a IllustrisTNG simulation based on mock observations of Sunyaev-Zel’dovich effects
1
School of Astronomy and Space Science, Nanjing University, Nanjing 210093, People’s Republic of China
2
Key Laboratory of Modern Astronomy and Astrophysics (Nanjing University), Ministry of Education, Nanjing 210093, People’s Republic of China
3
Shanghai Key Lab for Astrophysics, Shanghai Normal University, Shanghai 200234, China
⋆ Corresponding author; yong@nju.edu.cn
Received:
30
July
2024
Accepted:
24
January
2025
Galaxy clusters are the largest self-gravitational systems in the Universe. They are valuable probes of the structure growth of the Universe when we estimate their peculiar velocities with the kinematic Sunyaev-Zel’dovich (kSZ) effects. We investigate whether there is a systematic offset between the peculiar velocities (vz, kSZ) estimated with the kSZ effect and the true velocities of halos. We first created mocks of the 2D maps of SZ effects in seven frequency bands for galaxy clusters spanning a broad range of masses in TNG 300-3. We then derived the line-of-sight (LOS) peculiar velocities of galaxy clusters by applying an analytical formula to fit the spectra of the SZ effect. We find that the analytical formula-fitting method tends to overestimate the peculiar velocities of galaxy clusters, regardless of whether they approach us or recede from us. However, when galaxy clusters larger than 500 km/s were excluded, the slopes of the relations between vz, kSZ and the real LOS velocities of galaxy clusters (vz, halo) were consistent with unity within the errors. Additionally, we further accounted for observational noise for different-aperture telescopes under different precipitable water vapor. The slopes for most cases are consistent with unity within the errors, and the errors of the slopes slightly increase with higher observation noise. The |vz, kSZ − vz, halo| exhibits no obvious trend with increasing concentration and M200. It is noteworthy that the |vz, kSZ − vz, halo| is relatively high for galaxy clusters with strong active galactic nucleus feedback and high star formation rates. This indicates that these physical processes affect the cluster dynamic state and may reflect on the peculiar velocity estimated from the kSZ effect.
Key words: galaxies: clusters: general / galaxies: clusters: intracluster medium
© The Authors 2025
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1. Introduction
Galaxy clusters are the largest self-gravitationally bound systems in the Universe. They contain large quantities of gas, member galaxies, and dark matter (DM). In the hierarchical model of cosmic structure formation, galaxy clusters form and evolve through the merging of galaxy groups and clusters and continuous accretion of surrounding gas (Kravtsov & Borgani 2012). Galaxy clusters offer powerful probes of the large-scale structure and velocity field of the Universe (Carlstrom et al. 2002; Allen et al. 2011; Planck Collaboration XIII 2016) and provide laboratories for understanding the interaction of galaxies, active galactic nuclei (AGN), and star formation in different environments.
The traditional method for studying the intracluster medium (ICM) within galaxy clusters is through X-ray emissions from the hot gas. However, an X-ray detection becomes increasingly challenging at z ≳ 1 because the X-ray flux decreases rapidly with increasing z (van Marrewijk et al. 2024). Additionally, many existing X-ray observations are hampered by a low spectral resolution, which makes it difficult for X-ray observations to provide information about the motion of hot gas in clusters. In contrast, the Sunyaev-Zel’dovich (SZ) effect provides a unique and redshift-independent method of probing the ICM. The SZ effect contains the thermal-SZ (tSZ) effect (Sunyaev & Zeldovich 1972), characterized by the distortion of the cosmic microwave background (CMB) spectrum through inverse-Compton scattering of CMB photons by electrons in the ICM, and the kinematic-SZ (kSZ) effect (Sunyaev & Zeldovich 1980), arising from the scattering of CMB photons by the electrons in ICM due to the bulk motion of the cluster. The tSZ effect is proportional to the line-of-sight (LOS) integral of the electron pressure: 1, while the kSZ effect is related to the integrated LOS electron density and the gas velocity with respect to the CMB reference frame:
(Carlstrom et al. 2002; Planck Collaboration XIII 2016).
Compared to the tSZ effect, the kSZ effect is much weaker and difficult to detect unless the gas velocity reaches a few thousandths of the speed of light (Birkinshaw 1999). Therefore, the direct observation of the kSZ signal in individual clusters requires high-sensitivity observations. The kSZ signal was only detected for a few individual galaxy clusters (Sayers et al. 2013; Adam et al. 2017) or through statistical measurements (Hand et al. 2012; Planck Collaboration XXXVI 2016; Tanimura et al. 2021; Chen et al. 2022). This limits our understanding of the cosmic velocity field. Upcoming telescopes such as the Cerro-Chajnantor Atacama Telescope prime (CCAT-prime; Parshley et al. 2018; Chapman et al. 2022; Aravena 2023) and the Leighton Chajnantor Telescope (LCT; Yao et al. 2023) are expected to obtain measurements of the SZ effect in more than three bands and to detect the faint kSZ effect to trace the peculiar velocities of a large number of galaxy clusters (Terry et al. 2019). These measurements can provide us maps of the cosmic velocity field and powerful constraints on the structure growth of the Universe (Peebles 1980; Okumura & Taruya 2022; Alonso et al. 2016; Adams & Blake 2017, 2020; Turner et al. 2023; Lai et al. 2023; Shi et al. 2024). If clusters were ideally spherically symmetric systems with minimal internal motions, then the kSZ effect would provide a relatively straightforward method of measuring their peculiar velocities. Nonetheless, galaxy clusters are complex systems with significant internal flows (Ricker & Sarazin 2001; Sun & Murray 2002) driven by mergers and physical processes such as star formation and AGN feedback. Consequently, it becomes imperative to study the bias of the peculiar velocities that are estimated through the kSZ effect. The deviation of the slope of the vz, kSZ − vz, halo relation from unity determines how precisely the kSZ effect can estimate the true peculiar velocity and constrain the cosmological parameters.
Only a few previous studies have focused on how accurate the peculiar velocity from the kSZ effect is. Nagai et al. (2003) used a high-resolution simulation of a galaxy cluster with a complex internal structure and motion to investigate the errors in the peculiar velocity estimated through the kSZ effect. They found a dispersion in the peculiar velocity estimated through the kSZ effect of 50 − 100 km s−1, and the velocity averaged over the kSZ map over the cluster virial radius provided an unbiased estimate. However, Nagai et al. (2003) did not simulate the observations of the SZ effect at different frequencies, and they only studied one single galaxy cluster. Wang et al. (2021) estimated the peculiar velocities using deep neural networks based on Magneticum simulations and found an average uncertainty of the peculiar velocity estimated through the kSZ of about 25%. Although the sample of galaxy clusters in this work is large, the results were based on the temperature distortion maps of the SZ effect instead of mock observations that mimic real observations.
We carried out more realistic mocks of the SZ maps at seven frequencies for a sample of 155 galaxy clusters with masses higher than 1013.5 M⊙ and redshift ranges from 0.2 to 1.0. Based on the mock observations, we obtained the LOS peculiar velocity (vz, kSZ) via the kSZ effect for each cluster and compared it with the real peculiar velocity (vz, halo) to study the bias of vz, kSZ. We present the sample of galaxy clusters from the TNG300-3 simulation run of IllustrisTNG (Nelson et al. 2019) in Section 2.1 and the procedures of simulated observations of the SZ effect in Section 2.2. Section 3 shows the main results of this work and compares the peculiar velocities estimated via the kSZ effect with the real peculiar velocities of galaxy clusters. We also present the cases with observation noise. In Section 4, we discuss the impacts of the galaxy cluster properties on the vz, kSZ. Section 5 concludes our work. Throughout this work, we assume H0 = 70 km s−1 Mpc−1, Ω0 = 0.3, and ΩΛ = 0.7.
2. Method
2.1. Sample
We used data from IllustrisTNG (Nelson et al. 2019), a cosmological magnetohydrodynamical simulation that models the evolution of dark matter, gas, stars, and supermassive black holes (Springel 2010; Pakmor et al. 2016). We focused on the SZ effect of the whole galaxy cluster, whose typical size is a few megaparsec. We used the TNG300-3 simulation run with a large volume of 302.63 Mpc but a low baryonic and DM resolution of 1.1 × 107 M⊙ and 5.9 × 107 M⊙. We selected a large sample of subhalos with masses higher than 1013.5 M⊙ and primary_flag = 1, indicating that the halo is either the most massive or primary subhalo of the friends-of-friends (FoF) halo. The redshift range was set to be between 0.2 and 1. The sample contained 155 galaxy clusters in total, with 50 objects in each mass range of 1013.5 ∼ 1014.0 M⊙, 1014.0 ∼ 1014.5 M⊙, and 1014.5 ∼ 1015.0 M⊙, along with five objects with masses higher than 1015.0 M⊙. We only chose subhalos from snapshots of the TNG300-3 with full information (z = 0.2, 0.3, 0.4, 0.5, 0.7, 1.0). Our clusters followed the redshift distribution of the subhalos in each mass range, which is shown in Figure 1.
![]() |
Fig. 1. Redshift distributions of the chosen galaxy clusters for the mass ranges 1013.5 ∼ 1014.0 M⊙, 1014.0 ∼ 1014.5 M⊙, 1014.5 ∼ 1015.0 M⊙, and > 1015.0 M⊙. The whole sample includes 155 galaxy clusters, the first three mass ranges contain 50 sources each, and the most massive range only contains 5 sources. |
2.2. Mock observations of the SZ effect
To investigate whether there is a systematic error in the peculiar velocity estimated from the kSZ effect, we needed to simulate observations of the SZ effect of galaxy clusters at different frequencies. We used SZpack (Chluba et al. 2012, 2013), a numerical library allowing fast and precise computation of the SZ signals for each gas cell of the galaxy clusters in TNG300-3 and projected it onto the x-y plane.
To obtain the SZ signals at different frequencies, we needed several fundamental physical quantities of each gas particle as the input to SZpack, including the optical depth (τ), the gas temperatures (kT), and the velocities of gas particles (vgas). The gas velocities are directly given in the TNG300-3 dataset. The optical depths were derived from
where σT is the Thomson cross-section and ne is the electron number densities. ne is derived from the electron abundance xe, the hydrogen mass fraction XH = 0.76, and the density of gas particles following the equation
where nH = XH * ρ/mp and xe = ne/nH is given by PartType0 snapshot fields of the TNG dataset. The gas temperatures were derived from xe and the internal energy u following
where the adiabatic index γ is 5/3 and the Boltzmann constant kB is in centimeter-gram-second units. UnitMass and UnitEnergy are code units, and . u can also be found in PartType0 snapshot fields of the TNG dataset. μ is the mean molecular weight,
We converted all quantities from comoving units into physical units.
We simulated the SZ signal at seven frequently used bands, including 90 GHz (frequency band of MUSTANG-2 camera on the Green Bank Telescope, GBT; Dicker et al. 2014), 150 GHz, 260 GHz (frequency bands of the NIKA2 (New IRAM KID Arrays 2) instrument on the IRAM 30M; Adam et al. 2018; Catalano et al. 2018), and 230 GHz, 345 GHz, 460 GHz, and 650 GHz (frequency bands of the LCT; Yao et al. 2023).
We also studied the peculiar velocities obtained from the kSZ effect under certain observational noises. We estimated the observation noise for telescopes with assumed conditions, and we list them in Table 1. We considered one-hour on-source observations for each pointing of every galaxy cluster with 10-meter, 30-meter, and 100-meter radio telescopes at Chajnantor Plateau (CBI) under different precipitable water vapor (PWV). The observation noise was obtained following the equations in He et al. (2022),
PWV and the observation noise (rms) of a one-hour on-source observation for each pointing with 10-meter, 30-meter, and 100-meter radio telescopes at seven frequency bands, with an elevation of 45 degrees and under different PWV at CBI.
where Tsky is the sky temperature, Tatm is the representative atmosphere temperature, and τsky is the sky opacity. We adopted Tatm = 273 K, an average atmospheric temperature for the CBI (Cortés et al. 2020). In this equation, e−τsky = θtrans is the transmission, for which we adopted the transmission of the atmosphere at the CBI2. Furthermore, Tsys can be calculated following the equation in He et al. (2022),
in which Trx is the receiver temperature. It is typically 10–20 K, so we assumed Trx = 10 K here. The system equivalent flux density (SEFD) in Jy is
where G is the antenna gain, G = Aeff/2kB. Here, Aeff is the effective area of the antenna, and , where ηA ∼ 0.7 for a well-designed telescope, and D is the diameter of the telescope. Then, the rms noise in Jy/beam of the observation for a single-dish telescope is
where the c accounts for the quantization and correlator efficiencies, and we assumed c = 0.9. Δν is the frequency bandwidth and tint is the integration time of the observation. Then, we added Gaussian noise to each pixel on the projected maps of our mock observations with the rms given in Table 1. One caveat of our study is that we did not simulate the contamination of the cosmic infrared background (CIB), which includes the clustering of dusty star-forming galaxies (DSFGs) and the shot noise from DSFGs and radio sources. At higher frequencies, the CIB becomes significantly contaminated by the SZ signal (Hauser et al. 1998; Hauser & Dwek 2001). The clustering of DSFGs can be computed by using the CIB power spectrum derived from CIB anisotropies surveys (e.g. Planck Collaboration XVI 2014), which can then be extrapolated to the observation frequencies. For unresolved DSFGs and radio sources, the shot noise can be computed using the models from Béthermin et al. (2012) for DSFGs and Tucci et al. (2011) for radio sources. Additionally, contamination from compact radio sources and submillimeter point sources can, in principle, be subtracted using multiband and high-resolution observations. The CMB can also be modeled with the software HEALPix (Górski et al. 2005) to mitigate its contamination. Therefore, our analysis assumed that the CMB and CIB contaminations can be well separated.
3. Results
We projected the SZ signals of all gas particles onto the x–y plane for the galaxy clusters in our sample. An example of SZ effect maps at seven frequencies is shown in Figure 2a. After we obtained the projected maps of the SZ effect for each galaxy cluster, we measured the fluxes of the SZ effect within four different aperture sizes: 0.25R500, 0.5R500, R500, and R200 (corresponding to the circles in the SZ maps in Figure 2a). The four panels in Figure 2b depict examples of the derived broadband SED for these four apertures. We fit the data points with an analytic formula of the frequency for the SZ effect provided by Nozawa et al. (2006), which describes the change in photon number density over time due to Compton scattering with electrons with relativistic corrections for a moving system. The fitting procedure allowed us to derive the best-fit results for the optical depth (τ), gas temperatures (kT), and LOS (z-direction) peculiar velocities (vz, kSZ).
![]() |
Fig. 2. SZ effect of a galaxy cluster at various frequencies and the corresponding SZ spectra. Panel a: examples of the SZ effect of a galaxy cluster at various frequencies: 90 GHz, 150 GHz, 230 GHz, 260 GHz, 345 GHz, 460 GHz, and 650 GHz. The circles represent four different aperture sizes for photometry. The solid red circle corresponds to 0.25R500, the dotted cyan circle to 0.5R500, the dashed pink circle to R500, and the dash-dotted white circle to R200. Panel b: SZ spectra of the galaxy cluster in panel a measured at the cluster center using different aperture sizes of 0.25R500, 0.5R500, R500, and R200, corresponding to the circles in the SZ maps in panel a. The black dots represent the measured fluxes of the SZ effect at different frequencies. We fit the black dots with an analytic formula for the SZ effect (Nozawa et al. 2006). The dashed black lines show the results of fitting the total SZ signal, and the dashed red and blue lines represent the contributions from the tSZ and kSZ effects, respectively. The dotted green line corresponds to the best-fit result of the tSZ-only model. The residuals of the best-fit SZ model and the tSZ-only model with respect to the data are shown as black dots and green crosses in the bottom panels. |
The Compton y-parameter quantifies the photon energy change due to inverse-Compton scattering. The radial distribution of ytSZ shows the distribution of the electron pressure within the cluster, which is typically highest at the center and decreases with increasing radius. We show the mean radial profiles of the tSZ effect () for the four mass bins of our sample in Figure 3. In all mass bins, the ytSZ radial profiles decrease from the cluster center to the outer regions. More massive clusters typically exhibit higher ytSZ radial profiles that correspond to a higher electron pressure.
![]() |
Fig. 3. Radial profiles of the Compton ytSZ parameter for the four mass bins of our sample. The radii are normalized to R200 of the clusters. The lines represent the mean radial profile of ytSZ, and the shaded regions indicate the scatter. |
Our primary focus was the vz, kSZ estimated from the mock observations of the SZ effect. We first present the results without observational noise. Figure 4 presents one of the main results of our work, where we compare the vz, kSZ and the true LOS (z-direction) velocities of the entire galaxy cluster (vz, halo) based on the measurements with various aperture sizes. The vz, kSZ values were derived by fitting the SZ formula using a nonlinear least-squares function. As shown in the figure, there are outliers in all cases, and the slopes of the vz, kSZ − vz, halo relations for all clusters in our sample are greater than one, ranging from 1.07 to 1.12. Galaxy clusters with a high measured vz, kSZ are overestimated in LOS velocities and show large deviations from vz, halo, which is likely due to the discrepancy between vz, gas and vz, halo, as shown in Figure 5. In Figure 5, we show the relations between the vz, gas and vz, halo for all galaxy clusters. The rms scatter decreases with larger measuring aperture sizes. The slopes of the vz, gas − vz, halo relations of our selected clusters are slightly larger than one, which may also slightly contribute to the overestimation of the vz, kSZ, regardless of whether the galaxy clusters approach us or move away from us in Figure 4. We also checked the vz, gas − vz, halo relations for ∼3000 TNG halos with higher redshift and mass ranges. The slopes are within ∼3% of one. These overestimates affect the slopes of the linear-fitting results, which are significantly larger than one. Thus, we fit the vz, kSZ − vz, halo relation after excluding the clusters with |vz, kSZ|> 500 km s−1. Except for the case with a photometric aperture size of R200, the slopes are consistent with unity within the errors. There is a significant difference between the slope of the vz, kSZ − vz, halo relation with photometric aperture sizes of R200 (0.90) and unity, which is due to underestimates of vz, kSZ for many clusters. The underestimates are likely caused by the large fitting errors with very low average fluxes of the SZ effect out to R200. The root mean-square (rms) scatter around the vz, kSZ = vz, halo line at the given x-axis can reach approximately 132.9 km s−1 for the smallest aperture size of 0.25R500. The scatters roughly decrease with larger aperture sizes, and have values of ∼89.0 km s−1 and 77.8 km s−1 for apertures of 0.5R500 and R500, respectively. The rms scatter of the results obtained with apertures of R200 increases slightly to ∼89.8 km s−1.
![]() |
Fig. 4. Comparison between the LOS (z-direction) velocities of galaxy clusters (vz, halo) given by TNG300-3 and the velocities obtained from the SZ effect (vz, kSZ) measured with different aperture sizes of 0.25R500, 0.5R500, R500, and R200. The solid lines represent vz, halo = vz, kSZ, the colored regions represent the rms scatter of the vz, halo = vz, kSZ line at a given x-axis for all galaxy clusters in our sample. The dashed lines are the best-fit lines of the data points excluding |vz, kSZ|> 500 km s−1, and ‘m’ in the legends are the slopes of the best-fit lines. |
![]() |
Fig. 5. Comparison between mass-weighted gas velocities (vz, gas) and vz, halo within different aperture sizes of 0.25R500, 0.5R500, R500, and R200. The solid lines represent vz, gas = vz, halo, the colored regions represent the rms scatter of the vz, gas = vz, halo line at a given x-axis for all galaxy clusters in our sample. The dashed lines are the best-fit lines of the data points, and ‘m’ in the legends are the slopes of the best-fit lines. |
We also tested the case with measurements of an X-ray spectroscopic temperature TX by limiting kT to be within 0.65 ∼ 1.35 of the gas-mass-weighted temperature Tgas. This bound for kT is based on Nagai et al. (2007), which studied the relation between TX and Tgas. Figure 6 shows the vz, kSZ − vz, halo relations with vz, kSZ derived by fitting the SZ formula using the nonlinear least-squares function given the bounds of kT. We excluded the galaxy clusters with |vz, kSZ|> 500 km s−1 to carry out the linear fitting to the vz, kSZ − vz, halo relation. In this case, the rms scatter of the vz, kSZ − vz, halo relations is smaller than the results in Figure 4, and the number of outliers decreases as well. This suggests that measurements of TX are helpful when vz, kSZ is estimated from the kSZ effect, especially for galaxy clusters with faint SZ signals. The slopes, ranging from 0.94 to 0.99, are consistent with unity, which is also within the errors of the slopes. Our results differ from those of Wang et al. (2021), who identified slopes much smaller than unity. They estimated the vz, kSZ from temperature distortion maps of the SZ effect instead of simulated 2D observations and the derived broad SED at different frequencies as we did.
![]() |
Fig. 6. Comparison between the LOS (z-direction) velocities of the galaxy clusters (vz, halo) given by TNG300-3 and the velocities obtained with the SZ effect (vz, kSZ) measured with different aperture sizes of 0.25R500, 0.5R500, R500, and R200, given the bounds for kT. The solid lines represent vz, halo = vz, kSZ, and the colored regions represent the rms scatter of the vz, halo = vz, kSZ line at given x-axis for all galaxy clusters in our sample. The dashed lines are the best-fit lines of the data points excluding |vz, kSZ|> 500 km s−1, and ‘m’ in the legends are the slopes of the best-fit lines. |
In Figure 7, we present the median values of the differences between vz, kSZ and vz, halo, denoted as vz, kSZ − vz, halo, for the results obtained with different apertures under given bounds for kT. The error bars represent the standard deviations from the median of vz, kSZ − vz, halo. The median values of vz, kSZ − vz, halo with different aperture sizes are all negative, which is possibly because it is easier to detect galaxy clusters that approach us. When galaxy clusters approach us, the flux of the kSZ effect is positive and increases the peak value of the SZ effect, which is easier to detect. The median of vz, kSZ − vz, halo with larger apertures show larger deviations from 0 km s−1. We also compared the vz, kSZ with the true velocity of hot gas – vz, gas in Figure 8 and linearly fit the vz, kSZ − vz, gas relations for all galaxy clusters. The slopes are all larger than one, indicating that in general, the method that uses the analytic formula fitting (Nozawa et al. 2006) tends to overestimate the LOS peculiar velocities based on the kSZ effect, regardless of whether galaxy clusters approach us or move away from us.
![]() |
Fig. 7. Median values (dots with values) of the difference between vz, kSZ and vz, halo with measurements using different aperture sizes. The error bars show the standard deviation of (vz, kSZ − vz, halo). |
![]() |
Fig. 8. Comparison between vz, gas and vz, kSZ measured with different aperture sizes of 0.25R500, 0.5R500, R500, and R200, given the bounds for kT. The solid lines represent vz, kSZ = vz, gas, and the colored regions represent the rms scatter of the vz, kSZ = vz, gas line at given x-axis for all galaxy clusters in our sample. The dashed lines are the best-fit lines of the data points, and ‘m’ in the legends are the slopes of the best-fit lines. |
We further accounted for observational noises when simulating observations of the SZ effect. We adopted the noises observed with 10-meter, 30-meter, and 100-meter telescopes with an elevation of 45 degrees operating under various PWV conditions at the CBI as listed in Table 1. We conducted the measurements with aperture sizes of 0.25R500, 0.5R500, R500, and R200 separately. The results for the 30-meter telescope, with an aperture size of R200 for the SZ flux measurement, are visualized in Figure 9. We only included clusters with measurements with a signal-to-noise ratio (S/N) greater than 2 for at least three frequencies. The vz, kSZ values were derived by fitting the SZ signal with observational errors using a Markov chain Monte Carlo (MCMC) method. The linear fitting results of the vz, kSZ − vz, halo relation excluding the clusters with |vz, kSZ|> 500 km s−1 for all cases are listed in Table 2. As shown in Table 2, except for the 100-meter telescope, fewer sources can be detected with increasing PWV, and in general, the slope errors also increase slightly. For the 100-meter telescope, a decrease in PWV does not significantly improve the sensitivity and the accuracy of vz, kSZ. Instead, additional clusters with weak SZ signals will be detected under lower PWV, which introduces more uncertainties. Overall, in most cases, the slopes are consistent with one within the errors. However, the rms scatter of the vz, kSZ = vz, halo line at the given x-axis does not exhibit a clear decreasing trend with reduced observational noises. With lower observation noise, more faint sources are detected, but their average SZ fluxes are very low, leading to larger errors in estimating vz, kSZ and larger rms scatter of the vz, kSZ = vz, halo line at the given x-axis.
![]() |
Fig. 9. Comparison between vz, halo and vz, kSZ under different weather conditions (pwv = 0.5 mm, 1.0 mm, 1.5 mm, and 2.0 mm) with given bounds for kT for a 30-meter radio telescope. The results in these four figures are based on measurements with an aperture size of R200. The solid lines represent vz, kSZ = vz, halo, and the gray regions represent the rms scatters of the vz, kSZ = vz, halo lines at a given x-axis for all galaxy clusters in our sample. The dashed lines are the best-fit lines of the data points excluding |vz, kSZ|> 500 km s−1, and ‘m’ in the legends are the slopes of the best-fit lines. |
vz, kSZ − vz, halo relations for telescopes with different sizes under different PWV conditions.
4. Discussion
Any departure from spherical symmetry in the mass distribution of a galaxy cluster can introduce bias into the estimation of key physical properties such as mass, centroid, and peculiar velocity (Lee et al. 2018; Yuan & Han 2020). Consequently, a disturbed structure can have a significant impact on the cosmological constraints derived from various relations, including the YSZ − M relation and the concentration–mass relation (Macciò et al. 2007; Ludlow et al. 2012). The bias in estimating physical properties from SZ effect images depends on the dynamical state of the galaxy clusters (Yuan & Han 2020).
The concentration parameter, c200, represents the central density of the halo and is used to characterize the assembly history, which reflects the time when it formed (Bullock et al. 2001; Eke et al. 2001; Zhao et al. 2003). We compare |vz, kSZ − vz, halo| and c200 in Figure 10. In TNG300-3 run, the DM particles have the same mass of ∼0.255 × 1010 M⊙/h given in MassTable in Header field, and the positions of every DM particle are given in PartType1 snapshot fields of the TNG dataset. We fit the DM mass density distribution with the Navarro–Frenk–White (NFW) profile (Navarro et al. 1996, 1997),
![]() |
Fig. 10. Distribution of |vz, kSZ − vz, halo| against the concentration (c200) of galaxy clusters in our sample for measurements with Raper = R200. The colors of the data points represent their redshift. |
Then, we obtained rs by fitting the NFW formula to the DM density profile. R200 is the radius within which the mean density of the halo is 200 times the critical density, and it can be derived based on the radial density profile of all particles (DM, gas, stars, and black holes). The concentration, c200 is r200/rs. We do not observe an obvious trend between |vz, kSZ − vz, halo| and c200. This differs from our expectation that higher c200 would lead to higher |vz, kSZ − vz, halo| because clusters showing signs of dynamical activity tend to exhibit high concentrations (Sereno & Covone 2013), and there exists an empirical relation M200 − c200 within galaxy clusters, where the concentration decreases as the halo mass increases (Buote et al. 2007; Comerford & Natarajan 2007; Ludlow et al. 2012; Du et al. 2015; Amodeo et al. 2016). In Figure 11, we compare M200 and |vz, kSZ − vz, halo| and find no obvious trend in the median value of |vz, kSZ − vz, halo| in the halo mass bins either.
![]() |
Fig. 11. Distribution of |vz, kSZ − vz, halo| against M200 of galaxy clusters in our sample for measurements with Raper = R200. The black dots with error bars are the medians of |vz, kSZ − vz, halo| and their standard deviations at different M200. |
Physical processes such as AGN and star formation activity can affect the dynamical state of galaxy clusters. We plot the AGN-feedback/M200 − |vz, kSZ − vz, halo| relation in Figure 12. In this figure, most clusters with high AGN-feedback/M200 exhibit relatively high |vz, kSZ − vz, halo|, which can reach tens to some hundred km s−1. The dashed line in Figure 12 represents the rms scatter (∼55.15 km s−1) of the vz, kSZ = vz, gas line at a given x-axis with a photometric aperture size of R200, which reflects the error introduced by fitting the SZ spectrum with an analytic formula. This suggests that AGN feedback indeed induces deviations from equilibrium in galaxy clusters, and the perturbation of the dynamical state of clusters caused by AGN feedback appears to outweigh the influence of the self-gravitational potential energy of the cluster. Most of the clusters with high AGN feedback are located around or largely above the value of the fitting error. Although |vz, kSZ − vz, halo| for clusters with strong AGN feedback are relatively high and may reflect the effect of AGN feedback on the cluster dynamical state, these velocity differences are likely to be obscured for measurement errors and contamination. Figure 13 shows the relation between the star formation rate (SFR) and |vz, kSZ − vz, halo| relation. Clusters with high SFR show a relatively high |vz, kSZ − vz, halo|, but their |vz, kSZ − vz, halo| are just about or lower than the rms scatter of the vz, kSZ = vz, gas line at the given x-axis.
![]() |
Fig. 12. Distribution of |vz, kSZ − vz, halo| against the ratios of total cumulative AGN feedback to M200 of galaxy clusters in our sample for measurements with Raper = R200. The dashed line represents the rms scatter (∼55.15 km s−1) of the vz, kSZ = vz, gas line at a given x-axis for all galaxy clusters based on measurements using apertures of R200, as shown in Figure 8. The color bar represents the AGN feedback of the galaxy clusters. |
![]() |
Fig. 13. Distribution of |vz, kSZ − vz, halo| against the SFR of galaxy clusters in our sample for measurements with Raper = R200. The dashed line represents the rms scatter (∼55.15 km s−1) of the vz, kSZ = vz, halo line at a given x-axis for all galaxy clusters based on measurements using apertures of R200, as shown in Figure 8. |
5. Conclusions
We selected a sample of galaxy clusters with varying redshift and total mass from TNG300-3 and generated mock 2D maps of the SZ effect at seven frequencies. We measured the fluxes of the SZ effect and fit the spectra with an analytic formula (Nozawa et al. 2006), and we obtained the best-fit parameters including the peculiar velocity vz, kSZ. We compared the velocity vz, kSZ with the true LOS peculiar velocities vz, halo of the halos. The main conclusions are listed below.
(1) Overall, vz, kSZ represents the halo velocity well. For high peculiar velocities, they tend to be overestimated. When they are excluded, the slopes of the vz, kSZ − vz, halo relations are consistent with unity within the errors.
(2) The median values of (vz, kSZ − vz, halo) are lower than 0 km s−1, indicating that it is easier to detect the kSZ effect of galaxy clusters that approach us.
(3) When we considered the observational noise, most of the slopes of vz, kSZ − vz, halo relations are consistent with one within the errors. In general, with lower observation noise, more sources can be detected, and the slope errors decrease slightly.
(4) The |vz, kSZ − vz, halo| exhibits no discernible trend with increasing c200 parameter, which is associated with the halo assembly history. In addition, the |vz, kSZ − vz, halo| shows no obvious trend with increasing M200.
(5) Physical processes such as AGN feedback and star formation can influence the dynamical states of hot gas in galaxy clusters. The clusters with high AGN feedback/M200 and SFR exhibit relatively large discrepancies between vz, kSZ and vz, halo that reach several tens to some hundred km s−1. These variations may influence the peculiar velocity estimates derived with the kSZ effect. However, many of these values are comparable to or lower than the average fitting error.
Acknowledgments
This work acknowledges the support from the National Key R&D Program of China No. 2023YFA1608204, the National Natural Science Foundation of China (NSFC grants 12141301, 12121003, 12333002).
References
- Adam, R., Bartalucci, I., Pratt, G. W., et al. 2017, A&A, 598, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Adam, R., Adane, A., Ade, P. A. R., et al. 2018, A&A, 609, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Adams, C., & Blake, C. 2017, MNRAS, 471, 839 [NASA ADS] [CrossRef] [Google Scholar]
- Adams, C., & Blake, C. 2020, MNRAS, 494, 3275 [NASA ADS] [CrossRef] [Google Scholar]
- Allen, S. W., Evrard, A. E., & Mantz, A. B. 2011, ARA&A, 49, 409 [Google Scholar]
- Alonso, D., Louis, T., Bull, P., & Ferreira, P. G. 2016, Phys. Rev. D, 94, 043522 [NASA ADS] [CrossRef] [Google Scholar]
- Amodeo, S., Ettori, S., Capasso, R., & Sereno, M. 2016, A&A, 590, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Béthermin, M., Daddi, E., Magdis, G., et al. 2012, ApJ, 757, L23 [Google Scholar]
- Birkinshaw, M. 1999, Phys. Rep., 310, 97 [Google Scholar]
- Bullock, J. S., Kolatt, T. S., Sigad, Y., et al. 2001, MNRAS, 321, 559 [Google Scholar]
- Buote, D. A., Gastaldello, F., Humphrey, P. J., et al. 2007, ApJ, 664, 123 [NASA ADS] [CrossRef] [Google Scholar]
- Carlstrom, J. E., Holder, G. P., & Reese, E. D. 2002, ARA&A, 40, 643 [Google Scholar]
- Catalano, A., Adam, R., Ade, P. A. R., et al. 2018, Journal of Low Temperature Physics, 193, 916 [NASA ADS] [CrossRef] [Google Scholar]
- CCAT-Prime Collaboration (Aravena, M., et al.) 2023, ApJS, 264, 7 [NASA ADS] [CrossRef] [Google Scholar]
- Chapman, S. C., Huber, A. I., Sinclair, A. K., et al. 2022, in Millimeter, Submillimeter and Far-Infrared Detectors and Instrumentation for Astronomy XI, eds. J. Zmuidzinas, & J. R. Gao, SPIE Conf. Ser., 12190, 1219005 [NASA ADS] [Google Scholar]
- Chen, Z., Zhang, P., Yang, X., & Zheng, Y. 2022, MNRAS, 510, 5916 [NASA ADS] [CrossRef] [Google Scholar]
- Chluba, J., Nagai, D., Sazonov, S., & Nelson, K. 2012, MNRAS, 426, 510 [Google Scholar]
- Chluba, J., Switzer, E., Nelson, K., & Nagai, D. 2013, MNRAS, 430, 3054 [NASA ADS] [CrossRef] [Google Scholar]
- Comerford, J. M., & Natarajan, P. 2007, MNRAS, 379, 190 [NASA ADS] [CrossRef] [Google Scholar]
- Cortés, F., Cortés, K., Reeves, R., Bustos, R., & Radford, S. 2020, A&A, 640, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dicker, S. R., Ade, P. A. R., Aguirre, J., et al. 2014, in Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy VII, eds. W. S. Holland, & J. Zmuidzinas, SPIE Conf. Ser., 9153, 91530J [NASA ADS] [CrossRef] [Google Scholar]
- Du, W., Fan, Z., Shan, H., et al. 2015, ApJ, 814, 120 [NASA ADS] [CrossRef] [Google Scholar]
- Eke, V. R., Navarro, J. F., & Steinmetz, M. 2001, ApJ, 554, 114 [NASA ADS] [CrossRef] [Google Scholar]
- Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 [Google Scholar]
- Hand, N., Addison, G. E., Aubourg, E., et al. 2012, Phys. Rev. Lett., 109, 041101 [NASA ADS] [CrossRef] [Google Scholar]
- Hauser, M. G., & Dwek, E. 2001, ARA&A, 39, 249 [Google Scholar]
- Hauser, M. G., Arendt, R. G., Kelsall, T., et al. 1998, ApJ, 508, 25 [Google Scholar]
- He, H., Dent, W. R. F., & Wilson, C. 2022, PASP, 134, 125001 [NASA ADS] [CrossRef] [Google Scholar]
- Kravtsov, A. V., & Borgani, S. 2012, ARA&A, 50, 353 [Google Scholar]
- Lai, Y., Howlett, C., & Davis, T. M. 2023, MNRAS, 518, 1840 [Google Scholar]
- Lee, B. E., Le Brun, A. M. C., Haq, M. E., et al. 2018, MNRAS, 479, 890 [Google Scholar]
- Ludlow, A. D., Navarro, J. F., Li, M., et al. 2012, MNRAS, 427, 1322 [NASA ADS] [CrossRef] [Google Scholar]
- Macciò, A. V., Dutton, A. A., van den Bosch, F. C., et al. 2007, MNRAS, 378, 55 [Google Scholar]
- Nagai, D., Kravtsov, A. V., & Kosowsky, A. 2003, ApJ, 587, 524 [NASA ADS] [CrossRef] [Google Scholar]
- Nagai, D., Vikhlinin, A., & Kravtsov, A. V. 2007, ApJ, 655, 98 [Google Scholar]
- Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 563 [Google Scholar]
- Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 [Google Scholar]
- Nelson, D., Springel, V., Pillepich, A., et al. 2019, Computational Astrophysics and Cosmology, 6, 2 [CrossRef] [Google Scholar]
- Nozawa, S., Itoh, N., Suda, Y., & Ohhata, Y. 2006, Nuovo Cimento B Serie, 121, 487 [NASA ADS] [Google Scholar]
- Okumura, T., & Taruya, A. 2022, Phys. Rev. D, 106, 043523 [NASA ADS] [CrossRef] [Google Scholar]
- Pakmor, R., Springel, V., Bauer, A., et al. 2016, MNRAS, 455, 1134 [Google Scholar]
- Parshley, S. C., Kronshage, J., Blair, J., et al. 2018, in Ground-based and Airborne Telescopes VII, eds. H. K. Marshall, & J. Spyromilio, SPIE Conf. Ser., 10700, 107005X [NASA ADS] [Google Scholar]
- Peebles, P. J. E. 1980, The Large-scale Structure of the Universe (Princeton: Princeton University Press) [Google Scholar]
- Planck Collaboration XVI. 2014, A&A, 571, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration XXXVI. 2016, A&A, 586, A140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration XIII. 2016, A&A, 594, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ricker, P. M., & Sarazin, C. L. 2001, ApJ, 561, 621 [NASA ADS] [CrossRef] [Google Scholar]
- Sayers, J., Czakon, N. G., Mantz, A., et al. 2013, ApJ, 768, 177 [Google Scholar]
- Sereno, M., & Covone, G. 2013, MNRAS, 434, 878 [Google Scholar]
- Shi, Y., Zhang, P., Mao, S., & Gu, Q. 2024, MNRAS, 528, 4922 [NASA ADS] [CrossRef] [Google Scholar]
- Springel, V. 2010, MNRAS, 401, 791 [Google Scholar]
- Sun, M., & Murray, S. S. 2002, ApJ, 576, 708 [NASA ADS] [CrossRef] [Google Scholar]
- Sunyaev, R. A., & Zeldovich, Y. B. 1972, Comments on Astrophysics and Space Physics, 4, 173 [NASA ADS] [Google Scholar]
- Sunyaev, R. A., & Zeldovich, Y. B. 1980, MNRAS, 190, 413 [NASA ADS] [CrossRef] [Google Scholar]
- Tanimura, H., Zaroubi, S., & Aghanim, N. 2021, A&A, 645, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Terry, H., Battaglia, N., Basu, K., et al. 2019, Bulletin of the American Astronomical Society, 51, 213 [NASA ADS] [Google Scholar]
- Tucci, M., Toffolatti, L., de Zotti, G., & Martínez-González, E. 2011, A&A, 533, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Turner, R. J., Blake, C., & Ruggeri, R. 2023, MNRAS, 518, 2436 [Google Scholar]
- van Marrewijk, J., Di Mascolo, L., Gill, A. S., et al. 2024, A&A, 689, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wang, Y., Ramachandra, N., Salazar-Canizales, E. M., et al. 2021, MNRAS, 506, 1427 [NASA ADS] [CrossRef] [Google Scholar]
- Yao, Y.-W., Chen, W.-R., & Wang, Z. 2023, Research in Astronomy and Astrophysics, 23, 045013 [NASA ADS] [CrossRef] [Google Scholar]
- Yuan, Z. S., & Han, J. L. 2020, MNRAS, 497, 5485 [Google Scholar]
- Zhao, D. H., Jing, Y. P., Mo, H. J., & Börner, G. 2003, ApJ, 597, L9 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
PWV and the observation noise (rms) of a one-hour on-source observation for each pointing with 10-meter, 30-meter, and 100-meter radio telescopes at seven frequency bands, with an elevation of 45 degrees and under different PWV at CBI.
vz, kSZ − vz, halo relations for telescopes with different sizes under different PWV conditions.
All Figures
![]() |
Fig. 1. Redshift distributions of the chosen galaxy clusters for the mass ranges 1013.5 ∼ 1014.0 M⊙, 1014.0 ∼ 1014.5 M⊙, 1014.5 ∼ 1015.0 M⊙, and > 1015.0 M⊙. The whole sample includes 155 galaxy clusters, the first three mass ranges contain 50 sources each, and the most massive range only contains 5 sources. |
In the text |
![]() |
Fig. 2. SZ effect of a galaxy cluster at various frequencies and the corresponding SZ spectra. Panel a: examples of the SZ effect of a galaxy cluster at various frequencies: 90 GHz, 150 GHz, 230 GHz, 260 GHz, 345 GHz, 460 GHz, and 650 GHz. The circles represent four different aperture sizes for photometry. The solid red circle corresponds to 0.25R500, the dotted cyan circle to 0.5R500, the dashed pink circle to R500, and the dash-dotted white circle to R200. Panel b: SZ spectra of the galaxy cluster in panel a measured at the cluster center using different aperture sizes of 0.25R500, 0.5R500, R500, and R200, corresponding to the circles in the SZ maps in panel a. The black dots represent the measured fluxes of the SZ effect at different frequencies. We fit the black dots with an analytic formula for the SZ effect (Nozawa et al. 2006). The dashed black lines show the results of fitting the total SZ signal, and the dashed red and blue lines represent the contributions from the tSZ and kSZ effects, respectively. The dotted green line corresponds to the best-fit result of the tSZ-only model. The residuals of the best-fit SZ model and the tSZ-only model with respect to the data are shown as black dots and green crosses in the bottom panels. |
In the text |
![]() |
Fig. 3. Radial profiles of the Compton ytSZ parameter for the four mass bins of our sample. The radii are normalized to R200 of the clusters. The lines represent the mean radial profile of ytSZ, and the shaded regions indicate the scatter. |
In the text |
![]() |
Fig. 4. Comparison between the LOS (z-direction) velocities of galaxy clusters (vz, halo) given by TNG300-3 and the velocities obtained from the SZ effect (vz, kSZ) measured with different aperture sizes of 0.25R500, 0.5R500, R500, and R200. The solid lines represent vz, halo = vz, kSZ, the colored regions represent the rms scatter of the vz, halo = vz, kSZ line at a given x-axis for all galaxy clusters in our sample. The dashed lines are the best-fit lines of the data points excluding |vz, kSZ|> 500 km s−1, and ‘m’ in the legends are the slopes of the best-fit lines. |
In the text |
![]() |
Fig. 5. Comparison between mass-weighted gas velocities (vz, gas) and vz, halo within different aperture sizes of 0.25R500, 0.5R500, R500, and R200. The solid lines represent vz, gas = vz, halo, the colored regions represent the rms scatter of the vz, gas = vz, halo line at a given x-axis for all galaxy clusters in our sample. The dashed lines are the best-fit lines of the data points, and ‘m’ in the legends are the slopes of the best-fit lines. |
In the text |
![]() |
Fig. 6. Comparison between the LOS (z-direction) velocities of the galaxy clusters (vz, halo) given by TNG300-3 and the velocities obtained with the SZ effect (vz, kSZ) measured with different aperture sizes of 0.25R500, 0.5R500, R500, and R200, given the bounds for kT. The solid lines represent vz, halo = vz, kSZ, and the colored regions represent the rms scatter of the vz, halo = vz, kSZ line at given x-axis for all galaxy clusters in our sample. The dashed lines are the best-fit lines of the data points excluding |vz, kSZ|> 500 km s−1, and ‘m’ in the legends are the slopes of the best-fit lines. |
In the text |
![]() |
Fig. 7. Median values (dots with values) of the difference between vz, kSZ and vz, halo with measurements using different aperture sizes. The error bars show the standard deviation of (vz, kSZ − vz, halo). |
In the text |
![]() |
Fig. 8. Comparison between vz, gas and vz, kSZ measured with different aperture sizes of 0.25R500, 0.5R500, R500, and R200, given the bounds for kT. The solid lines represent vz, kSZ = vz, gas, and the colored regions represent the rms scatter of the vz, kSZ = vz, gas line at given x-axis for all galaxy clusters in our sample. The dashed lines are the best-fit lines of the data points, and ‘m’ in the legends are the slopes of the best-fit lines. |
In the text |
![]() |
Fig. 9. Comparison between vz, halo and vz, kSZ under different weather conditions (pwv = 0.5 mm, 1.0 mm, 1.5 mm, and 2.0 mm) with given bounds for kT for a 30-meter radio telescope. The results in these four figures are based on measurements with an aperture size of R200. The solid lines represent vz, kSZ = vz, halo, and the gray regions represent the rms scatters of the vz, kSZ = vz, halo lines at a given x-axis for all galaxy clusters in our sample. The dashed lines are the best-fit lines of the data points excluding |vz, kSZ|> 500 km s−1, and ‘m’ in the legends are the slopes of the best-fit lines. |
In the text |
![]() |
Fig. 10. Distribution of |vz, kSZ − vz, halo| against the concentration (c200) of galaxy clusters in our sample for measurements with Raper = R200. The colors of the data points represent their redshift. |
In the text |
![]() |
Fig. 11. Distribution of |vz, kSZ − vz, halo| against M200 of galaxy clusters in our sample for measurements with Raper = R200. The black dots with error bars are the medians of |vz, kSZ − vz, halo| and their standard deviations at different M200. |
In the text |
![]() |
Fig. 12. Distribution of |vz, kSZ − vz, halo| against the ratios of total cumulative AGN feedback to M200 of galaxy clusters in our sample for measurements with Raper = R200. The dashed line represents the rms scatter (∼55.15 km s−1) of the vz, kSZ = vz, gas line at a given x-axis for all galaxy clusters based on measurements using apertures of R200, as shown in Figure 8. The color bar represents the AGN feedback of the galaxy clusters. |
In the text |
![]() |
Fig. 13. Distribution of |vz, kSZ − vz, halo| against the SFR of galaxy clusters in our sample for measurements with Raper = R200. The dashed line represents the rms scatter (∼55.15 km s−1) of the vz, kSZ = vz, halo line at a given x-axis for all galaxy clusters based on measurements using apertures of R200, as shown in Figure 8. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.