© ESO, 2015

1. Introduction

Galaxy clusters offer several classes of observables reflecting their overall internal constitution: gravitational lensing effects in their weak and strong variants, X-ray emission, the thermal Sunyaev-Zel’dovich (tSZ) effect, and the kinematics of their member galaxies. We neglect radio emission, turbulence, and metal abundance in the intracluster gas and the population statistics of member galaxies because these effects are locally driven.

Gravitational lensing measures the gravitational tidal field, projected along the line of sight (l.o.s.) and thus directly probes the projected gravitational potential more precisely its curvature. The X-ray emission and the thermal Sunyaev-Zel’dovich effect depend on powers of the density and temperature of the intracluster medium. If equilibrium assumptions hold, hydrostatic and virial equilibrium foremost, the X-ray emission, the thermal Sunyaev-Zel’dovich effect and galaxy kinematics are also determined by the gravitational potential. Current discussions of the validity of these equilibrium assumptions can be found e.g. in Martino et al. (2014), von der Linden et al. (2014), Planck Collaboration Int. III (2013), Planck Collaboration XX (2014).

What is the gravitational potential that agrees best with all cluster observables? This question is relevant for different reasons. First, different observables trace the gravitational potential at different scales. The cluster core can be probed by the stellar kinematics of the brightest cluster galaxy, strong lensing and X-ray emission probe the innermost regions, weak lensing and galaxy kinematics probe large scales, and the tSZ effect falls in between (see Umetsu et al. 2014; Merten et al. 2015; Newman et al. 2013). Aiming at a reliable reconstruction of cluster density profiles from their cores to their outskirts, combining all observables into a unique potential reconstruction offers the advantage of covering all relevant scales in a single step. Second, comparing cluster potential reconstructions based on lensing, on the one hand, and based on the rest of the observables, on the other hand, allows us to test the equilibrium assumptions or possible deviations therefrom. Third, lensing and the other observables do not necessarily see the same gravitational potential. While lensing is sensitive to the sum of the Bardeen potentials, the other observables only probe the spatial potential. In general relativity, the two Bardeen potentials agree in case of negligible anisotropic stress. Differences in potential reconstructions based on lensing compared to other observables may also hint at deviations from relativity (see also Sartoris et al. 2014; Barreira et al. 2015, in this context).

The gravitational potential has the major and important advantage that it is a locally measurable quantity, which is directly related to the observables listed above (see also Angrick & Bartelmann 2009, 2012).

We have recently developed methods for reconstructing the projected gravitational potential of galaxy clusters from their X-ray emission, their tSZ effect and the kinematics of their member galaxies (Sarli et al. 2014; Konrad et al. 2013; Majer et al. 2013). These methods operate similarly, but with important differences in detail. That is, an observable is deprojected by means of the Richardson-Lucy algorithm, requiring symmetry assumptions. The deprojected quantities are related to the three-dimensional gravitational potential by relations derived from justifiable equilibrium assumptions. The gravitational potential can then be projected along the line of sight. These methods complement our techniques for joint cluster reconstruction from weak and strong gravitational lensing (Bartelmann et al. 1996; Cacciato et al. 2006; Merten et al. 2009; see also Coe et al. 2012; Merten et al. 2011; Merten 2014 for examples).

In this paper, we reconstruct the projected gravitational potential of the galaxy cluster MACS J1206.2−0847 based on galaxy kinematics, applying the technique developed in Sarli et al. (2014). This work is structured as follows: Sect. 2 briefly reviews the reconstruction method. In Sect. 3, we describe the selection and preparation of the data. Section 4 presents the results and compares the projected gravitational potential to that obtained from gravitational lensing. We summarise in Sect. 5 and discuss other methods in the Appendix.

2. Reconstruction method

The basic assumptions of our reconstruction method are the following: we treat the galaxy cluster as a spherically symmetric gas cloud of collisionless, pointlike test particles of mass m, i.e. galaxies, moving in the gravitational potential of their common dark matter halo. The system can then be described by the Jeans equation relating the radial velocity dispersion weighted by the galaxy number density ρ_gal to the gravitational potential φ: = Φ /m, (1)where the anisotropy parameter quantifies the ratio between the tangential and the radial velocity dispersions.

Expanding on the formal analogy with gas dynamics, we define an effective galaxy pressure . Using this definition and setting β to zero, (1) would turn into the equation of hydrostatic equilibrium for a gas. In addition, we introduce a polytropic relation between the effective galaxy pressure and the matter density, (2)This was justified for a simulated cluster in Sarli et al. (2014). In fact, the authors showed that the polytropic assumption is reasonable for several density profiles. Additionally, we shall see in our reconstruction results that the exact choice of the polytropic index only affects the final results very mildly.

The reconstruction algorithm described in detail in Sarli et al. (2014) proceeds along the following steps (see also Fig. 1):

• 1.

  In order to obtain the effective pressure P, we first have to deproject the actual observable, i.e. the line-of-sight projected velocity dispersion weighted by the galaxy number density, . This is accomplished via the Richardson-Lucy deconvolution or deprojection, see Lucy (1974, 1994).

• 2.

  Furthermore, we make use of the polytropic relation (2) to rewrite the Jeans Eq. (1) in terms of the effective pressure. This leads to a Volterra integral equation of the second kind for the gravitational potential Φ. After fixing the shape of the anisotropy profile, we can solve this equation in a quickly converging iteration process.

• 3.

  Finally, we project Φ along the line of sight to find the projected gravitational potential Ψ.

This entire algorithm depends on three parameters and a function, which is carefully adapted during the reconstruction. Noise suppression within the Richardson-Lucy deprojection requires regularisation, controlled by a smoothing scale L and an amplitude α. Furthermore, we have to choose the polytropic index γ in (2), usually being of order unity, as cluster analyses in Sarli et al. (2014) suggest. This also implies that the galaxy fluid can approximately be treated as an isothermal gas. The importance of the polytropic index is investigated in more detail in Sect. 4.1.

As we shall demonstrate later, the dominant parameter is the anisotropy profile β(r). However, if kinematic data are used exclusively for cluster reconstruction, there is in principle a degeneracy between the anisotropy parameter and the gravitational potential. We resolve this well-known degeneracy by fixing the β-profile to obtain a non-parametric estimate for the gravitational potential. To avoid an arbitrary, unmotivated guess for β(r), we choose it such that the reconstructed potential agrees best with a reconstruction based on gravitational-lensing data. Alternative methods for breaking the anisotropy-mass degeneracy and for recovering the gravitational potential with galaxy kinematics are discussed in the last section.

3. Data selection and preparation

The observational input data for the relaxed and massive galaxy cluster MACS J1206, Ebeling et al. (2001, 2009), (M_lens ≈ 5 × 10¹⁴M_⊙h^-1, Umetsu et al. 2012) investigated here consist of a catalogue of member galaxies (Biviano et al. 2013) observed in the context of the CLASH project (Postman et al. 2012) as part of a large spectroscopic campaign carried out with the Very Large Telescope (VLT; CLASH-VLT Large Programme; Rosati et al. 2014). According to Biviano et al. (2013) and Annunziatella et al. (2014), the spatial incompleteness of the spectroscopic sample varies with position in the cluster by less than 20%.

	Fig. 1 Scheme of the reconstruction algorithm. The observable is converted to a three-dimensional quantity by Richardson-Lucy deprojection. This is then turned into the three-dimensional gravitational potential, solving the Jeans equation. The resulting potential is finally projected along the line of sight.
Open with DEXTER

The following analysis is based on 592 member galaxies in total, selected via the method presented in Biviano et al. (2013). In order to arrive at the line-of-sight projected velocity dispersion profile, we assign to every cluster member a projected radius from the cluster centre, which we take to be marked by the brightest cluster galaxy. This is achieved by multiplying the angular separation with the angular-diameter distance, calculated with a standard ΛCDM-cosmology with recent cosmological parameters from Planck Collaboration XVI (2014).

We determine the line-of-sight projected velocities from the measured redshifts as described in Harrison (1974). In our case, we can use the non-relativistic relation between redshift z and line-of-sight velocity v_los, since all involved velocities (see below) are clearly in the non-relativistic regime. Thus, (3)As we are only interested in the galaxy velocities with respect to the cluster’s centre-of-mass, we have to take into account that the observed redshift includes four main contributions: the motion of the observer with respect to a local comoving observer, the motion of the cluster’s centre-of-mass with respect to a local comoving observer, the motion of each galaxy with respect to the cluster’s centre-of-mass, and finally the Hubble expansion. Harrison (1974) has shown that neglecting the motion of the cluster’s centre of mass leads to a multiplication of redshifts, (4)where z₀ is the redshift due to the observer’s motion with respect to a local comoving observer, z_cosm represents the Hubble expansion, and z_G the motion of each galaxy with respect to a local comoving observer. To leading order, we can calculate the cluster average of Eq. (4) and establish that the average redshift of all member galaxies with respect to the local comoving frame vanishes because we assume an isotropic distribution of velocities along the l.o.s., (5)Furthermore, we can safely assume that the contribution by the Hubble expansion is the same for all cluster members, (6)Thus, we end up with (7)Inserting Eq. (7) into Eq. (4), we find the following expression for v_los, see also Fig. 2: (8)To arrive at the density-weighted, projected velocity dispersion profile , we adapt the bin width to that chosen for the potential from lensing for better comparability, i.e. a constant bin width of 0.1 Mpc. Next, we calculate the velocity dispersion within each bin and reconstruct the projected gravitational potential as described above by solving the radial Jeans equation.

	Fig. 2 Velocities projected along the l.o.s. of cluster members.
Open with DEXTER

Because of the gauge freedom of the lensing potential, one is in particular allowed to add an arbitrary constant. By convention, we scale the potential such that Ψ(0) = 1 and use the normalisation Ψ(R_cut) = 0 for a given cut-off radius R_cut = 3 Mpc. This is even larger than the virial radius of MACS J1206, which is approximately 2 Mpc (Umetsu et al. 2012).

We perform a bootstrap analysis, i.e. for each bootstrap sample, we draw as many times with replacement from the original dataset as there are member galaxies to obtain error bars and to suppress the effect of outliers. Then, we reconstruct the gravitational potential for each individual bootstrap sample. This procedure is repeated 300 times. Finally, we calculate the mean of all potentials and their standard deviation. The bootstrap analysis also allows us to assess the effects of incompleteness and uncertain membership assignment by varying the sample of galaxies entering into the potential reconstruction.

4. Results

4.1. Effect of reconstruction parameters

As indicated above, we have to fix three reconstruction parameters to perform the reconstruction: the regularisation amplitude α, the smoothing scale L, and the polytropic index γ. Additionally, we have to model the anisotropy profile β(r).

Figure 3 shows how the variation of these reconstruction parameters affects the resulting potential Ψ. The amplitude of the anisotropy parameter β clearly has the dominant effect. In particular, an isotropic velocity distribution corresponding to β → 0 seems very implausible because of the strongly varying shape of the potential resulting from this assumption.

Although the assumption of a polytropic relation (2) may appear very bold, Fig. 3c shows that the particular choice of the polytropic index does not really matter. Also, the effects of the regularisation parameters α and L remain within the uncertainties, which is of similar order as shown in Fig. 6. When they are not varied, the parameters are set as follows: α = 0.3, L = 0.6, γ = 1.1, and β = 0.6. The value chosen for the polytropic index γ is motivated by the result of Sarli et al. (2014) that γ would usually be around unity. The regularisation parameters are chosen such that the agreement with the lensing reconstruction is best.

	Fig. 3 Reconstructed radial profiles of the projected gravitational potential of MACS J1206, obtained with different parameter sets for the smoothing amplitude α (panel a)), smoothing scale L (panel b)), the polytropic index γ (panel c)), and the anisotropy parameter β (panel d)).
Open with DEXTER

4.2. Variable anisotropy profile

We also systematically investigate the effect of an anisotropy parameter varying with radius. To accomplish this, we compare radial profiles of gravitational potentials obtained with an anisotropy parameter increasing or decreasing linearly in discrete steps every 0.5 Mpc. In case of the decreasing anisotropy-parameter profile, we choose β = 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, and reverse for the increasing case.

Figure 4 shows both cases together with the potential reconstruction from weak and strong gravitational lensing and a reconstruction with constant β = 0.6. First, one can observe that the qualitative form of the anisotropy profile affects the curvature of the potential: An increasing β profile leads to a potential with negative curvature. In contrast, a decreasing β-profile leads to a positive curvature and strong variations in the potential at radii beyond 1 Mpc. The reconstruction using a constant β leads to an almost straight line for the radial potential profile.

These results motivate the shape of the anisotropy profile used later in Eq. (9). Since the potential from lensing first shows a negative curvature at small radii and a linear evolution later, we choose an anisotropy profile increasing for small radii and turning constant afterwards to reproduce the linear behaviour.

	Fig. 4 Radial profiles of the projected gravitational potential, reconstructed from galaxy kinematics with different profiles of the anisotropy parameter β. The other reconstruction parameters remain fixed at α = 0.3, L = 0.6 and γ = 1.1.
Open with DEXTER

4.3. Comparison with lensing and X-ray data

As a first test, we compare the projected gravitational potential reconstructed from galaxy velocities with a reconstruction based on weak and strong gravitational lensing data of Merten et al. (2015) following constraints from Zitrin et al. (2012). After converting the 2D-lensing potential map (see Fig. 5) into a radial profile and using the same normalisation process as above, we can compare the result from gravitational lensing to the potential obtained from cluster kinematics (Fig. 6).

	Fig. 5 Two-dimensional map of the lensing potential of MACS1206 based on weak and strong lensing.
Open with DEXTER

The reconstruction parameters are chosen as listed above. For the radial profile of the anisotropy parameter β, we assume the following radially dependent profile, which is motivated in Sect. 4.2. This profile turns out to be in qualitative agreement with that in Biviano et al. (2013), i.e. (9)Additionally, we compare it with a reconstruction based on X-ray emission of the intracluster gas based on Donahue et al. (2014) using Chandra data up to 0.72 Mpc, which is a typical range of validity using this method for galaxy clusters.

The growing difference between the lensing and kinematic profile outside ~ 2 Mpc can be explained by the limited range of validity of each reconstruction method. The small difference at the innermost radii is due to increasing baryonic effects causing the equilibrium assumption ultimately to break down. The X-ray data agree well within the error boundaries.

	Fig. 6 Radial profiles of the gravitational potential of MACS J1206, reconstructed from combined strong and weak gravitational lensing (blue), X-ray emission (yellow), and from galaxy kinematics (red) with the parameters α = 0.3, L = 0.6, γ = 1.1 and β as in (9).
Open with DEXTER

4.4. Comparison with common lensing potentials

Furthermore, we can compare the radial profile of the projected gravitational potential obtained from galaxy kinematics with radial profiles common in gravitational-lensing studies, such as the singular (SIS) and the non-singular (NIS) isothermal spheres and the NFW profile (Navarro et al. 1996; Narayan & Bartelmann 1996); see Fig. 7, i.e.: (10)with the function h(x), depending on the radius in units of the scale radius x: = r/r_s (see Golse & Kneib (2002), i.e. (11)where s_c, r_s and A...F are free parameters determined by fitting the models (10) to the potential obtained from galaxy kinematics. The parameter values are summarised in Table 1.

We can compare the results for the scales s_c in case of the softened isothermal sphere and r_s in case of NFW with Umetsu et al. (2012). They obtain The agreement is quantified by a goodness-of-fit parameter Ω, which we define in analogy to the χ²-function appropriate for uncorrelated measurements. Since our data points are not independent but correlated by the bootstrap method, i.e. some data points appear multiple times, we cannot interpret Ω as a χ²-function. Because we just want to single out the best model, we do not go into a more elaborate analysis for correlated measurements here. For the ith data point x_i with standard deviation σ_i and the corresponding model prediction f(x_i), we write (14)As one can already see by eye, the softened isothermal sphere and the NFW-model describe the data equally well. More quantitatively, the goodness-of-fit parameters for all lensing model are (15)So far, all analytic lensing potentials are compatible with the reconstruction from galaxy kinematics, but the error bars are substantial.

	Fig. 7 Comparison of three analytic profiles of gravitational-lensing potentials with the gravitational-potential profile of MACS J1206 reconstructed from galaxy kinematics using α = 0.3, L = 0.6, γ = 1.1 and β as in Eq. (9).
Open with DEXTER

5. Summary

We have recently developed a method for reconstructing the projected gravitational potential of galaxy clusters from the kinematics of their member galaxies. In our study, we applied this method for the first time to a real galaxy cluster (MACS J1206). The reconstruction assumes that the three-dimensional gravitational potential of this cluster is approximately spherically symmetric. The potential reconstruction is based on 592 measured galaxy velocities. Our results can be summarised as follows: 1) The reconstruction algorithm requires four parameters to be set: the anisotropy parameter β; the polytropic index γ of the effective galaxy pressure; the amplitude α of the regularisation term in the Richardson-Lucy deconvolution and the length scale L of the regularisation. We found that the reconstructed potential is almost insensitive to the so-called nuisance parameters γ, α, and L, while it does depend very sensitively on the anisotropy parameter β. 2) Within the bootstrap error bars, the reconstructed potential profile is indistinguishable from the potential profile obtained from the combination of weak and strong gravitational lensing and from X-ray analysis. Typical analytic lensing profiles, such as the singular and non-singular isothermal spheres and the NFW profile, also agree well with our reconstruction from galaxy kinematics. 3) Small anisotropy parameters β ≲ 0.2 or anisotropy-parameter profiles decreasing with radius yield potential profiles differing strongly from the profile obtained from gravitational lensing. Best agreement with gravitational lensing is achieved with an anisotropy parameter increasing to β ≈ 0.6 within 0.5 Mpc and staying constant further out.

Our results show that the method returns convincing results on one well-studied galaxy cluster, which we use here as a test case. Even though the uncertainties are still large, the anisotropy parameter can be rather well constrained by comparing our results with those obtained from gravitational lensing. We see three main future applications of this method: (1) in joint reconstructions of galaxy-cluster potentials compatible with all cluster observables; (2) in constraints of the anisotropy parameter; and (3) in tests of fundamental assumptions, such as hydrostatic or virial equilibrium, and possibly also of general relativity, if applied to large and well-measured cluster samples.

Acknowledgments

This work was supported in part by the Collaborative Research Center TR 33 “The Dark Universe” of the Deutsche Forschungsgemeinschaft. This work is based on data collected at the ESO VLT (prog. ID 186.A-0798).

References

Appendix A: Discussion of other methods

On top of the conditions leading to the Jeans equation, the method presented above implicitly rests on the assumption that all test particles (i.e. galaxies) have the same mass, that their number density is proportional to their mass density and that one can establish a polytropic relation between their mass density ρ and the effective galaxy pressure . Fixing the velocity anisotropy profile β(r) enables us to constrain the gravitational potential of the cluster via the radial Jeans equation using the observed density-weighted galaxy velocity dispersion. Our method does not assume the mass to follow light a priori, i.e. ρ_DM does not necessarily need to follow ρ_gal. Since the resulting potential is non-parametric, one can test the validity of particular parametric models as done above.

Another common deprojection method is the Abel inversion using derivatives of observables. However, because of the fluctuating nature of most observables, which cause strongly fluctuating derivatives, we decided to use the Richardson-Lucy deprojection instead. This deprojection only involves integrals.

There are many other methods to determine the gravitational potential of a cluster, or its mass, or to constrain the velocity anisotropy. For example, Binney & Mamon (1982) use the line-of-sight velocity dispersions as observables together with the surface brightness in spherical galaxies. Assuming a constant mass-to-light ratio, they can replace the unknown density in the Jeans equation by the luminosity, and can thus determine the velocity anisotropy profile β(r). However, in the context of galaxy clusters it is not clear that the same relation between mass and luminosity should hold as well.

Dejonghe & Merritt (1992) determine the potential Φ by taking into account higher than second-order velocity moments of the collision-less Boltzmann equation (CBE). Given the velocity dispersions, and making use of a linear relation between the observables and Φ, they expand the gravitational potential in terms of basis functions with unspecified coefficients. Using the positivity constraint of the distribution function and higher order moments of the CBE, they are able to determine the coefficients and thus the potential. However, from a practical point of view, it becomes increasingly difficult to constrain higher-order moments of the velocity distribution from their observable line-of-sight projections, as the projection integrals turn out to be substantially more complex with increasing order (see Richardson & Fairbairn 2013, for instance).

A method very similar to ours is proposed by Mamon & Boué (2010). Starting again from the Jeans equation and fixing β(r), they determine both the number density and the velocity dispersion observationally. Besides their using the Abel inversion, the only difference from our method is that they determine the number density of galaxies appearing in the Jeans equation by counting. Hence, they do not need our assumed polytropic relation between density and effective pressure. In future work, one could use this to test directly to what degree the polytropic assumption is appropriate.

Wolf et al. (2010) and van der Marel et al. (2000) are considering a general parametrisation of β(r) and of the mass M(r) and constrain all free parameters via a maximum-likelihood analysis of the velocity dispersion predicted by the Jeans equation compared to the observed dispersion. van der Marel et al. (2000) use a constant velocity anisotropy whereas Wolf et al. (2010) take a varying parametric profile.

Another way to break the anisotropy mass degeneracy is taking the fourth order velocity moments of the CBE into account (Łokas (2002); see also Richardson & Fairbairn (2013) in this context). By expressing the velocity distribution function in terms of two integrals of motion, the energy E and the angular momentum L, and further assuming that the distribution is separable in E and L with β = const., the two fourth-order moments are reduced to one equation involving the radial velocity to the fourth power, β, the radial velocity dispersion and Φ. Thus, by measuring the velocity dispersion and , one can infer the gravitational potential by solving this fourth-order equation together with the Jeans equation.

The MAMPOSSt method introduced by Mamon et al. (2013) breaks the degeneracy assuming parametric forms for the gravitational potential, the velocity anisotropy, and the distribution of 3D-velocities. By performing a maximum-likelihood fit of the galaxy distribution in the projected phase space, all parameters can be determined. However, it has recently been shown (Richardson & Fairbairn (2013) that their assumption of a Gaussian 3D-velocity distribution is incompatible with the equilibrium assumption underlying the CBE. One would thus have to choose a different 3D-velocity distribution.

All Tables

All Figures

	Fig. 1 Scheme of the reconstruction algorithm. The observable is converted to a three-dimensional quantity by Richardson-Lucy deprojection. This is then turned into the three-dimensional gravitational potential, solving the Jeans equation. The resulting potential is finally projected along the line of sight.
Open with DEXTER
In the text

	Fig. 2 Velocities projected along the l.o.s. of cluster members.
Open with DEXTER
In the text

	Fig. 3 Reconstructed radial profiles of the projected gravitational potential of MACS J1206, obtained with different parameter sets for the smoothing amplitude α (panel a)), smoothing scale L (panel b)), the polytropic index γ (panel c)), and the anisotropy parameter β (panel d)).
Open with DEXTER
In the text

	Fig. 4 Radial profiles of the projected gravitational potential, reconstructed from galaxy kinematics with different profiles of the anisotropy parameter β. The other reconstruction parameters remain fixed at α = 0.3, L = 0.6 and γ = 1.1.
Open with DEXTER
In the text

	Fig. 5 Two-dimensional map of the lensing potential of MACS1206 based on weak and strong lensing.
Open with DEXTER
In the text

	Fig. 6 Radial profiles of the gravitational potential of MACS J1206, reconstructed from combined strong and weak gravitational lensing (blue), X-ray emission (yellow), and from galaxy kinematics (red) with the parameters α = 0.3, L = 0.6, γ = 1.1 and β as in (9).
Open with DEXTER
In the text

	Fig. 7 Comparison of three analytic profiles of gravitational-lensing potentials with the gravitational-potential profile of MACS J1206 reconstructed from galaxy kinematics using α = 0.3, L = 0.6, γ = 1.1 and β as in Eq. (9).
Open with DEXTER
In the text